hdu 1576 A/B

最新推荐文章于 2021-02-24 00:30:01 发布

o詤o訁o

最新推荐文章于 2021-02-24 00:30:01 发布

阅读量726

点赞数

分类专栏： ∟__ACM_数论文章标签： acm

本文链接：https://blog.csdn.net/u010304217/article/details/21470573

版权

∟__ACM_数论专栏收录该内容

8 篇文章 0 订阅

订阅专栏

A/B

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1576

Problem Description

要求(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

由题意可知n=A%9973,所以n=A-(A/9973)*9973。令x=A/B,则A=Bx，n+(A/9973)*9973=Bx,n=Bx-(A/9973)*9973,此时把（A/9973）看为一个整体y，原式为n=Bx-9973*y；

所以可以用扩展欧几里德算法计算出x，但是因为9973为素数，所以gcd(B,9973)=1,此时计算出的是1=Bx-9973*y；所以等译n是的x的解为n*x；

代码

#include<stdio.h>

int f(int n,int m,int &x,int &y)
{
    int x0,x1,y0,y1;
    x0=1;x1=0;
    y0=0;y1=1;
    x=0;y=1;
    int r=m%n;
    int q=(m-r)/n;
    while(r)
    {
        x=x0-q*x1;
        y=y0-q*y1;
        x0=x1;
        y0=y1;
        x1=x;
        y1=y;
        m=n;
        n=r;
        r=m%n;
        q=(m-r)/n;
    }
    return n;
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        int n,B;
        int x,y,z;
        scanf("%d%d",&n,&B);
        z=f(9973,B,x,y);
        //printf("z=%d\nx=%d\ny=%d\n",z,x,y);
        x*=n;
        printf("%d\n",(x%9973+9973)%9973);
    }
    return 0;
}