Probability generating function of negative binomial distributions I have recently took a course on probability theory and learned negative binomial distribution. The probability generating function (pgf) for negative binomial distribution under the interpretation that the the coefficient of $z^k$ is the number of trials needed to obtain exactly n successes is $F(z)=(\frac{pz}{1-qz })^n=\sum_k\ \binom {k-1} {k-n} p^nq^{k-n}z^k$. But Donald Knuth gave another definition of the pgf of negative binomial distribution on Concrete Mathematics, which is $G(z)=(\frac{p}{1-qz})^n=\sum_k\binom {n+k-1} k p^nq^kz^k$, which's combinatorial interpretation is the distribution of k failures when the number of successes n is specified. So my question  is which of these formulae is actually the pgf of negative binomial distribution, if not, what's the exact name for that pgf? And why there's inconsistency on different definitions? I’m very confused about this question. Actually the definition of negative binomial distribution on Wikipedia is different from that of on Introduction to Probability too.
 A: There are indeed two different versions of the negative binomial distribution: (i) the distribution of the number (say $X$) of independent Bernoulli trials (with a given success probability $p$ in each trial) needed to have $k$ successes and (ii)  the distribution of the number (say $Y$) of failures in an infinite series of such trials before the $k$th success. Obviously, $X=Y+k$. Each of these two "real-world" interpretations may have a certain appeal. 
The existence of these two different (but closely and simply related) versions of the negative binomial distribution is perhaps unfortunate and confusing, but it is a fact. There is apparently nothing that can be done about it.
A: As Iosif Pinelis noted there are two versions if the geometric distribution (the special negative binomial distribution with parameter $n = 1$ successes in your notation), which are related by the formular $X = Y + 1$. These have the pdf $p_X = \frac{pz}{1-qz}$ and $p_Y = \frac{p}{1-qz}$. Since the times between successes are independent and identically distributed, you get $F$ resp. $G$ simply by taking of $p_X$ and $p_Y$ the $n$-th power, as in your formulas. 
