The OEIS is supported by the many generous donors to the OEIS Foundation.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A328919 Smallest index m such that from the m-th term on, the sequence {sigma_k(n) mod n: k >= 0} enters into a cycle. 2
 0, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 0, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 3, 1, 1, 1, 1, 5, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 3, 2, 2, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 2, 6, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 2, 1, 1, 1, 1, 3, 4, 1, 1, 1, 1, 1, 1 (list; graph; refs; listen; history; text; internal format)
 OFFSET 1,4 COMMENTS sigma_k(n) = Sum_{d|n} d^k. It is easy to see that {sigma_k(n) mod n: k >= A051903(n)} is purely periodic, and the period divides psi(n) = A002322(n) (the Carmichael lambda). So we have a(n) <= A051903(n). The equality seems to hold for most n: it holds for 7125 n's in the range [1, 10000] and 70287 n's in the range [1, 100000]. The numbers n such that a(n) < A051903(n) are listed in A328930. Conjecture: a(n) = 0 only for n = 1, 12. If {sigma_k(n) mod n: k >= 0} is purely periodic, then sigma_0(n) == sigma_psi(n)(n) (mod n). Let p be any prime divisor of n, write n = p^e*s, then sigma_0(n) = d(n) = (e+1)*d(s), sigma_psi(n)(n) == Sum_{d|n, gcd(d,p)=1} 1 = d(s) (mod p^e), d = A000005. So sigma_0(n) == sigma_psi(n)(n) (mod n) means that p^e divides e*d(s) for every p dividing n, which seems highly impossible for n > 12. Of course, this shows that a(n) cannot be 0 if n is squarefree, so a(n) = 1 for squarefree n. Note that sigma_m(n) == sigma_(m+psi(n))(n) (mod n) does NOT mean the sequence {sigma_k(n) mod n: k >= 0} enters into a cycle from the m-th term on: for n = 112, psi(112) = 12, sigma_1(112) == sigma_13(n) (mod 112) == 24 (mod 112), but sigma_2(112) == 26 (mod 112) while sigma_14(112) == 82 (mod 112). Actually, a(112) = 3. Every number occurs in this sequence, because a(p^e) = e for primes p. Conjecture: a(2^e) is always the earliest occurrence of e. LINKS Table of n, a(n) for n=1..87. EXAMPLE {sigma_k(8) mod 8: k >= 0}: 4, 7, 5, (1). {sigma_k(8) mod 8: k >= 0} enters into the cycle (1) from the 3rd term on, so a(8) = 3. {sigma_k(12) mod 12: k >= 0}: (6, 4). {sigma_k(12) mod 12: k >= 0} enters into the cycle (6, 4) from the very beginning, so a(12) = 0. {sigma_k(24) mod 24: k >= 0}: 8, (12, 10). {sigma_k(24) mod 24: k >= 0} enters into the cycle (12, 10) from the 1st term on, so a(24) = 1. {sigma_k(90) mod 90: k >= 0}: 12, (54, 40, 18, 76). {sigma_k(90) mod 90: k >= 0} enters into the cycle (54, 40, 18, 76) from the 1st term on, so a(90) = 1. PROG (PARI) sigmamod(k, n) = my(d=divisors(n)); lift((sum(i=1, #d, Mod(d[i], n)^k))%n) a(n) = forstep(k=A051903(n)-1, 0, -1, if(sigmamod(k, n)!=sigmamod(k+A002322(n), n), return(k+1))); return(0) \\ See A002322 and A051903 for their programs CROSSREFS Cf. A051903, A328930, A002322, A000005. Cf. A000203 (sigma), A013960 (sigma_12). Sequence in context: A037861 A145037 A267115 * A277647 A296134 A306694 Adjacent sequences: A328916 A328917 A328918 * A328920 A328921 A328922 KEYWORD nonn AUTHOR Jianing Song, Oct 31 2019 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified June 4 10:51 EDT 2023. Contains 363121 sequences. (Running on oeis4.)