求斐波那契数列模n的循环节长度

HDU 3977
一、

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cmath>
#include<cstring>
#define ll long long
using namespace std;
const int maxn=1000006;
int prime[maxn];
bool ha[maxn]={false};
int tot=0;
void Prime(void)
{
    ha[1]=true;
    for(int i=2;i<maxn;i++)
    {
        if(!ha[i]) prime[tot++]=i;
        for(int j=0;j<tot&&(ll)prime[j]*i<maxn;j++)
        {
            ha[i*prime[j]]=true;
            if(i%prime[j]==0) break;
        }
    }
    return ;
}
int p[20];
int e[20];
int cnt=0;
void only(int n)
{
    cnt=0;
    memset(p,0,sizeof(p));
    memset(e,0,sizeof(e));

    for(int i=0;i<tot&&prime[i]*prime[i]<=n;i++)
    {
        if(n%prime[i]) continue;
        p[++cnt]=prime[i];
        while(n%prime[i]==0)
        {
            e[cnt]++;
            n/=prime[i];
        }
    }
    if(n>1) p[++cnt]=n,e[cnt]=1;
    return ;
}

ll mypow(ll a,ll b,ll p)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%p;
        a=a*a%p;
        b>>=1;
    }
    return ans;
}
ll mypow(ll a,ll b)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=ans*a;
        a=a*a;
        b>>=1;
    }
    return ans;
}

ll gcd(ll a,ll b)
{
    if(b==0) return a;
    else  return gcd(b,a%b);
}

ll lcm(ll a,ll b)
{
    return a/gcd(a,b)*b;
}

ll G(int p)
{
    int a=1,b=2,c=(2)%p;
    int pos=3;
    while(c) c=(a+b)%p,a=b,b=c,pos++;
    ll aa=a;
    for(int i=1;i<p;i++)
    {
        if(mypow(aa,i,p)==1) return (ll)pos*i;
    }
    return 0;
}


ll minn(int pi)
{
    only(pi);
    ll ans=1;
    for(int i=1;i<=cnt;i++)
    {
        ll k=G(p[i])*mypow(p[i],e[i]-1);
        ans=lcm(ans,k);
    }
    return ans;
}


int main(void)
{
    int t;
    scanf("%d",&t);
    Prime();
    for(int k=1;k<=t;k++)
    {
        int p;
        scanf("%d",&p);
        ll m=minn(p);

        printf("Case #%d: %lld\n",k,m);
    }
    return 0;
}

二、

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string>
#include<cmath>
#include<cstring>
#define ll long long
using namespace std;
const int maxn=1000006;
int prime[maxn];
bool ha[maxn]={false};
int tot=0;
void Prime(void)
{
    ha[1]=true;
    for(int i=2;i<maxn;i++)
    {
        if(!ha[i]) prime[tot++]=i;
        for(int j=0;j<tot&&(ll)prime[j]*i<maxn;j++)
        {
            ha[i*prime[j]]=true;
            if(i%prime[j]==0) break;
        }
    }
    return ;
}
int p[20];
int e[20];
int cnt=0;
void only(int n)
{
    cnt=0;
    memset(p,0,sizeof(p));
    memset(e,0,sizeof(e));

    for(int i=0;i<tot&&prime[i]*prime[i]<=n;i++)
    {
        if(n%prime[i]) continue;
        p[++cnt]=prime[i];
        while(n%prime[i]==0)
        {
            e[cnt]++;
            n/=prime[i];
        }
    }
    if(n>1) p[++cnt]=n,e[cnt]=1;
    return ;
}

ll mypow(ll a,ll b,ll p)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%p;
        a=a*a%p;
        b>>=1;
    }
    return ans;
}

ll mypow(ll a,ll b)
{
    ll ans=1;
    while(b)
    {
        if(b&1) ans=ans*a;
        a=a*a;
        b>>=1;
    }
    return ans;
}

struct node
{
    ll a[2][2];
};

node mult(node a,node b,int p)
{
    node c;
    memset(&c,0,sizeof(c));
    for(int i=0;i<2;i++)
    {
        for(int j=0;j<2;j++)
        {
            for(int k=0;k<2;k++)
            {
                c.a[i][j]=(c.a[i][j]+a.a[i][k]*b.a[k][j]%p)%p;
            }
        }
    }
    return c;
}



ll getf(ll b,int p)
{
    node ans,res;
    ans.a[0][0]=1,ans.a[0][1]=0,ans.a[1][0]=0,ans.a[1][1]=1;
    res.a[0][0]=1,res.a[0][1]=1,res.a[1][0]=1,res.a[1][1]=0;
    b--;
    while(b)
    {
        if(b&1) ans=mult(ans,res,p);
        res=mult(res,res,p);
        b>>=1;
    }
    return ans.a[0][0]%p;
}

bool check(int a,int p)
{
    return mypow(a,(p-1)>>1,p)==1;
}

ll fac[2008];
int  _cnt=0;
void f(ll n)
{
    memset(fac,0,sizeof(fac));
    _cnt=0;
    for(int i=1;(ll)i*i<=n;i++)
    {
        if(n%i) continue;
        if(i*i!=n) fac[_cnt++]=i,fac[_cnt++]=n/i;
        else fac[_cnt++]=i;
    }
    return ;
}

ll gcd(ll a,ll b)
{
    if(b==0) return a;
    else  return gcd(b,a%b);
}

ll lcm(ll a,ll b)
{
    return a/gcd(a,b)*b;
}

ll G(int p)
{
    if(p==2) return 3;
    if(p==3) return 8;
    if(p==5) return 20;
    ll sum=0;

    if(check(5,p))
        sum=p-1;
    else sum=2*(p+1);

    f(sum);
    sort(fac,fac+_cnt);

    for(int i=0;i<_cnt;i++)
    {
        if(getf(fac[i],p)==0&&getf(fac[i]+1,p)==1)
            return fac[i];
    }
    return 0;
}


ll minn(int pi)
{
    if(pi==1) return 1;
    only(pi);
    ll ans=1;
    for(int i=1;i<=cnt;i++)
    {
        ll k=G(p[i])*mypow(p[i],e[i]-1);
        ans=lcm(ans,k);
    }
    return ans;
}


int main(void)
{
    int t;
    scanf("%d",&t);
    Prime();
    for(int k=1;k<=t;k++)
    {
        int p;
        scanf("%d",&p);
        ll m=minn(p);

        printf("Case #%d: %lld\n",k,m);
    }
    return 0;
}