```
#include <stdio.h>
#include<string.h>
int main(){
long long S=1400000000;
long long day=1;
long long I=41;
long long E=60;
long long R=0;
float d=1.0/5;
long long? d3=10;
float w=0.6;
float re=0.7;
float ri=0.7;
long long dS,dE,dI,dR;
float p=0.1;
float dp;
while(day<68){
day++;
if(day<=12){
dp=(1-p)*(1-p)/20.0;
}
if(day==12){
w=0;
}
if(day>12&&day<20){
dp=(1-p)*(1-p)/(15);
}
if(day>=20&&day<25)
{
dp=(1-p)*(1-p)/(10);
}
if(day>=25&&day<29)
{
dp=(1-p)*(1-p)/3;
}
if(day>=29)
{
dp=(1-p)*(1-p);
}
dS=-(w*I*ri+(1-p)*E*re);
dE=w*I*ri+(1-p)*E*re-E*d;
dI=d*E-I/d3;
dR=I/d3;
S=S+dS;
E=E+dE;
I=I+dI;
R=R+dR;
p=p+dp;
printf("day %lld 確診病人變化:%lld,潛伏期病人變化:%lld,退出者變化:%lld\n",day,dI,dE,dR);
printf("確診累計(jì)病人:%lld,潛伏期病人:%lld,退出者:%lld\n",I+R,E,R);
//printf("%lld\n",I+R+E);
}
return 0;
}
```
