OFFSET
1,2
COMMENTS
Let d be the vector of divisors of n. The sequence d^(2^k) mod n has some period p. Thus if n divides sigma_(2^k)(n) for one period, then n divides sigma_(2^k)(n) for all k. For these n, the first period ends for k < 158. Hence it is easy to verify divisibility for all k. - T. D. Noe, Apr 11 2006
EXAMPLE
n=84 is here because 84 divides each one of sigma_4(n)=53771172, sigma_8(n)=2488859101224132, sigma_16(n)=6144339637187846520573009496452, etc.
MATHEMATICA
t={}; Do[If[Mod[DivisorSigma[4, n], n]==0, AppendTo[t, n]], {n, 10^8}]; Do[t=Select[t, Mod[DivisorSigma[2^k, # ], # ]==0&], {k, 3, 20}]; t (* T. D. Noe, Apr 11 2006 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Labos Elemer, Dec 12 2001
EXTENSIONS
Edited by T. D. Noe, Apr 11 2006
STATUS
approved