A324060 Smallest positive index m such that A000009(m) is divisible by n.

1, 3, 5, 6, 7, 8, 33, 9, 13, 10, 14, 11, 23, 33, 12, 16, 31, 13, 17, 31, 40, 14, 18, 27, 65, 23, 15, 33, 52, 32, 42, 16, 26, 31, 41, 38, 28, 17, 32, 37, 112, 41, 149, 79, 35, 18, 121, 27, 203, 65, 49, 23, 40, 19, 26, 33, 87, 52, 50, 41, 24, 42, 41, 20, 32, 79

Offset: 1

Author

Robert G. Wilson v, Jul 31 2019

Keywords

• nonn

Comments

Analogous to A046641.

Links

• Alois P. Heinz, Table of n, a(n) for n = 1..20000

Crossrefs

Cf. A000009, A046641.

Programs

Maple

b:= proc(n) option remember; `if`(n=0, 1, add(b(n-j)*add(
     `if`(d::odd, d, 0), d=numtheory[divisors](j)), j=1..n)/n)
    end:
a:= proc(n) local k; for k while irem(b(k), n)>0 do od; k end:
seq(a(n), n=1..100);  # Alois P. Heinz, Aug 17 2019

Mathematica

t[] := 0; k = 1; f[n] := Block[{p = PartitionsQ@n, lst}, lst = Select[Range@100, Mod[p, #] == 0 &]; If[t[#] == 0, t[#] = n] & /@ lst]; While[k < 1001, f[k]; k++]; t@# & /@ Range@100