表示我这种蒟蒻面对这种递推第一思想显然是打表啊

先贴个用来打表的暴力:

#include <cstdio>
struct node
{
    int l,r;
}p[];
bool used[];
int f[];
int n,cnt,cct;
int findf(int x)
{
    if(x==f[x])
    {
        return x;
    }
    return f[x]=findf(f[x]);
}
bool check()
{
    for(int i=;i<=n+;i++)
    {
        f[i]=i;
    }
    for(int i=;i<=cct;i++)
    {
        if(used[i])
        {
            int f1=findf(p[i].l);
            int f2=findf(p[i].r);
            f[f2]=f1;
        }
    }
    int ff=;
    for(int i=;i<=n+;i++)
    {
        int f1=findf(i);
        if(!ff)
        {
            ff=f1;
        }else if(ff!=f1)
        {
            return ;
        }
    }
    return ;
}
void dfs(int dep,int tot)
{
    if(dep==cct+)
    {
        if(tot==n)
        {
            if(check())
            {
                cnt++;
            }
        }
        return;
    }
    used[dep]=;
    dfs(dep+,tot+);
    used[dep]=;
    dfs(dep+,tot);
}
int main()
{
    scanf("%d",&n);    
    for(int i=;i<=n+;i++)
    {
        p[++cct].l=;
        p[cct].r=i;    
    }
    for(int i=;i<=n+;i++)
    {
        if(i==n+)
        {
            p[++cct].l=i;
            p[cct].r=;
        }else
        {
            p[++cct].l=i;
            p[cct].r=i+;
        }
    }
    dfs(,);
    printf("%d\n",cnt);
    return ;
}
/*
1 1
2 5
3 16
4 45
5 121
6 320
*/

实测这个打表程序是正确的(可以获得30分)

接下来是本人心路历程:

观察一下:1-1,2-5,3-16,4-45...找一下前后项吧!

观察前后项的倍数关系应该在2~3之间,那先定一个基础表达式

f[i]=2f[i-1]+...或f[i]=3f[i-1]+...

如果系数用2,发现剩下的部分长这样啊...

f[3]=2f[2]+6

f[4]=2f[3]+13

f[5]=2f[4]+31

...

好像后面的没啥规律...

那换系数用3!

f[4]=3f[3]-3

f[5]=3f[4]-14

...

好像也没啥啊...

等一下!

-3=-5+2

-14=-16+2

如果下面再写下来,应该是-43=-45+2!

这不就找出来了吗!

f[i]=3f[i-1]-f[i-2]+2!

于是敲个高精度这题就结束了...

#include <cstdio>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
struct Bignum
{
int a[10005];
int ilen;
}f[105],zero;
int n;
Bignum add(Bignum x)
{
Bignum ret=x;
ret.a[1]+=2;
int i=1;
while(ret.a[i]>=10)
{
ret.a[i+1]+=ret.a[i]/10;
ret.a[i]%=10;
i++;
}
if(ret.a[ret.ilen+1])
{
ret.ilen++;
}
return ret;
}
Bignum mul(Bignum x)
{
Bignum ret=zero;
for(int i=1;i<=x.ilen;i++)
{
ret.a[i]+=3*x.a[i];
ret.a[i+1]+=ret.a[i]/10;
ret.a[i]%=10;
}
ret.ilen=x.ilen;
while(ret.a[ret.ilen+1])
{
ret.ilen++;
}
return ret;
}
Bignum sub(Bignum x,Bignum y)
{
Bignum ret=zero;
for(int i=1;i<=y.ilen;i++)
{
ret.a[i]+=x.a[i]-y.a[i];
if(ret.a[i]<0)
{
ret.a[i]+=10;
ret.a[i+1]--;
}
}
for(int i=y.ilen+1;i<=x.ilen;i++)
{
ret.a[i]+=x.a[i];
}
ret.ilen=x.ilen;
while(!ret.a[ret.ilen]&&ret.ilen>1)
{
ret.ilen--;
}
return ret;
}
int main()
{
scanf("%d",&n);
f[1].a[1]=1;
f[1].ilen=1;
f[2].a[1]=5;
f[2].ilen=1;
f[3].a[1]=6;
f[3].a[2]=1;
f[3].ilen=2;
for(int i=4;i<=n;i++)
{
f[i]=mul(f[i-1]);
f[i]=add(f[i]);
f[i]=sub(f[i],f[i-2]);
}
for(int i=f[n].ilen;i>=1;i--)
{
printf("%d",f[n].a[i]);
}
printf("\n");
return 0;
}

  

05-11 11:11