A010021 -id:A010021 - cp's OEIS frontend

A008527 Coordination sequence for body-centered tetragonal lattice.

1, 10, 34, 74, 130, 202, 290, 394, 514, 650, 802, 970, 1154, 1354, 1570, 1802, 2050, 2314, 2594, 2890, 3202, 3530, 3874, 4234, 4610, 5002, 5410, 5834, 6274, 6730, 7202, 7690, 8194, 8714, 9250, 9802, 10370, 10954, 11554, 12170, 12802, 13450, 14114, 14794, 15490, 16202, 16930, 17674

Offset: 0

N. J. A. Sloane

Also sequence found by reading the segment (1, 10) together with the line from 10, in the direction 10, 34, ..., in the square spiral whose vertices are the generalized hexagonal numbers A000217. - Omar E. Pol, Nov 02 2012

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
M. O'Keeffe, Coordination sequences for lattices, Zeit. f. Krist., 210 (1995), 905-908.
M. O'Keeffe, Coordination sequences for lattices, Zeit. f. Krist., 210 (1995), 905-908. [Annotated scanned copy]
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Apart from leading term, same as A108100.
Cf. A206399.
Cf. A016754 (SE), A054554 (NE), A054569 (SW), A053755 (NW), A033951 (S), A054552 (E), A054556 (N), A054567 (W) (Ulam spiral spokes).
A143839 (SSE) + A143855 (ESE) = A143838 (SSW) + A143856 (ENE) = A143854 (WSW) + A143861 (NNE) = A143859 (WNW) + A143860 (NNW) = even bisection = a(2n) = A010021(n).

Concatenation([1], List([1..40], n-> 2*(1+4*n^2) )); # G. C. Greubel, Nov 09 2019

[1] cat [2*(1 + 4*n^2): n in [1..50]]; // G. C. Greubel, Nov 09 2019

Maple
```
1, seq(8*k^2+2, k=1..50);
```

a[0]:= 1; a[n_]:= 8n^2 +2; Table[a[n], {n,0,50}] (* Alonso del Arte, Sep 06 2011 *)
LinearRecurrence[{3,-3,1},{1,10,34,74},50] (* Harvey P. Dale, Feb 13 2022 *)

vector(51, n, if(n==1,1, 2*(1+(2*n-2)^2)) ) \\ G. C. Greubel, Nov 09 2019

[1]+[2*(1+4*n^2) for n in (1..40)] # G. C. Greubel, Nov 09 2019

a(0) = 1; a(n) = 8*n^2+2 for n>0.
From Gary W. Adamson, Dec 27 2007: (Start)
a(n) = (2n+1)^2 + (2n-1)^2 for n>0.
Binomial transform of [1, 9, 15, 1, -1, 1, -1, 1, ...]. (End)
From Colin Barker, Apr 13 2012: (Start)
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3) for n>3.
G.f.: (1+x)*(1+6*x+x^2)/(1-x)^3. (End)
From Bruce J. Nicholson, Jul 31 2019: (Start) Assume n>0.
a(n) = A016754(n) + A016754(n-1).
a(n) = 2 * A053755(n).
a(n) = A054554(n+1) + A054569(n+1).
a(n) = A033951(n) + A054552(n).
a(n) = A054556(n+1) + A054567(n+1). (End)
E.g.f.: -1 + 2*exp(x)*(1 + 2*x)^2. - Stefano Spezia, Aug 02 2019
Sum_{n>=0} 1/a(n) = 3/4+1/8*Pi*coth(Pi/2) = 1.178172.... - R. J. Mathar, May 07 2024

A244082 a(n) = 32*n^2.

Original entry on oeis.org

0, 32, 128, 288, 512, 800, 1152, 1568, 2048, 2592, 3200, 3872, 4608, 5408, 6272, 7200, 8192, 9248, 10368, 11552, 12800, 14112, 15488, 16928, 18432, 20000, 21632, 23328, 25088, 26912, 28800, 30752, 32768, 34848, 36992, 39200, 41472, 43808, 46208, 48672, 51200

Offset: 0

Wesley Ivan Hurt, Jun 19 2014

Geometric connections of a(n) to the area and perimeter of a square.
Area:
. half the area of a square with side 8n (cf. A008590);
. area of a square with diagonal 8n (cf. A008590);
. twice the area of a square with side 4n (cf. A008586);
. four times the area of a square with diagonal 4n (cf. A008586);
. eight times the area of a square with side 2n (cf. A005843);
. sixteen times the area of a square with diagonal 2n (cf. A005843);
. thirty two times the area of a square with side n (cf. A001477);
. sixty four times the area of a square with diagonal n (cf. A001477).
Perimeter:
. perimeter of a square with side 8n^2 (cf. A139098);
. twice the perimeter of a square with side 4n^2 (cf. A016742);
. four times the perimeter of a square with side 2n^2 (cf. A001105);
. eight times the perimeter of a square with side n^2 (cf. A000290).
Sequence found by reading the line from 0, in the direction 0, 32, ..., in the square spiral whose vertices are the generalized 18-gonal numbers. - Omar E. Pol, May 10 2018

Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Pentasection of A077221, A181900.
Cf. A000290, A001105, A010021, A016742, A016802, A139098.
Cf. A001477, A005843, A008586, A008590, A139098, A174312.

Magma
```
[32*n^2 : n in [0..50]];
```

A244082:=n->32*n^2; seq(A244082(n), n=0..50);

32 Range[0, 50]^2 (* or *)
Table[32 n^2, {n, 0, 50}] (* or *)
CoefficientList[Series[32 x (1 + x)/(1 - x)^3, {x, 0, 30}], x]

a(n)=32*n^2 \\ Charles R Greathouse IV, Jun 17 2017

G.f.: 32*x*(1+x)/(1-x)^3.
a(n) =  2 * A016802(n).
a(n) =  4 * A139098(n).
a(n) =  8 * A016742(n).
a(n) = 16 * A001105(n).
a(n) = 32 * A000290(n).
a(n) = A010021(n) - 2 for n > 0. - Bruno Berselli, Jun 24 2014
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Wesley Ivan Hurt, Nov 19 2021
From Elmo R. Oliveira, Dec 02 2024: (Start)
E.g.f.: 32*x*(1 + x)*exp(x).
a(n) = n*A174312(n) = A139098(2*n). (End)

A158563 a(n) = 32*n^2 - 1.

Original entry on oeis.org

31, 127, 287, 511, 799, 1151, 1567, 2047, 2591, 3199, 3871, 4607, 5407, 6271, 7199, 8191, 9247, 10367, 11551, 12799, 14111, 15487, 16927, 18431, 19999, 21631, 23327, 25087, 26911, 28799, 30751, 32767, 34847, 36991, 39199, 41471, 43807, 46207, 48671, 51199, 53791

Offset: 1

Vincenzo Librandi, Mar 21 2009

The identity (32*n^2 - 1)^2 - (256*n^2 - 16)*(2*n)^2 = 1 can be written as a(n)^2 - A158562(n)*A005843(n)^2 = 1. [comment rewritten by R. J. Mathar, Oct 16 2009]
From Omar E. Pol, Apr 21 2021: (Start)
Sequence found by reading the line from 31, in the direction 31, 127, ..., in the rectangular spiral whose vertices are the generalized 18-gonal numbers A274979.
The spiral begins as follows:
       46_   _   _   _   _   _   _   _   _   _18
        |                                                        |
        |                            0                           |
        |                            |  _   _   _   _  |
        |                            1                           15
        |
       51
(End)

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Vincenzo Librandi, X^2-AY^2=1, Math Forum, 2007. [Wayback Machine link]
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A005843, A010021, A158562, A158575, A244082.

Cf. A274979 (generalized 18-gonal numbers).

[32*n^2-1: n in [1..40]]; // Vincenzo Librandi, Sep 11 2013

32 Range[40]^2 - 1  (* Harvey P. Dale, Mar 04 2011 *)
CoefficientList[Series[(- 31 - 34 x + x^2) / (x - 1)^3, {x, 0, 40}], x] (* Vincenzo Librandi, Sep 11 2013 *)

a(n)=32*n^2-1 \\ Charles R Greathouse IV, Jun 17 2017

G.f.: x*(-31-34*x+x^2)/(x-1)^3.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
a(n) = A244082(n) - 1. - Omar E. Pol, Apr 21 2021
From Amiram Eldar, Mar 09 2023: (Start)
Sum_{n>=1} 1/a(n) = (1 - cot(Pi/(4*sqrt(2)))*Pi/(4*sqrt(2)))/2.
Sum_{n>=1} (-1)^(n+1)/a(n) = (cosec(Pi/(4*sqrt(2)))*Pi/(4*sqrt(2)) - 1)/2. (End)
E.g.f.: 1 + exp(x)*(32*x^2 + 32*x - 1). - Elmo R. Oliveira, Jan 25 2025

A158575 a(n) = 32*n^2 + 1.

Original entry on oeis.org

1, 33, 129, 289, 513, 801, 1153, 1569, 2049, 2593, 3201, 3873, 4609, 5409, 6273, 7201, 8193, 9249, 10369, 11553, 12801, 14113, 15489, 16929, 18433, 20001, 21633, 23329, 25089, 26913, 28801, 30753, 32769, 34849, 36993, 39201, 41473, 43809, 46209, 48673, 51201, 53793

Offset: 0

Vincenzo Librandi, Mar 21 2009

The identity (32*n^2 + 1)^2 - (256*n^2 + 16)*(2*n)^2 = 1 can be written as a(n)^2-A158574(n)*A005843(n)^2 = 1. - Comment rewritten by R. J. Mathar, Oct 16 2009

Sequence found by reading the line segment from 1 to 33 together with the line from 33, in the direction 33, 129, ..., in the square spiral whose vertices are the generalized 18-gonal numbers A274979. - Omar E. Pol, Apr 21 2021

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
Vincenzo Librandi, X^2-AY^2=1, Math Forum, 2007. [Wayback Machine link]
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000217, A005843, A010021, A081585, A158563, A158574, A244082.

Cf. A274979 (generalized 18-gonal numbers).

I:=[1, 33, 129]; [n le 3 select I[n] else 3*Self(n-1)-3*Self(n-2)+1*Self(n-3): n in [1..50]]; // Vincenzo Librandi, Feb 15 2012

LinearRecurrence[{3, -3, 1}, {1, 33, 129}, 50] (* Vincenzo Librandi, Feb 15 2012 *)
32*Range[0,40]^2+1 (* Harvey P. Dale, Jul 20 2021 *)

for(n=0, 50, print1(32*n^2+1", ")); \\ Vincenzo Librandi, Feb 15 2012

G.f.: (1+30*x+33*x^2)/(1-x)^3.
a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3).
For n > 0 a(n) = sqrt(8*(A000217(4*n-1)^2 + A000217(4*n)^2) + 1). - J. M. Bergot, Sep 03 2015
a(n) = A244082(n) + 1. - Omar E. Pol, Apr 21 2021
From Amiram Eldar, Mar 09 2023: (Start)
Sum_{n>=0} 1/a(n) = (1 + coth(Pi/(4*sqrt(2)))*Pi/(4*sqrt(2)))/2.
Sum_{n>=0} (-1)^n/a(n) = (1 + cosech(Pi/(4*sqrt(2)))*Pi/(4*sqrt(2)))/2. (End)
From Elmo R. Oliveira, Jan 16 2025: (Start)
E.g.f.: exp(x)*(1 + 32*x + 32*x^2).
a(n) = A081585(2*n). (End)

a(0) added by R. J. Mathar, Oct 16 2009

cp's OEIS Frontend

A008527 Coordination sequence for body-centered tetragonal lattice.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A244082 a(n) = 32*n^2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A158563 a(n) = 32*n^2 - 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

A158575 a(n) = 32*n^2 + 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

Extensions