A069782 Numbers k such that gcd(d(k^3), d(k)) = 2^w for some w.
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107
Offset: 1
Examples
Below 100000 only 314 integers are missing, collected in A069781.
Links
- Andrew Howroyd, Table of n, a(n) for n = 1..10000
- Tanya Khovanova, Non Recursions
Programs
-
Mathematica
f[x_] := GCD[DivisorSigma[0, x^3], DivisorSigma[0, x]]; Do[s=f[n]; If[IntegerQ[Log[2, s]], Print[{n, s}]], {n, 1, 100000}] -
PARI
is(n)=my(f=factor(n)[, 2], g=gcd(prod(i=1, #f, 3*f[i]+1), prod(i=1, #f, f[i]+1))); g>>valuation(g, 2)==1 \\ Charles R Greathouse IV, Oct 16 2015
Comments