Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A006171
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
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).

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | The OEIS Foundation | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified February 9 11:24 EST 2010. Contains 172296 sequences.


AT&T Labs Research