有理数取余（快速幂，费马定理或扩展欧几里得算法）

最新推荐文章于 2024-05-23 20:19:57 发布

慧狗仔

最新推荐文章于 2024-05-23 20:19:57 发布

阅读量616

点赞数 21

文章标签： c++ 算法

本文链接：https://blog.csdn.net/2301_80065123/article/details/138663747

版权

【模板】有理数取余（快速幂，费马定理或扩展欧几里得算法）

题目描述

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

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

输入格式

一共两行。

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

输出格式

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

样例 #1

样例输入 #1

233
666

样例输出 #1

18595654

提示说明

在这里插入图片描述

特别注意： $b*x\equiv 1\pmod{p}$ (只对于p是质数)

由费马定理得： $b^{p-1}\equiv 1\pmod{p}$ (p为质数)
即： $b*b^{p-2}\equiv 1\pmod{p}$
得： $b^{p-2} \pmod{p}$ ;

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

代码内容

快速幂，费马定理

// #include <iostream>
// #include <algorithm>
// #include <cstring>
// #include <stack>//栈
// #include <deque>//队列
// #include <queue>//堆/优先队列
// #include <map>//映射
// #include <unordered_map>//哈希表
// #include <vector>//容器，存数组的数，表数组的长度
#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
const ll p=19260817;

//快读（分字符读入）
ll getint()
{
    ll res=0,ch=getchar();
    while(!isdigit(ch)&&ch!=EOF)
        ch=getchar();
    while(isdigit(ch))
    {
        res=(res<<3)+(res<<1)+(ch-'0');
        res%=p;//取模p;
        ch=getchar();
    }
    return res;
}

ll qmi(ll a,ll b,ll p)
{
    //当b = 0,p = 1时
    //①不 % p ,结果是 1
    //② % p ,结果是 0(正确答案)
    ll res=1%p;//注意:5 0 1
    while(b)
    {
        if(b&1) res=res*a%p;
        a=a*a%p;
        b>>=1;//右移 1 位
    }
    return res;
}

int main()
{
    ll a,b;//快读（分字符读入）
    a=getint();
    b=getint();
    
    if(!b)
    {
        cout<<""<<endl;
        return 0;
    }

    //b * x ≡ 1 (mod p)//乘法逆元
    //b^{p-1} ≡ 1 (mod p)//费马定理
    //b * b^{p-2} ≡ 1 (mod p)
    //(b*x)模 p 余 1
    //x = b^{p-2} (mod p);
    //res = x

    //(p 为质数,即 a 与模数 p 互质)
    ll res=a*qmi(b,p-2,p);
    //非负整数
    cout<<(res%p+p)%p<<endl;
    return 0;
}

扩展欧几里得算法

// #include <iostream>
// #include <algorithm>
// #include <cstring>
// #include <stack>//栈
// #include <deque>//队列
// #include <queue>//堆/优先队列
// #include <map>//映射
// #include <unordered_map>//哈希表
// #include <vector>//容器，存数组的数，表数组的长度
#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
const ll p=19260817;

//快读（分字符读入）
ll getint()
{
    ll res=0,ch=getchar();
    while(!isdigit(ch)&&ch!=EOF)
        ch=getchar();
    while(isdigit(ch))
    {
        res=(res<<3)+(res<<1)+(ch-'0');
        res%=p;//取模p;
        ch=getchar();
    }
    return res;
}

//扩展欧几里得算法(裴蜀定理)
//根据费马定理，任意正整数a, b都存在整数x, y
//使得ax + by = gcd(a, b)
void exgcd(ll a,ll b,ll &x,ll &y)
{
    if(!b)
    {
        x=1,y=0;
        return;
    }
    exgcd(b,a%b,y,x);
    y-=a/b*x;
}

int main()
{
    ll a,b;//快读（分字符读入）
    a=getint();
    b=getint();
    
    if(!b)
    {
        cout<<""<<endl;
        return 0;
    }

    ll x,y;
    //扩展欧几里得算法
    //b*x≡1(mod p)
    //b*x=p*y+1
    /*b*x+p*y=1*/
    exgcd(b,p,x,y);
    
    //(p 为质数,即 a 与模数 p 互质)
    //ax + by = gcd(a, b) = d = 1
    //ll res=a*x*(1/d)=a*x
    ll res=a*x;
    //非负整数
    cout<<(res%p+p)%p<<endl;
    return 0;
}