A332969 a(n) = [x^n] (Sum_{j>=0} A002193(1-j) * x^j)^2.

1, 8, 18, 16, 37, 26, 34, 52, 70, 90, 87, 116, 127, 112, 157, 212, 158, 192, 252, 252, 249, 272, 349, 276, 287, 478, 482, 334, 407, 478, 465, 488, 544, 698, 562, 504, 682, 698, 738, 736, 742, 880, 907, 826, 944, 848, 998, 1110, 976, 1106, 1217, 1112, 1060

Offset

0,2

Links

Andrew Howroyd, Table of n, a(n) for n = 0..10000

Formula

G.f.: (Sum_{j>=0} A002193(1-j) * x^j)^2.

Sum_{k>=0} a(k)/10^k = 2.

a(n) = Sum_{j=0..n} A002193(1-j)*A002193(j-n+1).

Example

a(1) = 8 because the coefficient of x^1 in (1 + 4x + ... )^2 is 8.

Prog

(PARI) seq(n)={Vec(Ser(digits(sqrtint(2*100^n)))^2)} \\ Andrew Howroyd, Mar 04 2020

Crossrefs

Cf. A002193.

Sequence in context: A333828 A133202 A217424 * A386598 A322031 A196207

Adjacent sequences: A332966 A332967 A332968 * A332970 A332971 A332972

Keyword

nonn,base,easy

Author

Andrew Slattery, Mar 04 2020

Status

approved