A099434 a(n) = A066421(A099433(n)). The value of A066421 on its record-holders.
1, 2, 5, 9, 12, 13, 14, 30, 37, 291, 292
Offset: 1
Extensions
Offset corrected and a(10)-a(11) added by Amiram Eldar, Apr 11 2025
This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
sigma(sigma(sigma(8))) + 1 = sigma(sigma(15)) + 1 = sigma(24) + 1 = 60 + 1 = 61, a prime; hence a(8) = 3.
A066421[n_]:=Length[NestWhileList[DivisorSigma[1,#]&,DivisorSigma[1,n],!PrimeQ[#+1]&]];Array[A066421,100] (* Paolo Xausa, Oct 16 2023 *)
A066421(n) = { my(k=1,s=sigma(n)); while(!isprime(1+s),k++;s = sigma(s)); k; }; \\ Antti Karttunen, Nov 07 2017
Comments