「BZOJ1002」[FJOI2007] 轮状病毒
给定n(N<=100),编程计算有多少个不同的n轮状病毒。
Input
第一行有1个正整数n。
Output
将编程计算出的不同的n轮状病毒数输出
Sample Input
3
Sample Output
16
代码
基尔霍夫矩阵(我也不知道是什么)推出f[i]=(f[i-1]*3-f[i-2]+2)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 |
#include<iostream> #include<cstdio> using namespace std; struct data{ int a[101],len; }; int n; data mul(data a,int k) { for(int i=1;i<=a.len;i++) a.a[i]*=k; for(int i=1;i<=a.len;i++) { a.a[i+1]+=a.a[i]/10; a.a[i]%=10; } if(a.a[a.len+1]!=0)a.len++; return a; } data sub(data a,data b) { a.a[1]+=2; int j=1; while(a.a[j]>=10){a.a[j]%=10;a.a[j+1]++;j++;} for(int i=1;i<=a.len;i++) { a.a[i]-=b.a[i]; if(a.a[i]<0){a.a[i]+=10;a.a[i+1]--;} } while(a.a[a.len]==0)a.len--; return a; } int main() { data f[101];f[1].a[1]=1;f[2].a[1]=5; f[1].len=f[2].len=1; scanf("%d",&n); for(int i=3;i<=n;i++) f[i]=sub(mul(f[i-1],3),f[i-2]); for(int i=f[n].len;i>0;i--) printf("%d",f[n].a[i]); return 0; } |
是f :=(f[i-1]*3-f[i-2]+1)吧
+2…