hdu 5690 多种方法实现

太菜了 比赛的时候下面罗列的方法大致方向都想对了 但都差了一点 没敢写 不够自信

第一种方法：找循环节法 注意循环节不一定就从第一个开始

#include<cstdio>
#include<cstring>
#define ll __int64
using namespace std;
int vis[10005];
int num[10005];
int main()
{
    int T,kase=1,tmp,cnt,i,fi;
    ll x,m,k,c;
    scanf("%d",&T);
    while(T--)
    {
        memset(vis,0,sizeof(vis));
        scanf("%I64d%I64d%I64d%I64d",&x,&m,&k,&c);
        tmp=0;cnt=0;
        for(i=1;i<=m;i++)
        {
            tmp=(tmp*10+x)%k;
            if(i==m) break;
            if(!vis[tmp]) {num[++cnt]=tmp;vis[tmp]=1;}
            else break;
        }
        if(i==m)
        {
            if(tmp==c) printf("Case #%d:\nYes\n",kase++);
            else printf("Case #%d:\nNo\n",kase++);
            continue;
        }
        if(cnt==1)
        {
            if(x==c) printf("Case #%d:\nYes\n",kase++);
            else printf("Case #%d:\nNo\n",kase++);
            continue;
        }
        for(i=1;i<=cnt;i++)
            if(num[i]==tmp) break;
        fi=i;
        cnt=cnt-fi+1;
        m=(m-fi+1)%cnt;
        if(m==0)m=cnt;
        m=fi+m-1;
        x=num[m];
        if(x==c) printf("Case #%d:\nYes\n",kase++);
        else printf("Case #%d:\nNo\n",kase++);
    }
    return 0;
}

第二种方法：暴力算

#include<set>
#include<map>
#include<cmath>
#include<queue>
#include<vector>
#include<cstdio>
#include<algorithm>
#define ll __int64
using namespace std;
const int N=1e6;
int main()
{
    int T,kase=0;
    ll x,m,k,c;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%I64d%I64d%I64d%I64d",&x,&m,&k,&c);
        ll tmp=0,tmp2=1,ans=0;
        for(int i=0;i<10000;i++)
        {
            tmp=(tmp*10+x)%k;
            tmp2=(tmp2*10)%k;
        }
        for(int i=0;i<m/10000;i++)
            ans=(ans*tmp2+tmp)%k;
        for(int i=0;i<m%10000;i++)
            ans=(ans*10+x)%k;
        if(ans==c) printf("Case #%d:\nYes\n",++kase);
        else printf("Case #%d:\nNo\n",++kase);
    }
}

第三种方法：快速幂

#include<cstdio>
using namespace std;
typedef long long int ll;
ll quickmod(ll a,ll b,ll c)
{
    ll temp=1;
    while(b>0)
    {
        if(b%2==1)
            temp=temp*a%c;
        b/=2;
        a=a*a%c;
    }
    return temp%c;
}
int calc(ll x,ll m,ll k){
    if(m==1LL) return x;
    int ret=calc(x,m>>1LL,k);
    ret=(ret*quickmod(10,m>>1,k)+ret)%k;
    return (m&1LL)?((ret*10+x)%k):ret;
}
int main(){
    int T,kase=0;
    ll x,m,k,c;
    scanf("%d", &T);
    while(T--)
    {
        scanf("%I64d%I64d%I64d%I64d",&x,&m,&k,&c);
        printf("Case #%d:\n%s\n",++ kase,calc(x,m,k)==c?"Yes":"No");
    }
    return 0;
}

本条转自：http://blog.csdn.net/Miracle_ma/article/details/51472075

第四种方法：矩阵快速幂 其实类似上面的快速幂

#include<cstdio>
int T,kase,x,p,c;
long long m;
struct mat
{
    int a[2][2];
    mat()
    {
        for(int i=0;i<2;i++)
            for(int j=0;j<2;j++)
                a[i][j]=0;
    }
    mat operator*(const mat&b)
    {
        mat 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]+=a[i][k]*b.a[k][j])%=p;
        return c;
    }
}A,B,C;
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%I64d%d%d",&x,&m,&p,&c);
        A=B=C=mat();
        A.a[0][0]=10;A.a[0][1]=1;
        A.a[1][0]=0;A.a[1][1]=1;
        B.a[1][0]=x;
        C.a[0][0]=C.a[1][1]=1;
        for(;m;m>>=1,A=A*A)if(m&1)C=C*A;
        C=C*B;
        int tmp=C.a[0][0]%p;
        printf("Case #%d:\n%s\n",++kase,tmp==c?"Yes":"No");
    }
}