|
%I #7 Dec 05 2013 19:55:31
%S 1,5,11,19,29,41,59,79,101,127,157,191,227,269,313,359,409,461,521,
%T 587,659,733,809,887,967,1049,1151,1259,1373,1489,1607,1733,1861,1993,
%U 2129,2267,2411,2557,2707,2861,3019,3181,3347,3517,3691,3877,4073,4271,4481
%N a(0) = 1, c(0) = 1, a(n) for n > 0 is the smallest prime a(n-1) + c(n), where c(n) is composite and larger than c(n-1).
%C Same as A070865 after first term. - _David Wasserman_, Jun 27 2005
%e a(6) = 59 since 59 = a(5) + c(6) = 41 + 18, 18 is composite and larger than c(5) = 12 and the composite numbers 14, 15, 16 do not lead to a prime when added to 41.
%t a[n_] := a[n] = NextPrim[a[n - 1] + a[n - 1] - a[n - 2]]; a[0] = 1; a[1] = 5; Table[ a[n], {n, 0, 50}]
%o (PARI) {print1(a=1,","); c=1; for(n=1,48,c++; while(isprime(c)||!isprime(a+c),c++); a=a+c; print1(a,","))}
%K nonn
%O 0,2
%A _Amarnath Murthy_, Aug 14 2002
%E Edited and extended by _Klaus Brockhaus_ and _Robert G. Wilson v_, Aug 14 2002
|
