记录编号 74935 评测结果 AAAAAAAAAA
题目名称 drei 最终得分 100
用户昵称 Gravatarcstdio 是否通过 通过
代码语言 C++ 运行时间 0.022 s
提交时间 2013-10-26 18:55:12 内存使用 0.28 MiB
显示代码纯文本
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<vector>
#include<cmath>
using namespace std;
typedef long long ll;
#define INF 0x7fffffff
vector<int> prime;
void getprime(void){
	int l[40001]={0};
	int i,j;
	for(i=2;i<=40000;i++){
		if(!l[i]){
			prime.push_back(i);
			for(j=i*i;j<=40000;j+=i) l[j]=true;
		}
	}
}
ll quickpow(ll a,ll x,ll P){//a^x mod P
	ll s=1;
	while(x){
		if(x&1) s=(s*a)%P;
		a=(a*a)%P;
		x>>=1;
	}
	return s;
}
ll ord(ll m,ll a){//使得a^x=1(mod m)的最小x,其中m是素数
	ll i;
	ll ans=INF;
	ll high=sqrt((double)m-1);
	for(i=1;i<=high;i++){
		if((m-1)%i==0){
			if(quickpow(a,i,m)==1) ans=min(ans,i);
			if(quickpow(a,(m-1)/i,m)==1) ans=min(ans,(m-1)/i);
		}
	}
	return ans;
}
ll ord(ll m,ll mk,ll a){//使得a^x=1(mod mk)的最小x,其中m是素数,mk是m的幂
	ll astep=ord(m,a);
	ll tstep=quickpow(a,astep,mk);
	ll ans=astep;
	ll temp=tstep;
	while(temp!=1){
		ans+=astep;
		temp*=tstep;temp%=mk;
	};
	return ans;
}
ll gcd(ll a,ll b){
	return b==0?a:gcd(b,a%b);
}
ll lcm(ll a,ll b){
	return a*b/gcd(a,b);
}
void work(){
	ll a,p;
	scanf("%lld%lld",&a,&p);
	if(gcd(a,p)!=1||p==1){
		printf("-1\n");
		return;
	}
	ll i,j;
	ll ans=1;
	ll temp,temk;
	for(i=0;i<prime.size();i++){
		if(prime[i]>p) break;
		temp=1;temk=0;
		while(p%prime[i]==0){
			temk++,temp*=prime[i];
			p/=prime[i];
		}
		if(temk){
			ans=lcm(ans,ord(prime[i],temp,a));
		}
	}
	if(p>1) ans=lcm(ans,ord(p,a));
	printf("%lld\n",ans);
}
int main(){
	freopen("drei.in","r",stdin);
	freopen("drei.out","w",stdout);
	getprime();
	ll T;
	scanf("%lld",&T);
	for(ll i=1;i<=T;i++) work();
	return 0;
}