cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-10 of 10 results.

A056059 GCD of largest square and squarefree part of central binomial coefficients.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 3, 6, 2, 1, 1, 1, 3, 3, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 6, 3, 1, 1, 1, 2, 3, 6, 2, 1, 1, 2, 2, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 2, 1, 1, 1, 1, 2, 3, 6, 6, 3, 1, 2, 2, 1, 2, 1, 3, 6, 1, 1, 1, 2, 1, 2
Offset: 1

Views

Author

Labos Elemer, Jul 26 2000

Keywords

Examples

			n=14, binomial(14,7) = 3432 = 8*3*11*13. The largest square divisor is 4, squarefree part is 858. So a(14) = gcd(4,858) = 2.

Links

Crossrefs

Cf. A000188, A001405, A007913, A008833, A034974, A046098, A055229, A055231.
Cf. A056056, A056057, A056058, A056060, A056061.
A056201 is the cube of this sequence.

Programs

  • Mathematica
    Table[GCD[First@ Select[Reverse@ Divisors@ #, IntegerQ@ Sqrt@ # &], Times @@ Power @@@ Map[{#1, Mod[#2, 2]} & @@ # &, FactorInteger@ #]] &@ Binomial[n, Floor[n/2]], {n, 80}] (* Michael De Vlieger, Feb 18 2017, after Zak Seidov at A007913 *)
  • PARI
    A001405(n) = binomial(n, n\2);
A055229(n) = { my(c=core(n)); gcd(c, n/c); } \\ Charles R Greathouse IV, Nov 20 2012
A056059(n) = A055229(A001405(n)); \\ Antti Karttunen, Jul 20 2017
  • Python
    from sympy import binomial, gcd
from sympy.ntheory.factor_ import core
def a001405(n): return binomial(n, n//2)
def a055229(n):
    c = core(n)
    return gcd(c, n//c)
def a(n): return a055229(a001405(n))
print([a(n) for n in range(1, 151)]) # Indranil Ghosh, Jul 20 2017

Formula

a(n) = A055229(A001405(n)), where A055229(n) = gcd(A008833(n), A007913(n)).

Extensions

Formula clarified by Antti Karttunen, Jul 20 2017

A056057 The largest square which divides n-th central binomial coefficient.

Original entry on oeis.org

1, 1, 1, 1, 1, 4, 1, 1, 9, 36, 1, 4, 4, 4, 9, 9, 1, 4, 1, 4, 4, 4, 1, 4, 100, 100, 900, 900, 36, 144, 9, 9, 9, 36, 25, 100, 100, 100, 9, 36, 4, 4, 4, 4, 900, 3600, 225, 900, 1764, 1764, 1764, 1764, 196, 784, 4, 4, 4, 16, 4, 16, 16, 16, 441, 441, 49, 196, 49, 196, 36, 36, 1, 4
Offset: 1

Views

Author

Labos Elemer, Jul 26 2000

Keywords

Links

Crossrefs

Cf. A001405, A000188, A008833, A007913, A055229, A055231, A046098, A034974.
Cf. A056056, A056058, A056059, A056060, A056061.

Programs

  • Mathematica
    Table[First@ Select[Reverse@ Divisors@ Binomial[n, Floor[n/2]], IntegerQ@ Sqrt@ # &], {n, 72}] (* Michael De Vlieger, Feb 18 2017 *)
a[n_] := Times @@ (First[#]^(2*Floor[Last[#]/2]) & /@ FactorInteger[Binomial[n, Floor[n/2]]]); Array[a, 100] (* Amiram Eldar, Sep 06 2020 *)

Formula

a(n) = A008833(A001405(n)).
a(A046098(n)) = 1.

A056061 Number of square divisors of central binomial coefficients.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 2, 4, 1, 2, 2, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 4, 4, 8, 8, 4, 6, 2, 2, 2, 4, 2, 4, 4, 4, 2, 4, 2, 2, 2, 2, 8, 12, 4, 8, 8, 8, 8, 8, 4, 6, 2, 2, 2, 3, 2, 3, 3, 3, 4, 4, 2, 4, 2, 4, 4, 4, 1, 2, 2, 2, 4, 4, 8, 12, 2, 4, 12, 12, 4, 4, 8, 12, 12, 12, 4, 6, 8, 12, 12, 12, 8, 16, 8, 8, 6
Offset: 1

Views

Author

Labos Elemer Jul 26 2000

Keywords

Examples

			n=27: binomial(27,13) = 20058300, its largest square-divisor is 900=30^2 so a(27) = tau(30) = 8.

Links

Crossrefs

Cf. A001405, A046951, A000005, A000188, A056056, A056057, A056058, A056059, A056060.

Programs

  • Mathematica
    Table[Count[Divisors@ Binomial[n, Floor[n/2]], d_ /; IntegerQ@ Sqrt@ d], {n, 0, 84}] (* Michael De Vlieger, Feb 18 2017 *)
  • PARI
    a(n) = sumdiv(binomial(n, n\2), d, issquare(d)); \\ Michel Marcus, Feb 19 2017

Formula

a(n) = A046951(A001405(n)) = A000005(A000188(A001405(n))).

A056060 The powerfree part of the central binomial coefficients.

Original entry on oeis.org

1, 2, 3, 6, 10, 5, 35, 70, 14, 7, 462, 231, 429, 429, 715, 1430, 24310, 12155, 92378, 46189, 88179, 88179, 1352078, 676039, 52003, 52003, 7429, 7429, 1077205, 1077205, 33393355, 66786710, 43214930, 21607465, 181502706, 90751353, 176726319, 176726319, 7658140490
Offset: 1

Views

Author

Labos Elemer, Jul 26 2000

Keywords

Examples

			n=14, binomial(14,7) = 3432 = 8*3*11*13. The largest square divisor is 4, and the squarefree part is 858. So GCD(4,858) = 2 and a(14) = 858/2 = 429.

Links

Crossrefs

Cf. A001405, A000188, A008833, A007913, A055229, A055231, A046098, A034974, A056056, A056057, A056058, A056059, A056060, A056061.

Programs

  • Mathematica
    a[n_] := Denominator[(b = Binomial[n, Floor[n/2]])/(Times @@ First /@ FactorInteger[b])^2]; Array[a, 36] (* Amiram Eldar, Sep 05 2020 *)

Formula

a(n) = A055231(A001405(n)).

Extensions

New name and more terms from Amiram Eldar, Sep 05 2020

A056201 Characteristic cube divisor (A056191) of central binomial coefficient (A001405).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 8, 27, 216, 8, 1, 1, 1, 27, 27, 1, 1, 1, 8, 1, 1, 1, 8, 1, 8, 216, 27, 1, 1, 1, 8, 27, 216, 8, 1, 1, 8, 8, 1, 8, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 8, 1, 8, 8, 1, 1, 1, 1, 8, 27, 216, 216, 27, 1, 8, 8, 1, 8, 1, 27, 216
Offset: 1

Views

Author

Labos Elemer, Aug 02 2000

Keywords

Examples

			n=14, binomial(14,8) = 3432 = 2*2*2*3*11*13 so g(3432)=2, thus a(14)=8.

Crossrefs

Cf. A056056-A056060, A056202.
Equals A056059^3.

Formula

a(n) = A056059(n)^3 = g^3 and binomial(n, floor(n/2)) = a(n)^3 * L^2 * A056060(n), where L = A056056(n)/A056059(n).

A056647 a(n) = A056623(A001405(n)).

Original entry on oeis.org

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, 4, 441, 441, 49, 196, 49, 196, 36, 9, 1, 4, 4, 1, 100, 25, 1225
Offset: 1

Views

Author

Labos Elemer, Aug 09 2000

Keywords

Comments

Previous name "Largest unitary square divisor of central binomial coefficient" was incorrect. See A376553 for the correct sequence with this name. - Amiram Eldar, Sep 28 2024

Examples

			a(28) = A056623(binomial(28,14)) = A056623(40116600) = 25.

Crossrefs

Cf. A001405, A008833, A055229, A056056, A056057, A056059, A056623, A376553.

Formula

a(n) = A008833(A001405(n))/A055229(A001405(n))^2 = A056057(n)/A056059(n)^2.

Extensions

Incorrect name replaced with a formula by Amiram Eldar, Sep 28 2024

A056202 Central binomial coefficient A001405(n) divided by its characteristic cube divisor A056201(n).

Original entry on oeis.org

1, 2, 3, 6, 10, 20, 35, 70, 126, 252, 462, 924, 1716, 429, 6435, 12870, 24310, 48620, 92378, 184756, 352716, 88179, 1352078, 2704156, 5200300, 1300075, 742900, 185725, 9694845, 155117520, 300540195, 601080390, 43214930, 86429860
Offset: 1

Views

Author

Labos Elemer, Aug 02 2000

Keywords

Examples

			n=14, binomial(14,7) = 3432 and A056059(14) = 2, thus a(14) = 3432/(2*2*2) = 429.

Crossrefs

Cf. A001405, A056056-A056060, A056201.

Formula

a(n) = A001405(n)/A056059(n)^3 = binomial(n, floor(n/2))/A056059(n)^3 = A001405(n)/A056201(n).

A056646 a(n) = A056622(A001405(n)).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 3, 6, 1, 2, 2, 1, 3, 3, 1, 2, 1, 2, 2, 1, 1, 2, 10, 5, 10, 5, 3, 12, 3, 3, 1, 2, 5, 10, 10, 5, 3, 6, 2, 1, 2, 1, 5, 20, 15, 30, 42, 21, 14, 7, 7, 28, 2, 1, 1, 4, 1, 4, 4, 2, 21, 21, 7, 14, 7, 14, 6, 3, 1, 2, 2, 1, 10, 5, 35, 140, 7, 14, 126, 63, 2, 1, 5, 20, 90, 45, 3, 12
Offset: 1

Views

Author

Labos Elemer, Aug 09 2000

Keywords

Comments

Previous name "Square root of largest unitary square divisor of central binomial coefficient" was incorrect. See A376554 for the correct sequence with this name. - Amiram Eldar, Sep 28 2024

Examples

			a(28) = A056622(binomial(28,14)) = A056622(40116600) = 5.

Crossrefs

Cf. A000188, A001405, A055229, A056056, A056059, A056622, A376554.

Formula

a(n) = A000188(A001405(n))/A055229(A001405(n)) = A056056(n)/A056059(n).

Extensions

Incorrect name replaced with a formula by Amiram Eldar, Sep 28 2024

A056648 a(n) = A034444(A056647(n)).

Original entry on oeis.org

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

Views

Author

Labos Elemer, Aug 09 2000

Keywords

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

Examples

			a(28) = A034444(A056647(28)) = A034444(25) = 2.

Crossrefs

Cf. A000188, A001221, A001405, A008833, A034444, A055229, A056056, A056057, A056059, A056061, A376555.

Programs

  • 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 *)

Formula

a(n) = 2^r, where r = A001221(A000188(A001405(n))/A055229(A001405(n))).

Extensions

Incorrect name replaced with a formula by Amiram Eldar, Sep 28 2024

A056649 a(n) = A056061(n) - A034444(A056647(n)).

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 4, 6, 2, 2, 0, 0, 1, 2, 0, 0, 0, 2, 0, 0, 0, 1, 0, 1, 6, 8, 0, 0, 0, 4, 4, 6, 2, 2, 0, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 1, 0, 2, 4, 4, 0, 0, 4, 8, 2, 3, 6, 8, 4, 8, 2, 2, 4, 4, 8, 8, 0, 0, 0, 4, 2, 4, 3, 4, 2, 3, 4
Offset: 1

Views

Author

Labos Elemer, Aug 09 2000

Keywords

Comments

Previous name, "Number of non-unitary square divisors of central binomial coefficient", was incorrect. See A376556 for the correct sequence with this name. - Amiram Eldar, Sep 28 2024

Examples

			a(28) = A056061(28) - A034444(A056647(28)) = A056061(28) - A034444(25) = 8 - 2 = 6.

Crossrefs

Cf. A000188, A001221, A001405, A008833, A034444, A055229, A056056, A056057, A056059, A056061, A376556.

Programs

  • Mathematica
    A056061[n_] := Count[Divisors@Binomial[n, Floor[n/2]], d_ /; IntegerQ@Sqrt@d]; 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[A056061[n] - 2^(PrimeNu[Sqrt[A008833[Binomial[n, Floor[n/2]]]]/ A055229[Binomial[n, Floor[n/2]]]]), {n, 1, 15}] (* G. C. Greubel, May 20 2017 *)

Formula

a(n) = A056061(n) - 2^r, where r = A001221(A000188(A001405(n))/A055229(A001405(n))).

Extensions

Incorrect name replaced with a formula by Amiram Eldar, Sep 28 2024
Showing 1-10 of 10 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.