bzoj1002 [FJOI2007]轮状病毒

题目

先运用DFS把规律找出来，我们可以发现
f[n]=3*f[n-1]-f[n-2]+2

之后又是高精度了。。。

#include<bits/stdc++.h>
using namespace std;
int first[101],nex[201],to[201],size=0;
int n,Ans;
int A[100001],B[100001],C[100001];
int la,lb,lc;
long double f[101];
bool E[201],vis[201];
void add(int x,int y)
{   
    nex[size]=first[x];
    first[x]=size;
    to[size]=y;
    size++;
}
void DFS2(int x)
{
    vis[x]=true;
    for(int i=first[x];i!=-1;i=nex[i])
    {
        if(E[i])
        {
            int u=to[i];
            if(!vis[u])DFS2(u);
        }
    }
}
bool check()
{
    memset(vis,false,sizeof(vis));
    DFS2(0);
    for(int i=0;i<=n;i++)
        if(!vis[i])return false;
    return true;
}
void DFS1(int st,int num)
{
    if(st>size+1)return;
    if(num==n)
    {
        if(check())Ans++;
        return;
    }
    DFS1(st+2,num);
    E[st]=E[st^1]=true;
    DFS1(st+2,num+1);
    E[st]=E[st^1]=false;
}
int main()
{
//  for(n=1;n<=15;n++)
//  {
//      //scanf("%d",&n);
//      memset(first,-1,sizeof(first));
//      size=0;
//      for(int i=1;i<=n;i++)    
//          add(0,i),add(i,0);
//      for(int i=1;i<n;i++) 
//          add(i,i+1),add(i+1,i);
//      add(1,n),add(n,1);
//      size--;
//      Ans=0;
//      DFS1(0,0);
//      cout<<Ans<<" ";
//  }
    /*--1 5 16 45 121 320 841 2205 5776 15125 39601 103680 271441 710645 1860496--*/
    /*--------f[n]=3*f[n-1]-f[n-2]+2;f[1]=1,f[2]=5------*/
    cin>>n;
    if(n==1){
        cout<<1;
        return 0;
    }
    if(n==2){
        cout<<5;
        return 0;
    }
    la=lb=1;
    A[1]=1,B[1]=5;
    for(int j=3;j<=n;j++)
    {
        int x=0;
        for(int i=1;i<=lb;i++)
        {
            x=x+B[i]*3;
            C[i]=x%10;
            x=x/10;
        }
        lc=lb;
        if(x)lc++,C[lc]=x;
//      for(int i=lc;i>=1;i--)
//          cout<<C[i];
//      cout<<endl;
        for(int i=1;i<=lc;i++)
        {
            x=10+C[i]-A[i];
            C[i]=x%10;
            C[i+1]=C[i+1]-1+x/10;
        }
//      for(int i=lc;i>=1;i--)
//          cout<<C[i];
//      cout<<endl;
        while(C[lc]==0)lc--;
        x=C[1]+2;
        for(int i=1;i<=lc;i++)
        {       
            C[i]=x%10;
            C[i+1]=C[i+1]+x/10;
            x=C[i+1];
        }   
        memcpy(A,B,sizeof(B)),la=lb;
        memcpy(B,C,sizeof(C)),lb=lc;
    }
    for(int i=lc;i>=1;i--)
        cout<<C[i];
//  f[1]=1,f[2]=5;
//  for(int i=3;i<=n;i++)
//      f[i]=3*f[i-1]-f[i-2]+2;
//  cout<<endl;
//  printf("%.f",(double)f[n]);
    return 0;
}


附DFS代码