login
A069975
a(n) = n*(16*n^2 - 1).
2
15, 126, 429, 1020, 1995, 3450, 5481, 8184, 11655, 15990, 21285, 27636, 35139, 43890, 53985, 65520, 78591, 93294, 109725, 127980, 148155, 170346, 194649, 221160, 249975, 281190, 314901, 351204, 390195, 431970, 476625, 524256, 574959, 628830, 685965, 746460
OFFSET
1,1
COMMENTS
Sequence could or should start with a(0) = 0. - M. F. Hasler, Sep 06 2025
FORMULA
Sum_{n>=1} 1/a(n) = 3*log(2) - 2 = A016631 - 2. (Ramanujan)
Sum_{n>=1} (-1)^(n+1)/a(n) = 2 - log(2) + sqrt(2)*log(sqrt(2)-1). - Amiram Eldar, Jun 24 2022
From Elmo R. Oliveira, Sep 05 2025: (Start)
G.f.: 3*x*(5 + 22*x + 5*x^2)/(x-1)^4.
E.g.f.: x*(15 + 48*x + 16*x^2)*exp(x).
a(n) = 4*a(n-1) - 6*a(n-2) + 4*a(n-3) - a(n-4) for n > 4.
a(n) = A069140(n)/4. (End)
MATHEMATICA
Table[n(16n^2-1), {n, 40}] (* Harvey P. Dale, Dec 17 2018 *)
PROG
(PARI) a(n) = n*(16*n^2-1); \\ Michel Marcus, Nov 25 2013
(PARI) my(x='x+O('x^37)); Vec(3*x*(5+22*x+5*x^2)/(1-x)^4) \\ Elmo R. Oliveira, Sep 05 2025
(PARI) apply( {A069975(n)=n^3<<4-n}, [0..44]) \\ M. F. Hasler, Sep 06 2025
(Python) def A069975(n): return (n**3<<4)-n # M. F. Hasler, Sep 06 2025
CROSSREFS
Sequence in context: A171220 A071080 A193365 * A027779 A337958 A337957
KEYWORD
easy,nonn
AUTHOR
Benoit Cloitre, Apr 30 2002
STATUS
approved