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

A003983 Array read by antidiagonals with T(n,k) = min(n,k).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 2, 3, 2, 1, 1, 2, 3, 3, 2, 1, 1, 2, 3, 4, 3, 2, 1, 1, 2, 3, 4, 4, 3, 2, 1, 1, 2, 3, 4, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 6, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 7, 6, 5, 4, 3, 2, 1, 1, 2, 3, 4, 5, 6, 7, 7, 6, 5, 4, 3, 2, 1
Offset: 1

Views

Author

Keywords

Comments

Also, "correlation triangle" for the constant sequence 1. - Paul Barry, Jan 16 2006
Antidiagonal sums are in A002620.
As a triangle, row sums are A002620. T(2n,n)=n+1. Diagonal sums are A001399. Construction: Take antidiagonal triangle of MM^T where M is the sequence array for the constant sequence 1 (lower triangular matrix with all 1's). - Paul Barry, Jan 16 2006
From Franklin T. Adams-Watters, Sep 25 2011: (Start)
As a triangle, count up to ceiling(n/2) and back down again (repeating the central term when n is even).
When the first two instances of each number are removed from the sequence, the original sequence is recovered.
(End)

Examples

			Triangle version begins
  1;
  1, 1;
  1, 2, 1;
  1, 2, 2, 1;
  1, 2, 3, 2, 1;
  1, 2, 3, 3, 2, 1;
  1, 2, 3, 4, 3, 2, 1;
  1, 2, 3, 4, 4, 3, 2, 1;
  1, 2, 3, 4, 5, 4, 3, 2, 1;
  ...
		

Crossrefs

Programs

  • Haskell
    a003983 n k = a003983_tabl !! (n-1) !! (k-1)
    a003983_tabl = map a003983_row [1..]
    a003983_row n = hs ++ drop m (reverse hs)
       where hs = [1..n' + m]
             (n',m) = divMod n 2
    -- Reinhard Zumkeller, Aug 14 2011
    
  • Maple
    a(n) = min(floor(1/2 + sqrt(2*n)) - (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2+1, (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2) # Leonid Bedratyuk, Dec 13 2009
  • Mathematica
    Flatten[Table[Min[n-k+1, k], {n, 1, 14}, {k, 1, n}]] (* Jean-François Alcover, Feb 23 2012 *)
  • PARI
    T(n,k) = min(n,k) \\ Charles R Greathouse IV, Feb 06 2017
    
  • Python
    from math import isqrt
    def A003983(n):
        a = (m:=isqrt(k:=n<<1))+(k>m*(m+1))
        x = n-(a*(a-1)>>1)
        return min(x,a-x+1) # Chai Wah Wu, Jun 14 2025

Formula

Number triangle T(n, k) = Sum_{j=0..n} [j<=k][j<=n-k]. - Paul Barry, Jan 16 2006
G.f.: 1/((1-x)*(1-x*y)*(1-x^2*y)). - Christian G. Bower, Jan 17 2006
a(n) = min(floor( 1/2 + sqrt(2*n)) - (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2+1, (2*n + round(sqrt(2*n)) - round(sqrt(2*n))^2)/2). - Leonid Bedratyuk, Dec 13 2009

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Nov 08 2000
Entry revised by N. J. A. Sloane, Dec 05 2006

A173987 a(n) = denominator of ((Zeta(0,2,2/3) - Zeta(0,2,n+2/3))/9), where Zeta is the Hurwitz Zeta function.

Original entry on oeis.org

1, 4, 100, 1600, 193600, 9486400, 2741569600, 2741569600, 1450290318400, 245099063809600, 206128312663873600, 3298053002621977600, 3298053002621977600, 1190597133946533913600, 2001393782164123508761600
Offset: 0

Views

Author

Artur Jasinski, Mar 04 2010

Keywords

Crossrefs

Programs

  • Magma
    [1] cat [Denominator((&+[9/(3*k+2)^2: k in [0..n-1]])): n in [1..20]]; // G. C. Greubel, Aug 23 2018
  • Maple
    a := n -> (Zeta(0,2,2/3) - Zeta(0,2,n+2/3))/9:
    seq(denom(a(n)), n=0..14); # Peter Luschny, Nov 14 2017
  • Mathematica
    Table[FunctionExpand[(1/9)*(4*(Pi^2)/3 - Zeta[2, 1/3] - Zeta[2, (3*n + 2)/3])], {n, 0, 20}] // Denominator (* Vaclav Kotesovec, Nov 13 2017 *)
    Denominator[Table[Sum[9/(3*k + 2)^2, {k, 0, n - 1}], {n, 0, 20}]] (* G. C. Greubel, Aug 23 2018 *)
  • PARI
    for(n=0,20, print1(denominator(9*sum(k=0,n-1, 1/(3*k+2)^2)), ", ")) \\ G. C. Greubel, Aug 23 2018
    

Formula

a(n) = denominator of 2*(Pi^2)/3 - J - Zeta(2,(3*n+2)/3), where Zeta is the Hurwitz Zeta function and J is the constant A173973.
a(n) = denominator of Sum_{k=0..(n-1)} 9/(3*k+2)^2. - G. C. Greubel, Aug 23 2018

Extensions

Name simplified by Peter Luschny, Nov 14 2017

A173984 a(n) is the denominator of (Zeta(0,2,1/3) - Zeta(0,2,n+1/3)) where Zeta is the Hurwitz Zeta function.

Original entry on oeis.org

1, 1, 16, 784, 19600, 3312400, 52998400, 19132422400, 2315023110400, 57875577760000, 57875577760000, 55618430227360000, 16073726335707040000, 22004931353582937760000, 22004931353582937760000
Offset: 0

Views

Author

Artur Jasinski, Mar 04 2010

Keywords

Crossrefs

Programs

  • Magma
    [1,1] cat [Denominator((&+[9/(3*k+1)^2: k in [1..n-1]])): n in [2..20]]; // G. C. Greubel, Aug 24 2018
  • Maple
    a := n -> Zeta(0,2,1/3) - Zeta(0,2,n+1/3):
    seq(denom(a(n)), n=0..14); # Peter Luschny, Nov 14 2017
  • Mathematica
    Table[FunctionExpand[-Zeta[2, (3*n + 1)/3] + Zeta[2, 1/3]], {n, 0, 20}] // Denominator (* Vaclav Kotesovec, Nov 13 2017 *)
    Denominator[Table[Sum[9/(3*k + 1)^2, {k, 1, n - 1}], {n, 0, 30}]] (* G. C. Greubel, Aug 24 2018 *)
  • PARI
    for(n=0,20, print1(denominator(sum(k=1,n-1, 9/(3*k+1)^2)), ", ")) \\ G. C. Greubel, Aug 24 2018
    

Formula

a(n) = denominator of 2*(Pi^2)/3 + J - Zeta(2,(3*n+1)/3), where Zeta is the Hurwitz Zeta function and the constant J is A173973.
a(n) = denominator of Sum_{k=1..(n-1)} 9/(3*k+1)^2. - G. C. Greubel, Aug 24 2018

Extensions

Name simplified by Peter Luschny, Nov 14 2017

A173985 a(n) = numerator of (Zeta(0,2,2/3) - Zeta(0,2,n+2/3)), where Zeta is the Hurwitz Zeta function.

Original entry on oeis.org

0, 9, 261, 4401, 546921, 27234729, 7956214281, 8017899597, 4266143013213, 724241322449397, 611292843754229277, 9809672294036025657, 9833902887524676921, 3557459561652307818081, 5990804897343048247416561
Offset: 0

Views

Author

Artur Jasinski, Mar 04 2010

Keywords

Comments

All numbers in this sequence are divisible by 9.

Crossrefs

Programs

  • Magma
    [0] cat [Numerator((&+[9/(3*k+1)^2: k in [0..n-1]])): n in [1..20]]; // G. C. Greubel, Aug 23 2018
  • Maple
    A173985 := proc(n) add( 1/(2/3+i)^2,i=0..n-1) ; numer(%) ; end proc: seq(A173985(n),n=0..20) ; # R. J. Mathar, Apr 22 2010
  • Mathematica
    Table[FunctionExpand[4*(Pi^2)/3 - Zeta[2, 1/3] - Zeta[2, (3*n + 2)/3]], {n, 0, 20}] // Numerator (* Vaclav Kotesovec, Nov 13 2017 *)
    Numerator[Table[Sum[9/(3*k + 1)^2, {k, 0, n - 1}], {n, 0, 20}]] (* G. C. Greubel, Aug 23 2018 *)
  • PARI
    for(n=0,20, print1(numerator(9*sum(k=0,n-1, 1/(3*k+1)^2)), ", ")) \\ G. C. Greubel, Aug 23 2018
    

Formula

a(n) = numerator of 2*(Pi^2)/3 - J - Zeta(2, (3*n+2)/3), where Zeta is the Hurwitz Zeta function and the constant J is A173973.
a(n)/A173987(n) = sum_{i=0..n-1} 1/(i+2/3)^2 = psi'(2/3)-psi'(2/3+n). - R. J. Mathar, Apr 22 2010
a(n) = numerator of Sum_{k=0..(n-1)} 9/(3*k+1)^2. - G. C. Greubel, Aug 23 2018

Extensions

Name simplified by Peter Luschny, Nov 14 2017

A294967 Decimal expansion of (1/9)*Hurwitz Zeta(2, 2/3) = (1/9)*Psi(1, 2/3), with the Polygamma function Psi.

Original entry on oeis.org

3, 4, 0, 4, 3, 0, 6, 0, 1, 0, 3, 9, 8, 5, 7, 4, 8, 9, 9, 9, 8, 5, 9, 0, 8, 0, 3, 6, 9, 7, 2, 9, 8, 3, 5, 0, 3, 5, 9, 1, 8, 8, 3, 4, 3, 3, 7, 4, 8, 2, 3, 2, 6, 2, 2, 1, 5, 8, 6, 4, 7, 3, 7, 1, 2, 5, 4, 4, 8, 7, 1, 6, 7, 4, 2, 2, 8, 0, 1, 6, 8, 2, 1, 2, 9, 5, 5, 8, 3, 7, 0, 8, 1, 5, 6, 5, 6, 0, 5, 1, 8, 0, 1, 4, 7, 4, 1, 1, 0, 7, 7, 2, 2, 8, 6, 7, 7, 9, 7, 3, 1, 7, 9, 8, 3, 1
Offset: 0

Views

Author

Wolfdieter Lang, Nov 12 2017

Keywords

Comments

This is the value of the series Sum_{n=0..infinity} 1/(3*n+2)^2. For (3*n+2)^2 see A016790.
For the partial sums see A173986(n+1)/A173987(n+1), n >= 0.

Examples

			0.340430601039857489998590803697298350359188343374823262215864737125448716...
		

Crossrefs

Programs

Formula

From Amiram Eldar, Oct 02 2020: (Start)
Equals Integral_{1..oo} log(x)/(x^3-1) dx.
Equals 4*Pi^2/27 - A214550. (End)

Extensions

Data corrected by Amiram Eldar, Oct 02 2020
Showing 1-5 of 5 results.