A268497 Records in A268630 (= a(n)^2+a(n+1) is prime).
0, 2, 4, 8, 10, 16, 18, 23, 28, 32, 37, 48, 63, 66, 73, 76, 79, 82, 99, 100, 102, 106, 109, 126, 142, 147, 169, 172, 177, 178, 189, 190, 198, 201, 205, 208, 217, 225, 226, 228, 232, 235, 244, 253, 258, 264, 268, 278, 288, 306, 333, 336, 343, 348, 352, 358, 361, 370, 379, 382
Offset: 0
Keywords
Programs
-
PARI
show(n, f="/tmp/b268497.txt", c=0, a=0, u=[a-1])={for(n=1, n, a>=u[#u]&&if(f, write(f, c++-1, " "a), print1(a", ")); u=setunion(u, [a]); while(#u>1&&u[2]==u[1]+1, u=u[^1]); for(k=u[1]+1, 9e9, !setsearch(u, k) && isprime(a*a+k) && (a=k) && break))}
Comments