斐波那契数列找循环节

首先贴一个大佬的广义方法，，数学真的diao。。

广义斐波那契数列的循环节

暴力求法， 

一般情况下这样求long long 范围内的也就够了。

#include<bits/stdc++.h>

#define ll long long

using namespace std;

const int maxn=1e5+10;// 自调

ll f[maxn];

int main()

{

    int mod;

    cin>>mod;

    f[0]=0;

    f[1]=1;

    f[2]=1;

    for(int i=3;i<=maxn;i++)

    {

        f[i]=(f[i-1]+f[i-2])%mod;

        if(f[i]==1&&f[i-1]==0)

        {

            printf("循环节为 %d\n",i-1);// 循环节 i-1;

            break;

        }

    }

    return 0;

}

来一道题：

问题描述

　　某天，HWD老师开始研究兔子，因为他是个土豪

　　，所以他居然一下子买了一个可以容纳10^18代兔子的巨大笼子（好像比我们伟大的地球母亲大一点点？），并开始研究它们，之后，他了解了兔子的繁衍规律：即fibonacci数列。

　　兔子繁殖了n代后，HWD老师很开心。

　　但是，HWD老师有密集恐惧症，所以，他只能去卖了兔子，他找到了一个好的雇主，但是这个雇主有强迫症，他只每次收购1007只兔子，HWD老师为了避免自己的密集恐惧症，要尽量多的卖了兔子。

　　但是即便是密集恐惧症，也打击不了HWD老师研究兔子的决心，他数着数着自己剩下的兔子……

输入格式

　　HWD老师让兔子繁衍了几代（一个整数，没有其他字符）。

输出格式

　　HWD老师剩余（残余？）的兔子（一个整数，忽略行尾回车及空格）。

样例输入

1

样例输出

1

数据规模和约定

　　兔子的总量最大时小于HWD老师笼子的大小。

　　f[1]=1，f[2]=1，f[3]=2 ……

因为10^18，所以随便开了一个循环节。

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int main()
{
    ios::sync_with_stdio(false);
    ll n;
    cin>>n;
    if(n==1||n==2)
    {
        puts("1");
        return 0;
    }
    int mod=1007;
    int q=108;
    ll a,b,c;
    n%=q;
    for(int i=1;i<=n;i++)
    {
        if(i==1)
        {
            a=1;
            continue;
        }
        if(i==2)
        {
            b=1;
            continue;
        }
        c=(a+b)%mod;
        a=b;
        b=c;
    }
    cout<<b<<endl;
    return 0;
}