HDU——1576 A/B

　　　　　　　　　　　　　　A/B

　　　　　　　　　　Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
　　　　　　　　　　Total Submission(s): 6412    Accepted Submission(s): 5069

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

7922

6060

Author

xhd

Source

HDU 2007-1 Programming Contest 

 

思路：

(a/b)mod p=?
定义 c为b在mod p意义下的逆元
=a*c mod p = （a mod p * c mod p）mod p

这样我们就有把这道题转换成裸地扩展欧几里得了。。。

代码：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int n,m,t,a,b,x,y,gcd,ans;
int read()
{
    int x=0,f=1; char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
int exgcd(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1,y=0;
        return a;
    }
    int r=exgcd(b,a%b,x,y),tmp;
    tmp=x,x=y,y=tmp-a/b*y;
    return r;
}
int main()
{
    t=read();
    while(t--)
    {
        n=read(),m=read();
        a=m;b=9973;
        gcd=exgcd(a,b,x,y);
        while(x<=0) x+=b/gcd; 
        ans=(n%b*x%b)%b;
        printf("%d\n",ans);
    }
    return 0;
}

 

转载于:https://www.cnblogs.com/z360/p/7327322.html