怪物牛繁殖问题与常系数线性齐次递推关系式的通解

假设初始有一对牛，每个月生出一对小牛，新生的小牛过三个月成长为大牛，并开始繁殖，在假设所有牛不死亡的情况下，问第 $n$ 个月有多少对牛？ (假设第 $n$ 个月有 $a[n]$ 对牛,那么 $a[1]=2,a[2]=3,a[3]=4,a[4]=5,a[5]=\cdots$依次类推 )

因为新生的小牛经过三个月可以成为大牛并开始繁殖，所以如果假设第 $n$ 个月出生 $b[n]$ 对小牛，那么其等于上个月的大牛对数 $b[n-1]$ 和 第 $n-4$ 个月(经过三个月)出生的小牛对数 $b[n-4]$，所以有
$$b[n]=\{\begin{array}{ll}1& n=1,2,3\\ b[n-1]+b[n-4]& n\ge 4\end{array}$$

再加上有一对初始的大牛，所以第 $n$ 个月的牛对数为 $$a[n]=1+\sum _{i=1}^{n}b[i]$$

所以求出 $b[n]$，即可得出 $a[n]$，求 $b[n]$ 有两种方法：

方法一：
参考Linear recurrence relations with constant coefficients 由递推关系可知
$$\left(\begin{array}{cccc}1& 0& 0& 1\\ 1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\end{array}\right)\left(\begin{array}{c}b[n-1]\\ b[n-2]\\ b[n-3]\\ b[n-4]\end{array}\right)=\left(\begin{array}{c}b[n]\\ b[n-1]\\ b[n-2]\\ b[n-3]\end{array}\right)$$

令 $$M=\left(\begin{array}{cccc}1& 0& 0& 1\\ 1& 0& 0& 0\\ 0& 1& 0& 0\\ 0& 0& 1& 0\end{array}\right)$$

由于 $det(M)=-1\ne 0$，所以 $M$ 可以对角化为 $M=VD{V}^{-1}$，所以有 $V{D}^n{V}^{-1}(b[4],b[3],b[2],b[1]{)}^T=(b[n],b[n-1],b[n-2],b[n-3])$

方法二：
参考Solving homogeneous linear recurrence relations with constant coefficients，假设 $q_1.q_2,q_3,q_4$ 是方程 $x^4=x^3+1$ 的四个解，则显然有 $q^{n-4}\cdot q^4=q^{n-4}\cdot q^3+q^{n-4}$，所以递推关系$b[n]=b[n-1]+b[n-4]$ 的解是
$$b[n]=k_1{q}_1^n+k_2{q}_2^n+k_3{q}_3^n+k_4{q}_4^n$$ ，将初始解
$$\{\begin{array}{ll}b[1]=1& \\ b[2]=1& \\ b[3]=1& \\ b[4]=1& \end{array}$$
代入其中可解出系数$k_1,k_2,k_3,k_4$
所以可以直接推出
$$a[n]=k_1\frac{q_1^{n+1}-q_1}{q_1-1}+k_2\frac{q_2^{n+1}-q_2}{q_2-1}+k_1\frac{q_3^{n+1}-q_3}{q_3-1}+k_1\frac{q_4^{n+1}-q_4}{q_4-1}+1$$
其中 $q$ 为 $x^4-x^3-1=0$ 的解。

其中的一个 $q$ 解为
$q=-\frac{\sqrt{-\frac{\sqrt{3}{\left(\frac{\sqrt{283}}{23^{\frac{3}{2}}}-\frac{1}{2}\right)}^{\frac{1}{6}}}{2\sqrt{12{\left(\frac{\sqrt{283}}{23^{\frac{3}{2}}}-\frac{1}{2}\right)}^{\frac{2}{3}}+3{\left(\frac{\sqrt{283}}{23^{\frac{3}{2}}}-\frac{1}{2}\right)}^{\frac{1}{3}}-16}}-{\left(\frac{\sqrt{283}}{23^{\frac{3}{2}}}-\frac{1}{2}\right)}^{\frac{1}{3}}+\frac{4}{3{\left(\frac{\sqrt{283}}{23^{\frac{3}{2}}}-\frac{1}{2}\right)}^{\frac{1}{3}}}+\frac{1}{2}}}{2}-\frac{\sqrt{12{\left(\frac{\sqrt{283}}{23^{\frac{3}{2}}}-\frac{1}{2}\right)}^{\frac{2}{3}}+3{\left(\frac{\sqrt{283}}{23^{\frac{3}{2}}}-\frac{1}{2}\right)}^{\frac{1}{3}}-16}}{4\sqrt{3}{\left(\frac{\sqrt{283}}{23^{\frac{3}{2}}}-\frac{1}{2}\right)}^{\frac{1}{6}}}+\frac{1}{4}$

常系数线性齐次递推关系式的通解求法：

1. Linear recurrence relations with constant coefficients
2. Solving homogeneous linear recurrence relations with constant coefficients