How much weight is on each person in a human pyramid?

For a generating function approach, Let the coefficient of $x^k$ in $P_n(x)$ be the weight on the shoulders of person $k$ (starting at $k=0$ on the left) in row $n$ (starting at $n=0$ at the top). To compute $P_n(x)$ from $P_{n-1}(x)$, we first add the weight of the people in row $n-1$ to the weights on their shoulders: $$ P_{n-1}(x)+\frac{1-x^n}{1-x}\tag{1} $$ Those people distribute half their weight to each person below them, so we get $$ P_n(x)=\frac{1+x}2\left(P_{n-1}(x)+\frac{1-x^n}{1-x}\right)\tag{2} $$ Solving this recurrence yields $$ P_n(x)=\frac{1+x}{(1-x)^2}\left(1-2\left(\frac{1+x}{2}\right)^{n+1}+x^{n+1}\right)\tag{3} $$ $(3)$ can be checked by noting that $P_0(x)=0$ and then by checking that $(3)$ satisfies $(2)$.

Using the binomial theorem to multiply the terms in $(3)$, we get that $w_{m+1,n+1}$, the coefficient of $x^m$ in $P_n(x)$, is given by $$ 2m+1-\frac1{2^n}\sum_{k=0}^m(k+1)\binom{n+2}{m-k}\tag{4} $$ There might be some special function of which I am unaware (a hypergeometric function perhaps), but in terms of common functions, $(4)$ is as close to a closed form as I think there is.


Solving the recurrence

Let $Q_n(x)=\left(\frac2{1+x}\right)^nP_n(x)$, then multiplying $(2)$ by $\left(\frac2{1+x}\right)^n$ yields $$ \begin{align} Q_n(x) &=Q_{n-1}(x)+\frac{1-x^n}{1-x}\left(\frac2{1+x}\right)^{n-1}\\ &=Q_{n-1}(x)+\frac12\frac{1+x}{1-x}\left(1-x^n\right)\left(\frac2{1+x}\right)^n\\ &=Q_{n-1}(x)+\frac12\frac{1+x}{1-x}\left(\left(\frac2{1+x}\right)^n-\left(\frac{2x}{1+x}\right)^n\right)\tag{5} \end{align} $$ Summing the geometric series gets a bit messy, $$ \begin{align} Q_n(x) &=\frac12\frac{1+x}{1-x}\sum_{k=1}^n\left(\left(\frac2{1+x}\right)^k-\left(\frac{2x}{1+x}\right)^k\right)\\ &=\frac12\frac{1+x}{1-x}\left(\frac2{1+x}\frac{\left(\frac2{1+x}\right)^n-1}{\frac2{1+x}-1}-\frac{2x}{1+x}\frac{\left(\frac{2x}{1+x}\right)^n-1}{\frac{2x}{1+x}-1}\right)\\ &=\frac{1+x}{1-x}\left(\frac{\left(\frac2{1+x}\right)^n-1}{1-x}+x\frac{\left(\frac{2x}{1+x}\right)^n-1}{1-x}\right)\\ &=\frac{1+x}{(1-x)^2}\left(\left(\frac2{1+x}\right)^n-1+x\left(\frac{2x}{1+x}\right)^n-x\right)\\ &=\frac{1+x}{(1-x)^2}\left(\left(\frac2{1+x}\right)^n\left(1+x^{n+1}\right)-(1+x)\right)\tag{6} \end{align} $$ Multiplying by $\left(\frac{1+x}2\right)^n$, we arrive at $(3)$: $$ P_n(x)=\frac{1+x}{(1-x)^2}\left(1+x^{n+1}-2\left(\frac{1+x}2\right)^{n+1}\right)\tag{7} $$


Examples: $$ \begin{align} P_0(x)&=0\\ P_1(x)&=\frac12+\frac12x\\ P_2(x)&=\frac34+\frac32x+\frac34x^2\\ P_3(x)&=\frac78+\frac{17}{8}x+\frac{17}{8}x^2+\frac78x^3\\ P_4(x)&=\frac{15}{16}+\frac52x+\frac{25}{8}x^2+\frac52x^3+\frac{15}{16}x^4 \end{align} $$ The bottom-middle person on a $5$-tier pyramid is supporting more than $3$ people's weight on their back.


(I've deleted the previous answer which was wrong.)

Now I have, without proof, by heuristics a direct formula for the entries in the weight matrix if we consider its entries as $$ W=\small \begin{bmatrix} 0 & . & . & . & . & . \\ 1/2 & 1/2 & . & . & . & . \\ 3/4 & 3/2 & 3/4 & . & . & . \\ 7/8 & 17/8 & 17/8 & 7/8 & . & . \\ 15/16 & 5/2 & 25/8 & 5/2 & 15/16 & . \\ 31/32 & 87/32 & 61/16 & 61/16 & 87/32 & 31/32 \end{bmatrix}$$ as it was also shown by Rob Johnson.

It was easier to find the recurrences and then the direct formula when I included the unit-weight per person, such that I have: $$ X= \small \begin{bmatrix} 1 & . & . & . & . & . \\ 3/2 & 3/2 & . & . & . & . \\ 7/4 & 5/2 & 7/4 & . & . & . \\ 15/8 & 25/8 & 25/8 & 15/8 & . & . \\ 31/16 & 7/2 & 33/8 & 7/2 & 31/16 & . \\ 63/32 & 119/32 & 77/16 & 77/16 & 119/32 & 63/32 \end{bmatrix} $$ or, expanded to make it better visible, rowscaled by powers of 2: $$ Y =\small \begin{bmatrix} 1 & . & . & . & . & . \\ 3 & 3 & . & . & . & . \\ 7 & 10 & 7 & . & . & . \\ 15 & 25 & 25 & 15 & . & . \\ 31 & 56 & 66 & 56 & 31 & . \\ 63 & 119 & 154 & 154 & 119 & 63 \end{bmatrix} $$

Let $r$ denote the rowindex, $c$ denote the column-index, both beginning at zero, $c \le r$ then we have $$ w_{r,c} = 1+2c - {\sum_{k=0}^c \binom{r+2}{k}\cdot(1+c-k)\over 2^r}$$

In terms of the $X$-matrix this was: $$ x_{r,c} = 2\cdot(1+c) - {\sum_{k=0}^c \binom{r+2}{k}\cdot(1+c-k)\over 2^r}$$ and in terms of the $Y$-matrix this was $$ y_{r,c} = 2^{r+1}\cdot(1+c) - \sum_{k=0}^c \binom{r+2}{k}\cdot(1+c-k)$$


[added] I found it interesting to consider a more general aspect. The idea is to assume an infinitely wide pyramid of "zero-humans" or "ghosts" which have weight $o$ (which for the final evaluation is then to set to zero). This is the infinitely wide grey shaded area in the following picture. The structure for this area is very simple and homogenuous: in each row $r$ the weight in one entry is $r \cdot o$.
Now one of the top-most ghosts "gets a soul", becomes human and thus gravitation gives him a pressing weight "w". All ghosts which are touched by his feet become then human too, and their weight changes to "w" as well, and this iterates downwards.

The picture shows, how this distributes to the bottom: enter image description here

Note, how the coefficients at a w and at the o in the same cell add to make the equilibrium-coefficient of the same row, if we would set $w=o$.

It is in general the description of a certain disturbance inserted into a equilibrium, and I'm curious, whether this can even be made continuous (like the generalization of the binomial-coefficients using the Gamma-function)


[added2] With iterations of this type it is often useful to express the problem in a matrix-notation and convert iteration into matrix-powers. Then sometimes one can find a closed form in terms of fixed matrices and the iteration-parameter is then only in the exponent of a matrix. (Clearly that does not finally solve the problem of the finite sums, sometimes recursive, in the original, non-matrix, representation but moves it into the internal dotproducts for the matrix-powers. But for the reader of the formula this might be convenient).

First let's define a rowvector $U=[1,0,0,...]$ of some dimension n which represents the weight of the person on the top where n is the final width of the pyramide.

Then one iteration should have the effect, that $U_0$ changes to $V_1=[1/2,1/2,0,0,...]$ which is the received weight for the lower layer and then $U_1$ must also have the unit-weights for the persons to become $U_1=[3/2,3/2,0,0,0...]$ .

A matrix, which performs $U_0 \to V_1$ is $M= \frac12 \small\begin{bmatrix}1&1&0&0&0&0&...\\0&1&1&0&0&0&...\\0&0&1&1&0&0&... \\ \vdots & & & & & &\ddots \end{bmatrix}$ such that $$U_0 \cdot M = V_1 = \left[ \frac 12,\frac12, 0,0,0,... \right] $$.
After that we need a matrix, which transfers two ones for the unit weights into $V_1$ and we can do this simply by addition of the unscaled $M$ such that $$ U_0 \cdot 2M = [1,1,0,0,...]$$ and we have $$ U_0 \cdot M + U_0 \cdot 2M = U_1 = [3/2,3/2,0,0,...]$$

But now, if the number of persons increase by the iterations, the second summand is not sufficient; we need so many unit-weigths as we have persons in one iteration. So we need to find a matrix, whose k't power gives k leading ones in the top row. A simple candidate is the matrix of the upper unit-subdiagonal, let's call it $L$: $$ L = \small \begin{bmatrix}0&1&0&0&0&0&...\\0&0&1&0&0&0&...\\0&0&0&1&0&0&... \\ \vdots & & & & & &\ddots \end{bmatrix} $$ Then the k'th power of shifts the diagonal to the right to the k'th position, and the sum of $L^0 + L^1 + L^2 + ... + L^K$ gives $k+1$ ones. Then for the first iteration we can write $$ U_0 \cdot M + U_0 \cdot (L^0 + L^1) = U_1 $$ or $$ U_0 \cdot (M + (L^0 + L^1)) = U_1 $$

We can then simply proceed by $$ U_1 \cdot (M + (L^0 + L^1 + L^2)) = U_2 $$ which is also $$ U_0 \cdot (M^2 + (L^0 + L^1)M+ (L^0 + L^1 + L^2)) = U_2 $$ and it is obvious how this iterates.

But we see now, that we have partial geometric series depending on the iteration-parameter, and they have a closed form, namely the cyclotomic polynomial of the k'th order.
For the cyclotomic polynomial we need to write $$ L^0+L^1+... + L^k = (I - L^{k+1}) \cdot (I - L)^{-1} $$ Let's denote that inverse as $R$ . We find that it is just the upper triangular unit-matrix (as expected by the series representation!) $$R = \small \begin{bmatrix}1&1&1&1&1&1&...\\0&1&1&1&1&1&...\\0&0&1&1&1&1&... \\ \vdots & & & & & &\ddots \end{bmatrix} $$ Then we have $$ U_0 \cdot ((I-L^1)R M^k + (I-L^2)R M^{k-1}+ (I-L^3)R M^{k-2} +...+(I-L^{k+1})R) = U_k $$ (where I smoothed the formula at the leading $M^k$

Because it is also $M=I + L$ the matrices $M$ and $L$ and $R$ commute all and we can rotate the coefficients for any convenience. We can isolate a geometric sum of powers of $M$:

$$ U_0 \cdot ( R( M^k + M^{k-1}+ M^{k-2} +...I) \\ -R(L^1 M^k + L^2 M^{k-1}+ L^3 M^{k-2} +...+L^{k+1})) = U_k $$ Now the powers of $M$ are also partial geometric series, and we need the reciprocal, let's denote it as $S$: $$ S= (I - M)^-1 =\small \begin{matrix} 2&2&2&... \\ 0&2&2& ...\\ 0&0&2&...\\ ... \end{matrix}$$ which happens to be just $S=2R$ and is also commuting with every matrix in the game so far.

Then we have $$ U_0 \cdot \left( R( I - M^{k+1})S \\ -R(L^1 M^k + L^2 M^{k-1}+ L^3 M^{k-2} +...+L^{k+1}) \right) = U_k \\ $$
and the last parenthese of variable length can be reduced similarly because $M$ is invertible and with a matrix $W=L \cdot M^{-1}$ and $T=(I -W)^{-1} $we have $$ U_0 \cdot \left( R( I - M^{k+1})S -R L M^k (I + W+ W^2 +...+W^k) \right) \\ = U_0 \cdot \left( R( I - M^{k+1})S -R L M^k (I -W^{k+1}) T) \right) \\ \qquad \qquad = U_k \\ $$

This formula has a fixed length for any iteration parameter. Only the variable matrix powers contain that parameter. (Possibly from this one can invent fractional iterates, but that's not the focus here...)


Setup and change of indices

Put $W=1$. Also change the indices of $w_{i,j}$ so that $(i,j)$ points relatively to the top person at (0,0) by moving $i$ persons down-right and then $j$ persons down-left. So (i,j) points to the $i$'th person in the $(i+j)$'th row counting from zero.

The weight distribution of each person

Let us first consider how the top persons weight is distributed down through the pyramid. Then we have the recurrence $$ w_{i,j}=\frac{w_{i-1,j}+w_{i,j-1}}{2} $$ where $w_{0,0}=1$ and $w_{i,j}=0$ outside of the pyramid (that is $i<0$ or $j<0$). This can easily be solved to get $w_{i,j}=\frac{1}{2^{i+j}}\binom{i+j}{i}$. So this is just a variation over Pascal's Triangle where the $n$'th row is divided by $\frac{1}{2^n}$. Let VPT denote this Variation over Pascal's Triangle.

Now move on to consider how the weight of the person at position $(i,j)$ is distributed. The pattern is basically the same! Isolating each person like this, we can see how the distibution of a single persons weight forms a copy of VPT having its top corner at $(i,j)$.

Analysis of the general case

Consider the person at $(i,j)$. This person is included in all VPT's with top corners in $P=\{0,...,i\}\times\{0,...,j\}$ and carries the weight distributed according to all these VPT's. Only we have added 1 too much since the person is included as top corner in its own VPT but does not carry its own weight in the scope of this question.

Now consider $(s,t)\in P$ as one such top corner. Then the position of $(i,j)$ relative to this will be $(x,y)=(i-s,j-t)\in P$. Note here how the map $(s,t)\mapsto (i-s,j-t)$ maps $P$ to itself bijectively. Given this observation one may deduce that the distributed weight contributions form a sum that can be expressed as: $$ w_{i,j}+1=\sum_{x=0}^i\sum_{y=0}^j\frac{1}{2^{x+y}}\binom{x+y}{x} $$ Though not closed form, I believe this to be a very nice and symmetric expression for $w_{i,j}$. Still, I consider the answer given by robjohn just brilliant. BTW my formula yields the same values as robjohn got, and the next example of his should read:

Example (continued from robjohn)

$P_n(x)=\tfrac{31}{32}+\tfrac{87}{32}x+\tfrac{61}{16}x^2+\tfrac{61}{16}x^3+\tfrac{87}{32}x^4+\tfrac{31}{32}x^5$

Also note that the sum of coefficeints of each row in the pyramid should be a triangular number since the persons in the $n$'th row carries the weight of $T_n$ people (the $n$'th triangular number).