bzoj千题计划213：bzoj2660: [Beijing wc2012]最多的方案

http://www.lydsy.com/JudgeOnline/problem.php?id=2660

很容易想到是先把n表示成最大的两个斐波那契数相加，然后再拆分这两个斐波那契数

把数表示成斐波那契进制的形式，第i位表示有没有第i个斐波那契数

比如16=13+3     001001

那么拆分一个数就是把一个1变成0，左边的两个0变成1

前面的1不影响后面

后面1拆出的两个1不能拆到前面1的前面

所以b[i] 表示n的第i个1是第几项斐波那契数

所以dp[i][0/1] 表示b中的i所在位（n的第b[i]个1）是0/1的方案数

如果这个位是1，dp[i][0]=dp[i-1][0]+dp[i-1][1]

如果这个位是0，即这个1被拆了，他能拆的次数是 与前面的1之间的0的个数/2

所以若i-1是1，两个1之间有 b[i]-b[i-1]-1个0

若i-1是0，两个1之间有b[i]-b[i-1]个0

dp[i][1]=dp[i-1][1]*(b[i]-b[i-1]-1)/2+dp[i-1][0]*(b[i]-b[i-1])/2

#include<cstdio>
#include<algorithm>

typedef long long LL;

using namespace  std;

LL f[];

int b[];

LL dp[][]; 

int main()
{
    LL n;
    scanf("%lld",&n);
    f[]=; f[]=;
    int t;
    for(t=;f[t-]+f[t-]<=n;++t) f[t]=f[t-]+f[t-];
    int m=;
    for(int i=t-;i;--i)
        if(n>=f[i]) b[++m]=i,n-=f[i];
    reverse(b+,b+m+);
    dp[][]=;
    dp[][]=b[]->>;
    for(int i=;i<=m;++i)
    {
        dp[i][]=dp[i-][]+dp[i-][];
        dp[i][]=dp[i-][]*(b[i]-b[i-]->>)+dp[i-][]*(b[i]-b[i-]>>);
    }
    printf("%lld",dp[m][]+dp[m][]);
}