Codeforces 126D

感谢wuyiqi博客中提供的题，这道题我觉得还不错，难度适中~

首先有一点就是如何去找任意一组解，其实贪心就可以，从大到小来，只要n>=f[i],n-=f[i].为什吗这样的贪心是可行的，因为当i逐渐减小时，相邻的f[i]差值递减，并且最终等于0，而且f[1]=1，大概有一种极限的思想。。。

想明白这个就好办了，把n的某一种表达法换成01序列的模式，比如14=13+1,可以写成100001（1,2,3,5,8,13）

然后就dp了，由于每个001可以拆成110，进而还可以再拆，知道遇到1为止。附代码：

#include <iostream>
#include <cstdio>
#include <algorithm>

#define N 100
#define ll __int64
#define ss(a) scanf("%I64d",&a)

using namespace std;

ll inf,f[N],dp[N][2],b[N];

void init(int &m)
{
    int i;
    inf=1;
    for (i=1;i<=18;i++) inf*=10;
    i=1;
    f[i]=1;f[++i]=2;
    while (f[i]<=inf) f[++i]=f[i-1]+f[i-2];
    m=i-1;
}

int main()
{
    int T,i,m;
    ll n;
    cin>>T;
    init(m);
    while (T--)
    {
        ss(n);
        int t=0;
        for (i=m;i>=1;i--) 
            if (n>=f[i]) 
            {
                n-=f[i];
                b[++t]=i;    
            }     
        sort(b+1,b+t+1);
        b[0]=0;dp[0][1]=1;
        dp[0][0]=0;
        for (i=1;i<=t;i++) 
        {
            dp[i][1]=dp[i-1][1]+dp[i-1][0];
            dp[i][0]=dp[i-1][1]*((b[i]-b[i-1]-1)/2)+dp[i-1][0]*((b[i]-b[i-1])/2);
        }
        ll res=dp[t][1]+dp[t][0];
        printf("%I64d\n",res);
    }
    return 0;
}