[luogu4133 BJOI2012] 最多的方案 (计数dp)

题目描述

第二关和很出名的斐波那契数列有关，地球上的OIer都知道：F1=1, F2=2, Fi = Fi-1 + Fi-2，每一项都可以称为斐波那契数。现在给一个正整数N，它可以写成一些斐波那契数的和的形式。如果我们要求不同的方案中不能有相同的斐波那契数，那么对一个N最多可以写出多少种方案呢？

输入输出格式

输入格式：

只有一个整数N。

输出格式：

一个方案数

输入输出样例

输入样例#1： 复制

16

输出样例#1： 复制

4

说明

Hint：16=3+13=3+5+8=1+2+13=1+2+5+8

对于30%的数据，n<=256

对于100%的数据，n<=10^18

首先需要知道的是任何一个自然数可以像这样被斐波那契数分解

我们又知道n>=3时斐波那契数可以被前两个数的和替换

那么我们可以先贪心找到一种每个数最大的方案

然后尝试将其中的数替换，使用dp计数

f[i][1/0] 表示最大分解中的第i个是否被替换（0表示被替换）

得到动规方程：

f[i][1]=f[i-1][1]+f[i-1][0];
f[i][0]=f[i-1][0]*(pos[i]-pos[i-1])/2 +
            f[i-1][1]*(pos[i]-pos[i-1]-1)/2;
    //ps:当该数(pos[i])要分解时考虑pos[i-1]能否使用
//即若pos[i-1]被替换(f[i-1][0])则可使用，若没替换则不行（-1即可）

code:

//By Menteur_Hxy
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <string>
#include <map>
#include <vector>
#include <queue>
#include <set>
#include <ctime>
#define M(a,b) memset(a,(b),sizeof(a))
#define F(i,a,b) for(register int i=(a);i<=(b);i++)
#define LL long long
using namespace std;

inline LL rd() {
	LL x=0,fla=1; char c=' ';
	while(c>'9'|| c<'0') {if(c=='-') fla=-fla; c=getchar();}
	while(c<='9' && c>='0') x=x*10+c-'0',c=getchar();
	return x*fla;
}

inline void out(LL x){
    int a[25],wei=0;
    if(x<0) putchar('-'),x=-x;
    for(;x;x/=10) a[++wei]=x%10;
    if(wei==0){ puts("0"); return;}
    for(int j=wei;j>=1;--j) putchar('0'+a[j]);
    putchar('\n');
}

const int MAX=1010;
const int INF=0x3f3f3f3f;
LL n,cnt,scnt;
LL f[MAX][2],fib[MAX],pos[MAX];

int main() {
	n=rd();
	fib[1]=1,fib[2]=2;
	for(cnt=3;fib[cnt-1]<=n;cnt++) fib[cnt]=fib[cnt-1]+fib[cnt-2];
	cnt--;
	for(int i=cnt;i>=1;i--) if(n>=fib[i]) {
		n-=fib[i]; pos[++scnt]=i;
	}
	sort(pos+1,pos+scnt+1);
	f[1][1]=1; f[1][0]=((pos[1]-1)>>1);
	F(i,2,scnt) {
		f[i][1]=f[i-1][0]+f[i-1][1];
		f[i][0]=f[i-1][0]*(pos[i]-pos[i-1]>>1) + f[i-1][1]*(pos[i]-pos[i-1]-1>>1);
	}
	out(f[scnt][0]+f[scnt][1]);
	return 0;
}