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[101];

int b[101];

LL dp[101][2]; 

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