OFFSET
1,2
COMMENTS
Indices of records in A073802.
This sequence is infinite since A073802 is unbounded. For example, for any odd number m we have A073802(2^(m-1)*(2^m-1)) >= A000005(m) and the number of divisors of odd numbers is unbounded.
The corresponding record values are 1, 4, 6, 16, 24, 40, 60, 96, 144, 216, 240, 336, ... .
a(26) <= 799959888000.
MATHEMATICA
seq[kmax_] := Module[{d, dmax = 0, s = {}}, Do[d = DivisorSigma[0, GCD[k, DivisorSigma[1, k]]]; If[d > dmax, dmax = d; AppendTo[s, k]], {k, 1, kmax}]; s]; seq[10^6]
PROG
(PARI) lista(kmax) = {my(d, dmax = 0); for(k = 1, kmax, d = numdiv(gcd(k, sigma(k))); if(d > dmax, dmax = d; print1(k, ", "))); }
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Nov 21 2024
STATUS
approved