Amoeba细菌繁殖



问题

 

Amoeba

 

农部区域正在研究新发现的阿米巴虫。繁殖力很强的阿米巴虫有如下一样的特征。

• 它做无性生殖。

• 出生之后的第 a天会成为成体。

• 成为成体的那天开始每天都会造出一只新的个体。它成为成体后立马会造出第一个个体，之后每天造出新的个体。新的个体也一样，出生之后第a天开始成为成体，然后造出新的个体。

• 到了出生之后的第 b天开始，不能再造出新的个体了。出生之后的第 a天开始到 b天的前一天可以造出新的个体，所以一生总共能造出 (b−a)只个体。

• 到了出生之后的第d天的那瞬间，他会死亡。

  如下是当 a=2,b=4,d=6时，水槽里放入一只新出生的阿米巴虫后，每天进行观察并记录的结果。括号里面的数字整数表示水槽里面的每个阿米巴虫出生后，经过的天数。

• 出生的天: (0) – 放入新的个体

• 第1天: (1) – 阿米巴虫生长

• 第2天: (2, 0) – 阿米巴虫出生第2天会成为成体，所以会造出新的个体

• 第3天: (3, 1, 0) – 成为成体第2天的阿米巴虫，今天也造出了一个新的个体。

• 第4天: (4, 2, 1, 0) – 第2天造出的阿米巴虫造出了新的个体 (第一次放入的阿米巴虫不能再造出新的个体)

• 第5天: (5, 3, 2, 1, 0, 0)

• 第6天: (4, 3, 2, 1, 1, 0, 0) – 第一次放入的阿米巴虫会死。

  第6天，水槽里面活下来的阿米巴虫总共为7只。

  利用阿米巴虫的繁殖信息a, b, d, 请编写输出当新出生的阿米巴虫放到水槽里后，第N天活着的阿米巴虫个数除以1000后的余数的程序。

   

  输入

   

  第一行上，表示 a, b, d, N的四个整数以空格划分给出。(0<a<b<d≤10,000,1≤N≤1000000)

   

  输出

  输出水槽里放入一只阿米巴虫后，第 N天水槽里活着的阿米巴虫数除以1000的余数。

   

  案例输入 1

2 4 6 6

案例输出 1

7

案例输入 2

3 5 7 20000

案例输出 2

609

 

#include <cstdio>
#include<cstring>

using namespace std;

const int maxn=1000002;
int dp[maxn], sumi[maxn];
int a, b, d, n;

int main()
{
	scanf("%d%d%d%d", &a, &b, &d, &n);
	// init
	memset(dp, 0, sizeof(dp));
	memset(sumi, 0, sizeof(sumi));

	dp[0] = 1; sumi[0] = 1;
	// calc dp for every day
	for (int i = 1; i <= n; i++)
	{
		// add all new babies 
		if (i >= a)	dp[i] += sumi[i - a];
		// minus all babies who mother can not make baby,fix above line
		if (i >= b) dp[i] -= sumi[i - b];
		dp[i] += 2000; // if dp[i] is minus value,1000 can fix it ;
		dp[i] %= 1000;

		sumi[i] = sumi[i - 1] + dp[i];
		sumi[i] %= 1000;
	}
	int ans = 0;
	// last day,all living AMEOBAs are borned in [n, n-1,n-2,...,n-d+1] days.
	for (int i = n; i > n - d && i >= 0; i--)
	{
		ans += dp[i];
		ans %= 1000;
	}
	printf("%d\n", ans);

	return 0;
}