|
%I #5 Jul 07 2021 08:20:32
%S 9,341,1105,1541,1729,1891,2465,29341,162401,252601,294409,334153,
%T 340561,399001,410041,488881,530881,636641,954271,1024651,1152271,
%U 1193221,1461241,1615681,1857241
%N Records in A098650.
%t PrimeFactors[n_] := Flatten[ Table[ # [[1]], {1} ] & /@ FactorInteger[ n ]]; f[n_] := Block[{k = n + 1}, If[EvenQ[k], k++ ]; While[ PrimeQ[k] || Union[ PowerMod[ PrimeFactors[n], k - 1, k]] != {1}, k += 2]; k]; a = {1}; b = {9}; Do[ If[ SquareFreeQ[n], c = f[n]; If[c > b[[ -1]], AppendTo[a, n]; AppendTo[b, c]; Print[n]]], {n, 2, 145000}]; b
%Y Cf. A098650, A098651.
%K nonn,more
%O 1,1
%A _Robert G. Wilson v_, Sep 18 2004
%E a(13)-a(25) from _Amiram Eldar_, Jul 07 2021
|
