A227000 - OEIS

%I #7 Jun 26 2013 17:19:35

%S 1,2,3,6,9,11,14,23,30,41,48,63,74,81,86,90,95,105,119,125,128,140,

%T 153,156,158,165,179,186,191,209,216,219,224,233,245,251,296,303,308,

%U 315,321,350,354,359,375,398,405,419,429,441,443,468,485,506,524,531,543,546,576

%N Numbers n such that (n+1)^2-n^2 and (n+1)^3-n^3 are both prime.

%e n=23; n+1=24; 24^2-23^2=47 and 24^3-23^3=1657.

%t Select[Range[576], PrimeQ[(# + 1)^2 - #^2] && PrimeQ[(# + 1)^3 - #^3] &] (* _T. D. Noe_, Jun 26 2013 *)

%o (PARI) forprime(p=3,1e3,n=p\2;if(isprime(3*n*(n+1)+1),print1(n", "))) \\ _Charles R Greathouse IV_, Jun 26 2013

%K nonn

%O 1,2

%A _César Aguilera_, Jun 26 2013