login
A346182
a(n) = A346463(n) / radical(A346463(n)).
1
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,9
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}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Jul 20 2021
STATUS
approved