Vector Partition Functions (Mattias Beck) Question(s) On page 35 of Beck's exposition on Vector Partition Functions (found here: http://math.sfsu.edu/beck/papers/vectorpar.slides.pdf), Beck introduces the notion of the quasi-polynomial
$$\phi _{\textbf{A}} ( \textbf{b}) = \text{const} \left( \frac{1}{(1- \textbf{z} ^{\textbf{c}_1})(1-\textbf{z}^{\textbf{c}_2}) (1- \textbf{z} ^{\textbf{c}_3}) (1 - \textbf{z}^{\textbf{c}_4}) \textbf{z} ^{ \textbf{b}}} \right)$$
Which he then expands into partial fractions, treating $z_1$ as a variable and the rest of $z_k$ as constants to attain the form
$$\text{const} \left( \frac{1}{(1- \textbf{z} ^{\textbf{c}_1})(1-\textbf{z}^{\textbf{c}_2}) (1- \textbf{z} ^{\textbf{c}_3}) (1 - \textbf{z}^{\textbf{c}_4}) \textbf{z} ^{ \textbf{b}}} \right) = \frac{1}{z_2^{b_2} \cdots z_m ^{b_m}} \left( \sum ^d _{k = 1} \frac{A_k(\textbf{z},b_1)}{1-\textbf{z}^{c_k}} + \sum ^{b_1} _{j=1} \frac{B_j(\textbf{z})}{z_1^j} \right)$$
where "$A_k$ and $B_j$ are polynomials in $z_1$, rational functions in $z_2, ..., z_m$ and exponential in $b_1$." First off, what does this even mean? How do you get such a form, and how would you find these polynomials for a given example?
Then, how does Beck reduce this to
$$\text{const} \frac{1}{z_2^{b_2} \cdots z_m ^{b_m}} \sum ^d _{k=1} \frac{A_k (0, z_2, ..., z_m, b_1)}{1-(0,z_2,...,z_m)^{c_k}} $$
???
And finally, how in the world is this applicable to my (very concrete) example at Finding Constant Term in Product of Series ?
I would greatly appreciate any help with understanding what is going on here and especially how I can apply it. I'm sorry, I don't really have "progress" to speak about because I'm honestly very lost. 
 A: General remark: the slides by Matthias Beck are not a tutorial, so it was quite challenging to extract what was intended from just the slides. Trying to do examples was my best shot at getting to understand them.
So I'll do another example that probably better shows what Beck is up to. This example is $2$-dimensional, so we'll see some $z_i$ being used formally in rational expressions rather than as indeterminate. Also, importantly, the set of vectors $\def\c{\mathbf c}\c_i$ for which the partition function is computed will avoid having some being rational multiples of others (as was inevitable in dimension$~1$), which causes the polynomials $1-z^{\c_i}$ to be mutually relatively prime; one can then do partial fraction decomposition while leaving these denominator factors intact.
Concretely I'll take the collection $\c_1=(1,0)$, $\c_2=(0,1)$ and $\c_3=(1,1)$, which is about the smallest interesting case possible (isomorphic to the $A_2$ system of positive roots, a case Beck mentions on a previous slide). So we have as matrix
$$ A=\pmatrix{1&0&1\\0&1&1}.$$
As target vector I'll first take a vector of the form $\def\b{\mathbf b}\b=(b_1,b_2)$ with $b_1=3$, and then see what happens when $b_1$ is varied (the $b_2$ is separated out in the formula, so it really does not matter and can be left to vary right from the start. For readability call $z_1=X$ (which is being used as indeterminate, whence the capital letter) and $z_2=y$ (which is used as parameter whence the lower case). 
Now $\phi_A(\b)$ is the coefficient of $X^3y^{b_2}$ in the expansion of $\frac1{(1-X)(1-y)(1-Xy)}$ as a formal power series in $X$. (Since we are not doing any series in $y$, it does not really make sense to mention $y^{b_2}$, and I only mention it because Beck does, but we'll soon get rid of it anyway). This question is translated into finding the constant term in the Laurent series expansion of
$$
  \frac1{(1-X)(1-y)(1-Xy)X^3y^{b_2}}.
$$
To do so, Beck asks for a partial fraction decomposition of the expression using denominators $1-z^{\c_i}$, so concretely $1-X$, $1-y$, and $1-Xy$, as well as $X$, $X^2$ and $X^3$. A first thing to note is that $1-y$ is just a scalar in the field $\def\Q{\Bbb Q}\Q(y)$ we are going to work over, so it makes no sense to include a term with denominator $1-y$, and I will just drop it (though Beck mentions nothing about such cases). Also I will write $1-yX$ for $1-Xy$ so my brain has less difficulty of considering $y$ as a scalar. The factor $\frac1{1-y}$ is going to be present as a factor in all terms of our result, so it simplifies matters to single it out explicitly. Grouping together for now the terms with powers of $X$ as denominator, we are looking for polynomials $a,b,C\in\Q(y)[X]$ such that
$$
  \frac1{(1-X)(1-y)(1-Xy)X^3y^{b_2}}=
  \frac1{(1-y)y^{b_2}}\left(\frac a{1-X}+\frac b{1-yX}+\frac C{X^3}\right),\tag1
$$
and with numerators of degree less than their denominators, the is $\deg_X(a)<1$, $\deg_X(b)<1$ and $\deg(C)<3$ (so $a$ and $b$ are just scalars in $\Q(y)$, which is why I chose lower case for them). We turn equation $(1)$ into a polynomial one:
$$
  1 = a(1-yX)X^3+b(1-X)X^3+C(1-X)(1-yX). \tag2
$$
The trick to avoid solving a linear system (over $\Q(y)$) is to use special evaluations for $a,b$, and some modular reduction for $C$. Concretely, setting $X:=1$ kills the last two factors and the remaining equation $1=a(1-y)$ gives $a=\frac1{1-y}$; similarly setting $X=\frac1y$ kills the outer two terms, and $1=b(1-y^{-1})y^{-3}$ solves to $b=\frac{y^3}{1-y^{-1}}=-\frac{y^4}{1-y}$. Finally $(2)$ clearly implies $C(1-X)(1-yX)\equiv1\pmod{X^3}$, which is most easily solved from the series expansion of $\frac1{(1-X)(1-yX)}$ to give (using the degree bound for $C$)
$$
  C=1+(1+y)X+(1+y+y^2)X^2.
$$
Now the term $\frac C{X^3}$ can be split as $\frac{1+y+y^2}X+\frac{1+y}{X^2}+\frac1{X^3}$. So to get back to Beck's terminology we have found (taking care to reinsert the factor $\frac1{1-y}$ that we had set aside): $A_1(X,y,3)=\frac1{(1-y)^2}$, $A_2$ is absent (or zero), $A_3(X,y,3)=\frac{-y^4}{(1-y)^2}$, as well as $B_1(X,y)=\frac{1+y+y^2}{1-y}$, $B_2(X,y)=\frac{1+y}{1-y}$, and $B_3(X,y)=\frac1{1-y}$. Note that the dependency of the $B_i$ on $X$ is bogus: they must not contain $X$ as terms that have $X$ could be transferred to a lower-indexed $B_{i'}$; I only wrote $X$ because Beck did (and he calls them polynomials, while they are just scalars).
Now what will change if we replace $b_1=3$ by a larger value$~n$? Obviously in our argument $X^3$ gets replaced by $X^n$. The derivations of $a,b$ are very similar: we get $a=\frac1{1-y}$ independently of$~n$, and $b=-\frac{y^{n+1}}{1-y}$. For $C$ our series argument continues to work, and gives $C=\sum_{k=0}^{n-1}(\sum_{i=0}^ky^i)X^k$. This gives some justification to the claim that "the $A_k$ are exponential in $b_1$" (i.e., in $n$), though I cannot really see what the same claim for the $B_j$ should be taken to mean (possibly if one writes $B_j=\frac{1-y^{n+1-j}}{(1-y)^2}$ it makes some kind of sense).
Back to the original problem, we need to take the constant term of our result as a series in$~X$, which results in a rational function of$~y$, and then take the constant term in$~y$ of the Laurent series expansion of that rational function. The first step completely discards the terms coming from $C$ (with the $B_j$), as these have purely negative powers of $X$, and in the terms with $a$ and $b$ we can replace the denominators by$~1$ (as here we are basically setting $X:=0$). So in the general case we are looking for the constant term in the Laurent series expansion of
$$\frac1{(1-y)y^{b_2}}\left(\frac1{1-y}-\frac{y^{n+1}}{1-y}\right)
=\frac{1-y^{n+1}}{(1-y)^2y^{b_2}}.
$$
This is the same a asking for the coefficient of $Y^{b_2}$ in the formal power series for $\frac{1-Y^{n+1}}{(1-Y)^2}$, which can be found to be $\min(b_2+1,n+1)$, and indeed this is the number of ways to write the vector $(n,b_2)$ as a non-negative integer linear combination of $(1,0)$, $(0,1)$, and $(1,1)$.
