login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A094003
a(1) = 2, a(n+1) = (a(n)^(n+1) - 1)/(a(n) - 1).
1
2, 3, 13, 2380, 32098914298781, 34076243757937210559857502330309024747435410425446137941014310689706
OFFSET
1,1
COMMENTS
a(7) has 406 digits and is too large to include.
LINKS
MAPLE
a[1]:=2: for n from 2 to 7 do a[n]:=(a[n-1]^n-1)/(a[n-1]-1) od:seq(a[n], n=1..7); # Emeric Deutsch, Apr 17 2005
MATHEMATICA
RecurrenceTable[{a[1] == 2, a[n] == (a[n-1]^n - 1)/(a[n-1] -1)}, a, {n, 6}] (* Vincenzo Librandi, Dec 20 2015 *)
PROG
(PARI) { a=2; write("b094003.txt", "1 2"); for (n=2, 7, a=(a^n-1)/(a-1); write("b094003.txt", n, " ", a); ) } \\ Harry J. Smith, Jun 16 2009
(Magma) [n le 1 select 2 else (Self(n-1)^n - 1)/(Self(n-1) - 1): n in [1..6]]; // Vincenzo Librandi, Dec 20 2015
CROSSREFS
Sequence in context: A258357 A160858 A062648 * A119987 A093553 A253575
KEYWORD
nonn
AUTHOR
Amarnath Murthy, May 29 2004
EXTENSIONS
More terms from Emeric Deutsch, Apr 17 2005
STATUS
approved