http://acm.hdu.edu.cn/showproblem.php?pid=1041
有一个初始只有一个1的串 每次都按①0 -> 10;②1 -> 01;这两条规则进行替换
形如:n = 1 1
n = 2 01
n = 3 1001
...
求经过n步替换之后 串中只含复数个0的连续子串(不难发现,这种子串只能是‘00’)的出现次数
因为0<n<=1000的限制 在最坏情况下(n==1000)串的长度将达到2^1000位 排除了直接模拟上述替换过程的可能
列出前几项的替换结果:
n = 0 1 ‘00’=0 ‘01’=0
n = 1 01 ‘00’=0 ‘01’=1
n = 2 1001 ‘00’=1 ‘01’=0
n = 3 01101001 ‘00’=1 ‘01’=2
n = 4 1001011001101001 ‘00’=3 ‘01’=2
n = 5 01101001100101101001011001101001 ‘00’=5 ‘01’=6
在上面的数据 不仅给出结果串 还统计出了串中‘00’对和‘01’对的个数【注:‘01’的个数是取出全部的‘00’对后才统计】
从n = 3开始 观察可以发现 01 -> 1001 即每个00对(1001)都是由01对转化而来
而每个00对(1001)又将产生一个 00对 和 两个 01对(见n=2->3)
由上述可得 n步时00对的个数 = n-1步的00对的个数 + n-1的01对的个数(即n-2步的00对的个数 * 2)
所以可以推得 f(n) = f(n - 1) + f(n - 2) * 2
如果细心的话,可以发现,当n大到一定程度时,由于上面给出的公式有着与斐波拉契数列相似的递增效应, 最终结果的数值会非常大大大大大,用__int64都远远无法满足
所以采用字符数组来模拟大数加法运算
# include <stdio.h>
# include <string.h>
# define MAX 1001
# define LEN 1001
char Number[MAX][LEN];
void BigNumPlus()
{
for(int i = 3; i < MAX; i++)
{
memset(Number[i], '0', LEN);//初始化
for(int j = 0; j < LEN; j++)//模拟公式
{
Number[i][j] += (Number[i - 1][j] - '0') + (Number[i - 2][j] - '0') + (Number[i - 2][j] - '0');
}
for(int j = 0; j < LEN; j++)//处理进位
{
if(Number[i][j] > '9')
{
Number[i][j + 1] += (Number[i][j] - '0') / 10;
Number[i][j] = (Number[i][j] - '0') % 10 + '0';
}
}
}
}
int main()
{
//初始值
memset(Number[1], '0', LEN);
memset(Number[2], '0', LEN);
Number[1][0] = '0';
Number[2][0] = '1';
BigNumPlus();
int n;
while(scanf("%d",&n) != EOF)
{
if(n == 1)//特殊处理
{
printf("0\n");
continue;
}
for(int i = LEN - 1; i >= 0; i--)
{
if(Number[n][i] != '0')//格式输出
{
while(i >= 0) printf("%c",Number[n][i--]);
printf("\n");
break;
}
}
}
return 0;
}