比赛 ZLXSCDay1 评测结果 AAAAAAAAAA
题目名称 多-有丝分裂 最终得分 100
用户昵称 frontier 运行时间 0.319 s
代码语言 C++ 内存使用 0.32 MiB
提交时间 2016-03-18 20:16:40
显示代码纯文本
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<map>
using namespace std;
typedef long long ll;
ll qmul(ll a,ll b,ll p){
	ll ans=1;
	for(;b;b>>=1,a=a*a%p)if(b&1)ans=ans*a%p;
	return ans;
}
ll baby_giant_step(ll a,ll b,ll p){
	ll m=sqrt(p+0.5);
	map<int,int>mp;
	for(int i=0;i<m;i++)
	mp[b]=i,b=b*a%p;
	ll x;
	for(int i=1;i<=m+1;i++){
		x=qmul(a,i*m,p);
		if(mp.count(x))return i*m-mp[x];
	}
	return -1;
}
ll gcd(ll a,ll b){
	return b?gcd(b,a%b):a;
}
ll exgcd(ll a,ll b,ll &d,ll &x,ll &y){
	if(b)exgcd(b,a%b,d,y,x),y-=x*(a/b);
	else d=a,x=1,y=0;
}
ll inv(ll a,ll p){
	ll d,x,y;
	exgcd(a,p,d,x,y);
	return d==1?(x+p)%p:-1;
}
ll log_mod(ll a,ll b,ll p){
	ll g=gcd(a,p);
	a%=p;b%=p;
	if(g==1)return baby_giant_step(a,b,p);
	else if(b%g)return -1;
	else{
		ll ans=log_mod(a,b/g*inv(a/g,p/g),p/g);
		if(ans<0)return -1;
		return ans+1;
	}
}
int main(){
	freopen("mitotic_division.in","r",stdin);
	freopen("mitotic_division.out","w",stdout);
	int T;scanf("%d",&T);
	ll a,b,c;
	while(T--){
		scanf("%lld%lld%lld",&a,&b,&c);
		ll ans=log_mod(a,b,c);
		if(ans<0)puts("no solution");
		else printf("%lld\n",ans);
	}
	return 0;
}