OFFSET
1,2
MATHEMATICA
Select[ Range[ 10^6 ], Mod[ PowerMod[ 7, #, # ] + PowerMod[ 6, #, # ] + PowerMod[ 5, #, # ] + PowerMod[ 4, #, # ] + PowerMod[ 3, #, # ], # ] == 0 & ]
With[{c=Mod[Total[Table[PowerMod[k, #1, #1], {k, 3, 7}]], #1]==0&}, Select[Range[10^6], c]] (* Harvey P. Dale, Mar 07 2026 *)
PROG
(Python)
def ok(k): return k and sum(pow(m, k, k) for m in range(3, 8))%k == 0
print([k for k in range(10**6) if ok(k)]) # Michael S. Branicky, Mar 08 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Robert G. Wilson v, Sep 21 2000
EXTENSIONS
a(29) onward from Michael S. Branicky, Mar 08 2026
STATUS
approved
