%I #17 Dec 07 2015 01:54:42
%S 73,2,3,439,7,5,6729871,103,23,92581,13,19,2453563465139998636061911,
%T 739,3167,47356379285063777,463,3673,2918137,41,17,
%U 2307841395358410336056217199460000033494100011180619106269258300884070797073446703818111
%N Euclid-Mullin sequence (A000945) with initial value a(1)=73 instead of a(1)=2.
%H Robert Price, <a href="/A051325/b051325.txt">Table of n, a(n) for n = 1..44</a>
%t a[1]=73; a[n_] := First[ Flatten[ FactorInteger[ 1+Product[ a[ j ], {j, 1, n-1} ] ] ] ]; Array[a, 10]
%o (PARI) spf(n)=my(f=factor(n)[1, 1]); f
%o first(m)=my(v=vector(m)); v[1]=73; for(i=2, m, v[i]=spf(1+prod(j=1, i-1, v[j]))); v \\ _Anders Hellström_, Dec 06 2015
%Y Cf. A000945, A000946, A005265, A005266.
%K easy,nonn
%O 1,1
%A _Labos Elemer_
|