```

#include <stdio.h>

#include<string.h>

int main(){

int S1=330000000;

int S=330000000;

int day=1;

int I=34;

int E=65;

int R=0;

float d=1.0/4;

float? d3=40;

float w=0;

float re=0.7;

float ri=0.7;

int dS,dE,dI,dR;

float p=0.001;

float dp;

while(day<1000&&S>0){

day++;

if(day>=40){

dp=0;

p=1-d/re;

}

else

dp=(1-p)*(1-p)/36;

if(day>=80&&d3>10){

d3=d3-0.2;

}

dS=-(w*I*ri+(1-p)*E*re)*(float)S/S1;

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 %d 確診病人變化:%d,潛伏期病人變化:%d,退出者變化:%d,普通人變化：%d ",day,dI,dE,dR,dS);

printf("確診累計(jì)病人:%lld,潛伏期病人:%lld,退出者:%lld，現(xiàn)存確診：%d\n",I+R,E,R,I);

// printf("%lld\n",I);

}

return 0;

}

```