%I #20 Nov 25 2015 20:48:14
%S 83,2,167,3,7,67,5,13,719,37,11,31,1367,31440216015620321911,
%T 988487183108868589955299792587646370011,19,499,937,23,29,863,
%U 1852812869,1157393039341097158113740816072403984079,445807,13833892959043,359,109,331
%N Euclid-Mullin sequence (A000945) with initial value a(1)=83 instead of a(1)=2.
%H Robert Price, <a href="/A051327/b051327.txt">Table of n, a(n) for n = 1..64</a>
%t a[1]=83; a[n_] := First[ Flatten[ FactorInteger[ 1+Product[ a[ j ], {j, 1, n-1} ] ] ] ]; Array[a, 10]
%o (PARI) gpf(n)=my(f=factor(n)[, 1]); f[#f];
%o first(m)=my(v=vector(m)); v[1]=83; for(i=2, m, v[i]=gpf(1+prod(j=1, i-1, v[j]))); v \\ _Anders Hellström_, Nov 25 2015
%Y Cf. A000945, A000946, A005265, A005266.
%K nonn
%O 1,1
%A _Labos Elemer_
%E More terms from _Sean A. Irvine_, Apr 12 2011
|