|
PROG
|
(PARI) is(p)={my(W, V1, V2, V3, V4, q1, q2, q3, q4, i1, i2, i3, i4, l1, l2, l3, l4); W=0; V1=factor(p+1); l1=length(V1[, 1]); for(i1=1, l1, q1=V1[i1, 1]; V2=factor(q1^2+q1+1); l2=length(V2[, 1]); for(i2=1, l2, q2=V2[i2, 1]; V3=factor(q2^2+q2+1); l3=length(V3[, 1]); for(i3=1, l3, q3=V3[i3, 1]; V4=factor(q3^2+q3+1); l4=length(V4[, 1]); for(i4=1, l4, q4=V4[i4, 1]; if(q4==p, W=[p, q1, q2, q3]))))); W}
|