%I #23 Mar 28 2019 13:13:34
%S 15,49,91,133,217,259,301,427,469,511,553,679,721,763,889,973,1015,
%T 1057,1099,1141,1267,1351,1393,1477,1561,1603,1687,1897,1939,1981,
%U 2107,2149,2191,2317,2359,2443,2569,2611,2653,2779,2863,2947,3031,3073,3199,3241,3409
%N 7-Knödel numbers.
%H Amiram Eldar, <a href="/A208155/b208155.txt">Table of n, a(n) for n = 1..10000</a>
%H Eric Weisstein, <a href="http://mathworld.wolfram.com/KnoedelNumbers.html">Knödel Numbers</a>.
%p with(numtheory);
%p knodel:= proc(i,k)
%p local a,n,ok;
%p for n from k+1 to i do
%p ok:=1;
%p for a from 1 to n do
%p if gcd(a,n)=1 then if (a^(n-k) mod n)<>1 then ok:=0; break; fi; fi;
%p od;
%p if ok=1 then print(n); fi;
%p od;
%p end:
%p knodel(10000,7);
%t (* First run program for A208154 to define knodelQ *) Select[Range[3500], knodelQ[#, 7] &] (* _Alonso del Arte_, Feb 24 2012 *)
%Y Cf. A002997, A050990, A033553, A050992, A050993, A208154, A208156-A208158.
%K nonn
%O 1,1
%A _Paolo P. Lava_, Feb 24 2012
|