A346182 a(n) = A346463(n) / radical(A346463(n)).
1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 7, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 9, 1, 1, 1, 1, 7, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 5, 1, 1, 7, 9, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 17, 1, 1, 1, 3, 1, 49, 1, 1, 1, 31, 1, 1, 27, 1, 1, 1, 1
Offset: 1
Keywords
Programs
-
Maple
with(NumberTheory): a := proc(n) local v, d; v := (4^n - 2)*(4^n - 1); for d in Divisors(2*n) do if isprime(d+1) then v := v/(d+1) fi od; v / Radical(v) end: seq(a(n), n = 1..85);
-
Mathematica
radical[n_] := Product[p, {p, FactorInteger[n][[All, 1]]}]; A346463[n_] := 6 QBinomial[2 n, 2, 2] / Denominator[BernoulliB[2 n]]; A346182[n_] := With[{r = A346463[n]}, r / radical[r]]; Table[A346182[n], {n, 1, 40}]