OFFSET
1,2
COMMENTS
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..52
EXAMPLE
MATHEMATICA
r[n_] := DivisorSigma[0, n]/DivisorSum[n, Boole[BitOr[#, n] == n] &];
seq[kmax_] := Module[{rm = 0, k = 1, s = {}, r1}, Do[r1 = r[k]; If[r1 > rm, rm = r1; AppendTo[s, k]], {k, 1 , kmax}]; s]; seq[10^6]
PROG
(PARI) r(n) = numdiv(n)/sumdiv(n, d, bitor(d, n) == n);
lista(kmax) = {my(rm = 0, r1); for(k = 1, kmax, r1 = r(k); if(r1 > rm, rm = r1; print1(k, ", "))); }
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Amiram Eldar, Mar 31 2023
STATUS
approved