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-2 of 2 results.

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

Formula

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

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.
		

Crossrefs

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
Showing 1-2 of 2 results.