|
PROG
|
(PARI) pol(n) = {if (n<2, return (n)); if (n%2, pol((n+1)/2) + pol((n-1)/2), x*pol(n/2)); } \\ A125184
ispow2(n) = if ((n==1) || (n==2), return (1)); my(p); isprimepower(n, &p) && (p==2);
istrivial(n) = ispow2(n+1) || ispow2(n+3);
isokrm(n, r, m) = {if ((n%2) && !istrivial(n), my(p=pol(n), d=poldegree(p)); Mod(p, m) == Mod(1+r*sum(i=1, d, x^i), m); ); }
lista(nn) = forstep(n=1, nn, 2, if (isokrm(n, 1, 3), print1(n, ", ")));
|