A071619 -id:A071619 - cp's OEIS frontend

A005893 Number of points on surface of tetrahedron; coordination sequence for sodalite net (equals 2*n^2+2 for n > 0).

1, 4, 10, 20, 34, 52, 74, 100, 130, 164, 202, 244, 290, 340, 394, 452, 514, 580, 650, 724, 802, 884, 970, 1060, 1154, 1252, 1354, 1460, 1570, 1684, 1802, 1924, 2050, 2180, 2314, 2452, 2594, 2740, 2890, 3044, 3202, 3364, 3530, 3700, 3874, 4052, 4234

Offset: 0

N. J. A. Sloane

Number of n-matchings of the wheel graph W_{2n} (n > 0). Example: a(2)=10 because in the wheel W_4 (rectangle ABCD and spokes OA,OB,OC,OD) we have the 2-matchings: (AB, OC), (AB, OD), (BC, OA), (BC,OD), (CD,OA), (CD,OB), (DA,OB), (DA,OC), (AB,CD) and (BC,DA). - Emeric Deutsch, Dec 25 2004
For n > 0 a(n) is the difference of two tetrahedral (or pyramidal) numbers: binomial(n+3, 3) = (n+1)(n+2)(n+3)/6. a(n) = A000292(n+1) - A000292(n-3) = (n+1)(n+2)(n+3)/6 - (n-3)(n-2)(n-1)/6. - Alexander Adamchuk, May 20 2006; updated by Peter Munn, Aug 25 2017 due to changed offset in A000292
Equals binomial transform of [1, 3, 3, 1, -1, 1, -1, 1, -1, 1, ...]. Binomial transform of A005893 = nonzero terms of A053545: (1, 5, 19, 63, 191, ...). - Gary W. Adamson, Apr 28 2008
Disregarding the terms < 10, the sums of four consecutive triangular numbers (A000217). - Rick L. Shepherd, Sep 30 2009
Use a set of n concentric circles where n >= 0 to divide the plane. a(n) is the maximal number of regions after the 2nd division. - Frank M Jackson, Sep 07 2011
Euler transform of length 4 sequence [4, 0, 0, -1]. - Michael Somos, May 14 2014
Also, growth series for affine Coxeter group (or affine Weyl group) A_3 or D_3. - N. J. A. Sloane, Jan 11 2016
For n > 2 the generalized Pell's equation x^2 - 2*(a(n) - 2)y^2 = (a(n) - 4)^2 has a finite number of positive integer solutions. - Muniru A Asiru, Apr 19 2016
Union of A188896, A277449, {1,4}. - Muniru A Asiru, Nov 25 2016
Interleaving of A008527 and A108099. - Bruce J. Nicholson, Oct 14 2019

			G.f. = 1 + 4*x + 10*x^2 + 20*x^3 + 34*x^4 + 52*x^5 + 74*x^6 + 100*x^7 + ...

N. Bourbaki, Groupes et Algèbres de Lie, Chap. 4, 5 and 6, Hermann, Paris, 1968. See Chap. VI, Section 4, Problem 10b, page 231, W_a(t).
H. S. M. Coxeter, "Polyhedral numbers," in R. S. Cohen et al., editors, For Dirk Struik. Reidel, Dordrecht, 1974, pp. 25-35.
B. Grünbaum, Uniform tilings of 3-space, Geombinatorics, 4 (1994), 49-56. See tiling #28.
R. W. Marks and R. B. Fuller, The Dymaxion World of Buckminster Fuller. Anchor, NY, 1973, p. 46.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Vincenzo Librandi, Table of n, a(n) for n = 0..10000
J. H. Conway and N. J. A. Sloane, Low-Dimensional Lattices VII: Coordination Sequences, Proc. Royal Soc. London, A453 (1997), 2369-2389 (pdf).
Steven Edwards and William Griffiths, On Generalized Delannoy Numbers, J. Int. Seq., Vol. 23 (2020), Article 20.3.6.
J. M. Grau, C. Miguel, and A. M. Oller-Marcén, Generalized Quaternion Rings over Z/nZ for an odd n, arXiv:1706.04760 [math.RA], 2017. See Theorem 1, p. 10.
Guo-Niu Han, Enumeration of Standard Puzzles. [Cached copy.]
Milan Janjić, Hessenberg Matrices and Integer Sequences, J. Int. Seq. 13 (2010), Article 10.7.8.
M. O'Keeffe, N-dimensional diamond, sodalite and rare sphere packings, Acta Cryst., A 47 (1991), 749-753.
Simon Plouffe, Approximations de séries génératrices et quelques conjectures, Dissertation, Université du Québec à Montréal, 1992; arXiv:0911.4975 [math.NT], 2009.
Simon Plouffe, 1031 Generating Functions, Appendix to Thesis, Montreal, 1992.
Reticular Chemistry Structure Resource, sod.
Aditya Sivakumar and Dmitri Tymoczko, Intuitive Musical Homotopy, 2018.
B. K. Teo and N. J. A. Sloane, Magic numbers in polygonal and polyhedral clusters, Inorgan. Chem. 24 (1985),4545-4558. DOI: 10.1021/ic00220a025.
Index entries for linear recurrences with constant coefficients, signature (3,-3,1).

Cf. A000217, A000292, A053545, A206399.
Cf. similar sequences listed in A255843.
The growth series for the affine Coxeter groups D_3 through D_12 are A005893 and A266759-A266767.
For partial sums see A005894.
The 28 uniform 3D tilings: cab: A299266, A299267; crs: A299268, A299269; fcu: A005901, A005902; fee: A299259, A299265; flu-e: A299272, A299273; fst: A299258, A299264; hal: A299274, A299275; hcp: A007899, A007202; hex: A005897, A005898; kag: A299256, A299262; lta: A008137, A299276; pcu: A005899, A001845; pcu-i: A299277, A299278; reo: A299279, A299280; reo-e:  A299281, A299282; rho: A008137, A299276; sod: A005893, A005894; sve: A299255, A299261; svh: A299283, A299284; svj: A299254, A299260; svk: A010001, A063489; tca: A299285, A299286; tcd: A299287, A299288; tfs: A005899, A001845; tsi: A299289, A299290; ttw: A299257, A299263; ubt: A299291, A299292; bnn: A007899, A007202. See the Proserpio link in A299266 for overview.

[2*n^2-0^n+2: n in [0..60]]; // Vincenzo Librandi, Sep 27 2011

A005893:=-(z+1)*(1+z^2)/(z-1)^3; # Simon Plouffe in his 1992 dissertation

Join[{1}, Table[2*(n + 1)^2 + 2, {n, 0, 200}]] (* Vladimir Joseph Stephan Orlovsky, Jul 10 2011 *)
Join[{1},LinearRecurrence[{3,-3,1},{4,10,20},50]] (* Harvey P. Dale, Feb 26 2012 *)
a[ n_] := SeriesCoefficient[ (1 - x^4) / (1 - x)^4, {x, 0, Abs@n}]; (* Michael Somos, May 14 2014 *)
a[ n_] := 2 n^2 + 2 - Boole[n == 0]; (* Michael Somos, May 14 2014 *)

a(n)=2*n^2-0^n+2 \\ Charles R Greathouse IV, Sep 24 2015

G.f.: (1 - x^4)/(1-x)^4.
a(n) = A071619(n-1) + A071619(n) + A071619(n+1), n > 0. - Ralf Stephan, Apr 26 2003
a(n) = binomial(n+3, 3) - binomial(n-1, 3) for n >= 1. - Mitch Harris, Jan 08 2008
a(n) = (n+1)^2 + (n-1)^2. - Benjamin Abramowitz, Apr 14 2009
a(n) = A000217(n-2) + A000217(n-1) + A000217(n) + A000217(n+1) for n >= 2. - Rick L. Shepherd, Sep 30 2009
a(n) = 2*n^2 - 0^n + 2. - Vincenzo Librandi, Sep 27 2011
a(0)=1, a(1)=4, a(2)=10, a(3)=20, a(n) = 3*a(n-1) - 3*a(n-2) + a(n-3). - Harvey P. Dale, Feb 26 2012
a(n) = A228643(n+1,2) for n > 0. - Reinhard Zumkeller, Aug 29 2013
a(n) = a(-n) for all n in Z. - Michael Somos, May 14 2014
For n >= 2: a(n) = a(n-1) + 4*n - 2. - Bob Selcoe, Mar 22 2016
E.g.f.: -1 + 2*(1 + x + x^2)*exp(x). - Ilya Gutkovskiy, Apr 19 2016
a(n) = 2*A002522(n), n>0. - R. J. Mathar, May 30 2022
From Amiram Eldar, Sep 16 2022: (Start)
Sum_{n>=0} 1/a(n) = (coth(Pi)*Pi + 3)/4.
Sum_{n>=0} (-1)^n/a(n) = (cosech(Pi)*Pi + 3)/4. (End)
Empirical: Integral_{u=-oo..+oo} sigmoid(u)*log(sigmoid(n * u)) du = -Pi^2*a(n) / (24*n), where sigmoid(x) = 1/(1+exp(-x)). Also works for non-integer n>0. - Carlo Wood, Dec 04 2023
Let P(k,n) be the n-th k-gonal number. Then P(a(k),n) = (k*n-k+1)^2 + (k-1)^2*(n-1). - Charlie Marion, May 15 2024

A042968 Numbers not divisible by 4.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 18, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 45, 46, 47, 49, 50, 51, 53, 54, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 81, 82, 83, 85, 86, 87, 89, 90, 91, 93, 94, 95, 97, 98, 99, 101, 102

Offset: 1

N. J. A. Sloane, Dec 11 1999

Equivalently, numbers whose square part is odd. Cf. A028982. - Peter Munn, Jul 14 2020
More generally the sequence of numbers not divisible by some fixed integer m >= 2 is given by a(n,m) = 1 + n + floor(n/(m-1)). - Benoit Cloitre, Jul 11 2009
Also a(n,m) = floor((m*n-1)/(m-1)) [with offset 1]. - Gary Detlefs, May 14 2011
Numbers not having more even than odd divisors: A048272(a(n)) >= 0. - Reinhard Zumkeller, Jan 21 2012
Extending the comments of Benoit Cloitre (Jul 11 2009) and Gary Detlefs (May 14 2011), the g.f. is A(m,x) = (1-x^m) / ((1-x^(m-1))*(1-x)^2) where m >= 2 is fixed. - Werner Schulte, Apr 26 2018

			G.f. = 1 + 2*x + 3*x^2 + 5*x^3 + 6*x^4 + 7*x^5 + 9*x^6 + 10*x^7 + 11*x^8 + ... - _Michael Somos_, Jun 17 2018

G. C. Greubel, Table of n, a(n) for n = 1..10000
Index entries for linear recurrences with constant coefficients, signature (1,0,1,-1).

Cf. A001651, A001935, A028982, A070048, A042965, A064680.
Cf. A071619 (partial sums); A008586 (complement).
Numbers that are congruent to {k0,k1,k2} mod 4: A004772, A004773, A042965, a(n).

a042968 = (`div` 3) . (subtract 1) . (* 4)
a042968_list = filter ((/= 0) . (`mod` 4)) [1..]
-- Reinhard Zumkeller, Sep 02 2012

[n+1+Floor(n/3): n in [0..80]]; // Vincenzo Librandi, Aug 03 2015

seq(n+floor((n-1)/3), n=1..80); # Muniru A Asiru, Feb 17 2019

Select[Table[n,{n,200}], Mod[#,4] != 0&] (* Vladimir Joseph Stephan Orlovsky, Feb 18 2011 *)
LinearRecurrence[{1,0,1,-1},{1,2,3,5},80]  (* or *) Drop[Range[110],{4,-1,4}] (* Harvey P. Dale, Jan 07 2023 *)

PARI
```
{a(n) = 1 + n + n\3};
```

def A042968(n): return n+(n-1)//3 # Chai Wah Wu, Apr 13 2025

[1+n+floor(n/3) for n in (0..80)] # G. C. Greubel, Feb 17 2019

a(n) = a(n-1) + a(n-3) - a(n-4).
a(n) = a(n-3) + 4, with a(1) = 1.
G.f.: x * (1+x) * (1+x^2) / ( (1+x+x^2)*(1-x)^2 ). - Michael Somos, Jan 12 2000
A064680(A064680(a(n))) = a(n). - Reinhard Zumkeller, Oct 19 2001
Nearest integer to (Sum_{k>n} 1/k^4)/(Sum_{k>n} 1/k^5). - Benoit Cloitre, Jun 12 2003
a(n) = n + 1 + floor(n/3). - Benoit Cloitre, Jul 11 2009
a(n) = floor((4*n+3)/3). - Gary Detlefs, May 14 2011
A214546(a(n)) >= 0 for n > 0. - Reinhard Zumkeller, Jul 20 2012
a(n) = 2*n - ceiling(2*n/3) + 1. - Arkadiusz Wesolowski, Sep 21 2012
Sum_{k=0..n} a(n) = A071619(n+1). - L. Edson Jeffery, Jul 30 2014
The g.f. A(x) satisfies x*A(x)^2 = (B(x)/x)^2 + (B(x)/x), where B(x) is the o.g.f. of A042965. - Peter Bala, Apr 12 2017
a(n) = (12*n + 6 + 3*cos(2*n*Pi/3) + sqrt(3)*sin(2*n*Pi/3))/9. - Wesley Ivan Hurt, Sep 30 2017
Euler transform of length 4 sequence [2, 0, 1, -1]. - Michael Somos, Jun 17 2018
a(n) = -a(-1-n) for all n in Z. - Michael Somos, Jun 17 2018
E.g.f.: (2/3)*exp(x)*(1 + 2*x) + (1/9)*exp(-x/2)*(3*cos(sqrt(3)*x/2) + sqrt(3)*sin(sqrt(3)*x/2)). - Stefano Spezia, Nov 16 2019
a(n) = (12*n + 6 + w^(2*n)*(w + 2) - w^n*(w - 1))/9 where w = (-1 + sqrt(-3))/2. - Guenther Schrack, Jun 07 2021
Sum_{n>=1} (-1)^(n+1)/a(n) = (2*sqrt(2)-1)*Pi/8. - Amiram Eldar, Dec 05 2021

Edited by Peter Munn, Nov 16 2019

I restored my original (1999) definition and offset, which in the intervening 21 years had been lost. - N. J. A. Sloane, Jun 12 2021

A156040 Number of compositions (ordered partitions) of n into 3 parts (some of which may be zero), where the first is at least as great as each of the others.

Original entry on oeis.org

1, 1, 3, 4, 6, 8, 11, 13, 17, 20, 24, 28, 33, 37, 43, 48, 54, 60, 67, 73, 81, 88, 96, 104, 113, 121, 131, 140, 150, 160, 171, 181, 193, 204, 216, 228, 241, 253, 267, 280, 294, 308, 323, 337, 353, 368, 384, 400, 417, 433, 451, 468, 486, 504, 523, 541, 561, 580, 600

Offset: 0

Jack W Grahl, Feb 02 2009, Feb 11 2009

For n = 1, 2 these are just the triangular numbers. a(n) is always at least 1/3 of the corresponding triangular number, since each partition of this type gives up to three ordered partitions with the same cyclical order.
An alternative definition, which avoids using parts of size 0: a(n) is the third diagonal of A184957. - N. J. A. Sloane, Feb 27 2011
Diagonal sums of the triangle formed by rows T(2, k) k = 0, 1, ..., 2m of ascending m-nomial triangles (see A004737):
1
1 2 1
1 2 3 2 1
1 2 3 4 3 2 1
1 2 3 4 5 4 3 2 1
1 2 3 4 5 6 5 4 3 2 1
- Bob Selcoe, Feb 07 2014
Arrange A004396 in rows successively shifted to the right two spaces and sum the columns:
1  1  2  3  3  4  5  5  6 ...
      1  1  2  3  3  4  5 ...
            1  1  2  3  3 ...
                  1  1  2 ...
                        1 ...
------------------------------
1  1  3  4  6  8 11 13 17 ... - L. Edson Jeffery, Jul 30 2014
a(n) is the dimension of three-dimensional (2n + 2)-homogeneous polynomial vector fields with full tetrahedral symmetry (for a given orthogonal representation), and which are solenoidal. - Giedrius Alkauskas, Sep 30 2017
Also the number of compositions of n + 3 into three parts, the first at least as great as each of the other two. Also the number of compositions of n + 4 into three parts, the first strictly greater than each of the other two. - Gus Wiseman, Oct 09 2020

			G.f. = 1 + x + 3*x^2 + 4*x^3 + 6*x^4 + 8*x^5 + 11*x^6 + 13*x^7 + 17*x^8 + 20*x^9 + ...
The a(4) = 6 compositions of 4 are: (4 0 0), (3 1 0), (3 0 1), (2 2 0), (2 1 1), (2 0 2).
From _Gus Wiseman_, Oct 05 2020: (Start)
The a(0) = 1 through a(7) = 13 triples of nonnegative integers summing to n where the first is at least as great as each of the other two are:
  (000)  (100)  (101)  (111)  (202)  (212)  (222)  (313)
                (110)  (201)  (211)  (221)  (303)  (322)
                (200)  (210)  (220)  (302)  (312)  (331)
                       (300)  (301)  (311)  (321)  (403)
                              (310)  (320)  (330)  (412)
                              (400)  (401)  (402)  (421)
                                     (410)  (411)  (430)
                                     (500)  (420)  (502)
                                            (501)  (511)
                                            (510)  (520)
                                            (600)  (601)
                                                   (610)
                                                   (700)
(End)

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Giedrius Alkauskas, Projective and polynomial superflows. I, arxiv.org/1601.06570 [math.AG] (2017), Section 4.3.
Index entries for linear recurrences with constant coefficients, signature (1,1,0,-1,-1,1).

For compositions into 4 summands see A156039; also see A156041 and A156042.
Cf. A184957, A071619 (bisection).
A001399(n-2)*2 is the strict case.
A001840(n-2) is the version with opposite relations.
A001840(n-1) is the version with strict opposite relations.
A069905 is the case with strict relations.
A014311 ranks 3-part compositions, with strict case A337453.
A014612 ranks 3-part partitions, with strict case A007304.
Cf. A000217, A001523, A211540, A218004, A220377, A337483, A337484.

a:= proc(n) local m, r; m := iquo(n, 6, 'r'); (4 +6*m +2*r) *m + [1, 1, 3, 4, 6, 8][r+1] end: seq(a(n), n=0..60); # Alois P. Heinz, Jun 14 2009

nn = 58; CoefficientList[Series[x^3/(1 - x^2)^2/(1 - x^3) + 1/(1 - x^2)^2/(1 - x), {x, 0, nn}], x] (* Geoffrey Critzer, Jul 14 2013 *)
CoefficientList[Series[(1 + x^2)/((1 + x) * (1 + x + x^2) * (1 - x)^3), {x, 0, 58}], x] (* L. Edson Jeffery, Jul 29 2014 *)
LinearRecurrence[{1, 1, 0, -1, -1, 1}, {1, 1, 3, 4, 6, 8}, 60] (* Harvey P. Dale, May 28 2015 *)
Table[Length[Select[Join@@Permutations/@IntegerPartitions[n+3,{3}],#[[1]]>=#[[2]]&&#[[1]]>=#[[3]]&]],{n,0,15}] (* Gus Wiseman, Oct 05 2020*)

{a(n) = n*(n+4)\6 + 1}; /* Michael Somos, Mar 26 2017 */

G.f.: (x^2+1) / (1-x-x^2+x^4+x^5-x^6). - Alois P. Heinz, Jun 14 2009
Slightly nicer g.f.: (1+x^2)/((1-x)*(1-x^2)*(1-x^3)). - N. J. A. Sloane, Apr 29 2011
a(n) = A007590(n+2) - A000212(n+2). - Richard R. Forberg, Dec 08 2013
a(2*n) = A071619(n+1). - L. Edson Jeffery, Jul 29 2014
a(n) = a(n-1) + a(n-2) - a(n-4) - a(n-5) + a(n-6), with a(0) = 1, a(1) = 1, a(2) = 3, a(3) = 4, a(4) = 6, a(5) = 8. - Harvey P. Dale, May 28 2015
a(n) = (n^2 + 4*n + 3)/6 + IF(MOD(n, 2) = 0, 1/2) + IF(MOD(n, 3) = 1, -1/3). - Heinrich Ludwig, Mar 21 2017
a(n) = 1 + floor((n^2 + 4*n)/6). - Giovanni Resta, Mar 21 2017
Euler transform of length 4 sequence [1, 2, 1, -1]. - Michael Somos, Mar 26 2017
a(n) = a(-4 - n) for all n in Z. - Michael Somos, Mar 26 2017
0 = a(n)*(-1 + a(n) - 2*a(n+1) - 2*a(n+2) + 2*a(n+3)) + a(n+1)*(+1 + a(n+1) + 2*a(n+2) - 2*a(n+3)) + a(n+2)*(+1 + a(n+2) - 2*a(n+3)) + a(n+3)*(-1 + a(n+3)) for all n in Z. - Michael Somos, Mar 26 2017
a(n) = round((n+1)*(n+3)/6). - Bill McEachen, Feb 16 2021
Sum_{n>=0} 1/a(n) = 3/2 + Pi^2/36 + (tan(c1)-1)*c1 + 3*c2*sinh(c2)/(1+2*cosh(c2)), where c1 = Pi/(2*sqrt(3)) and c2 = Pi*sqrt(2)/3. - Amiram Eldar, Dec 10 2022
E.g.f.: ((16 + 15*x + 3*x^2)*cosh(x) + 2*exp(-x/2)*(cos(sqrt(3)*x/2) - sqrt(3)*sin(sqrt(3)*x/2)) + (7 + 15*x + 3*x^2)*sinh(x))/18. - Stefano Spezia, Apr 05 2023

More terms from Alois P. Heinz, Jun 14 2009

A008747 Expansion of (1+x^4)/((1-x)(1-x^2)(1-x^3)).

Original entry on oeis.org

1, 1, 2, 3, 5, 6, 9, 11, 14, 17, 21, 24, 29, 33, 38, 43, 49, 54, 61, 67, 74, 81, 89, 96, 105, 113, 122, 131, 141, 150, 161, 171, 182, 193, 205, 216, 229, 241, 254, 267, 281, 294, 309, 323, 338, 353, 369, 384, 401, 417, 434, 451, 469, 486, 505, 523, 542, 561

Offset: 0

N. J. A. Sloane

For n>=1, the set {A008747(6n+-1)} is the set of numbers of the form a^2 + 5*(a+1)^2 for -inf < a < inf. Furthermore the set A008747(6n) is A033581(n). - Kieren MacMillan, Dec 19 2007

For n>1, a(n-1) is the number of aperiodic necklaces (Lyndon words) with k<=3 black beads and n-k white beads. For n=4 we have for example a(3)=3 aperiodic necklaces: BWWW, BBWW and BBBW. BWBW is periodic and is not counted. - Herbert Kociemba, Oct 23 2016

			G.f. = 1 + x + 2*x^2 + 3*x^3 + 5*x^4 + 6*x^5 + 9*x^6 + 11*x^7 + 14*x^8 + ...

Seiichi Manyama, Table of n, a(n) for n = 0..10000
Robert Morris, Minimal percolating sets in bootstrap percolation, arXiv:math/0702370 [math.CO], 2007-2008.
Index entries for linear recurrences with constant coefficients, signature (1,1,0,-1,-1,1).

Cf. A008747, A033581, A051274, A061347, A071619, A077043.

a:=[1,1,2,3,5,6];; for n in [7..60] do a[n]:=a[n-1]+a[n-2]-a[n-4] -a[n-5]+a[n-6]; od; a; # G. C. Greubel, Aug 03 2019

R:=PowerSeriesRing(Integers(), 60); Coefficients(R!( (1+x^4)/((1-x)*(1-x^2)*(1-x^3)) )); // G. C. Greubel, Aug 03 2019

A008747:=n->ceil((n+1)^2/6): seq(A008747(n), n=0..100); # Wesley Ivan Hurt, Oct 25 2016

CoefficientList[Series[(1+x^4)/((1-x)(1-x^2)(1-x^3)),{x,0,60}],x] (* or *) LinearRecurrence[{1,1,0,-1,-1,1},{1,1,2,3,5,6},60] (* Harvey P. Dale, Sep 05 2012 *)

Vec((1+x^4)/((1-x)*(1-x^2)*(1-x^3))+O(x^60)) \\ Charles R Greathouse IV, Sep 25 2012

((1+x^4)/((1-x)*(1-x^2)*(1-x^3))).series(x, 60).coefficients(x, sparse=False) # G. C. Greubel, Aug 03 2019

G.f.: (1+x^4)/((1-x)*(1-x^2)*(1-x^3)).
a(n) = ceiling((n+1)^2/6).
a(n) = (12*n + 23 + 6*n^2 + 9*(-1)^n + 4*A061347(n))/36. - R. J. Mathar, Mar 15 2011
a(0)=1, a(1)=1, a(2)=2, a(3)=3, a(4)=5, a(5)=6, a(n) = a(n-1) + a(n-2) - a(n-4) - a(n-5) + a(n-6) for n > 5. - Harvey P. Dale, Sep 05 2012
From Michael Somos, Oct 25 2016: (Start)
Euler transform of length 8 sequence [ 1, 1, 1, 1, 0, 0, 0, -1].
a(n) = a(-2-n) for all n in Z.
a(2*n-1) = A071619(n).
a(3*n-1) = 2*A077043(n).
a(n) - a(n-1) = A051274(n). (End)

A071318 Lesser of 2 consecutive numbers which are cubefree and not squarefree, i.e., numbers k such that both k and k+1 are in A067259.

Original entry on oeis.org

44, 49, 75, 98, 99, 116, 147, 171, 244, 260, 275, 315, 332, 363, 387, 475, 476, 507, 524, 531, 548, 549, 603, 604, 636, 692, 724, 725, 747, 764, 774, 819, 844, 845, 846, 867, 908, 924, 931, 963, 980, 1035, 1075, 1083, 1179, 1196, 1251, 1274, 1275, 1324

Offset: 1

Labos Elemer, May 29 2002

nonn

The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 0, 5, 41, 407, 4125, 41215, 412331, 4123625, 41236308, ... . Apparently, the asymptotic density of this sequence exists and equals 0.041236... . - Amiram Eldar, Jan 18 2023

The asymptotic density of this sequence is Product_{p prime} (1 - 2/p^3) - 2 * Product_{p prime} (1 - 1/p^2 - 1/p^3) + Product_{p prime} (1 - 2/p^2) = 0.041236147082334172926... . - Amiram Eldar, Jan 05 2024

			75 is a term since 75 = 3*5^2 and 76 = 2^2*19.

Amiram Eldar, Table of n, a(n) for n = 1..10000 (terms 1..1000 from T. D. Noe)

Subsequence of A067259, A340152.
Cf. A007674, A063528, A071619, A071620, A071621, A071622, A071623.
Cf. A008966, A051903, A212793.

a071318 n = a071318_list !! (n-1)
a071318_list = [x | x <- [1..],  a212793 x == 1, a008966 x == 0,
                    let y = x+1, a212793 y == 1, a008966 y == 0]
-- Reinhard Zumkeller, May 27 2012

With[{s = Select[Range[1350], And[MemberQ[#, 2], FreeQ[#, k_ /; k > 2]] &@ FactorInteger[#][[All, -1]] &]}, Function[t, Part[s, #] &@ Position[t, 1][[All, 1]]]@ Differences@ s] (* Michael De Vlieger, Jul 30 2017 *)

isok(n) = (n>1) && (vecmax(factor(n)[, 2])==2) && (vecmax(factor(n+1)[, 2])==2); \\ Michel Marcus, Aug 02 2017

A051903(k) = A051903(k+1) = 2 when k is a term.

A143976 Rectangular array R by antidiagonals: label each unit square in the first quadrant lattice by its northeast vertex (x,y) and mark squares having x + y == 1 (mod 3); then R(m,n) is the number of UNmarked squares in the rectangle [0,m] X [0,n].

Original entry on oeis.org

1, 2, 2, 2, 3, 2, 3, 4, 4, 3, 4, 6, 6, 6, 4, 4, 7, 8, 8, 7, 4, 5, 8, 10, 11, 10, 8, 5, 6, 10, 12, 14, 14, 12, 10, 6, 6, 11, 14, 16, 17, 16, 14, 11, 6, 7, 12, 16, 19, 20, 20, 19, 16, 12, 7, 8, 14, 18, 22, 24, 24, 24, 22, 18, 14, 8, 8, 15, 20, 24, 27, 28, 28, 27, 24, 20, 15, 8

Offset: 1

Clark Kimberling, Sep 06 2008

Rows numbered 3,6,9,12,15,... are, except for initial terms, multiples of (1,2,3,4,5,6,7,...) = A000027.

			Northwest corner:
  1  2  2  3  4  4  5
  2  3  4  6  7  8 10
  2  4  6  8 10 12 14
  3  6  8 11 14 16 18
  4  7 10 14 17 20 24
See A143974.

Stefano Spezia, First 140 antidiagonals of the array, flattened

Rows: A004523, A004772, A005843, A047399, et al.
Main diagonal: A071619.
Cf. A143974, A143977, A143979.

T[m_,n_]:=m*n-Floor[m*n/3]; Flatten[Table[T[n-k+1,k],{n,12},{k,n}]] (* Stefano Spezia, Oct 25 2022 *)

R(m,n) = m*n - floor(m*n/3).

A338432 Triangle read by rows: T(n, k) = (n - k + 1)^2 + 2*k^2, for n >= 1, and k = 1, 2, ..., n.

Original entry on oeis.org

3, 6, 9, 11, 12, 19, 18, 17, 22, 33, 27, 24, 27, 36, 51, 38, 33, 34, 41, 54, 73, 51, 44, 43, 48, 59, 76, 99, 66, 57, 54, 57, 66, 81, 102, 129, 83, 72, 67, 68, 75, 88, 107, 132, 163, 102, 89, 82, 81, 86, 97, 114, 137, 166, 201

Offset: 1

Wolfdieter Lang, Dec 09 2020

This triangle is obtained from the array A(m, k) = m^2 + 2*k^2, for k and m >= 1, read by upwards antidiagonals. This array A is of interest for representing numbers as a sum of three non-vanishing squares with two squares coinciding.
For the numbers represented this way, see A154777. To find the actual values for m and k (taken positive), given a representable number from A154777, one can also use the number triangle T(n, k) = A(n-k+1, k).
To find the number of representations of value N (from A154777), it is sufficient to consider the rows n >= 1 not exceeding n_{max} = Floor(N, Min), where the sequence Min gives the minima of the numbers in each row: Min = {min(n)}_{n>=1} with min(n) = min(T(n, 1), T(n, 2), ..., T(n, n)) and Floor(N, Min) is the greatest member of Min not exceeding N.
Conjecture: min(n) = T(n, ceiling(n/3)), n >= 1. This is the sequence (n+1)^2 - ceiling(n/3)*(2*(n+1) - 3*ceiling(n/3)) = A071619(n+1) = ceiling((2/3)*(n+1)^2) = (n+1)^2 - floor((1/3)*(n+1)^2) = 3, 6, 11, 17, 24, 33, 43, .... (Proof of these identities by considering the three n (mod 3) cases.)
For the multiplicities of the representable values A154777(n), see A339047.
The author met this representation problem in connection with special triples of integer curvatures in the Descartes-Steiner five circle problem.

			The triangle T(n, k) begins:
n \ k  1   2   3   4   5   6   7   8   9  10  11  12 ...
1:     3
2:     6   9
3:    11  12  19
4:    18  17  22  33
5:    27  24  27  36  51
6:    38  33  34  41  54  73
7:    51  44  43  48  59  76  99
8:    66  57  54  57  66  81 102 129
9:    83  72  67  68  75  88 107 132 163
10:  102  89  82  81  86  97 114 137 166 201
11:  123 108  99  96  99 108 123 144 171 204 243
12:  146 129 118 113 114 121 134 153 178 209 246 289
...
----------------------------------------------------
T(5, 1) = 5^2 + 2*1^2 = 27 = T(5, 3) = 3^2 + 2*3^2. A338433(11) = 2 for A154777(11) = 27.
T(4, 4) = 1^2 + 2*4^2 = 33 = T(6, 2) = 5^2 + 2*2^2. A338433(12) = 2 for A154777(12) = 33.
T(5, 5) = 1^2 + 2*5^2 = 51 = T(7, 1) = 7^2 + 2*1^2. A338433(20) = 2 for A154777(20) = 51.
T(7, 7) = 1^1 - 2*7^2 = 99 = T(11, 3) = 9^2 + 2*3^2 = 99 = T(11, 5) = 7^2 + 2*5^2. A338433(39) = 3 for A154777(39) = 99.
The first multiplicity 4 appears for 297.

Cf. A154777, A339047.
Cf. Columns k = 1..3: A059100, A189833, A241848.
Cf. Diagonals m = 1..4: A058331, A255843, A339048, A255847.

T(n, k) = A(n - k + 1, k), with the array A(m, k) = m^2 + 2*k^2, for n >= 1 and k = 1, 2, ..., n, and 0 otherwise.
G.f. of T and A column k (offset 0): G(k, x) = (1 + x + 2*(1 - x)^2*k^2)/(1-x)^3, for k >= 1.
G.f. of T diagonal m (A row m) (offset 0): D(m, x) = ((2*(1+x) + (1-x)^2*m^2)/(1-x)^3), for m >= 1.
G.f. of row polynomials in x (that is, g.f. of the triangle): G(z,x) = (3 - 3*z + (2 - 6*x + x^2)*z^2 + (2 + x)*x*z^3)*x*z / ((1 - z)*(1 - x*z))^3.

A358159 a(n) is the permanent of the n X n matrix M(n) that is defined by M[i,j] = ij - floor(ij/3).

Original entry on oeis.org

1, 1, 7, 102, 4396, 374216, 49857920, 11344877568, 3879729283968, 1804571320405248, 1195546731955854336, 1058730877124859138048, 1184751018265831288602624, 1725335046543668616765112320, 3147123030650561978295975936000, 6934187745940804400441946931200000, 18840570649600136750602236509552640000

Offset: 0

Stefano Spezia, Nov 01 2022

nonn

The matrix M(n) is the n-th principal submatrix of the rectangular array A143976 and it is singular for n > 3.

			a(5) = 374216:
    1   2   2   3   4
    2   3   4   6   7
    2   4   6   8  10
    3   6   8  11  14
    4   7  10  14  17

Cf. A143976.

Cf. A071619 (matrix element M[n,n]), A358042 (trace of M(n)), A358160 (hafnian of M(2*n)).

Join[{1},Table[Permanent[Table[i*j-Floor[i*j/3],{i,n},{j,n}]],{n,17}]]

from sympy import Matrix
def A358159(n): return Matrix(n,n,[i*j-i*j//3 for i in range(1,n+1) for j in range(1,n+1)]).per() if n else 1 # Chai Wah Wu, Nov 02 2022

A147623 The 3rd Witt transform of A040000.

Original entry on oeis.org

0, 2, 6, 12, 22, 34, 48, 66, 86, 108, 134, 162, 192, 226, 262, 300, 342, 386, 432, 482, 534, 588, 646, 706, 768, 834, 902, 972, 1046, 1122, 1200, 1282, 1366, 1452, 1542, 1634, 1728, 1826, 1926, 2028, 2134, 2242, 2352, 2466, 2582, 2700, 2822, 2946, 3072

Offset: 0

R. J. Mathar, Nov 08 2008

The 2nd Witt transform of A040000 is represented by A042964.

Vincenzo Librandi, Table of n, a(n) for n = 0..1000
Pieter Moree, The formal series Witt transform, Discr. Math. no. 295 vol. 1-3 (2005) 143-160.
Index entries for linear recurrences with constant coefficients, signature (2,-1,1,-2,1).

Cf. A040000, A042964, A049347, A071619.

[n le 2 select 1+(-1)^n else 4*(1+(n-2)^2) - Self(n-1) - Self(n-2): n in [1..30]]; // G. C. Greubel, Oct 24 2022

CoefficientList[Series[2x(1+x)(1 +x^2)/((1-x)^3 (1+x+x^2)), {x,0,40}], x] (* Vincenzo Librandi, Dec 14 2012 *)
LinearRecurrence[{2,-1,1,-2,1},{0,2,6,12,22},50] (* Harvey P. Dale, Jul 04 2021 *)

[2*(2*(1+3*n^2) -(2*chebyshev_U(n, -1/2) +chebyshev_U(n-1, -1/2)))/9 for n in range(41)] # G. C. Greubel, Oct 24 2022

G.f.: 2*x*(1+x)*(1+x^2)/((1-x)^3*(1+x+x^2)).
a(n) = 2*A071619(n).
From G. C. Greubel, Oct 24 2022: (Start)
a(n) = 4*(2 - 2*n + n^2) - a(n-1) - a(n-2).
a(n) = 2*(2*(1 + 3*n^2) - (2*A049347(n) + A049347(n-1)))/9. (End)

cp's OEIS Frontend

A358042 Partial sums of A071619.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

Formula

A005893 Number of points on surface of tetrahedron; coordination sequence for sodalite net (equals 2*n^2+2 for n > 0).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A042968 Numbers not divisible by 4.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Sage

Formula

Extensions

A156040 Number of compositions (ordered partitions) of n into 3 parts (some of which may be zero), where the first is at least as great as each of the others.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

Extensions

A008747 Expansion of (1+x^4)/((1-x)*(1-x^2)*(1-x^3)).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

GAP

Magma

Maple

Mathematica

PARI

Sage

Formula

A071318 Lesser of 2 consecutive numbers which are cubefree and not squarefree, i.e., numbers k such that both k and k+1 are in A067259.

Views

Author

Keywords

Comments

Examples

Links

A008747 Expansion of (1+x^4)/((1-x)(1-x^2)(1-x^3)).

A358159 a(n) is the permanent of the n X n matrix M(n) that is defined by M[i,j] = ij - floor(ij/3).