首先,我们简单列一下前几组数据,寻找下规律
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);
}
}