login
A056648
a(n) = A034444(A056647(n)).
0
1, 1, 1, 1, 1, 2, 1, 1, 2, 4, 1, 2, 2, 1, 2, 2, 1, 2, 1, 2, 2, 1, 1, 2, 4, 2, 4, 2, 2, 4, 2, 2, 1, 2, 2, 4, 4, 2, 2, 4, 2, 1, 2, 1, 2, 4, 4, 8, 8, 4, 4, 2, 2, 4, 2, 1, 1, 2, 1, 2, 2, 2, 4, 4, 2, 4, 2, 4, 4, 2, 1, 2, 2, 1, 4, 2, 4, 8, 2, 4, 8, 4, 2, 1, 2, 4, 8, 4, 2, 4, 4, 8, 4, 4, 8, 16, 8, 4, 4, 2, 1, 2, 2, 1
OFFSET
1,6
COMMENTS
Previous name, "Number of unitary square divisors of central binomial coefficient", was incorrect. See A376555 for the correct sequence with this name. - Amiram Eldar, Sep 28 2024
FORMULA
a(n) = 2^r, where r = A001221(A000188(A001405(n))/A055229(A001405(n))).
EXAMPLE
a(28) = A034444(A056647(28)) = A034444(25) = 2.
MATHEMATICA
A008833[n_] := First[Select[Reverse[Divisors[n]], IntegerQ[Sqrt[#]] &, 1]]; A055229[n_] := With[{sf = Times @@ Power @@@ ({#[[1]], Mod[#[[2]], 2]} & /@ FactorInteger[n])}, GCD[sf, n/sf]]; Table[2^(PrimeNu[ Sqrt[A008833[Binomial[n, Floor[n/2]]]]/A055229[Binomial[n, Floor[n/2]]]]), {n, 1, 25}] (* G. C. Greubel, May 20 2017 *)
KEYWORD
nonn
AUTHOR
Labos Elemer, Aug 09 2000
EXTENSIONS
Incorrect name replaced with a formula by Amiram Eldar, Sep 28 2024
STATUS
approved