51nod 1350 斐波那契表示

首先，我们简单列一下前几组数据，寻找下规律

x     1  2  3  4  5  6  7  8  9  10  11  12  13
     *  *  *     *       *                *
F(x)   1  1  1  2  1  2  2  1  2   2   2   3   1  
G(x)   1  2  3     6       11               21

打星表示该位的x为斐波那契数

以 x 在9~12的范围为例，其 F(x) 的值等于 x 在1~4范围的F(x)+1。

为什么？很显然，以9为例，减去不大于9的最大的斐波那契数，即8，得到1。所以 F(9)=1+F(1)

即 F(n)=1+F(m)

（m 为不大于 n 的最大的斐波那契数）

如果将9~13作为一整段进行考虑，会发现

F(9)+F(10)+F(11)+F(12)+F(13)= F(1)+F(2)+F(3)+F(4)+F(5) +5-1

可以得到广义递推式： G(n)=G(m)+G(n-m)+(n-m)

（m 为不大于 n 的最大的斐波那契数）

特殊的，当n为斐波那契数时，可以得到

G( fib(x) ) - G( fib(x-1) )  =G( fib(x-2) )+fib(x-2)-1

移项得到

G( fib(x) )=G( fib(x-1) )+G( fib(x-2) )+fib(x-2)-1

设 g_fib(x) =G( fib(x) ) ，我们可以利用递推式，打表求出 g_fib 数组。

那么做法就可以相出：对于 n ，二分查找 m，将 G(m) + (n-m) 累加到答案上（实际算法中应该累加 G(m) + (n-m)*cnt，cnt为循环层数），并将n-m赋值给n，cnt++。直到n为0为止。

#include<bits/stdc++.h>
using namespace std;

const long long MAXM=1e17;
const int MAXN=100;
long long fib[MAXN],g_fib[MAXN];

int main()
{
	long long i,fsum,T,n,ans,cnt,index;
	fib[0]=fib[1]=1;
	g_fib[0]=g_fib[1]=1;
	for(i=2;i<MAXN;i++)
	{
		fib[i]=fib[i-1]+fib[i-2];
		g_fib[i]=g_fib[i-1]+g_fib[i-2]+fib[i-2]-1;
		if(fib[i]>MAXM)
		{
			fsum=i+1;
			break;
		} 
	}
	scanf("%lld",&T);
	while(T--)
	{
		scanf("%lld",&n);
		ans=0;
		cnt=0;
		while(n)
		{
			index=upper_bound(fib,fib+fsum,n)-fib-1;
			ans+=g_fib[index]+cnt*fib[index];
			cnt++;
			n-=fib[index];
		}
		printf("%lld\n",ans);
	}
}