A362805 Primitive terms of A362804: terms k of A362804 such that k/2 is not a term of A362804.
1, 6, 28, 30, 45, 496, 8128, 16380, 57720, 65472, 235246, 683520, 33550336, 50426880, 60945408, 105553910, 131297280, 3052879872, 8589869056
Offset: 1
Programs
-
Mathematica
q[n_] := IntegerQ[HarmonicMean[Select[Divisors[n], BitAnd[n, #] == # &]]]; Select[Range[10^6], q[#] && (OddQ[#] || ! q[#/2]) &]
-
PARI
div(n) = select(x->(bitor(x, n) == n), divisors(n)); is1(n) = {my(d = div(n)); denominator(#d/sum(i = 1, #d ,1/d[i])) == 1;} is(n) = is1(n) && (n%2 || !is1(n/2));
Comments