[算法设计与分析]3.4.5特征根求解递推方程(粒子分裂)

[算法设计与分析]3.4.5特征根求解递推方程(粒子分裂)

#include<stdio.h>
#include<iostream>
#include<math.h>using namespace std;void ParticleSplitting1(int t);
void ParticleSplitting2(int t);int main ()
{int t = 3;ParticleSplitting1(t);ParticleSplitting2(t);
}void ParticleSplitting1(int t)
{int a[100], b[100];//代表α和β粒子a[0] = 1;b[0] = 0;for(int i = 1; i <= t; i++){a[i] = b[i - 1];//因此a[1]=0因为1s之后α粒子变成了3个β粒子 此时原α粒子消失b[i] = 3 * a[i - 1] + 2 * b[i - 1];//当前时刻的β粒子是由前一时刻的α粒子和β粒子分裂而来的}cout << t << "s之后" << "α粒子的个数是：" << a[t] << "ntβ粒子的数目是：" << b[t];
}
//数学模型见书P116
void ParticleSplitting2(int t)
{//等同于exp(n)等同于e^nint a = int(exp(t * log(3)));//log(x) 的作用就是ln(x)int b = int(exp((t + 1) * log(3)));if(t %