A000980 -id:A000980 - cp's OEIS frontend

A004171 a(n) = 2^(2n+1).

2, 8, 32, 128, 512, 2048, 8192, 32768, 131072, 524288, 2097152, 8388608, 33554432, 134217728, 536870912, 2147483648, 8589934592, 34359738368, 137438953472, 549755813888, 2199023255552, 8796093022208, 35184372088832, 140737488355328, 562949953421312

Offset: 0

N. J. A. Sloane

Same as Pisot sequences E(2, 8), L(2, 8), P(2, 8), T(2, 8). See A008776 for definitions of Pisot sequences.
In the Chebyshev polynomial of degree 2n, a(n) is the coefficient of x^2n. - Benoit Cloitre, Mar 13 2002
1/2 - 1/8 + 1/32 - 1/128 + ... = 2/5. - Gary W. Adamson, Mar 03 2009
From Adi Dani, May 15 2011: (Start)
Number of ways of placing an even number of indistinguishable objects in n + 1 distinguishable boxes with at most 3 objects in box.
Number of compositions of even natural numbers into n + 1 parts less than or equal to 3 (0 is counted as part). (End)
Also the number of maximal cliques in the (n+1)-Sierpinski tetrahedron graph for n > 0. - Eric W. Weisstein, Dec 01 2017
Assuming the Collatz conjecture is true, any starting number eventually leads to a power of 2. A number in this sequence can never be the first power of 2 in a Collatz sequence except of course for the Collatz sequence starting with that number. For example, except for 8, 4, 2, 1, any Collatz sequence that includes 8 must also include 16 (e.g., 5, 16, 8, 4, 2, 1). - Alonso del Arte, Oct 01 2019
First differences of A020988, and thus the "wavelengths" of the local maxima in A020986. See the Brillhart and Morton link, pp. 855-856. - John Keith, Mar 04 2021

			G.f. = 2 + 8*x + 32*x^2 + 128*x^3 + 512*x^4 + 2048*x^5 + 8192*x^6 + 32768*x^7 + ...
From _Adi Dani_, May 15 2011: (Start)
a(1) = 8 because all compositions of even natural numbers into 2 parts less than or equal to 3 are:
  for 0: (0, 0)
  for 2: (0, 2), (2, 0), (1, 1)
  for 4: (1, 3), (3, 1), (2, 2)
  for 6: (3, 3).
a(2) = 32 because all compositions of even natural numbers into 3 parts less than or equal to 3 are:
  for 0: (0, 0, 0)
  for 2: (0, 0, 2), (0, 2, 0), (2, 0, 0), (0, 1, 1), (1, 0, 1) , (1, 1, 0)
  for 4: (0, 1, 3), (0, 3, 1), (1, 0, 3), (1, 3, 0), (3, 0, 1), (3, 1, 0), (0, 2, 2), (2, 0, 2), (2, 2, 0), (1, 1, 2), (1, 2, 1), (2, 1, 1)
  for 6: (0, 3, 3), (3, 0, 3), (3, 3, 0), (1, 2, 3), (1, 3, 2), (2, 1, 3), (2, 3, 1), (3, 1, 2), (3, 2, 1), (2, 2, 2)
  for 8: (2, 3, 3), (3, 2, 3), (3, 3, 2).
(End)

Adi Dani, Quasicompositions of natural numbers, Proceedings of III congress of mathematicians of Macedonia, Struga Macedonia 29 IX -2 X 2005 pages 225-238.

Vincenzo Librandi, Table of n, a(n) for n = 0..200
J. Brillhart and P. Morton, A case study in mathematical research: the Golay-Rudin-Shapiro sequence, Amer. Math. Monthly, 103 (1996) 854-869.
Ling Gao, Graph assembly for spider and tadpole graphs, Master's Thesis, Cal. State Poly. Univ. (2023).
Milan Janjic, Enumerative Formulas for Some Functions on Finite Sets.
Tanya Khovanova, Recursive Sequences.
Mitchell Paukner, Lucy Pepin, Manda Riehl, and Jarred Wieser, Pattern Avoidance in Task-Precedence Posets, arXiv:1511.00080 [math.CO], 2015-2016.
Eric Weisstein's World of Mathematics, Maximal Clique.
Eric Weisstein's World of Mathematics, Sierpinski Tetrahedron Graph.
Index entries for linear recurrences with constant coefficients, signature (4).
Index to divisibility sequences.

Cf. A013708-A013729.
Absolute value of A009117. Essentially the same as A081294.
Cf. A132020, A164632. Equals A000980(n) + 2*A181765(n). Cf. A013776.

List([0..30],n->2^(2*n+1)); # Muniru A Asiru, Mar 12 2019

a004171 = (* 2) . a000302
a004171_list = iterate (* 4) 2  -- Reinhard Zumkeller, Jan 09 2013

[2^(2*n+1): n in [0..30]]; // Vincenzo Librandi, May 16 2011

seq(2^(2*n+1),n=0..24); # Nathaniel Johnston, Jun 25 2011

Table[2^(2 n + 1), {n, 0, 24}]
2^(2 Range[20] - 1) (* Eric W. Weisstein, Dec 01 2017 *)
LinearRecurrence[{4}, {2}, 20] (* Eric W. Weisstein, Dec 01 2017 *)
CoefficientList[Series[2/(1 - 4 x), {x, 0, 20}], x] (* Eric W. Weisstein, Dec 01 2017 *)

a(n)=2<<(2*n) \\ Charles R Greathouse IV, Apr 07 2012

a(n) = 2^(2*n+1) \\ Michel Marcus, Aug 12 2014

[2**(2*n+1) for n in range(0,25)] # Stefano Spezia, Jul 23 2025

((List.fill(20)(4: BigInt)).scanLeft(1: BigInt)( * )).map(2 * ) // _Alonso del Arte, Sep 12 2019

a(n) = 2*4^n.
a(n) = 4*a(n-1).
1 = 1/2 + Sum_{n >= 1} 3/a(n) = 3/6 + 3/8 + 3/32 + 3/128 + 3/512 + 3/2048 + ...; with partial sums: 1/2, 31/32, 127/128, 511/512, 2047/2048, ... - Gary W. Adamson, Jun 16 2003
From Philippe Deléham, Nov 23 2008: (Start)
a(n) = 2*A000302(n).
G.f.: 2/(1-4*x). (End)
a(n) = A081294(n+1) = A028403(n+1) - A000079(n+1) for n >= 1. a(n-1) = A028403(n) - A000079(n). - Jaroslav Krizek, Jul 27 2009
E.g.f.: 2*exp(4*x). - Ilya Gutkovskiy, Nov 01 2016
a(n) = A002063(n)/3 - A000302(n). - Zhandos Mambetaliyev, Nov 19 2016
a(n) = Sum_{k = 0..2*n} (-1)^(k+n)*binomial(4*n + 2, 2*k + 1); a(2*n) = Sum_{k = 0..2*n} binomial(4*n + 2, 2*k + 1) = A013776(n). - Peter Bala, Nov 25 2016
Product_{n>=0} (1 - 1/a(n)) = A132020. - Amiram Eldar, May 08 2023

A053632 Irregular triangle read by rows giving coefficients in expansion of Product_{k=1..n} (1 + x^k).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1, 2, 2, 3, 3, 3, 3, 3, 3, 2, 2, 1, 1, 1, 1, 1, 1, 2, 2, 3, 4, 4, 4, 5, 5, 5, 5, 4, 4, 4, 3, 2, 2, 1, 1, 1, 1, 1, 1, 2, 2, 3, 4, 5, 5, 6, 7, 7, 8, 8, 8, 8, 8, 7, 7, 6, 5, 5, 4, 3, 2, 2, 1, 1, 1, 1, 1, 1, 2, 2, 3, 4

Offset: 0

N. J. A. Sloane, Mar 22 2000

Or, triangle T(n,k) read by rows, giving number of subsets of {1,2,...,n} with sum k. - Roger CUCULIERE (cuculier(AT)imaginet.fr), Nov 19 2000
Row n consists of A000124(n) terms. These are also the successive vectors (their nonzero elements) when one starts with the infinite vector (of zeros) with 1 inserted somewhere and then shifts it one step (right or left) and adds to the original, then shifts the result two steps and adds, three steps and adds, etc. - Antti Karttunen, Feb 13 2002
T(n,k) = number of partitions of k into distinct parts <= n. Triangle of distribution of Wilcoxon's signed rank statistic. - Mitch Harris, Mar 23 2006
T(n,k) = number of binary words of length n in which the sum of the positions of the 0's is k. Example: T(4,5)=2 because we have 0110 (sum of the positions of the 0's is 1+4=5) and 1001 (sum of the positions of the 0's is 2+3=5). - Emeric Deutsch, Jul 23 2006
A fair coin is flipped n times. You receive i dollars for a "success" on the i-th flip, 1<=i<=n. T(n,k)/2^n is the probability that you will receive exactly k dollars. Your expectation is n(n+1)/4 dollars. - Geoffrey Critzer, May 16 2010
From Gus Wiseman, Jan 02 2023: (Start)
With offset 1, also the number of integer compositions of n whose partial sums add up to k for k = n..n(n+1)/2. For example, row n = 6 counts the following compositions:
  6 15 24 33  42  51  141  231  321  411  1311 2211  3111  12111 21111 111111
          114 123 132 222  312  1131 1221 2121 11121 11211
                  213 1113 1122 1212 2112 1111
(End)

			Triangle begins:
  1;
  1, 1;
  1, 1, 1, 1;
  1, 1, 1, 2, 1, 1, 1;
  1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1;
  1, 1, 1, 2, 2, 3, 3, 3, 3, 3, 3, 2, 2, 1, 1, 1;
  1, 1, 1, 2, 2, 3, 4, 4, 4, 5, 5, 5, 5, 4, 4, 4, 3, 2, 2, 1, 1, 1;
  ...
Row n = 4 counts the following binary words, where k = sum of positions of zeros:
  1111  0111  1011  0011  0101  0110  0001  0010  0100  1000  0000
                    1101  1110  1001  1010  1100
Row n = 5 counts the following strict partitions of k with all parts <= n (0 is the empty partition):
  0  1  2  3  4  5  42  43  53  54  532  542  543  5431 5432 54321
           21 31 32 51  52  431 432 541  5321 5421
                 41 321 421 521 531 4321

A. V. Yurkin, New binomial and new view on light theory, (book), 2013, 78 pages, no publisher listed.

Alois P. Heinz, Rows n = 0..40, flattened
S. R. Finch, Signum equations and extremal coefficients.
Steven R. Finch, Signum equations and extremal coefficients, February 7, 2009. [Cached copy, with permission of the author]
FindStat - Combinatorial Statistic Finder, The major index of an integer composition
Alexander Rosa and Štefan Znám, A combinatorial problem in the theory of congruences. (Russian), Mat.-Fys. Casopis Sloven. Akad. Vied 15 1965 49-59. [Annotated scanned copy] See Table 1.
F. Wilcoxon, Individual Comparisons by Ranking Methods, Biometrics Bulletin, v. 1, no. 6 (1945), pp. 80-83.
A. V. Yurkin, On similarity of systems of geometrical and arithmetic triangles, in Mathematics, Computing, Education Conference XIX, 2012.
A. V. Yurkin, New view on the diffraction discovered by Grimaldi and Gaussian beams, arXiv preprint arXiv:1302.6287 [physics.optics], 2013.
A. V. Yurkin, About the evident description of distribution of beams and "wavy geometrical trajectories" in long thin pipes, 2014 (original in Russian).
A. V. Yurkin, Symmetric triangle of Pascal and non-linear arithmetic parallelepiped, Book Manuscript, Research Gate 2015.

Cf. A053633, A068009.
Rows reduced modulo 2 and interpreted as binary numbers: A068052, A068053. Rows converge towards A000009.
Row sums give A000079.
Cf. A001788, A028362.
Cf. A285101 (multiplicative encoding of each row), A285103 (number of odd terms on row n), A285105 (number of even terms).
Row lengths are A000124.
A reciprocal version is (A033999, A219977, A291983, A291984, A291985, ...).
A negative version is A231599.
A version for partitions is A358194, reversed partitions A264034.
Cf. A025591, A029931, A063865, A152947, A318283, A359042.

with(gfun,seriestolist); map(op,[seq(seriestolist(series(mul(1+(z^i), i=1..n),z,binomial(n+1,2)+1)), n=0..10)]); # Antti Karttunen, Feb 13 2002
# second Maple program:
g:= proc(n) g(n):= `if`(n=0, 1, expand(g(n-1)*(1+x^n))) end:
T:= n-> seq(coeff(g(n), x, k), k=0..degree(g(n))):
seq(T(n), n=0..10);  # Alois P. Heinz, Nov 19 2012

Table[CoefficientList[ Series[Product[(1 + t^i), {i, 1, n}], {t, 0, 100}], t], {n, 0, 8}] // Grid (* Geoffrey Critzer, May 16 2010 *)

From Mitch Harris, Mar 23 2006: (Start)
T(n,k) = T(n-1, k) + T(n-1, k-n), T(0,0)=1, T(0,k) = 0, T(n,k) = 0 if k < 0 or k > (n+1 choose 2).
G.f.: (1+x)*(1+x^2)*...*(1+x^n). (End)
Sum_{k>=0} k * T(n,k) = A001788(n). - Alois P. Heinz, Feb 09 2017
max_{k>=0} T(n,k) = A025591(n). - Alois P. Heinz, Jan 20 2023

A025591 Maximal coefficient of Product_{k<=n} (1 + x^k). Number of solutions to +- 1 +- 2 +- 3 +- ... +- n = 0 or 1.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 5, 8, 14, 23, 40, 70, 124, 221, 397, 722, 1314, 2410, 4441, 8220, 15272, 28460, 53222, 99820, 187692, 353743, 668273, 1265204, 2399784, 4559828, 8679280, 16547220, 31592878, 60400688, 115633260, 221653776, 425363952, 817175698

Offset: 0

David W. Wilson

If k is allowed to approach infinity, this gives the partition numbers A000009.

a(n) is the maximal number of subsets of {1,2,...,n} that share the same sum.

T. D. Noe, Alois P. Heinz and Ray Chandler, Table of n, a(n) for n = 0..3339 (terms < 10^1000, first 201 terms from T. D. Noe, next 200 terms from Alois P. Heinz)
Dorin Andrica and Ioan Tomescu, On an Integer Sequence Related to a Product of Trigonometric Functions, and Its Combinatorial Relevance , Journal of Integer Sequences, Vol. 5 (2002), Article 02.2.4.
Vlad-Florin Dragoi and Valeriu Beiu, Fast Reliability Ranking of Matchstick Minimal Networks, arXiv:1911.01153 [cs.DM], 2019.
Steven R. Finch, Signum equations and extremal coefficients, February 7, 2009. [Cached copy, with permission of the author]
Erich Friedman and Mike Keith, Magic Carpets, J. Int Sequences, 3 (2000), Article 00.2.5.
Susumu Kubo, Partially Ordered Sets Corresponding to the Partition Problem, arXiv:2405.05544 [cs.DM], 2024. See pp. 2, 16.
Marco Mondelli, S. Hamed Hassani, and Rüdiger Urbanke, Construction of Polar Codes with Sublinear Complexity, arXiv preprint arXiv:1612.05295 [cs.IT], 2016-2017. See Sect. I.
Robert A. Proctor, Solution of two difficult combinatorial problems with linear algebra, American Mathematical Monthly 89, 721-734.
Blair D. Sullivan, On a conjecture of Adrica and Tomescu, J. Int. Sequences 16 (2013), Article 13.3.1.

Cf. A039828, A063865, A069918, A063866, A063867, A083309, A083527, A086376.

Cf. A053632, A160235, A359319, A359320.

b:= proc(n, i) option remember; `if`(n>i*(i+1)/2, 0,
      `if`(i=0, 1, b(n+i, i-1)+b(abs(n-i), i-1)))
    end:
a:=n-> b(0, n)+b(1, n):
seq(a(n), n=0..40);  # Alois P. Heinz, Mar 10 2014

f[n_, s_] := f[n, s]=Which[n==0, If[s==0, 1, 0], Abs[s]>(n*(n+1))/2, 0, True, f[n-1, s-n]+f[n-1, s+n]]; Table[Which[Mod[n, 4]==0||Mod[n, 4]==3, f[n, 0], Mod[n, 4]==1||Mod[n, 4]==2, f[n, 1]], {n, 0, 40}]
(* Second program: *)
p = 1; Flatten[{1, Table[p = Expand[p*(1 + x^n)]; Max[CoefficientList[p, x]], {n, 1, 50}]}] (* Vaclav Kotesovec, May 04 2018 *)
b[n_, i_] := b[n, i] = If[n > i(i+1)/2, 0, If[i == 0, 1, b[n+i, i-1] + b[Abs[n-i], i-1]]];
a[n_] := b[0, n] + b[1, n]; a /@ Range[0, 40] (* Jean-François Alcover, Feb 17 2020, after Alois P. Heinz *)

a(n)=if(n<0,0,polcoeff(prod(k=1,n,1+x^k),n*(n+1)\4))

from collections import Counter
def A025591(n):
    c = {0:1,1:1}
    for i in range(2,n+1):
        d = Counter(c)
        for k in c:
            d[k+i] += c[k]
        c = d
    return max(c.values()) # Chai Wah Wu, Jan 31 2024

a(n) = A063865(n) + A063866(n).
a(n) ~ sqrt(6/Pi) * 2^n / n^(3/2) [conjectured by Andrica and Tomescu (2002) and proved by Sullivan (2013)]. - Vaclav Kotesovec, Mar 17 2020
More precise asymptotics: a(n) ~ sqrt(6/Pi) * 2^n / n^(3/2) * (1 - 6/(5*n) + 589/(560*n^2) - 39/(50*n^3) + ...). - Vaclav Kotesovec, Dec 30 2022
a(n) = max_{k>=0} A053632(n,k). - Alois P. Heinz, Jan 20 2023

A063865 Number of solutions to +- 1 +- 2 +- 3 +- ... +- n = 0.

Original entry on oeis.org

1, 0, 0, 2, 2, 0, 0, 8, 14, 0, 0, 70, 124, 0, 0, 722, 1314, 0, 0, 8220, 15272, 0, 0, 99820, 187692, 0, 0, 1265204, 2399784, 0, 0, 16547220, 31592878, 0, 0, 221653776, 425363952, 0, 0, 3025553180, 5830034720, 0, 0, 41931984034, 81072032060, 0, 0

Offset: 0

N. J. A. Sloane, suggested by J. H. Conway, Aug 27 2001

Number of sum partitions of the half of the n-th-triangular number by distinct numbers in the range 1 to n. Example: a(7)=8 since triangular(7)=28 and 14 = 2+3+4+5 = 1+3+4+6 = 1+2+5+6 = 3+5+6 = 7+1+2+4 = 7+3+4 = 7+2+5 = 7+1+6. - Hieronymus Fischer, Oct 20 2010
The asymptotic formula below was stated as a conjecture by Andrica & Tomescu in 2002 and proved by B. D. Sullivan in 2013. See his paper and H.-K. Hwang's review MR 2003j:05005 of the JIS paper. - Jonathan Sondow, Nov 11 2013
a(n) is the number of subsets of {1..n} whose sum is equal to the sum of their complement. See example below. - Gus Wiseman, Jul 04 2019

			From _Gus Wiseman_, Jul 04 2019: (Start)
For example, the a(0) = 1 through a(8) = 14 subsets (empty columns not shown) are:
  {}  {3}    {1,4}  {1,6,7}    {3,7,8}
      {1,2}  {2,3}  {2,5,7}    {4,6,8}
                    {3,4,7}    {5,6,7}
                    {3,5,6}    {1,2,7,8}
                    {1,2,4,7}  {1,3,6,8}
                    {1,2,5,6}  {1,4,5,8}
                    {1,3,4,6}  {1,4,6,7}
                    {2,3,4,5}  {2,3,5,8}
                               {2,3,6,7}
                               {2,4,5,7}
                               {3,4,5,6}
                               {1,2,3,4,8}
                               {1,2,3,5,7}
                               {1,2,4,5,6}
(End)

T. D. Noe, N. J. A. Sloane and Ray Chandler, Table of n, a(n) for n = 0..3339 (terms < 10^1000, first 101 terms from T. D. Noe, next 300 terms from N. J. A. Sloane)
Dorin Andrica and Ovidiu Bagdasar, On k-partitions of multisets with equal sums, The Ramanujan J. (2021) Vol. 55, 421-435.
Dorin Andrica, Ovidiu Bagdasar, and George Cătălin Ţurcaş, The Number of Partitions of a Set and Superelliptic Diophantine Equations, Disc. Math. and Applications, Springer, Cham (2020), 35-55.
D. Andrica and E. J. Ionascu, Variations on a result of Erdős and Surányi, INTEGERS 2013 slides.
D. Andrica and I. Tomescu, On an Integer Sequence Related to a Product of Trigonometric Functions, and Its Combinatorial Relevance, J. Integer Seq., 5 (2002), Article 02.2.4
Ovidiu Bagdasar and Dorin Andrica, New results and conjectures on 2-partitions of multisets, 2017 7th International Conference on Modeling, Simulation, and Applied Optimization (ICMSAO).
Steven R. Finch, Signum equations and extremal coefficients, February 7, 2009. [Cached copy, with permission of the author]
B. D. Sullivan, On a Conjecture of Andrica and Tomescu, J. Int. Sequences, 16 (2013), Article 13.3.1.
zbMATH, Review of Andrica and Tomescu

Cf. A000980, A025591, A058377, A063866, A063867, A113036, A113037, A141000.
"Decimations": A060468 = 2*A060005, A123117 = 2*A104456.
Analogous sequences for sums of squares and cubes are A158092, A158118, see also A019568. - Pietro Majer, Mar 15 2009
Cf. A053632, A059529, A326173, A326174, A326175.

M:=400; t1:=1; lprint(0,1); for n from 1 to M do t1:=expand(t1*(x^n+1/x^n)); lprint(n, coeff(t1,x,0)); od: # N. J. A. Sloane, Jul 07 2008

f[n_, s_] := f[n, s]=Which[n==0, If[s==0, 1, 0], Abs[s]>(n*(n+1))/2, 0, True, f[ n-1, s-n]+f[n-1, s+n]]; a[n_] := f[n, 0]
nmax = 50; d = {1}; a1 = {};
Do[
  i = Ceiling[Length[d]/2];
  AppendTo[a1, If[i > Length[d], 0, d[[i]]]];
  d = PadLeft[d, Length[d] + 2 n] + PadRight[d, Length[d] + 2 n];
  , {n, nmax}];
a1 (* Ray Chandler, Mar 13 2014 *)

a(n)=my(x='x); polcoeff(prod(k=1,n,x^k+x^-k)+O(x),0) \\ Charles R Greathouse IV, May 18 2015

a(n)=0^n+floor(prod(k=1,n,2^(n*k)+2^(-n*k)))%(2^n) \\ Tani Akinari, Mar 09 2016

Asymptotic formula: a(n) ~ sqrt(6/Pi)*n^(-3/2)*2^n for n = 0 or 3 (mod 4) as n approaches infinity.
a(n) = 0 unless n == 0 or 3 (mod 4).
a(n) = constant term in expansion of Product_{ k = 1..n } (x^k + 1/x^k). - N. J. A. Sloane, Jul 07 2008
If n = 0 or 3 (mod 4) then a(n) = coefficient of x^(n(n+1)/4) in Product_{k=1..n} (1+x^k). - D. Andrica and I. Tomescu.
a(n) = 2*A058377(n) for any n > 0. - Rémy Sigrist, Oct 11 2017

More terms from Dean Hickerson, Aug 28 2001

Corrected and edited by Steven Finch, Feb 01 2009

A047653 Constant term in expansion of (1/2) * Product_{k=-n..n} (1 + x^k).

Original entry on oeis.org

1, 2, 4, 10, 26, 76, 236, 760, 2522, 8556, 29504, 103130, 364548, 1300820, 4679472, 16952162, 61790442, 226451036, 833918840, 3084255128, 11451630044, 42669225172, 159497648600, 597950875256, 2247724108772, 8470205600640, 31991616634296, 121086752349064

Offset: 0

N. J. A. Sloane

nonn

Or, constant term in expansion of Product_{k=1..n} (x^k + 1/x^k)^2. - N. J. A. Sloane, Jul 09 2008
Or, maximal coefficient of the polynomial (1+x)^2 * (1+x^2)^2 *...* (1+x^n)^2.
a(n) = A000302(n) - A181765(n).
From Gus Wiseman, Apr 18 2023: (Start)
Also the number of subsets of {1..2n} that are empty or have mean n. The a(0) = 1 through a(3) = 10 subsets are:
  {}  {}   {}       {}
      {1}  {2}      {3}
           {1,3}    {1,5}
           {1,2,3}  {2,4}
                    {1,2,6}
                    {1,3,5}
                    {2,3,4}
                    {1,2,3,6}
                    {1,2,4,5}
                    {1,2,3,4,5}
Also the number of subsets of {-n..n} with no 0's but with sum 0. The a(0) = 1 through a(3) = 10 subsets are:
  {}  {}      {}           {}
      {-1,1}  {-1,1}       {-1,1}
              {-2,2}       {-2,2}
              {-2,-1,1,2}  {-3,3}
                           {-3,1,2}
                           {-2,-1,3}
                           {-2,-1,1,2}
                           {-3,-1,1,3}
                           {-3,-2,2,3}
                           {-3,-2,-1,1,2,3}
(End)

T. D. Noe, Alois P. Heinz and Ray Chandler, Table of n, a(n) for n = 0..1669 (terms < 10^1000, first 201 terms from T. D. Noe, next 200 terms from Alois P. Heinz)
Ovidiu Bagdasar and Dorin Andrica, New results and conjectures on 2-partitions of multisets, 2017 7th International Conference on Modeling, Simulation, and Applied Optimization (ICMSAO).
Dorin Andrica and Ovidiu Bagdasar, On k-partitions of multisets with equal sums, The Ramanujan J. (2021) Vol. 55, 421-435.
R. C. Entringer, Representation of m as Sum_{k=-n..n} epsilon_k k, Canad. Math. Bull., 11 (1968), 289-293.
Steven R. Finch, Signum equations and extremal coefficients, February 7, 2009. [Cached copy, with permission of the author]
R. P. Stanley, Weyl groups, the hard Lefschetz theorem and the Sperner property, SIAM J. Algebraic and Discrete Methods 1 (1980), 168-184.

Cf. A025591.
Cf. A053632; variant: A127728.
For median instead of mean we have A079309(n) + 1.
Odd bisection of A133406.
A000980 counts nonempty subsets of {1..2n-1} with mean n.
A007318 counts subsets by length, A327481 by mean.
Cf. A024718, A047997, A057552, A070925, A212352, A326512, A326513, A327475, A361866, A362046.

f:=n->coeff( expand( mul((x^k+1/x^k)^2,k=1..n) ),x,0);
# second Maple program:
b:= proc(n, i) option remember; `if`(n>i*(i+1)/2, 0,
      `if`(i=0, 1, 2*b(n, i-1)+b(n+i, i-1)+b(abs(n-i), i-1)))
    end:
a:=n-> b(0, n):
seq(a(n), n=0..40);  # Alois P. Heinz, Mar 10 2014

b[n_, i_] := b[n, i] = If[n>i*(i+1)/2, 0, If[i == 0, 1, 2*b[n, i-1]+b[n+i, i-1]+b[Abs[n-i], i-1]]]; a[n_] := b[0, n]; Table[a[n], {n, 0, 40}] (* Jean-François Alcover, Mar 10 2014, after Alois P. Heinz *)
nmax = 26; d = {1}; a1 = {};
Do[
  i = Ceiling[Length[d]/2];
  AppendTo[a1, If[i > Length[d], 0, d[[i]]]];
  d = PadLeft[d, Length[d] + 2 n] + PadRight[d, Length[d] + 2 n] +
    2 PadLeft[PadRight[d, Length[d] + n], Length[d] + 2 n];
, {n, nmax}];
a1 (* Ray Chandler, Mar 15 2014 *)
Table[Length[Select[Subsets[Range[2n]],Length[#]==0||Mean[#]==n&]],{n,0,6}] (* Gus Wiseman, Apr 18 2023 *)

PARI
```
a(n)=polcoeff(prod(k=-n,n,1+x^k),0)/2
```

{a(n)=sum(k=0,n*(n+1)/2,polcoeff(prod(m=1,n,1+x^m+x*O(x^k)),k)^2)} \\ Paul D. Hanna, Nov 30 2010

Sum of squares of coefficients in Product_{k=1..n} (1+x^k):
a(n) = Sum_{k=0..n(n+1)/2} A053632(n,k)^2. - Paul D. Hanna, Nov 30 2010
a(n) = A000980(n)/2.
a(n) ~ sqrt(3) * 4^n / (sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Sep 11 2014
From Gus Wiseman, Apr 18 2023 (Start)
a(n) = A133406(2n+1).
a(n) = A212352(n) + 1.
a(n) = A362046(2n) + 1.
(End)

More terms from Michael Somos, Jun 10 2000

A007576 Number of solutions to k_1 + 2k_2 + ... + nk_n = 0, where k_i are from {-1,0,1}, i=1..n.

Original entry on oeis.org

1, 1, 1, 3, 7, 15, 35, 87, 217, 547, 1417, 3735, 9911, 26513, 71581, 194681, 532481, 1464029, 4045117, 11225159, 31268577, 87404465, 245101771, 689323849, 1943817227, 5494808425, 15568077235, 44200775239, 125739619467

Offset: 0

Simon Plouffe, Robert G. Wilson v and Vladeta Jovovic

Also, number of maximally stable towers of 2 X 2 LEGO blocks.

			For n=4 there are 7 solutions: (-1,-1,1,0), (-1,0,-1,1), (-1,1,1,-1), (0,0,0,0), (1,-1,-1,1), (1,0,1,-1), (1,1,-1,0).

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
P. J. S. Watson, On "LEGO" towers, J. Rec. Math., 12 (No. 1, 1979-1980), 24-27.

Ray Chandler, Table of n, a(n) for n = 0..2106 (terms < 10^1000; first 101 terms from T. D. Noe)
D. Andrica and O. Bagdasar, Some remarks on 3-partitions of multisets, Electron. Notes Discrete Math., TCDM'18 (2018).
Dorin Andrica and Ovidiu Bagdasar, On k-partitions of multisets with equal sums, The Ramanujan J. (2021) Vol. 55, 421-435.
Andrei Asinowski, Axel Bacher, Cyril Banderier, and Bernhard Gittenberger, Analytic combinatorics of lattice paths with forbidden patterns, the vectorial kernel method, and generating functions for pushdown automata, Laboratoire d'Informatique de Paris Nord (LIPN 2019).
Steven R. Finch, Signum equations and extremal coefficients [Broken link]
Steven R. Finch, Signum equations and extremal coefficients, February 7, 2009. [Cached copy, with permission of the author]
P. J. S. Watson, On "LEGO" towers, J. Rec. Math., 12 (No. 1, 1979-1980), 24-27. (Annotated scanned copy)
Index entry for sequences related to LEGO blocks

Cf. A007575, A063865, A039826.

f[0] = 1; f[n_] := Coefficient[Expand@ Product[1 + x^k + x^(2k), {k, n}], x^(n(n + 1)/2)]; Table[f@n, {n, 0, 28}] (* Robert G. Wilson v, Nov 10 2006 *)

a(n):=coeff(expand(product(1+x^k+x^(2*k),k,1,n)),x,binomial(n+1,2));
makelist(a(n),n,0,24);

Coefficient of x^(n*(n+1)/2) in Product_{k=1..n} (1+x^k+x^(2*k)).
Equivalently, the coefficient of x^0 in Product_{k=1..n} (1/x^k + 1 + x^k). - Paul D. Hanna, Jul 10 2018
a(n) ~ 3^(n + 1) / (2 * sqrt(Pi) * n^(3/2)). - Vaclav Kotesovec, Jul 11 2018
a(n) = (1/(2*Pi))*Integral_{t=0..2*Pi} ( Product_{k=1..n} (1+2*cos(k*t)) ) dt. - Ovidiu Bagdasar, Aug 08 2018

More terms from David Wasserman, Mar 29 2005

Edited by N. J. A. Sloane, Nov 07 2006. This is a merging of two sequences which, thanks to the work of Søren Eilers, we now know are identical.

A086376 Maximal coefficient of the polynomial (1-x)(1-x^2)...*(1-x^n).

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 2, 2, 2, 3, 2, 4, 3, 3, 4, 6, 5, 6, 7, 8, 8, 10, 11, 16, 16, 18, 21, 28, 29, 34, 41, 50, 56, 66, 80, 100, 114, 131, 158, 196, 225, 263, 320, 388, 455, 532, 644, 786, 921, 1083, 1321, 1600, 1891, 2218, 2711, 3280, 3895, 4588, 5591, 6780, 8051, 9519, 11624

Offset: 0

Yuval Dekel (dekelyuval(AT)hotmail.com), Sep 07 2003

Alois P. Heinz, Table of n, a(n) for n = 0..1000
S. R. Finch, Signum equations and extremal coefficients.
Steven R. Finch, Signum equations and extremal coefficients, February 7, 2009. [Cached copy, with permission of the author]
J. W. Meijer and M. Nepveu, Euler's ship on the Pentagonal Sea, Acta Nova, Volume 4, No.1, December 2008. pp. 176-187. [From _Johannes W. Meijer_, Jun 21 2010]
E. M. Wright, A closer estimate for a restricted partition function, Q. J. Math. 15 (1964) 283-287.

Cf. A025591, A086781, A160089.

A086376 := proc(n)
        g := expand(mul( 1-x^k,k=1..n) );
        convert(PolynomialTools[CoefficientVector](g, x), list):
        max(%);
end proc:
seq(A086376(n),n=0..60) ; # R. J. Mathar, Jun 01 2011

b[0] = 1; b[n_] := b[n] = b[n-1]*(1-x^n) // Expand;
a[n_] := CoefficientList[b[n], x] // Max;
Table[a[n], {n, 0, 100}] (* Jean-François Alcover, Apr 13 2017 *)

a(n)=vecmax(Vec(prod(k=1,n,1-x^k)));
vector(100,n,a(n-1)) \\ Joerg Arndt, Jan 31 2024

More terms from Sascha Kurz, Sep 22 2003

a(0)=1 prepended by Alois P. Heinz, Apr 12 2017

A070925 Number of subsets of A = {1,2,...,n} that have the same center of gravity as A, i.e., (n+1)/2.

Original entry on oeis.org

1, 1, 3, 3, 7, 7, 19, 17, 51, 47, 151, 137, 471, 427, 1519, 1391, 5043, 4651, 17111, 15883, 59007, 55123, 206259, 193723, 729095, 688007, 2601639, 2465133, 9358943, 8899699, 33904323, 32342235, 123580883, 118215779, 452902071, 434314137, 1667837679, 1602935103

Offset: 1

Sharon Sela (sharonsela(AT)hotmail.com), May 20 2002

nonn

From Gus Wiseman, Apr 15 2023: (Start)
Also the number of nonempty subsets of {0..n} with mean n/2. The a(0) = 1 through a(5) = 7 subsets are:
  {0}  {0,1}  {1}      {0,3}      {2}          {0,5}
              {0,2}    {1,2}      {0,4}        {1,4}
              {0,1,2}  {0,1,2,3}  {1,3}        {2,3}
                                  {0,2,4}      {0,1,4,5}
                                  {1,2,3}      {0,2,3,5}
                                  {0,1,3,4}    {1,2,3,4}
                                  {0,1,2,3,4}  {0,1,2,3,4,5}
(End)

			Of the 32 (2^5) sets which can be constructed from the set A = {1,2,3,4,5} only the sets {3}, {2, 3, 4}, {2, 4}, {1, 2, 4, 5}, {1, 2, 3, 4, 5}, {1, 3, 5}, {1, 5} give an average of 3.

Fausto A. C. Cariboni, Table of n, a(n) for n = 1..40

The odd bisection is A000980(n) - 1 = 2*A047653(n) - 1.
For median instead of mean we have A100066, bisection A006134.
Including the empty set gives A222955.
The one-based version is A362046, even bisection A047653(n) - 1.
A007318 counts subsets by length.
A067538 counts partitions with integer mean, strict A102627.
A231147 counts subsets by median.
A327481 counts subsets by integer mean.
Cf. A024718, A057552, A079309, A133406, A326512, A326513, A326537, A327475, A349156, A359893.

Needs["DiscreteMath`Combinatorica`"]; f[n_] := Block[{s = Subsets[n], c = 0, k = 2}, While[k < 2^n + 1, If[ (Plus @@ s[[k]]) / Length[s[[k]]] == (n + 1)/2, c++ ]; k++ ]; c]; Table[ f[n], {n, 1, 20}]
(* second program *)
Table[Length[Select[Subsets[Range[0,n]],Mean[#]==n/2&]],{n,0,10}] (* Gus Wiseman, Apr 15 2023 *)

From Gus Wiseman, Apr 18 2023: (Start)
a(2n+1) = A000980(n) - 1.
a(n) = A222955(n) - 1.
a(n) = 2*A362046(n) + 1.
(End)

Edited by Robert G. Wilson v and John W. Layman, May 25 2002

a(34)-a(38) from Fausto A. C. Cariboni, Oct 08 2020

A160235 The maximal coefficient of (1+x)(1+x^4)(1+x^9)...(1+x^(n^2)).

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 3, 5, 6, 9, 14, 21, 32, 54, 87, 144, 230, 383, 671, 1158, 1981, 3408, 6246, 10925, 19463, 34624, 63941, 114954, 208429, 380130, 707194, 1298600, 2379842, 4398644, 8253618, 15303057, 28453809, 53091455, 100061278, 187446097

Offset: 0

Theodore Kolokolnikov, May 05 2009

nonn

Vaclav Kotesovec, Table of n, a(n) for n = 0..1000 (terms 0..200 from Seiichi Manyama)
Steven R. Finch, Signum equations and extremal coefficients, February 7, 2009. [Cached copy, with permission of the author]

Cf. A025591, A033461, A359319, A359320, A369728.

for N from 1 to 40 do
p := expand(product(1+x^(n^2), n=1..N)):
L:=convert(PolynomialTools[CoefficientVector](p, x), list):
mmax := max(op(map(abs, L)));
lprint(mmax):
end:

p = 1; Table[p = Expand[p*(1 + x^(n^2))]; Max[CoefficientList[p, x]], {n, 1, 50}] (* Vaclav Kotesovec, May 04 2018 *)
nmax = 100; poly = ConstantArray[0, nmax*(nmax+1)*(2*nmax+1)/6 + 1]; poly[[1]] = 1; poly[[2]] = 1; Do[Do[poly[[j + 1]] += poly[[j - k^2 + 1]], {j, k*(k+1)*(2*k+1)/6, k^2, -1}]; Print[k, " ", Max[poly]], {k, 2, nmax}]; (* Vaclav Kotesovec, Dec 30 2022 *)

An asymptotic formula is a(n) ~ sqrt(10/Pi) * n^(-5/2) * 2^n. See for example the reference by Finch.

More precise asymptotics: a(n) ~ sqrt(10/Pi) * 2^n / n^(5/2) * (1 - 35/(18*n) + ...). - Vaclav Kotesovec, Dec 30 2022

a(0)=1 prepended by Seiichi Manyama, Dec 26 2022

A060005 Number of ways of partitioning the integers {1,2,..,4n} into two (unordered) sets such that the sums of parts are equal in both sets (parts in either set will add up to (4n)*(4n+1)/4). Number of solutions to {1 +- 2 +- 3 +- ... +- 4n=0}.

Original entry on oeis.org

1, 1, 7, 62, 657, 7636, 93846, 1199892, 15796439, 212681976, 2915017360, 40536016030, 570497115729, 8110661588734, 116307527411482, 1680341334827514, 24435006625667338, 357366669614512168, 5253165510907071170

Offset: 0

Roland Bacher, Mar 15 2001

nonn

			a(1)=1 since there is only one way of partitioning {1,2,3,4} into two sets of equal sum, namely {1,4}, {2,3}.

Alois P. Heinz and Ray Chandler, Table of n, a(n) for n = 0..835 (terms < 10^1000, first 251 terms from Alois P. Heinz)
Steven R. Finch, Signum equations and extremal coefficients, February 7, 2009. [Cached copy, with permission of the author]
L. Sallows, M. Gardner, R. K. Guy and D. E. Knuth, Serial isogons of 90 degrees, Math. Mag. 64 (1991), 315-324.

Cf. A060468, A007219, A107350, a(n)=A058377(4n), A227850.

b:= proc(n, i) option remember; local m; m:= i*(i+1)/2;
      `if`(n>m, 0, `if`(n=m, 1, b(abs(n-i), i-1) +b(n+i, i-1)))
    end:
a:= n-> b(4*n, 4*n-1):
seq(a(n), n=0..30);  # Alois P. Heinz, Oct 30 2011

b[n_, i_] := b[n, i] = Module[{m = i*(i+1)/2}, If[n > m, 0, If[n == m, 1, b[Abs[n-i], i-1] + b[n+i, i-1]]]]; a[n_] := b[4*n, 4*n-1]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Sep 26 2013, translated from Alois P. Heinz's Maple program *)

a(0)=1 and a(n) is half the coefficient of q^0 in product((q^(-k)+q^k), k=1..4*n) for n >= 1.

For n>=1, a(n) = (1/Pi)*16^n*J(4n) where J(n) = integral(t=0, Pi/2, cos(t)cos(2t)...cos(nt)dt). - Benoit Cloitre, Sep 24 2006

More terms from Alois P. Heinz, Oct 30 2011

cp's OEIS Frontend

A004171 a(n) = 2^(2n+1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

PARI

SageMath

Scala

Formula

A053632 Irregular triangle read by rows giving coefficients in expansion of Product_{k=1..n} (1 + x^k).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A025591 Maximal coefficient of Product_{k<=n} (1 + x^k). Number of solutions to +- 1 +- 2 +- 3 +- ... +- n = 0 or 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

A063865 Number of solutions to +- 1 +- 2 +- 3 +- ... +- n = 0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

Extensions

A047653 Constant term in expansion of (1/2) * Product_{k=-n..n} (1 + x^k).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Formula

Extensions

A007576 Number of solutions to k_1 + 2*k_2 + ... + n*k_n = 0, where k_i are from {-1,0,1}, i=1..n.

Views

Author

Keywords

A007576 Number of solutions to k_1 + 2k_2 + ... + nk_n = 0, where k_i are from {-1,0,1}, i=1..n.

A086376 Maximal coefficient of the polynomial (1-x)(1-x^2)...*(1-x^n).

A160235 The maximal coefficient of (1+x)(1+x^4)(1+x^9)...(1+x^(n^2)).