OFFSET
1,1
COMMENTS
LINKS
Michael De Vlieger, Table of n, a(n) for n = 1..10000
MATHEMATICA
rad[x_] := Times @@ FactorInteger[x][[All, 1]];
a332785 = Select[Range[270], Nor[SquareFreeQ[#], Divisible[#, rad[#]^2]] &];
Complement[a332785, Select[a332785, Function[k, r = rad[k];
s = Select[Range[k], CoprimeQ[#, k] &];
t = Union@ Flatten@ Mod[
TensorProduct @@
Map[(p = #; NestWhileList[Mod[p*#, k] &, 1, UnsameQ, All]) &,
FactorInteger[k][[All, 1]]], k];
Count[Complement[Range[0, k - 1],
Union[s, t]], _?(Divisible[#, r] &)] > 0] ] ]
CROSSREFS
KEYWORD
nonn
AUTHOR
Michael De Vlieger, Jun 14 2026
STATUS
approved
