P2613 【模板】有理数取余 扩展欧几里得+同余定理

题目描述
给出一个有理数 c=a/b,求 c mod 19260817 的值。

输入格式
一共两行。

第一行,一个整数 a。
第二行,一个整数 b。

输出格式
一个整数,代表求余后的结果。如果无解,输出 Angry!。
说明/提示
对于所有数据,保证 0<=a ,b<=10 ^ 10001 ,且 a,b 不同时是 1926081719260817 的倍数。
乍一看题,这不是求逆元嘛,美滋滋。再看数据范围,一声晴天霹雳啊
仔细想想c是要取余的,b的逆也是要取余的,那就直接把a,b先取余吧,大数取余当然要用同余定理啦,求逆元得用拓展欧几里得,题目说存在无解的情况。费马小定理成立的前提是mod是质数,它是判断不了是否无解的。

#include<iostream>
using namespace std;
typedef long long ll;
const ll mod=19260817;
ll a,b,c,x,y;
string A,B;
ll Mod(string s)//同余定理
{
	ll res=0;
	for(int i=0;i<s.size();i++){
		res=(res*10+s[i]-'0')%mod;
	}
	return res;
}

ll exgcd(ll a,ll b,ll &x,ll &y)//拓展欧几里得
{
	if(!b) 
	{
		x=1,y=0;
		return a;
	 } 
	ll g=exgcd(b,a%b,y,x);
	y-=a/b*x;
	return g;
}
int main()
{
	cin>>A>>B;
	a=Mod(A);b=Mod(B);
	ll g=exgcd(b,mod,x,y);
	if(g==1) cout<<(a*x%mod+mod)%mod;//只有当b与mod的gcd等1时才会存在逆元
	else cout<<"Angry!";
	return 0;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章