【模板】有理数取余（exgcd）

最新推荐文章于 2024-09-25 13:54:13 发布

green qwq

最新推荐文章于 2024-09-25 13:54:13 发布

阅读量580

点赞数 28

分类专栏：好题文章标签：算法

本文链接：https://blog.csdn.net/m0_60738889/article/details/138684520

版权

好题专栏收录该内容

235 篇文章 0 订阅

订阅专栏

【模板】有理数取余

题目描述

给出一个有理数 $c=\frac{a}{b}$ ，求 $\bmod 19260817$ 的值。

这个值被定义为 $bx\equiv a\pmod{19260817}$ 的解。

输入格式

一共两行。

第一行，一个整数 $a$ 。
第二行，一个整数 $b$ 。

输出格式

一个整数，代表求余后的结果。如果无解，输出 Angry!。

样例 #1

样例输入 #1

233
666

样例输出 #1

18595654

提示

对于所有数据，保证 $0\leq a \leq 10^{10001}$ ， $\leq b \leq 10^{10001}$ ，且 $a, b$ 不同时是 $19260817$ 的倍数。

思路

首先，这个首先读入就是个麻烦事情，我们可以在读入的时候先求模，肯定不会对结果造成影响。
其次，这道题要求 $bx\equiv a\pmod{19260817}$ ，这个我们可以转化成 $b x + 19260817 y = a$ 其中 $g c d (b, 19260817) = 1$ ，如果不等于，则无解。
最后，剩下的都是扩展欧几里得的内容了（说实话，因为19260817是素数，这道题还可以用费马定理（求乘法逆元的））。

代码

//可以在读入的时候%mod
#include<iostream>

using namespace std;

const int mod = 19260817;

typedef long long LL;

LL a,b;

LL getint(){
    string aa;
    cin>>aa;
    LL ans=0;
    for(int i=0;i<aa.size();i++){
        LL t=aa[i]-'0';
        ans=(ans*10+t)%mod;
    }
    return ans;
    
}

int exgcd(int a,int b,LL &x,LL &y){
    if(!b){
        x=1,y=0;
        return a;
    }
    
    int d=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}

int main(){
    a=getint();
    b=getint();
    LL x,y;
    int d=exgcd(b,mod,x,y);
    //bx-mody=a
    
    if(a%d){
        puts("Angry!");
        return 0;
    }
    x*=a/d;
    int t=mod/d;
    cout<<(x%t+t)%t;
    return 0;
}