OFFSET
1,1
COMMENTS
The corresponding record values are 2, 4, 6, 8, 9, 10, 12, 15, 16, 18, 20, 24, 32, 40, 48, ... .
LINKS
David A. Corneth, Table of n, a(n) for n = 1..160 (first 56 terms from Amiram Eldar)
FORMULA
a(n) <= 2*a(n-1) for n >= 2. - David A. Corneth, Jun 18 2023
MATHEMATICA
seq[kmax_] := Module[{s = {}, dm = 0, d1}, Do[d1 = DivisorSigma[0, k]; If[d1 > dm && DivisorSum[k, Boole[BitOr[#, k] == k] &] == 2, dm = d1; AppendTo[s, k]], {k, 1, kmax}]; s]; seq[10^5]
PROG
(PARI) lista(kmax) = {my(dm = 0, d1); for(k = 1, kmax, d1 = numdiv(k); if(d1 > dm && sumdiv(k, d, bitor(d, k) == k) == 2, dm = d1; print1(k, ", "))); }
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Amiram Eldar, Jun 16 2023
STATUS
approved