A376553 Largest unitary square divisor of binomial(n, floor(n/2)).
1, 1, 1, 1, 1, 1, 4, 1, 1, 9, 36, 1, 4, 4, 1, 9, 9, 1, 4, 1, 4, 4, 1, 1, 4, 100, 25, 100, 25, 9, 144, 9, 9, 1, 4, 25, 100, 100, 25, 9, 36, 4, 1, 4, 1, 25, 400, 225, 900, 1764, 441, 196, 49, 49, 784, 4, 1, 1, 16, 1, 16, 16, 1, 441, 441, 49, 196, 49, 196, 36, 9
Offset: 0
Keywords
Links
- Amiram Eldar, Table of n, a(n) for n = 0..1000
Programs
-
Mathematica
f[p_, e_] := If[EvenQ[e], p^e, 1]; a[0] = a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[Binomial[n, Floor[n/2]]]; Array[a, 100, 0]
-
PARI
a(n) = {my(f = factor(binomial(n, n\2))); prod(i = 1, #f~, if(f[i, 2]%2, 1, f[i, 1]^f[i, 2]));}
Comments