|
|
A268496
|
|
Indices of records in A268630 (= a(n)^2+a(n+1) is prime).
|
|
3
|
|
|
0, 1, 3, 5, 7, 9, 15, 16, 17, 23, 24, 29, 40, 49, 52, 63, 72, 73, 74, 77, 83, 89, 92, 93, 105, 110, 130, 145, 154, 161, 162, 165, 175, 182, 184, 189, 190, 198, 205, 213, 217, 218, 219, 222, 237, 239, 247, 251, 253, 255, 278, 293, 308, 311, 313, 331, 336, 339, 346, 351, 353, 362, 366
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
See A268497 for the actual record values.
|
|
LINKS
|
|
|
PROG
|
(PARI) show(n, f="/tmp/b268496.txt", c=0, a=0, u=[a-1])={for(n=1, n, a>=u[#u]&&if(f, write(f, c++-1, " "n-1), print1(n-1", ")); 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))}
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|