OFFSET
1,1
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Knödel Numbers.
MAPLE
with(numtheory);
knodel:= proc(i, k)
local a, n, ok;
for n from k+1 to i do
ok:=1;
for a from 1 to n do
if gcd(a, n)=1 then if (a^(n-k) mod n)<>1 then ok:=0; break; fi; fi;
od;
if ok=1 then print(n); fi;
od;
end:
knodel(10000, 6);
MATHEMATICA
knodelQ[m_Integer?PrimeQ, n_Integer] := False; knodelQ[m_Integer, n_Integer] := Module[{i = n + 1}, While[i < m && (GCD[i, m] > 1 || Mod[i^(m - n), m] == 1), i++]; (i == m)]; Select[Range[1000], knodelQ[#, 6] &] (* Alonso del Arte, Feb 24 2012 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Paolo P. Lava, Feb 24 2012
STATUS
approved