|
|
|
|
3, 3, 307, 307, 11, 3, 23, 3, 3, 467, 3, 467, 467, 37, 2, 5, 31, 3, 127, 127, 127, 127, 127, 1667, 73, 19, 19, 19, 19, 19, 31, 31, 19, 31, 19, 19, 31, 31, 31, 31, 31, 31, 19, 31, 3, 3, 31, 31, 3, 3, 3, 3, 19, 19, 31, 31, 73, 73, 3, 3, 3, 3, 3, 3, 3, 19, 19, 31
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
PROG
|
(PARI) addtovec(vec) = my(w=[], vmax=0); for(t=1, #vec, if(vecmax(vec[t]) > vmax, vmax=vecmax(vec[t]))); for(k=1, #vec, forprime(q=1, vmax, if(Mod(vec[k][#vec[k]], q^2)^(q-1)==1, w=concat(w, [0]); w[#w]=concat(vec[k], [q])))); w
removefromvec(vec) = my(w=[]); for(k=1, #vec, if(vecsort(vec[k])==vecsort(vec[k], , 8), w=concat(w, [0]); w[#w]=vec[k])); w
printfromvec(vec) = for(k=1, #vec, if(vec[k][1]==vec[k][#vec[k]], print1(vec[k][2], ", ")))
forprime(p=1, , my(v=[[p]]); while(#v > 0, v=addtovec(v); printfromvec(v); v=removefromvec(v)))
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|