A359083 Numbers k such that A246600(k) = A000005(k) and A000005(k) sets a new record.
1, 3, 15, 63, 255, 891, 4095, 262143, 1048575, 16777215, 68719476735
Offset: 1
Links
Programs
-
Mathematica
s[n_] := DivisorSum[n, 1 &, BitAnd[n, #] == # &]; seq={}; dm = 0; Do[d = DivisorSigma[0, n]; If[d > dm && d == s[n], dm = d; AppendTo[seq, n]], {n, 1, 2*10^7}]; seq
-
PARI
lista(nmax) = {my(list = List(), ndmax = 0, d, s); for(n = 1, nmax, nd = numdiv(n); if(nd > ndmax && sumdiv(n, d, bitand(d, n)==d) == nd, ndmax = nd; listput(list, n))); Vec(list)};
Comments