A/B 拓展欧几里得算法

最新推荐文章于 2022-04-29 17:35:33 发布

quest_ni

最新推荐文章于 2022-04-29 17:35:33 发布

阅读量235

点赞数 1

分类专栏： c语言文章标签：数论 c/c

本文链接：https://blog.csdn.net/quest_ni/article/details/81105225

版权

c语言专栏收录该内容

8 篇文章 0 订阅

订阅专栏

要求(A/B)%9973，但由于A很大，我们只给出n(n=A%9973)(我们给定的A必能被B整除，且gcd(B,9973) = 1)。

Input

数据的第一行是一个T，表示有T组数据。
每组数据有两个数n(0 <= n < 9973)和B(1 <= B <= 10^9)。

Output

对应每组数据输出(A/B)%9973。

Sample Input

2
1000 53
87 123456789

Sample Output

7922
6060

思路：因为n=A%9973，n=A-A%9973*9973；而且x=A/B，A=Bx，所以n=Bx-A/9973*9973，Bx-9973y=n。要求（A/B）%9973，只需要求出x%9973就好。

利用扩展欧几里德算法，可求出gcd(B,9973)=Bx1+9973y1=1的x1。由题知gcd(B,9973)=1；所以等式两边同乘以n，得B(n*x1)-9973(-n*y1)=n。可得n*x1就是B*x-9973*y=n的解，即x=n*x1。（注意x可能为负值，需要进行x=(x%9973+9973)%9973）。

拓展欧几里得算法模板：

ll x,y;
ll gcd(ll a,ll b)
{
    if(b==0)
    {
        x=1;y=0;return a;
    }
    ll h=gcd(b,a%b);
    ll t=x;
    x=y;
    y=t-(a/b)*y;
    return h;
}

代码如下：

#include<bits/stdc++.h>
using namespace std;
long long int x,y;
void gcd(int a,int b)
{
    int m;
    if(b==0)
    {
        x=1;y=0;return;
    }
    else
    {
        gcd(b,a%b);
        m=x;x=y;
        y=m-(a/b)*y;
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        long long int n,b;
        scanf("%lld %lld",&n,&b);
        gcd(b,9973);
        if(x<0)x+=9973;
        x*=n;x%=9973;
        printf("%lld\n",x);
    }
    return 0;
}