OFFSET
1,2
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..1688 (terms up to 10^10)
Jose María Grau and Antonio M. Oller-Marcen, Power sums over commutative and unitary rings, arXiv:1603.05787 [math.NT], 2016.
EXAMPLE
8000 is a term since 8000 = 2^6 * 5^3 and 8000 == -2^6 (mod 2^7) and 8000 == -5^3 (mod 5^4).
MATHEMATICA
fa = FactorInteger; mas[1]=True; mas[n_] := Union@Table[Mod[n + fa[n][[i, 1]]^ fa[n][[i, 2]], fa[n][[i, 1]]^(fa[n][[i, 2]] + 1)], {i, Length[fa[n]]}] == {0}; Select[Range[10000], mas ]
PROG
(PARI) is(k) = {my(f = factor(k)); for(i = 1, #f~, if((k / f[i, 1]^f[i, 2] + 1) % f[i, 1], return(0))); 1; } \\ Amiram Eldar, Jul 24 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
José María Grau Ribas, Mar 12 2016
STATUS
approved