题目描述
给出一个有理数 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;
}