【Fibonacci 数列前四位&后四位】hdu 3117

最新推荐文章于 2020-10-07 16:37:01 发布

leolin_

最新推荐文章于 2020-10-07 16:37:01 发布

阅读量884

点赞数

分类专栏：数论文章标签：百度

本文链接：https://blog.csdn.net/leolin_/article/details/6896764

版权

数论专栏收录该内容

11 篇文章 0 订阅

订阅专栏

这题的前四位是参考别人的思路的，我发现求前几位通常做法是取对数！

后四位不用说，当n>=40时，位数>8。另外，后四位的周期是15000，就这样打表求出后四位。

至于前四位，就要用到Fibonacci 数列的通项公式，我看了百度百科后才知道这个公式==，另外，有趣的是当n无穷大，前一项与后一项之比逼近黄金分割0.618。。。orz。话说回来，an可以表示为t*10^k（t为>1的浮点数），取对数log10(an)= log10(t) + k，其中0<log10(t) <1，所以取完对数减去整数部分就是log10(t) ，然后pow(10,log10(t) )还原t。注意，这里有一个技巧，就是当n>=40时通式中括号里后一半趋向于无穷小，小到可以忽略。

上代码：

#include <map>
#include <set>
#include <list>
#include <queue>
#include <deque>
#include <stack>
#include <string>
#include <time.h>
#include <cstdio>
#include <math.h>
#include <iomanip>
#include <cstdlib>
#include <limits.h>
#include <string.h>
#include <iostream>
#include <fstream>
#include <algorithm>
using namespace std;

#define LL long long
#define MIN INT_MIN
#define MAX INT_MAX
#define PI acos(-1.0)
#define N 2
#define FRE freopen("input.txt","r",stdin)
#define FF freopen("output.txt","w",stdout);
#define MOD 10000
int Fib[45] = {0,1,1,2,3};
int last[20000] = {0,1,1,2,3};
double gao = (1 + sqrt(5.0)) / 2.0;
int main(){
    int i,j;
    for(i = 5; i < 15000; i++){
        last[i] = (last[i-1] + last[i-2]) % MOD;
    }
    for(i = 5; i < 40; i++)Fib[i] = Fib[i-1] + Fib[i-2];
    int n;
    while(scanf("%d",&n)!=EOF){
        int i,j;
        if(n<40){
            printf("%d\n",Fib[n]);
        }
        else{
            double ans = n * log10(gao) - log10(5.0) / 2.0;
            ans -= (int)ans;
            ans = pow(10,ans);
            while(ans < 1000){
                ans *= 1000;
            }
            printf("%d...%4.4d\n",(int)ans,last[n%15000]);
        }
    }
    return 0;
}