洛谷P2613 【模板】有理数取余 题解
题目链接:P2613 【模板】有理数取余
题意:给出 c = a b c=\frac{a}{b} c=ba ,求 a b m o d 19260817 \frac{a}{b} \bmod 19260817 bamod19260817
0 ≤ a ≤ 1 0 10001 , 1 ≤ b ≤ 1 0 10001 0 \le a \le 10^{10001},1\le b \le 10^{10001} 0≤a≤1010001,1≤b≤1010001
令 m = 19260817 m=19260817 m=19260817
注意这里模意义下的有理数其实就是
a
b
(
m
o
d
m
)
=
a
b
−
1
(
m
o
d
m
)
\frac{a}{b} \pmod m = ab^{-1} \pmod m
ba(modm)=ab−1(modm)
简单推推柿子吧
显然有
b
c
≡
a
(
m
o
d
m
)
bc \equiv a \pmod{m}
bc≡a(modm)
令
b
0
=
b
m
o
d
m
,
a
0
=
a
m
o
d
m
b_0=b \bmod m,a_0 = a \bmod m
b0=bmodm,a0=amodm
当
b
0
≠
0
b_0 \ne 0
b0=0 即
m
∤
b
m \nmid b
m∤b 时,有
c
≡
a
0
b
0
−
1
(
m
o
d
m
)
c \equiv a_0 b_0^{-1} \pmod{m}
c≡a0b0−1(modm)
那么exgcd一下就出来了
当 b 0 = 0 b_0=0 b0=0 即 m ∣ b m \mid b m∣b 时,分类讨论
- 当 m ∣ a m \mid a m∣a 时, ∀ c ∈ Z , b c ≡ a ( m o d m ) \forall c \in \Z,bc \equiv a \pmod{m} ∀c∈Z,bc≡a(modm) 恒成立,题目保证了没有这种情况
- 当 m ∤ a m\nmid a m∤a 时, b c ≡ a ( m o d m ) bc \equiv a \pmod{m} bc≡a(modm) 无解
因此特判 b 0 = 0 b_0 = 0 b0=0 的情况即可
注意读入 a , b a,b a,b 的时候要一边读一边取模
代码如下
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define INF 0x3f3f3f3f3f3f3f3f
#define gc() getchar()
const int mod=19260817;
int readint()
{
char ch=gc();int res=0;
while(!isdigit(ch))ch=gc();
while(isdigit(ch))
{
res=(res<<1)+(res<<3)+(ch^48);
res%=mod;
ch=gc();
}
return res;
}
int exgcd(int a,int b,int &x,int &y)
{
int d=a;
if(!b)x=1,y=0;
else d=exgcd(b,a%b,y,x),y-=a/b*x;
return d;
}
signed main()
{
ios::sync_with_stdio(0);
cin.tie(0);cout.tie(0);
// freopen("check.in","r",stdin);
// freopen("check.out","w",stdout);
int a,b;
a=readint();b=readint();
if(!b)return cout << "Angry!" << endl,0;
int x,y;
exgcd(b,mod,x,y);
x=(x%mod+mod)%mod;
cout << a*x%mod << endl;
return 0;
}
转载请说明出处