打卡信奥刷题(777)用Scratch图形化工具信奥P6056[普及组/提高] [加油武汉] SIR 模型

[加油武汉] SIR 模型

题目背景

SIR 模型将总人口分为以下三类:

  • 易感者(susceptibles),其数量记为 s(t)s(t)s(t) ,表示 ttt 时刻未染病但有可能被该类疾病传染的人数;
  • 染病者(infectives),其数量记为 i(t)i(t)i(t),表示 ttt 时刻已被感染成为病人而且具有传染力的人数;
  • 恢复者(recovered),其数量记为 r(t)r(t)r(t),表示 ttt 时刻已从染病者中移出的人数。

设总人口为N(t)N(t)N(t),则有N(t)=s(t)+i(t)+r(t)N(t)=s(t)+i(t)+r(t)N(t)=s(t)+i(t)+r(t)

SIR模型的建立基于以下三个假设:

  1. 不考虑人口的出生、死亡、流动等种群动力因素。人口始终保持一个常数,即 $N(t) \equiv K $。
  2. 一个病人一旦与易感者接触就必然具有一定的传染力。假设 ttt 时刻单位时间内,一个病人能传染的易感者数目与此环境内易感者总数 s(t)s(t)s(t) 成正比,比例系数为 β\betaβ,从而在t时刻单位时间内被所有病人传染的人数为 βs(t)i(t)\beta s(t)i(t)βs(t)i(t)
  3. ttt 时刻,单位时间内从染病者中移出的人数与病人数量成正比,比例系数为γ\gammaγ,单位时间内移出者的数量为 γi(t)\gamma i(t)γi(t)

题目描述

我们将这个模型简化一下,初始有感染者 III 人和易感者 SSS 人,对于每一天当前有 IiI_iIi 个感染者,SiS_iSi 个易感者,RiR_iRi 个恢复者,则每天会有 ⌈βSiIi⌉\lceil \beta S_iI_i \rceilβSiIi 人被感染(由易感者变成感染者),有 ⌈γIi⌉\lceil \gamma I_i \rceilγIi 人被治愈(由感染者变成恢复者) 。

其中 β\betaβ 为感染系数 γ\gammaγ 为恢复系数 ⌈⌉\lceil \rceil 为上取整符号。

nnn 天后,有多少易感者 SSS,感染者 III,和恢复者 RRR

注: 感染者和恢复者都是每天结算的,结算的结果只和当天开始的时候的值有关,即感染者当天恢复不影响他当天感染别人。

若计算被感染人数超过易感者人数则全员被感染。

输入格式

第一行三个正整数,分别表示第 000 天易感者人数 S0S_0S0 和感染者人数I0I_0I0,以及天数 nnn(刚开始恢复者数 R0=0R_0=0R0=0)。

第二行两个浮点数,分别表示感染系数 β\betaβ 和恢复系数 γ\gammaγ

输出格式

一行三个整数,分别表示 nnn 天后的易感者人数 SSS 、感染者人数 III 和恢复者RRR

样例 #1

样例输入 #1

980 20 2
0.0005 0.00001

样例输出 #1

955 43 2

样例 #2

样例输入 #2

1400000000 1 10
0.000000003 0.001

样例输出 #2

1386791252 13205592 3157

样例 #3

样例输入 #3

1919 810 1
0.00001 0.1

样例输出 #3

1903 745 81

提示

对于 30%30\%30% 的数据,n=1n=1n=1
对于另外 30%30\%30% 的数据,S0,R0≤104S_0, R_0\le 10^4S0,R0104
对于 100%100\%100% 的数据,1≤S0+R0≤2×109,0<β,γ<1,1≤n≤1001 \le S_0+R_0\le 2\times 10^9, 0 < \beta, \gamma < 1, 1 \le n \le 1001S0+R02×109,0<β,γ<1,1n100

Scratch实现

在这里插入图片描述

后续

接下来我会不断用scratch来实现信奥比赛中的算法题、Scratch考级编程题实现、白名单赛事考题实现,感兴趣的请关注,我后续将继续分享相关内容

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值