UVA 11582

思路:直接暴力肯定会超时,先用快速幂计算出a^b%n的值,因为mod为n时最多有n种余数,而且f[i]由前两项决定,所以周期至多为n*n,当f[i]==f[1]&&f[i-1]==f[0]时可以找到周期。

注意:我在每一次进行时用memset初始化f数组结果TLE,去掉后就AC了。。

#include<bits/stdc++.h>
using namespace std;
#define inf 0x3f3f3f3f
#define ull unsigned long long
#define fo freopen("in.txt","r",stdin)
#define fc fclose(stdin)
#define fu0(i,n) for(i=0;i<n;i++)
#define fu1(i,n) for(i=1;i<=n;i++)
#define fd0(i,n) for(i=n-1;i>=0;i--)
#define fd1(i,n) for(i=n;i>0;i--)
#define mst(a,b) memset(a,b,sizeof(a))
#define sd(n) scanf("%d",&n)
#define sdd(n,m) scanf("%d %d",&n,&m)
#define ss(s) scanf("%s",s)
#define sddd(n,m,k) scanf("%d %d %d",&n,&m,&k)
#define pans() printf("%d\n",ans)
#define aull(a) a.begin(),a.end()
#define sc(c) scanf("%c",&c)
#define we(a) while(scanf("%d",&a)!=EOF)
const int maxn=200005;
ull a,b;
int n;
int f1(ull a,ull b,int T)
{
    int t=1;
    while(b)
    {
        if(b & 1)
        {
            t=t*a%T;
            b--;
        }
        a=a*a%T;
        b>>=1;
    }
    return t;
}
int f[1000010];
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int T=1;
        //memset(f,0,sizeof(f));
        cin>>a>>b>>n;
        if(n==1||a==0)
        {
            cout<<0<<endl;
            continue;
        }
        f[0]=0;
        f[1]=1%n;
        for(int i=2; i<=n*n+100; i++)
        {
            f[i]=(f[i-1]+f[i-2])%n;
            if(f[i]==f[1]&&f[i-1]==f[0])
            {
                T=i-1;
                break;
            }
        }
        int gg=f1(a%T,b,T);
        cout<<f[gg]<<endl;
    }
    return 0;
}

阅读更多
个人分类: 递推找规律 数论
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭