A050872 a(n) = (1/2)*A050871 (row sums of array T in A050870, periodic binary words).

0, 1, 2, 5, 8, 17, 38, 65, 128, 284, 518, 1025, 2168, 4097, 8198, 16907, 32768, 65537, 133088, 262145, 524408, 1056731, 2097158, 4194305, 8421248, 16777712, 33554438, 67239680, 134217848, 268435457, 537396698, 1073741825, 2147483648

Offset

0,3

Links

Table of n, a(n) for n=0..32.

Prog

(PARI) T(n, k) = binomial(n, k) - sumdiv(gcd(n+!n, k), d, moebius(d)*binomial(n/d, k/d)); \\ A050870
row(n) = vector(n+1, k, k--; T(n, k));
a(n) = n*=2; vecsum(row(n))/2; \_ Michel Marcus, Aug 20 2021

Crossrefs

Cf. A050870, A050871.

Sequence in context: A181586 A112361 A343129 * A086324 A293830 A073708

Adjacent sequences: A050869 A050870 A050871 * A050873 A050874 A050875

Keyword

nonn

Author

Clark Kimberling

Extensions

a(29) onward corrected by Sean A. Irvine, Aug 20 2021

Status

approved