%I #40 Aug 23 2014 13:06:34
%S 2,6,8,12,16,10,4,30,26,22,24,14,50,42,18,64,46,60,32,36,20,34,28,108,
%T 48,44,68,282,90,54,76,62,180,66,132,86,74,38,58,106,120,52,244,94,
%U 100,82,138,156,98,72,172,150,248,154,166,114,162,126,124,208,222,324,212
%N a(1)=2; for n >= 2, a(n) = smallest even number such that a(1)*a(2)*...*a(n) + 1 is prime.
%C Is this a permutation of the even numbers?
%C For any even positive integers a_1, a_2, ..., a_n, there are infinitely many even positive integers t such that a_1 a_2 ... a_n t + 1 is prime: this follows from Dirichlet's theorem on primes in arithmetic progressions. As far as I know there is no guarantee that the sequence defined here leads to a permutation of the even numbers, i.e. there might be some even integer that never appears in the sequence. However, if the partial products a_1 ... a_n grow like 2^n n!, heuristically the probability of a_1 ... a_n t + 1 being prime is on the order of 1/log(a_1 ... a_n) ~ 1/(n log n), and since sum_n 1/(n log n) diverges we might expect that there should be infinitely many n for which some a_1 ... a_n t + 1 is prime, and thus every even integer should occur. - _Robert Israel_, Dec 20 2012
%H Robert Israel, <a href="/A083769/b083769.txt">Table of n, a(n) for n = 1..200</a>
%e 2+1=3, 2*6+1=13, 2*6*8+1=97, 2*6*8*12+1=1153, etc. are primes.
%e After 200 terms the prime is
%e 224198929826405912196464851358435330956778558123234657623126\
%e 069546460095464785674042966210907411841359152393200850271694\
%e 899718487202330385432243578646330245831108247815285116235792\
%e 875886417750289946171599027675234787802312202111702704952223\
%e 563058999855839876391430601719636148884060097930252529666254\
%e 756431522481046758186320659298713737639441014068272279177710\
%e 551232067814381240340990584869121776471244800000000000000000\
%e 00000000000000000000000000000 (449 digits). - _Robert Israel_, Dec 21 2012
%p N := 200: # number of terms desired
%p P := 2:
%p a[1] := 2:
%p C := {seq(2*j, j = 2 .. 10)}:
%p Cmax := 20:
%p for n from 2 to N do
%p for t in C do
%p if isprime(t*P+1) then
%p a[n]:= t;
%p P:= t*P;
%p C:= C minus {t};
%p break;
%p end if;
%p end do;
%p while not assigned(a[n]) do
%p t0:= Cmax+2;
%p Cmax:= 2*Cmax;
%p C:= C union {seq(j, j=t0 .. Cmax, 2)};
%p for t from t0 to Cmax by 2 do
%p if isprime(t*P+1) then
%p a[n]:= t;
%p P:= t*P;
%p C:= C minus {t};
%p break;
%p end if
%p end do;
%p end do;
%p end do;
%p [seq(a[n],n=1..N)];
%t f[s_List] := Block[{k = 2, p = Times @@ s}, While[ MemberQ[s, k] || !PrimeQ[k*p + 1], k += 2]; Append[s, k]]; Nest[f, {2}, 62] (* _Robert G. Wilson v_, Dec 24 2012 *)
%Y Cf. A036013, A046966, A046972, A051957, A073673, A073674, A083769, A083770, A083771, A084401, A084402, A084724, A087338.
%K nonn
%O 1,1
%A _Amarnath Murthy_ and Meenakshi Srikanth (menakan_s(AT)yahoo.com), May 06 2003
%E More terms from _David Wasserman_, Nov 23 2004
%E Edited by _N. J. A. Sloane_, Dec 20 2012
%E Comment edited, Maple code and additional terms by _Robert Israel_, Dec 20 2012