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.

Previous Showing 41-50 of 121 results. Next

A085582 The number of rectangles (orthogonal or not) with corners on an n X n grid of points.

Original entry on oeis.org

0, 1, 10, 44, 130, 313, 640, 1192, 2044, 3305, 5078, 7524, 10750, 14993, 20388, 27128, 35448, 45665, 57922, 72636, 89970, 110297, 133976, 161440, 192860, 228857, 269758, 316012, 367974, 426417, 491468, 564120, 644640, 733633, 831674, 939292
Offset: 1

Views

Author

Yuval Dekel (dekelyuval(AT)hotmail.com), Jul 06 2003

Keywords

Examples

			a(3) = 10 because on the 3 X 3 grid there are four 1 X 1 rectangles, two 1 X 2s, two 2 X 1's, one 2 X 2 and one 45-degree rectangle, sqrt(2) X sqrt(2).
		

Crossrefs

Cf. A000537, A002415, A113751 (diagonal rectangles on an n X n grid).

Formula

a(n) = A000537(n-1) + A113751(n). - T. D. Noe, Nov 09 2005 [corrected by David Radcliffe, Feb 06 2020]
a(n) = n*(n-1)^2*(2n-1)/6 + 2*Sum_{a,b>0, 0David Radcliffe, Feb 06 2020

Extensions

Edited by Don Reble, Nov 05 2005

A156308 Inverse of triangle S(n,m) defined by sequence A156290, n >= 1, 1 <= m <= n.

Original entry on oeis.org

1, 4, 1, 9, 6, 1, 16, 20, 8, 1, 25, 50, 35, 10, 1, 36, 105, 112, 54, 12, 1, 49, 196, 294, 210, 77, 14, 1, 64, 336, 672, 660, 352, 104, 16, 1, 81, 540, 1386, 1782, 1287, 546, 135, 18, 1, 100, 825, 2640, 4290, 4004, 2275, 800, 170, 20, 1
Offset: 1

Views

Author

Hartmut F. W. Hoft, Feb 07 2009

Keywords

Comments

From Wolfdieter Lang, Jun 26 2011: (Start)
This triangle S(n,m) appears as U_m(n) in the Knuth reference on p. 285. It is related to the Riordan triangle T_m(n) = A111125(n,m) by S(n,m) = A111125(n,m) - A111125(n-1,m), n >= m >= 1 (identity on p. 286).
Also, S(n,m)-S(n-1,m) = A111125(n-1,m-1), n >= 2, m >= 1 (identity on p. 286). (End)
These polynomials may be expressed in terms of the Faber polynomials of A263916 and are embedded in A127677 and A208513. - Tom Copeland, Nov 06 2015

Examples

			Triangle starts:
  n=1:  1;
  n=2:  4,  1;
  n=3:  9,  6,  1;
  n=4: 16, 20,  8,  1;
  ...
		

Crossrefs

Same as triangle A208513 with the first column truncated.
Columns: A000290 (m=1), A002415 (m=2), A040977 (m=3), A053347 (m=4), A054334 (m=5).

Programs

  • Magma
    [(n/k)*Binomial(n+k-1, 2*k-1): k in [1..n], n in [1..12]]; // G. C. Greubel, Feb 01 2022
  • Mathematica
    S[m_] := Flatten[Table[k/j Binomial[k + j - 1, 2 j - 1], {k, 1, m}, {j, 1, k}]]
  • Sage
    flatten([[(n/k)*binomial(n+k-1, 2*k-1) for k in (1..n)] for n in (1..12)]) # G. C. Greubel, Feb 01 2022
    

Formula

S(n, m) = (n/m) * binomial(n + m - 1, 2*m - 1).
From Peter Bala, May 01 2014: (Start)
The n-th row o.g.f. is polynomial R(n,x) = 2/x*( T(n,(x + 2)/2) - 1 ), where T(n,x) is Chebyshev polynomial of the first kind. They form a divisibility sequence: if n divides m then R(n,x) divides R(m,x) in the ring Z[x].
R(2*n,x) = (x + 4)*U(n-1,(x + 2)/2)^2;
R(2*n + 1,x) = ( U(n,(x + 2)/2) + U(n-1,(x + 2)/2) )^2.
O.g.f.: Sum_{n >= 0} R(n,x)*z^n = z*(1 + z)/( (1 - z)*(1 - (x + 2)*z + z^2) ). (End)
The polynomial R(n,x) defined above satisfies (x + 1/x - 2) * R(n, x + 1/x - 2) = x^n + 1/x^n - 2. - Alexander Burstein, May 23 2021

Extensions

Edited by Max Alekseyev, Mar 05 2018

A107985 Triangle read by rows: T(n,k) = (k+1)*(n+2)*(n-k+1)/2 for 0 <= k <= n.

Original entry on oeis.org

1, 3, 3, 6, 8, 6, 10, 15, 15, 10, 15, 24, 27, 24, 15, 21, 35, 42, 42, 35, 21, 28, 48, 60, 64, 60, 48, 28, 36, 63, 81, 90, 90, 81, 63, 36, 45, 80, 105, 120, 125, 120, 105, 80, 45, 55, 99, 132, 154, 165, 165, 154, 132, 99, 55, 66, 120, 162, 192, 210, 216, 210, 192, 162, 120, 66
Offset: 0

Views

Author

Emeric Deutsch, Jun 12 2005

Keywords

Comments

Kekulé numbers for certain benzenoids.
T(n,k) is the number of Dyck (n+3)-paths with 3 peaks (UDs) and last descent of length k+1. For example, T(1,1)=3 counts UUDUDUDD, UDUUDUDD, UDUDUUDD. The number of Dyck n-paths containing k peaks and with last descent of length j is (j/n)*binomial(n,k-1)*binomial(n-j-1,k-2) (where as usual binomial(a,b)=0 for b < 0 except that binomial(-1,-1):=1). - David Callan, Jun 26 2006
As a rectangular array, this is the accumulation array (cf. A144112) of the rectangular array W given by w(i,j)=i+j-1; i.e., W=A002024 as a rectangular array. - Clark Kimberling, Sep 16 2008
T(n,k) gives the dimension of an irreducible representation of SU(3) whose Young diagram (n,k) has two rows of length n and k, respectively. - Dimitris Cardaris, May 10 2025

Examples

			Triangle begins:
   1;
   3,  3;
   6,  8,  6;
  10, 15, 15, 10;
  15, 24, 27, 24, 15;
  ...
		

References

  • S. J. Cyvin and I. Gutman, Kekulé structures in benzenoid hydrocarbons, Lecture Notes in Chemistry, No. 46, Springer, New York, 1988 (p. 237, K{B(n,2,-l)}).

Crossrefs

Cf. A000217 (column 0 and main diagonal), A002024, A002415 (row sums), A098737, A144112.

Programs

  • Maple
    T:=proc(n,k) if k<=n then (k+1)*(n+2)*(n-k+1)/2 else 0 fi end: for n from 0 to 11 do seq(T(n,k),k=0..n) od; # yields sequence in triangular form
  • Mathematica
    T[n_,k_]:= (k+1)(n+2)(n-k+1)/2; Table[T[n,k],{n,0,10},{k,0,n}]//Flatten (* Stefano Spezia, Jan 06 2025 *)
  • Python
    from math import isqrt, comb
    def A107985(n):
        a = (m:=isqrt(k:=n+1<<1))+(k>m*(m+1))
        b = n-comb(a,2)
        return (b+1)*(a+1)*(a-b)>>1 # Chai Wah Wu, Jun 14 2025

Formula

T(n,n-k) = T(n,k); T(2n,n) = (n+1)^3.
G.f.: (1 - x^2*y)/((1 - x)^3*(1 - x*y)^3). - Stefano Spezia, Oct 01 2023

A101104 a(1)=1, a(2)=12, a(3)=23, and a(n)=24 for n>=4.

Original entry on oeis.org

1, 12, 23, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24, 24
Offset: 1

Views

Author

Cecilia Rossiter, Dec 15 2004

Keywords

Comments

Original name: The first summation of row 4 of Euler's triangle - a row that will recursively accumulate to the power of 4.

Crossrefs

For other sequences based upon MagicNKZ(n,k,z):
..... | n = 1 | n = 2 | n = 3 | n = 4 | n = 5 | n = 6 | n = 7
---------------------------------------------------------------------------
z = 0 | A000007 | A019590 | .......MagicNKZ(n,k,0) = A008292(n,k+1) .......
z = 1 | A000012 | A040000 | A101101 | thisSeq | A101100 | ....... | .......
z = 2 | A000027 | A005408 | A008458 | A101103 | A101095 | ....... | .......
z = 3 | A000217 | A000290 | A003215 | A005914 | A101096 | ....... | .......
z = 4 | A000292 | A000330 | A000578 | A005917 | A101098 | ....... | .......
z = 5 | A000332 | A002415 | A000537 | A000583 | A022521 | ....... | A255181
Cf. A101095 for an expanded table and more about MagicNKZ.

Programs

  • Mathematica
    MagicNKZ = Sum[(-1)^j*Binomial[n+1-z, j]*(k-j+1)^n, {j, 0, k+1}];Table[MagicNKZ, {n, 4, 4}, {z, 1, 1}, {k, 0, 34}]
    Join[{1, 12, 23},LinearRecurrence[{1},{24},56]] (* Ray Chandler, Sep 23 2015 *)

Formula

a(k) = MagicNKZ(4,k,1) where MagicNKZ(n,k,z) = Sum_{j=0..k+1} (-1)^j*binomial(n+1-z,j)*(k-j+1)^n (cf. A101095). That is, a(k) = Sum_{j=0..k+1} (-1)^j*binomial(4, j)*(k-j+1)^4.
a(1)=1, a(2)=12, a(3)=23, and a(n)=24 for n>=4. - Joerg Arndt, Nov 30 2014
G.f.: x*(1+11*x+11*x^2+x^3)/(1-x). - Colin Barker, Apr 16 2012

Extensions

New name from Joerg Arndt, Nov 30 2014
Original Formula edited and Crossrefs table added by Danny Rorabaugh, Apr 22 2015

A112742 a(n) = n^2*(n^2 - 1)/3.

Original entry on oeis.org

0, 0, 4, 24, 80, 200, 420, 784, 1344, 2160, 3300, 4840, 6864, 9464, 12740, 16800, 21760, 27744, 34884, 43320, 53200, 64680, 77924, 93104, 110400, 130000, 152100, 176904, 204624, 235480, 269700, 307520, 349184, 394944, 445060, 499800, 559440
Offset: 0

Views

Author

Matthew T. Cornick (maruth(AT)gmail.com), Sep 16 2005

Keywords

Comments

Second derivative of the n-th Chebyshev polynomial (of the first kind) evaluated at x=1.
The second derivative at x=-1 is just (-1)^n * a(n).
The difference between two consecutive terms generates the sequence a(n+1) - a(n) = A002492(n).
Consider the partitions of 2n into two parts (p,q) where p <= q. Then a(n) is the total volume of the family of rectangular prisms with dimensions p, |q-p| and |q-p|. - Wesley Ivan Hurt, Apr 15 2018

Examples

			a(4)=80 because
C_4(x) = 1 - 8x^2 + 8x^4,
C'_4(x) = -16x + 32x^3,
C''_4(x) = -16 + 96x^2,
C''_4(1) = -16 + 96 = 80.
		

Crossrefs

Programs

  • Mathematica
    Table[D[ChebyshevT[n, x], {x, 2}], {n, 0, 100}] /. x -> 1
  • PARI
    a(n)=n^2*(n^2-1)/3 \\ Charles R Greathouse IV, Oct 07 2015

Formula

a(n) = (n-1)*n^2*(n+1)/3 = 4*A002415(n).
a(n) = 2*( A000914(n-1) + C(n+1,4) ). - David Scambler, Nov 27 2006
From Colin Barker, Jan 26 2012: (Start)
a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5).
G.f.: 4*x^2*(1+x)/(1-x)^5. (End)
E.g.f.: exp(x)*x^2*(6 + 6*x + x^2)/3. - Stefano Spezia, Dec 11 2021
a(n) = A053126(n+2) - A006324(n-1). - Yasser Arath Chavez Reyes, Feb 22 2024

A108279 a(n) = number of squares with corners on an n X n grid, distinct up to congruence.

Original entry on oeis.org

0, 1, 3, 5, 8, 11, 15, 18, 23, 28, 33, 38, 45, 51, 58, 65, 73, 80, 89, 97, 107, 116, 126, 134, 146, 158, 169, 180, 192, 204, 218, 228, 243, 257, 270, 285, 302, 316, 331, 346, 364, 379, 397, 414, 433, 451, 468, 484, 505, 523, 544, 563, 584, 603, 625
Offset: 1

Views

Author

Hugo Pfoertner, Jun 05 2005

Keywords

Comments

Number of different sizes occurring among the A002415(n) = n^2*(n^2-1)/12 squares that can be drawn using points of an n X n square array as corners.
a(n) is also the number of rectangular isosceles triangles, distinct up to congruence, on an n X n grid (or geoboard). - Martin Renner, May 03 2011

Examples

			a(3)=3 because the 6 different squares that can be drawn on a 3 X 3 square lattice come in 3 sizes:
  4 squares of side length 1:
  x.x.o    o.x.x    o.o.o    o.o.o
  x.x.o    o.x.x    x.x.o    o.x.x
  o.o.o    o.o.o    x.x.o    o.x.x
  1 square of side length sqrt(2):
  o.x.o
  x.o.x
  o.x.o
  1 square of side length 2:
  x.o.x
  o.o.o
  x.o.x
.
a(4)=5 because there are 5 different sizes of squares that can be drawn using the points of a 4 X 4 square lattice:
  x.x.o.o    o.x.o.o    x.o.x.o    o.x.o.o    x.o.o.x
  x.x.o.o    x.o.x.o    o.o.o.o    o.o.o.x    o.o.o.o
  o.o.o.o    o.x.o.o    x.o.x.o    x.o.o.o    o.o.o.o
  o.o.o.o    o.o.o.o    o.o.o.o    o.o.x.o    x.o.o.x
		

Crossrefs

Programs

  • Mathematica
    a[n_] := Module[{v = Table[0, (n - 1)^2]}, Do[v[[k^2 + (w - k)^2]] = 1, {w, 1, n - 1}, {k, 0, w - 1}]; Total[v]]; Array[a, 55](* Jean-François Alcover, Oct 08 2017, after Andrew Howroyd *)
  • PARI
    a(n) = my(v=vector((n-1)^2)); for(w=1, n-1, for(k=0, w-1, v[k^2+(w-k)^2]=1)); vecsum(v); \\ Andrew Howroyd, Sep 17 2017

Extensions

More terms from David W. Wilson, Jun 07 2005

A131198 Triangle T(n,k), 0 <= k <= n, read by rows, given by [1,0,1,0,1,0,1,0,...] DELTA [0,1,0,1,0,1,0,1,...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 3, 1, 0, 1, 6, 6, 1, 0, 1, 10, 20, 10, 1, 0, 1, 15, 50, 50, 15, 1, 0, 1, 21, 105, 175, 105, 21, 1, 0, 1, 28, 196, 490, 490, 196, 28, 1, 0, 1, 36, 336, 1176, 1764, 1176, 336, 36, 1, 0, 1, 45, 540, 2520, 5292, 5292, 2520, 540, 45, 1, 0
Offset: 0

Views

Author

Philippe Deléham, Oct 20 2007

Keywords

Comments

Mirror image of triangle A090181, another version of triangle of Narayana (A001263).
Equals A133336*A130595 as infinite lower triangular matrices. - Philippe Deléham, Oct 23 2007

Examples

			Triangle begins:
  1;
  1,  0;
  1,  1,   0;
  1,  3,   1,   0;
  1,  6,   6,   1,   0;
  1, 10,  20,  10,   1,   0;
  1, 15,  50,  50,  15,   1,  0;
  1, 21, 105, 175, 105,  21,  1, 0;
  1, 28, 196, 490, 490, 196, 28, 1, 0; ...
		

Crossrefs

Programs

  • Magma
    [[n le 0 select 1 else (n-k)*Binomial(n,k)^2/(n*(k+1)): k in [0..n]]: n in [0..10]]; // G. C. Greubel, Feb 06 2018
  • Maple
    T := (n,k) -> `if`(n=0, 0^n, binomial(n,k)^2*(n-k)/(n*(k+1)));
    seq(print(seq(T(n,k), k=0..n)), n=0..5); # Peter Luschny, Jun 08 2014
    R := n -> simplify(hypergeom([1 - n, -n], [2], x)):
    Trow := n -> seq(coeff(R(n, x), x, k), k = 0..n):
    seq(print(Trow(n)), n = 0..9); # Peter Luschny, Apr 26 2022
  • Mathematica
    Table[If[n == 0, 1, (n-k)*Binomial[n,k]^2/(n*(k+1))], {n,0,10}, {k,0,n}] //Flatten (* G. C. Greubel, Feb 06 2018 *)
  • PARI
    for(n=0,10, for(k=0,n, print1(if(n==0,1, (n-k)*binomial(n,k)^2/(n* (k+1))), ", "))) \\ G. C. Greubel, Feb 06 2018
    

Formula

Sum_{k=0..n} T(n,k)*x^k = A000012(n), A000108(n), A001003(n), A007564(n), A059231(n), A078009(n), A078018(n), A081178(n), A082147(n), A082181(n), A082148(n), A082173(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 respectively.
Sum_{k=0..n} T(n,k)*x^(n-k) = A000007(n), A000108(n), A006318(n), A047891(n+1), A082298(n), A082301(n), A082302(n), A082305(n), A082366(n), A082367(n), for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Oct 23 2007
Sum_{k=0..floor(n/2)} T(n-k,k) = A004148(n). - Philippe Deléham, Nov 06 2007
T(2*n,n) = A125558(n). - Philippe Deléham, Nov 16 2011
T(n, k) = [x^k] hypergeom([1 - n, -n], [2], x). - Peter Luschny, Apr 26 2022

A210664 Square array read by upwards antidiagonals: T(m,n) is the number of simple 3-connected triangulations of a closed region in the plane with m+3 given external edges and 3n+m internal edges, m>=0, n>=1.

Original entry on oeis.org

1, 1, 0, 1, 2, 1, 1, 5, 6, 3, 1, 9, 20, 22, 12, 1, 14, 50, 85, 91, 52, 1, 20, 105, 254, 385, 408, 241, 1, 27, 196, 644, 1287, 1836, 1938, 1173, 1, 35, 336, 1448, 3696, 6630, 9120, 9614, 5929, 1, 44, 540, 2967, 9468, 20790, 34846, 46805, 49335, 30880
Offset: 0

Views

Author

N. J. A. Sloane, Mar 28 2012

Keywords

Comments

A triangulation is simple if it contains no separating 3-cycle. There are n interior nodes and m+3 nodes on the boundary. - Andrew Howroyd, Feb 24 2021

Examples

			Array begins:
  1,  0,   1,    3,   12, ... (A000256)
  1,  2,   6,   22,   91, ...
  1,  5,  20,   85,  385, ...
  1,  9,  50,  254, 1287, ...
  1, 14, 105,  644, 3696, ...
  1, 20, 196, 1448, 9468, ...
  ...
From _Andrew Howroyd_, Feb 24 2021: (Start)
The array transposed for comparability with A341856 begins:
==================================================
n\m |   0    1    2     3      4      5      6
----+---------------------------------------------
  1 |   1    1    1     1      1      1      1 ...
  2 |   0    2    5     9     14     20     27 ...
  3 |   1    6   20    50    105    196    336 ...
  4 |   3   22   85   254    644   1448   2967 ...
  5 |  12   91  385  1287   3696   9468  22131 ...
  6 |  52  408 1836  6630  20790  58564 151146 ...
  7 | 241 1938 9120 34846 116641 353056 983664 ...
(End)
		

Crossrefs

Rows m=0..3 are A000256, A000139, A341920, A341921.
Columns are A000012, A000096, A002415, A004305.
Antidiagonal sums give A341922.
Cf. A341856.

Programs

  • PARI
    \\ here H is A000256 as g.f., U(n,m) is A341856 for m > 0.
    H(n)={my(g=1+serreverse(x/(1+x)^4 + O(x*x^n) )); 2 - sqrt(serreverse(x*(2-g)^2*g^4)/x)}
    U(n, m)={(3*(m+2)!*(m-1)!/(3*n+3*m+3)!)*sum(j=0, min(m, n-1), (4*n+3*m-j+1)!*(m+j+2)*(m-3*j)/(j!*(j+1)!*(m-j)!*(m-j+2)!*(n-j-1)!))}
    R(N, m)={my(g=2-H(N)); Vec(if(m==0, 1-g, g^(m+1)*subst(O(x*x^N) + sum(n=1, N, U(n,m)*x^n), x, x*g^2)))}
    M(m, n=m)={Mat(vectorv(m+1, i, R(n,i-1)))}
    M(7) \\ Andrew Howroyd, Feb 23 2021

Formula

From Andrew Howroyd, Feb 24 2021: (Start)
G.f. of row m > 0: R(x) satisfies g(x^2)^(m+1)*R(x*g(x^2)) = B(x^2) where g(x) is the g.f. of column 0 of A341856 and B(x) is the g.f. of column m of A341856.
G.f. of row m > 0: h(x)^(m+1)*B(x*h(x)^2) where 2-h(x) is the g.f. of A000256 and B(x) is the g.f. of column m of A341856.
(End)

Extensions

Terms a(21) and beyond from Andrew Howroyd, Feb 23 2021

A051602 a(n) is the maximal number of squares that can be formed from n points in the plane.

Original entry on oeis.org

0, 0, 0, 0, 1, 1, 2, 3, 4, 6, 7, 8, 11, 13, 15, 17, 20, 22
Offset: 0

Views

Author

Keywords

Comments

Sascha Kurz has proved that one can assume that the points belong to the square grid Z X Z.
So we obtain the same values if we replace the definition by: a(n) is the maximal number of squares that can be formed from n points chosen from the infinite square grid.
In other words, take the infinite square grid. Pick a set S of n grid points, and let c(S) be the number of subsets of four points of S that form a square of any (nonzero) size. Then a(n) = maximum of c(S) over all choices of S.
The general problem of estimating the maximal number of similar figures within point-configurations was studied in [Elekes--Erdos]. References can be found in [Matousek, p. 47] and [AFKK]. Note that the present problem concerning squares shares several properties with the case of right isosceles triangles treated in [AFR].
The following remarks are out of date, and will be revised soon to reflect progress made in October 2021 by a number of members of the Sequence Fans Mailing List.
A more detailed account will be found in the report by Sascha Kurz et al. which is nearing completion.
[Comments revised to this point by N. J. A. Sloane, Nov 02 2021]
Values for n <= 9 [now 16] are exact and are the same for a and b (see proofs below by Hugo van der Sanden and Benoît Jubin for n <= 7 and Sascha Kurz for n <= 9, which furthermore classify all optimal configurations for these values). Values for n > 9 are conjectural since they were obtained by exhaustive search for grid points within a square of side ceil(sqrt(n)), which is a reasonable assumption. A proof that optimal configurations (with gcd of side lengths equal to 1) have a diameter at most f(n) with f of moderate growth would permit exact computation of values from exhaustive searches.
Asymptotic behavior:
One has a(n), b(n) = Theta(n^2):
Upper bound: Since two vertices determine squares, one has b(n) = O(n^2). More explicitly: a pair of points uniquely determines the square that has it as diagonal, and a square has two diagonals, so b(n) <= n(n-1)/4 ~ n^2/4.
Lower bound: When n = m^2, the set of all grid points in [0, m-1]^2 yields S = n(n-1)/12 squares. Indeed, for a in [0, m-1] and b in [1, m], the square formed on (a,0) and (0,b) (as its "lower-left side") has other vertices (a+b,b) and (a,a+b), so there are (m-(a+b))^2 translates of that square. Therefore, S = Sum_{a=0..m-1} Sum_{b=1..m} (m - (a + b))^2. By change of summation indices ((a,c) := (a,a+b)), expanding, using the sum of the first m integers, squares, cubes, and refactoring, one obtains S = n(n-1)/12. Since a is nondecreasing, one has a(n) >= (n-1)(n-2)/12 ~ n^2/12.
We actually have better lower and upper bounds:
0.09... = (1-2/Pi)/4 <= liminf a(n)/n^2 <= limsup b(n)/n^2 <= 1/6 = 0.16...
The upper bound is given in [AFR] (which counts isosceles right triangles, so their value has to be divided by four, the number of isosceles right triangles formed on a square). This gives b(n) <= (2/3*(n - 1)^2 - 5/3)/4.
Lower bound (due to Peter Munn, see SeqFan post in the links): For r >= 0, denote by D(r) the disc centered at the origin with radius r. If A is a point on the boundary of D(r), then the set of points B such that the square with diagonal AB is included in D(r) is a lens-shaped region of area (Pi-2)r^2 (as a proportion of the disc area, this is twice A258146). Therefore, the number S of grid-squares included in D(R) can be estimated as follows: since the set of squares with at least two vertices equidistant from the origin is negligible, we can assume that every square has a unique vertex furthest from the origin, say at distance r, which corresponds to the A above. Then the opposite vertex B is in the region computed above, and the L^1 (aka rectilinear, or Manhattan) distance between A and B is even (being opposite vertices, they are two sides apart), so we have to divide that area by 2. There are approximately 2*Pi*r grid points at a distance approximately r from the origin, so at first order, S ~= Integral_{r=0..R} Pi*r(Pi - 2)r^2 dr = Pi(Pi-2)/4 R^4. Since the disc D(R) contains approximately Pi*R^2 points, one obtains S ~= (1-2/Pi)/4 n^2.
Conjecture: the asymptotic density of the numbers n such that there is no maximal arrangement formed by all the grid points within a suitably chosen circle, is 0. - Peter Munn, Sep 30 2021
For the known values of a(n) (n <= 17), there is a maximal arrangement formed using circles as specified by the conjecture above. For n <= 100 no arrangement has yet been found to contain more squares than the best attained using a circle as specified by A348469. - Peter Munn, Nov 12 2021

Examples

			Lower bounds:
Computer searches, using glutton algorithms starting with all grid points in a convex polygon and adding successive points, have given the following current record configurations, which thus yield lower bounds for a(n).  They are due mainly for n <= 36 to Sean A. Irvine and for 37 <= n <= 50 to Sascha Kurz.  The representations below are given for ranges of n, and the integers indicate at what stage a given node is added (the letters A--Z encode the integers 10--36).
For instance, the first representation encodes the sequence of configurations
                  X
XX ; XXX ; XXX ; XXX ; etc.
XX ; XX  ; XXX ; XXX
----------
n = 4--11
435
0016
0027
----------
n = 12--19
..7
4003
00005
00006
1002
----------
n = 20--36
.GA78B
.93100C
E500000
F400000
.600000
.D2000
----------
n = 37--47
.60007
5000008
0000000
00000009
0000000A
4000003
.20001
----------
n = 48--50
..0000
.000000
20000000
00000000
00000000
.0000000
.000000
...001
----------
In particular, a(25)>=51, a(36)>=109, a(37)>=117, a(48)>=198, a(49)>=207, a(50)>=216.
----------
Another optimal configuration for a(8) = 4 due to Sascha Kurz:
.XX
XXXX
.XX
----------
Configurations and values for larger values of n can be found in the links below.
		

References

  • Elekes, Erdos, Similar configurations and pseudo grids, Coll. Math. Soc. Janos Bolyai 63 Intuitive Geometry, Budapest (Hungary), 1994.
  • J. Matousek, Lectures on Discrete Geometry, GTM 212, Springer, 2002.

Crossrefs

A348768 gives the number of inequivalent solutions.

Formula

If n = m^2, then a(n) >= m^2*(m^2-1)/12 (see A002415). If n = m^2-1, then a(n) >= (m-1)*(m-2)*(m^2+3*m+6)/12. - N. J. A. Sloane, Sep 28 2021

Extensions

a(13)-a(16) from Sean A. Irvine, Sep 23 2021
a(13)-a(16) confirmed and a(17) from Sascha Kurz, Oct 30 2021
Revised following a rich discussion on the seqfan mailing list, with contributions by the persons cited in the text, Allan Wechsler, Alex Meiburg, and Benoit Jubin. - Benoit Jubin, Oct 07 2021
Partially revised by N. J. A. Sloane, Nov 02 2021

A101095 Fourth difference of fifth powers (A000584).

Original entry on oeis.org

1, 28, 121, 240, 360, 480, 600, 720, 840, 960, 1080, 1200, 1320, 1440, 1560, 1680, 1800, 1920, 2040, 2160, 2280, 2400, 2520, 2640, 2760, 2880, 3000, 3120, 3240, 3360, 3480, 3600, 3720, 3840, 3960, 4080, 4200, 4320, 4440, 4560, 4680, 4800, 4920, 5040, 5160, 5280
Offset: 1

Views

Author

Cecilia Rossiter, Dec 15 2004

Keywords

Comments

Original Name: Shells (nexus numbers) of shells of shells of shells of the power of 5.
The (Worpitzky/Euler/Pascal Cube) "MagicNKZ" algorithm is: MagicNKZ(n,k,z) = Sum_{j=0..k+1} (-1)^j*binomial(n + 1 - z, j)*(k - j + 1)^n, with k>=0, n>=1, z>=0. MagicNKZ is used to generate the n-th accumulation sequence of the z-th row of the Euler Triangle (A008292). For example, MagicNKZ(3,k,0) is the 3rd row of the Euler Triangle (followed by zeros) and MagicNKZ(10,k,1) is the partial sums of the 10th row of the Euler Triangle. This sequence is MagicNKZ(5,k-1,2).

Crossrefs

Fourth differences of A000584, third differences of A022521, second differences of A101098, and first differences of A101096.
For other sequences based upon MagicNKZ(n,k,z):
...... | n = 1 | n = 2 | n = 3 | n = 4 | n = 5 | n = 6 | n = 7 | n = 8
--------------------------------------------------------------------------------------
z = 0 | A000007 | A019590 | ....... MagicNKZ(n,k,0) = T(n,k+1) from A008292 .......
z = 1 | A000012 | A040000 | A101101 | A101104 | A101100 | ....... | ....... | .......
z = 2 | A000027 | A005408 | A008458 | A101103 | thisSeq | ....... | ....... | .......
z = 3 | A000217 | A000290 | A003215 | A005914 | A101096 | ....... | ....... | .......
z = 4 | A000292 | A000330 | A000578 | A005917 | A101098 | ....... | ....... | .......
z = 5 | A000332 | A002415 | A000537 | A000583 | A022521 | ....... | A255181 | .......
z = 12 | A001288 | A057788 | ....... | A254870 | A254471 | A254683 | A254646 | A254642
z = 13 | A010965 | ....... | ....... | ....... | A254871 | A254472 | A254684 | A254647
z = 14 | A010966 | ....... | ....... | ....... | ....... | A254872 | ....... | .......
--------------------------------------------------------------------------------------
Cf. A047969.

Programs

  • Magma
    I:=[1,28,121,240,360]; [n le 5 select I[n] else 2*Self(n-1)-Self(n-2): n in [1..50]]; // Vincenzo Librandi, May 07 2015
    
  • Mathematica
    MagicNKZ=Sum[(-1)^j*Binomial[n+1-z, j]*(k-j+1)^n, {j, 0, k+1}];Table[MagicNKZ, {n, 5, 5}, {z, 2, 2}, {k, 0, 34}]
    CoefficientList[Series[(1 + 26 x + 66 x^2 + 26 x^3 + x^4)/(1 - x)^2, {x, 0, 50}], x] (* Vincenzo Librandi, May 07 2015 *)
    Join[{1,28,121,240},Differences[Range[50]^5,4]] (* or *) LinearRecurrence[{2,-1},{1,28,121,240,360},50] (* Harvey P. Dale, Jun 11 2016 *)
  • PARI
    a(n)=if(n>3, 120*n-240, 33*n^2-72*n+40) \\ Charles R Greathouse IV, Oct 11 2015
  • Sage
    [1,28,121]+[120*(k-2) for k in range(4,36)] # Danny Rorabaugh, Apr 23 2015
    

Formula

a(k+1) = Sum_{j=0..k+1} (-1)^j*binomial(n + 1 - z, j)*(k - j + 1)^n; n = 5, z = 2.
For k>3, a(k) = Sum_{j=0..4} (-1)^j*binomial(4, j)*(k - j)^5 = 120*(k - 2).
a(n) = 2*a(n-1) - a(n-2), n>5. G.f.: x*(1+26*x+66*x^2+26*x^3+x^4) / (1-x)^2. - Colin Barker, Mar 01 2012

Extensions

MagicNKZ material edited, Crossrefs table added, SeriesAtLevelR material removed by Danny Rorabaugh, Apr 23 2015
Name changed and keyword 'uned' removed by Danny Rorabaugh, May 06 2015
Previous Showing 41-50 of 121 results. Next