A119848 Positions where A119842 is not zero.

1, 2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 18, 19, 20, 23, 25, 27, 28, 29, 31, 32, 36, 37, 41, 43, 44, 45, 47, 48, 49, 50, 52, 53, 59, 61, 63, 64, 67, 68, 71, 72, 73, 75, 76, 79, 80, 81, 83, 89, 92, 97, 98, 99, 100, 101, 103, 107, 108, 109, 112, 113, 116, 117, 120, 121

Offset: 1

Author

Antti Karttunen, Jun 04 2006

Keywords

• nonn

Links

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

Crossrefs

Complement: A119847. Cf. A119850.

Programs

Mathematica

b[s_, t_] := b[s, t] = If[Length[s] < 1, 1, Sum[If[Mod[PrimeOmega[x], 2] == 1 - t && Length[Select[s, Mod[#, x] == 0 &]] == 1, b[s~Complement~{x}, 1 - t], 0], {x, s}]]; c[n_] := Module[{l, m}, l = Sort[FactorInteger[n], #1[[2]] > #2[[2]] &]; m = Product[Prime[i]^l[[i]][[2]], {i, 1, Length[l]} ]; b[Divisors[m][[2 ;; -2]], Mod[PrimeOmega[m], 2]]]; A119848 = Position[ Table[c[n], {n, 1, 200}], ?Positive] // Flatten (* _Jean-François Alcover, Mar 03 2016, after Alois P. Heinz *)