Permítanme asumir que el problema realmente es la siguiente:
Encontrar todos los de la familia de los números de $(P(n,r))_{n\geqslant1,r\geqslant0}$ tal que, para cada $n\geqslant1$, $m\geqslant1$ y $r\geqslant0$, $P(n,0)=P(1,r)=1$ y
$$
P(n+m,r)=\sum_{s=0}^rP(n,s)P(m,r-s).
$$
Un método para resolver este problema es introducir, para cada una de las $n\geqslant1$, el formal de la serie de $Q_n(x)$ definido por
$$
Q_n(x)=\sum_{i=0}^{+\infty}P(n,r)x^r,
$$
para algunos el parámetro $x$, y el uso de estos objetos de la siguiente manera:
- Identificar las $Q_1(x)$.
- Traducir la recursión anterior en una relación entre cada $Q_{n+m}(x)$, $Q_n(x)$ y $Q_m(x)$.
- Deducir, a partir de 1. y 2. una fórmula explícita para cada $Q_n(x)$.
- Deducir, a partir de 3. el valor de cada $P(n,r)$.
Ahora, disparar...
Edit: Para diferentes condiciones iniciales tales como $P(n,0)=1$$P(1,r)=\cos(r)$, el método permanece exactamente como descrito anteriormente, sólo el resultado del paso 1. los cambios y la viabilidad del paso 4. es disminuido. En este caso específico, se obtiene
$$
\sum\limits_{i=0}^{+\infty}P(n,r)x^r=\left(\frac{1-x\cos 1}{1-2x\cos 1+x^2}\right)^n.
$$
Para ver esto, observe que $P(1,r)$ es la parte real de la $\mathrm e^{\mathrm i r}$, por lo tanto $Q_1(x)$ es la parte real de
$$
\sum\limits_{i=0}^{+\infty}\mathrm e^{\mathrm i r}x^r=\frac1{1-x\mathrm e^{\mathrm i}}=\frac{1-x\mathrm e^{-\mathrm i}}{|1-x\mathrm e^{\mathrm i}|^2}=\frac{1-x\cos1+\mathrm ix\sin1}{1-2x\cos 1+x^2}.
$$