bzoj 1002
表示我这种蒟蒻面对这种递推第一思想显然是打表啊
先贴个用来打表的暴力：
#include <cstdio> struct node { int l,r; }p[100005]; bool used[100005]; int f[100005]; 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=1;i<=n+1;i++) { f[i]=i; } for(int i=1;i<=cct;i++) { if(used[i]) { int f1=findf(p[i].l); int f2=findf(p[i].r); f[f2]=f1; } } int ff=0; for(int i=1;i<=n+1;i++) { int f1=findf(i); if(!ff) { ff=f1; }else if(ff!=f1) { return 0; } } return 1; } void dfs(int dep,int tot) { if(dep==cct+1) { if(tot==n) { if(check()) { cnt++; } } return; } used[dep]=1; dfs(dep+1,tot+1); used[dep]=0; dfs(dep+1,tot); } int main() { scanf("%d",&n); for(int i=2;i<=n+1;i++) { p[++cct].l=1; p[cct].r=i; } for(int i=2;i<=n+1;i++) { if(i==n+1) { p[++cct].l=i; p[cct].r=2; }else { p[++cct].l=i; p[cct].r=i+1; } } dfs(1,0); printf("%d\n",cnt); return 0; } /* 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;
}
