A331338 E.g.f.: Sum_{k>=1} (1 - sech(x^k)) (even powers only).

1, 7, 421, 10375, 1864921, 177588115, 43788506581, 5364386887855, 3746015060385841, 743097477083711275, 562069712763196977901, 208700966175150043937095, 201649817635812110915892361, 89621739407126401163202051715, 158792595367325562826629059282821

Offset

1,2

Links

Table of n, a(n) for n=1..15.

Formula

a(n) = -(2*n)! * Sum_{d|n} A028296(d) / (2*d)!.

Mathematica

nmax = 15; Table[(CoefficientList[Series[Sum[1 - Sech[x^k], {k, 1, nmax}], {x, 0, 2 nmax}], x] Range[0, 2 nmax]!)[[n]], {n, 1, 2 nmax + 1, 2}] // Rest
Table[-(2 n)! DivisorSum[n, EulerE[2 #]/(2 #)! &], {n, 1, 15}]

Crossrefs

Cf. A000364, A028296, A176474, A330255, A330528.

Sequence in context: A287033 A367867 A362677 * A269555 A027876 A158599

Adjacent sequences: A331335 A331336 A331337 * A331339 A331340 A331341

Keyword

nonn

Author

Ilya Gutkovskiy, Jan 14 2020

Status

approved