Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Coefficients in expansion of product(1-x^a, a=1n)?

  1. Jan 22, 2010 #1
    I have an equation of the form [tex](1-x)(1-x^2)(1-x^3)\cdot\cdot\cdot(1-x^n)[/tex]. Or, in maple notation, product(1-x^a, a=1..n).

    I've been trying to find a way to express this as an expanded polynomial in sigma notation so that I can extract coefficients (for an enumerative generating function). If the exponents were all the same then the binomial theorem could be applied, but unfortunately it is not that simple.

    I have made little progress towards successfully expanding. I tried reforming it as a different counting problem and solving it indirectly, but still came to this same fundamental problem, just in a different form. However, in doing so, I think I've found a connection to a sum of sums; the answer may involve summing over i=0..n and j=i..n?

    Has anyone seen this problem before? Any helpful links?

    Thank you very much in advance,

    (P.S. This is not homework. This is extracurricular fun to find a more efficient way to solve a computer science problem I stumbled across a while ago.)
  2. jcsd
  3. Jan 22, 2010 #2
    Notice that in the expansion, the coefficient of [tex]x^k, 1\leq k\leq n[/tex], is the number of unordered integer partitions of [tex]k[/tex].
  4. Jan 22, 2010 #3
    That observation is related to the problem I was trying to solve! I realize now that a subproblem of what I'm doing is counting the number of ways to partition a natural number into any number of unordered positive parts. This equation was a piece of the resulting generating function (the part I had trouble expanding), and thus that's what the coefficients should be counting.

    You gave me the right words to look up. Found this page which looks like what I need:
    http://en.wikipedia.org/wiki/Partition_(number_theory [Broken])

    Thank you very much.

    However, there are two problems:

    (1) Maple seems to disagree with your post. Typing into Maple: sort(collect(expand(product(1-x^a,a=1..20)),x),x,ascending), does not appear to count partitions of [tex]k[/tex]. Edit: According to Wikipedia, [tex]\frac{1}{1-x^a}[/tex] will. Is there a way I can still use this fact?

    (2) There does not seem to be a usable expression for [tex]p(n)[/tex]. But, if I feel confident that my expansion equals a summation over [tex]p(n)x^n[/tex], then I can at least use that to algebraically solve for a final answer, and then implement [tex]p(n)[/tex] into my program at the end.
    Last edited by a moderator: May 4, 2017
  5. Jan 22, 2010 #4
    You are right, I forgot to include "without repetition" in my earlier post. The product you refer gives the number of unrestricted partitions, in which the same numerical factor may apper several times; in your product, it cannot.

    I am not sure, of the top of my head, if there is a relation between the two (like the classical one that states that the number of odd partitions equals the distinct factors ones), but I'll check it.

    As far as I know, there isn't simple closed expression for that product (as there is not one to the p(n) numbers); it's one of the more classical problems in Combinatorics. Maybe Macmahon polynomials will help but, again, I'll have to check it; I haven't looked at these problems lately.
  6. Jan 23, 2010 #5
    I just found an excellent resource.


    It covers both of the functions you mentioned, [tex]p(n)[/tex] (unordered partitions of n with repetition) and [tex]q(n)[/tex] (unordered partitions of n without repetition) concisely, in the perfect amount of detail, as well as provides the formula that connects the two. And it turns out that I will need both.

    There seems to be only one problem remaining since my product goes to [tex]n[/tex] rather than infinity. When [tex]k > n[/tex] (which it usually will be), the coefficient of [tex]x^k[/tex] is no longer [tex]p(n)[/tex] or [tex]q(n)[/tex]. What do I do then? It seems like the answer will either be very simple, or next to unreasonable.

    I have tried separating the infinite product to go from 1 .. n, and then n+1 .. infinity. Then I re-indexed the second product to make it look like what we already know. But that did not work. Any ideas? If it's a dead end, I can try to reformulate the problem yet again to have an infinite upper bound. Sigh, this is how I'm spending a Friday night.

    Edit: The original equation I posted is apparently the Euler function (http://en.wikipedia.org/wiki/Euler's_function), which seems to have extractable coefficients by the Pentagonal Number Theorem (http://en.wikipedia.org/wiki/Pentagonal_number_theorem)... if the upper bound was infinity.
    Last edited by a moderator: Apr 24, 2017
  7. Jan 24, 2010 #6
    I stumbled across an interesting connection between partitioning [tex]n[/tex] into [tex]k[/tex] parts from any finite subset of the natural numbers and the partition functions [tex]p(n)[/tex] and [tex]q(n)[/tex].

    To be precise, if [tex]n,k[/tex] are non-negative integers, [tex]S \subset \mathbb{N}[/tex], and [tex]S[/tex] has a maximum:

    Let [tex]p_{S-k}(n)[/tex] be the number of ways to partition [tex]n[/tex] into [tex]k[/tex] unordered parts, where each part is restricted to the elements in set [tex]S[/tex], with repetition.

    Let [tex]w = k + n(max(S)k+1)[/tex], [tex]P'(x) = 1 + x(max(S)k+1)[/tex], and [tex]Q'(x) = (k+1)P'(x)[/tex].

    Let [tex]p_S(n)[/tex] be the number of ways to partition [tex]n[/tex] into unordered parts, where each part is restricted to the elements in set [tex]S[/tex] under the mapping [tex]P'[/tex], with repetition.

    Let [tex]q_{S-even}(n)[/tex] be the number of ways to partition [tex]n[/tex] into an even number of unordered parts, where each part is restricted to the elements in set [tex]S[/tex] under the mapping [tex]Q'[/tex], without repetition.

    Let [tex]q_{S-odd}(n)[/tex] be the number of ways to partition [tex]n[/tex] into an odd number of unordered parts, where each part is restricted to the elements in set [tex]S[/tex] under the mapping [tex]Q'[/tex], without repetition.

    Then, [tex]p_{S-k}(n) = \sum_{i=0}^{w} p_S(i)(q_{S-even}(w-i) - q_{S-odd}(w-i))[/tex]​

    I have not found this result to be computationally useful for my initial problem. Nevertheless, I found it interesting. My proof is weak, but I have verified several cases computationally. Is this kind of result worth delving into to show similar connections and publishing? Or is it too trivial, obvious, or useless? Sorry if that question sounds naive.
    Last edited: Jan 24, 2010
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook

Similar Discussions: Coefficients in expansion of product(1-x^a, a=1n)?
  1. A/n = a^(1/n) (Replies: 7)