bzoj1002【fjoi2007】轮状病毒

1002: [FJOI2007]轮状病毒

Time Limit: 1 Sec  Memory Limit: 162 MB
Submit: 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;}