|
Search: id:A006171
|
|
|
| A006171 |
|
Number of factorization patterns of polynomials of degree n over integers. (Formerly M2479)
|
|
+0 13
|
|
| 1, 1, 3, 5, 11, 17, 34, 52, 94, 145, 244, 370, 603, 899, 1410, 2087, 3186, 4650, 6959, 10040, 14750, 21077, 30479, 43120, 61574, 86308, 121785, 169336, 236475, 326201, 451402, 618135, 848209, 1153733, 1571063, 2123325, 2871419, 3857569, 5182999
(list; graph; listen)
|
|
|
OFFSET
|
0,3
|
|
|
COMMENT
|
Number of partitions of n where there are unlimited distinguishable but unlabeled objects of each size. E.g., in splitting 2 into two parts of size 1, we distinguish whether the same object is used for each part. Also number of factorization patterns over rationals, or many other UFDs (but not over real or complex numbers). - Frank Adams-Watters (FrankTAW(AT)Netscape.net), Jun 19 2006
Contribution from Gary W. Adamson (qntmpkt(AT)yahoo.com), Jun 16 2009: (Start)
Equals the "aerate and convolve" convergent of A000041 =
(1, 1, 2, 3, 5, 7, 11,...) * (1, 0, 1, 0, 2, 0, 3, 0, 5,...) * (1, 0, 0, 1, 0, 0, 2, 0, 0, 3,...)
(End)
|
|
REFERENCES
|
R. A. Hultquist, G. L. Mullen and H. Niederreiter, Association schemes and derived PBIB designs of prime power order, Ars. Combin., 25 (1988), 65-82.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
T. D. Noe, Table of n, a(n) for n=0..1000
N. J. A. Sloane, Transforms
|
|
FORMULA
|
Euler transform of tau(n), tau(n)=the number of divisors of n, cf. A000005. G.f.: Product_{k=1..infinity} (1 - x^k)^(-tau(k)). a(n)=1/n*Sum_{k=1..n} a(n-k)*b(k), n>1, a(0)=1, b(k)=Sum_{d|k} d*tau(d), cf. A060640. - Vladeta Jovovic (vladeta(AT)eunet.rs), Apr 21 2001
a(n) = sum_{<b(i)^k(i)> partition of n} product p(k(i)), where p(n) is the partition function A000041. E.g., for the partition [4,2^3,1^4], the product is p(1)*p(3)*p(4) = 1*3*5 = 15. - Frank Adams-Watters (FrankTAW(AT)Netscape.net), Jun 19 2006
G.f.: A(x) = exp( Sum_{n>=1} sigma(n)*x^n/(1-x^n)/n ). [From Paul D. Hanna (pauldhanna(AT)juno.com), Mar 28 2009]
|
|
EXAMPLE
|
For n=3 we have 3=(3*1)=(1*3)=(2*1)+(1*1)=(1*2)+(1*1)=(1*1)+(1*1)+(1*1) so a(3)=5.
|
|
MAPLE
|
with (numtheory): etr:= proc(p) local b; b:=proc(n) option remember; local d, j; if n=0 then 1 else add (add (d*p(d), d=divisors(j)) *b(n-j), j=1..n)/n fi end end: a:=etr (tau): seq (a(n), n=0..38); [From Alois P. Heinz (heinz(AT)hs-heilbronn.de), Sep 08 2008]
|
|
PROGRAM
|
(PARI) a(n)=if(n<0, 0, polcoeff(1/prod(k=1, n, (1-x^k+x*O(x^n))^numdiv(k)), n))
(PARI) 1/prod(j=1, N, eta(x^j)); Vec(%) - Joerg Arndt (arndt(AT)jjj.de), May 03 2008
(PARI) {a(n)=if(n==0, 1, polcoeff(exp(sum(m=1, n, sigma(m)*x^m/(1-x^m+x*O(x^n))/m)), n))} [From Paul D. Hanna (pauldhanna(AT)juno.com), Mar 28 2009]
|
|
CROSSREFS
|
Cf. A000005, A060640, A061255, A061256, A001970, A061257.
Cf. A006167-A006170.
Cf. A000041, A115621.
Sequence in context: A091610 A006170 A147071 this_sequence A060647 A155989 A125557
Adjacent sequences: A006168 A006169 A006170 this_sequence A006172 A006173 A006174
|
|
KEYWORD
|
nonn,nice
|
|
AUTHOR
|
N. J. A. Sloane (njas(AT)research.att.com).
|
|
|
Search completed in 0.002 seconds
|