%I #9 Mar 29 2013 14:32:03
%S 1,2,3,5,11,19,37,67,127,223,383,659,1151,1987,3433,5923,10211,17579,
%T 30259,52103,89689,154387,265747,457421,787357,1355261,2332783,
%U 4015339,6911461,11896457,20476933,35246203,60668017,104425627,179743997
%N Let a(1)=1; for n>1, a(n)=nextprime( a(n-1)^(n/(n-1)) ).
%C Limit of a(n+1)/a(n) -> 1.7212631438796...
%t RecurrenceTable[{a[1]==1,a[n]==NextPrime[a[n-1]^(n/(n-1))]},a,{n,40}] (* _Harvey P. Dale_, Mar 29 2013 *)
%o (PARI) a(n)=if(n<2,1,nextprime(a(n-1)^(n/(n-1)))); for(n=1,50,print1(a(n),","))
%K nonn
%O 1,2
%A _Paul D. Hanna_, May 31 2003
%E Wrong formula removed by _Paul D. Hanna_, May 29 2010
|