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;
}

  

posted @ 2018-10-16 16:02  lleozhang  Views(162)  Comments(0Edit  收藏  举报
levels of contents