login
A125733
a(n) gives k numbers from A125732.
3
1, 2, 4, 8, 5, 16, 18, 32, 36, 11, 64, 72, 100, 128, 68, 144, 38, 162, 23, 196, 200, 256, 288, 324, 392, 400, 450, 512, 576, 648, 369, 374, 47, 784, 800, 882, 221, 900, 968, 1024, 1152, 1296, 1352, 1458, 1568, 1600, 1764, 1800, 1936, 2048, 536, 2178, 284, 2304
OFFSET
1,2
LINKS
FORMULA
a(n) = A125732(n)/A103977(A125732(n)). - Amiram Eldar, Dec 11 2019
MATHEMATICA
midi[n_] := Module[{d = Divisors[n], c, p, m}, c = CoefficientList[Product[1 + x^i, {i, d}], x]; p = -1 + Position[c, _?(# > 0 &)] // Flatten; m = Length[p]; If[OddQ[m], If[(d = p[[(m + 1)/2]] - p[[(m - 1)/2]]) == 1, 0, d], p[[m/2 + 1]] - p[[m/2]]]]; s = {}; Do[m = midi[n]; If[m > 0 && Divisible[n, m], AppendTo[s, n/m]], {n, 1, 2000}]; s (* Amiram Eldar, Dec 11 2019 *)
CROSSREFS
Sequence in context: A248573 A372129 A369414 * A333555 A352378 A280426
KEYWORD
nonn
AUTHOR
Yasutoshi Kohmoto, Jan 01 2007
EXTENSIONS
a(7) corrected and more terms added by Amiram Eldar, Dec 11 2019
STATUS
approved