dsy 1002 轮状病毒

基尔霍夫矩阵。。。虽然不知道是什么，只知道 f(n)=f(n-1)*3-f(n-2)+2,在直接用高精度就可以。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>

using namespace std;

struct node
{
    int a[101],len;
}f[110];
int n;
node Uzi(node 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;
} 
node OMG(node a,node b)
{
    a.a[1]+=2;
    int j=1;
    while(a.a[j]>=10)
    {
	    a.a[j]%=10;
		a.a[++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]-=1;
		}
    }
    while(a.a[a.len]==0)a.len--;
    return a;
}
int main()
{
    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]=OMG(Uzi(f[i-1],3),f[i-2]);
    }
	for(int i=f[n].len;i>0;i--) printf("%d",f[n].a[i]);
    return 0;
}