gp > forprime(p=3,30,print1(p," ");f=polcyclo(p^2);for(m=2,500,if(ispseudoprime(subst(f,x,m))==1,print1(m," ")));print)
3 2 3 8 11 20 21 26 30 50 51 56 60 78 98 102 117 129 134 146 159 171 186 189 191 198 200 209 210 212 222 240 249 267 269 278 279 299 300 333 344 363 383 390 398 399 425 429 438 444 450 458
5 22 33 39 43 62 74 134 142 167 212 238 287 313 335 369 414 415 418 432
7 2 5 62 67 104 123 138 151 186 194 233 236 261 269 306 359 416 425 439
11 43 174 236 261 293 300 306 312
13 24 120 205 262 408