A145738 a(n) = squarefree part of A145609(n).
3, 1, 1, 761, 61, 509, 1171733, 8431, 39541, 55835135, 36093, 1347822955, 34395742267, 375035183, 9682292227, 586061125622639, 54062195834749, 40030624861, 2053580969474233, 1236275063173, 6657281227331
Offset: 1
Keywords
Programs
-
Mathematica
m = 1; aa = {}; Do[k = 0; Do[k = k + m^(2 r + 1 - d)/d, {d, 1, 2 r}]; b = Sqrt[Numerator[k]] /. Sqrt[] -> 1; c = Sqrt[Numerator[k]]/b /. Sqrt[x] -> x; AppendTo[aa, c], {r, 1, 25}]; aa