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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值