OFFSET
0,4
COMMENTS
Pairwise sums of sequence A002865 (partitions in which the least part is at least 2).
Also number of partitions of n into parts with at most one 1. - Reinhard Zumkeller, Oct 25 2004
Also number of partitions of n into parts with at least half of the parts having size 1; equivalently (by duality) number of partitions of n where the large part is at least twice as big as the second largest part. - Franklin T. Adams-Watters, Jun 08 2005
Also number of 2-regular not necessarily connected graphs with loops allowed but no multiple edges. - Jason Kimberley, Jan 05 2011
LINKS
Kevin Beanland and Hung Viet Chu, On Schreier-type Sets, Partitions, and Compositions, arXiv:2311.01926 [math.CO], 2023.
P. Chinn and S. Heubach, Integer Sequences Related to Compositions without 2's, J. Integer Seqs., Vol. 6, 2003.
V. Jelinek, T. Mansour, and M. Shattuck, On multiple pattern avoiding set partitions, Advances in Applied Mathematics Volume 50, Issue 2, February 2013, Pages 292-326. - N. J. A. Sloane, Jan 01 2013
Jerome Kelleher and Barry O'Sullivan, Generating All Partitions: A Comparison Of Two Encodings, arXiv:0909.2331 [cs.DS], 2009-2014. [Peter Luschny, Oct 24 2010]
Krishna Menon and Anurag Singh, Pattern avoidance and dominating compositions, arXiv:2104.07274 [math.CO], 2021.
Mircea Merca, Fast algorithm for generating ascending compositions, arXiv:1903.10797 [math.CO], 2019.
FORMULA
G.f.: (1 - x^2)*Product_{m>=1} 1/(1 - x^m).
a(n) = p(n) - p(n-2) for n >= 2, where p(n) are the partition numbers (A000041); follows at once from the g.f. - Emeric Deutsch, Feb 18 2006
a(n) ~ exp(sqrt(2*n/3)*Pi)*Pi / (6*sqrt(2)*n^(3/2)) * (1 - (3*sqrt(3/2)/Pi + 25*Pi/(24*sqrt(6)))/sqrt(n) + (25/8 + 9/(2*Pi^2) + 817*Pi^2/6912)/n). - Vaclav Kotesovec, Nov 04 2016
MAPLE
with(combinat): a:=proc(n) if n=0 then 1 elif n=1 then 1 else numbpart(n)-numbpart(n-2) fi end: seq(a(n), n=0..49); # Emeric Deutsch, Feb 18 2006
MATHEMATICA
a[n_] = PartitionsP[n] - PartitionsP[n-2]; a /@ Range[0, 49] (* Jean-François Alcover, Jul 13 2011, after Emeric Deutsch *)
PROG
(PARI) a(n)=if(n<0, 0, polcoeff((1-x^2)/eta(x+x*O(x^n)), n))
(Magma) A41 := func<n|n ge 0 select NumberOfPartitions(n) else 0>;
[A41(n)-A41(n-2):n in [0..49]]; // Jason Kimberley, Jan 05 2011
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Benoit Cloitre, Dec 10 2002
STATUS
approved