|
|
A119843
|
|
A119842(n) sets a new record.
|
|
2
|
|
|
1, 36, 72, 120, 144, 180, 420, 720, 1080, 1680, 2520, 6300, 9240
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
a(9) is greater than 960.
|
|
LINKS
|
|
|
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}]]; a[n_] := a[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]]]; Reap[For[record = 0; k = 1, k < 10^4, k++, If[a[k] > record, record = a[k]; Print[k, " ", a[k]]; Sow[k]]]][[2, 1]] (* Jean-François Alcover, Mar 03 2016, after Alois P. Heinz *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|