bzoj1002: [FJOI2007]轮状病毒(基尔霍夫矩阵)

1002: [FJOI2007]轮状病毒

题目：传送门

 

题解：

　　 决定开始板刷的第一题...

　　 看到这题的时候想：这不就是求有多少种最小生成树的方式吗？

　　 不会啊！！！%题解。。。

　　 什么鬼？基尔霍夫矩阵？？？？OTZ...

　　 什么叫基尔霍夫矩阵就自己去学吧，博主太菜也不会啊...

　　 总之答案就是递归出来的：F(n)=3*F(n-1)-F(n-2)+2

　　 在码个高精度就莫名其妙的A了...真的是毒瘤...

代码：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<algorithm>
using namespace std;
struct node
{
    int a[110],len;
}f[110];
node chengfa(node x,int k)
{
    for(int i=1;i<=x.len;i++)x.a[i]*=k;
    for(int i=1;i<=x.len;i++)
    {
        x.a[i+1]+=x.a[i]/10;
        x.a[i]%=10;
    }
    while(x.a[x.len+1]!=0)x.len++;
    return x;
}
node jianfa(node x,node y)
{
    x.a[1]+=2;
    int j=1;while(x.a[j]>=10){x.a[j]%=10;x.a[++j]++;}
    for(int i=1;i<=x.len;i++)
    {
        x.a[i]-=y.a[i];
        if(x.a[i]<0)
        {
            x.a[i]+=10;
            x.a[i+1]--;
        }
    }
    while(x.a[x.len]==0)x.len--;
    return x;
}
int n;
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]=jianfa(chengfa(f[i-1],3),f[i-2]);
    for(int i=f[n].len;i>=1;i--)printf("%d",f[n].a[i]);
    return 0;
}