OFFSET
0,5
COMMENTS
For all p prime, a(p) == 0 (mod p*(p-1)).
FORMULA
E.g.f.: (-1 + 2*x - 2*x^2 + x^3 + (1 - x)*(log((1 - x)^(1 - 2*x)) - (log(1 - x))^2))*exp(x).
a(n) ~ 2 * exp(1) * log(n) * n! / n^2 * (1 + (gamma - 3/2)/log(n)), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Dec 09 2021
EXAMPLE
E.g.f.: -1 - 3*x^4/4! - 20*x^5/5! - 70*x^6/6! - 84*x^7/7! + 1267*x^8/8! + 18824*x^9/9! + ...
a(11) = Sum_{k=0..5} (-1)^(k-1)*(k-1)^2*A106828(11, k).
a(11) = (-1)*1*0 + (1)*0*3628800 + (-1)*1*6636960 + (1)*4*3678840 + (-1)*9*705320 + (1)*16*34650 = 2284920.
For k = 0, A106828(11,0) = 0.
For k = 1, (1-1)^2 = 0.
For 2 <= k <= 5, A106828(11, k) == 0 (mod 11*10).
Result a(11) == 0 (mod 11*10).
MAPLE
a := series((-1+2*x-2*x^2+x^3+(1-x)*(log((1-x)^(1-2*x))-(log(1-x))^2))*exp(x), x=0, 24):
seq(n!*coeff(a, x, n), n=0..23);
# second program:
a := n -> add((-1)^(k-1)*(k-1)^2*A106828(n, k), k=0..iquo(n, 2)):
seq(a(n), n=0..23);
MATHEMATICA
CoefficientList[Series[(-1+2*x-2*x^2+x^3+(1-x)*(Log[(1-x)^(1-2*x)]-(Log[1-x])^2))*Exp[x], {x, 0, 23}], x]*Range[0, 23]!
PROG
(PARI) my(x='x+O('x^30)); Vec(serlaplace((-1 + 2*x - 2*x^2 + x^3 + (1 - x)*(log((1 - x)^(1 - 2*x)) - (log(1 - x))^2))*exp(x))) \\ Michel Marcus, Oct 07 2021
CROSSREFS
KEYWORD
sign
AUTHOR
Mélika Tebni, Oct 07 2021
STATUS
approved