HAUT 1286: czl的复仇计划【fib公式演变】

xiaoxiao2021-02-28  77

1286: czl的复仇计划

时间限制: 1 秒   内存限制: 128 MB 提交: 107   解决: 20 提交  状态 

题目描述

还记得上次czl给ykc做苦力的事情了吗,czl越想越气,于是决定报复一下ykc,可他知道ykc是不会轻易上当的,呢该怎么办呢,咦……有了,“嗨,ykc,我最近爱上了编程,可是遇到了一个难题,想了一晚上都没有任何思路,是关于斐波那契的,能请教一下你吗?”,“啊,斐波那契呀,呢不是so easy嘛。”ykc自信的说道,czl内心已经开始窃喜了,“既然你这么自信,呢要是你不能帮我解决,就要给我买好吃的,我可想吃肉了。”czl咽了一下口水。“当然可以了,怎么可能有我解决不了的斐波那契”ykc不屑的说道,ykc果然上当了。。。。。czl面无表情的说道:“呢道题是这样描述的,给你一个整数n,让你求斐波那契数列的第n项的前5位和后5位。”“哈哈,看来你是吃不上肉了,递推一下,求一求啦。”ykc仍然自信的说道。“等等,我还没说完呢,n有100000000这么大”,czl补充道。“啊。。。。哦”,ykc感觉到自己被坑了,可是来不及了,眼看就要月底了,ykc早就没钱了,聪明的你能帮帮他吗? (注意:假如该项的斐波那契数列总位数小于等于10位,呢就直接输出就好了,否则按照 前5位***后5位  的格式输出)。

以上全是扯淡,只用求前5位就好了(毕竟czl也不想刁难ykc)

输入

多组测试实例,每行一个整数n(n<=100000000)

输出

要求输出的第n项

样例输入

0 1 2 3 4 5 35 36 37 38 39 40 64 65

样例输出

0 1 1 2 3 5 92274 14930 24157 39088 63245 10233 10610 17167 思路:(标程)

首先回顾一下对数的性质:loga(b^c)=c*loga(b),loga(b*c)=loga(b)+loga(c);

假设给出一个数10234432,那么log10(10234432)=log10(1.0234432*10^7)=log10(1.0234432)+7;

log10(1.0234432)就是log10(10234432)的小数部分.

log10(1.0234432)=0.010063744    10^0.010063744=1.023443198, 要求该数的前4位,则将1.023443198*1000即可,因此,pow(10.0,x的小数部分)即可方便求出x的前几位。 参考代码: [cpp] view plain copy int  main(void)   {       double x,tmp;       while(scanf("%lf",&x)!=EOF)       {           tmp=log(x)/log(10.0);           tmp=tmp-floor(tmp);           tmp=pow(10.0,tmp);           while(tmp<1000)               tmp*=10;           printf("%d\n",(int)tmp);       }       return 0;   }   然而本题要求的第n项斐波那契数列的话,要用到斐波那契数列的通项公式

取完对数

又因为当n很大时,log10(1-((1-√5)/(1+√5))^n)->0

故原始可化为log10(an)=-0.5*log10(5.0)+((double)n)*log(f)/log(10.0); 最后取小数部分即可

#include <cstdio> #include <cstring> #include <algorithm> #include <cmath> using namespace std; typedef long long LL; LL f[30]; int main() { f[0] = 0; f[1] = f[2] = 1; for(int i = 3; i <= 25; i++) { f[i] = f[i - 1] + f[i - 2]; // printf("%d\n", f[i]); 第25位正好是最后一个五位数 } LL n; while(scanf("%lld", &n) != EOF) { if(n <= 25) { printf("%lld\n", f[n]); continue; } double x, y, z; x = -0.5 * log(5.0) / log(10.0); y = x + double(n) * log((1.0 + sqrt(5.0)) / 2.0) / log(10.0); y -= floor(y); z = pow(10.0, y); while(z < 10000) { z *= 10.0; } printf("lld\n", (LL)z); } return 0; }

转载请注明原文地址: https://www.6miu.com/read-48515.html

最新回复(0)