太菜了 比赛的时候下面罗列的方法大致方向都想对了 但都差了一点 没敢写 不够自信
第一种方法:找循环节法 注意循环节不一定就从第一个开始
#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");
}
}