close
login
A258807
a(n) = n^5 - 1.
7
0, 31, 242, 1023, 3124, 7775, 16806, 32767, 59048, 99999, 161050, 248831, 371292, 537823, 759374, 1048575, 1419856, 1889567, 2476098, 3199999, 4084100, 5153631, 6436342, 7962623, 9765624, 11881375, 14348906, 17210367, 20511148, 24299999, 28629150, 33554431
OFFSET
1,2
FORMULA
G.f.: x^2*(31 + 56*x + 36*x^2 - 4*x^3 + x^4)/(1 - x)^6.
a(n) = 6*a(n-1) - 15*a(n-2) + 20*a(n-3) - 15*a(n-4) + 6*a(n-5) - a(n-6).
a(n) = -A024003(n). - Bruno Berselli, Jun 11 2015
Sum_{n>=2} 1/a(n) = Sum_{n>=1} (zeta(5*n) - 1) = 0.0379539032... - Amiram Eldar, Nov 06 2020
MAPLE
seq(n^5-1, n=1..35); # Muniru A Asiru, Oct 28 2018
MATHEMATICA
Table[n^5 - 1, {n, 1, 50}] (* or *) LinearRecurrence[{6, -15, 20, -15, 6, -1}, {0, 31, 242, 1023, 3124, 7775}, 50]
PROG
(Magma) [n^5-1: n in [1..50]];
(Magma) I:=[0, 31, 242, 1023, 3124, 7775]; [n le 6 select I[n] else 6*Self(n-1)-15*Self(n-2)+20*Self(n-3)-15*Self(n-4)+ 6*Self(n-5)-Self(n-6): n in [1..50]];
(SageMath) [n^5-1 for n in (1..50)] # Bruno Berselli, Jun 11 2015
(PARI) a(n)=n^5-1 \\ Charles R Greathouse IV, Jun 11 2015
(GAP) List([1..35], n->n^5-1); # Muniru A Asiru, Oct 28 2018
(Python) for n in range(1, 50): print(n**5 - 1, end=', ') # Stefano Spezia, Oct 28 2018
CROSSREFS
Subsequence of A181124.
Sequences of the type n^k-1: A132411 (k=2), A068601 (k=3), A123865 (k=4), this sequence (k=5), A123866 (k=6), A258808 (k=7), A258809 (k=8), A258810 (k=9), A123867 (k=10), A258812 (k=11), A123868 (k=12).
Sequence in context: A189923 A059378 A024003 * A358934 A221848 A344723
KEYWORD
nonn,easy
AUTHOR
Vincenzo Librandi, Jun 11 2015
STATUS
approved