1002: [FJOI2007]轮状病毒
Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 3153 Solved: 1749
[ Submit][ Status][ Discuss]
Description
给定n(N<=100),编程计算有多少个不同的n轮状病毒。
Input
第一行有1个正整数n。
Output
将编程计算出的不同的n轮状病毒数输出
Sample Input
3
Sample Output
16
根据一些推理可得到f[i]=3*f[i-1]-f[i-2]+2 (证明详见神犇vfk博客http://vfleaking.blog.163.com/blog/static/17480763420119685112649/)
然后用高精度就好了。
其实这道题我就是用来练习高精度了。。。
#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#define F(i,j,n) for(int i=j;i<=n;i++)
#define D(i,j,n) for(int i=j;i>=n;i--)
using namespace std;
int a[105][105];
int n;
int main()
{
scanf("%d",&n);
if (n==1)
{
cout<<1<<endl;
return 0;
}
if (n==2)
{
cout<<5<<endl;
return 0;
}
memset(a,0,sizeof(a));
a[1][0]=1;
a[1][1]=1;
a[2][0]=1;
a[2][1]=5;
F(i,3,n)
{
a[i][1]+=2;
F(j,1,a[i-1][0])
{
a[i][j]=a[i][j]+3*a[i-1][j]-a[i-2][j];
while (a[i][j]<0)
{
a[i][j]+=10;
a[i][j+1]--;
}
a[i][j+1]+=a[i][j]/10;
a[i][j]%=10;
}
int t=a[i-1][0]+1;
if (a[i][t]==0) t--;
while (a[i][t]>10)
{
t++;
a[i][t]+=a[i][t-1]/10;
a[i][t-1]%=10;
}
a[i][0]=t;
}
D(i,a[n][0],1) cout<<a[n][i];
cout<<endl;
}