HDU 1576 A/B（数论-简单逆元）

A/B

 

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2432    Accepted Submission(s): 1761
点击打开题目链接

 

 

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、则称作a为b模m的逆元，b为a模m的逆元；逆元求解方法：1、如果m为素数，且a不是m的倍数，那么就简单了，由费马小定理可知,,即为a模m的逆元；2、如果m不是素数，则可以借用欧几里得算法，在求出最大公约数的同时求出x，y分别是a模b的逆元和b模a的逆元；

 

逆元详见：http://blog.csdn.net/acdreamers/article/details/8220787

                   http://blog.csdn.net/cqlf__/article/details/7953039

考虑（a/b）mod m；

 

代码：

 

#include <iostream>
#include <vector>
#include <queue>
#include <set>
#include <map>
#include <stack>
#include <stdlib.h>
#include <string.h>
#include <stdio.h>
#define max(a,b) (a>b?a:b)
#define min(a,b) (a<b?a:b)
#define mod 9973
using namespace std;
long long exc_gcd(long long a,long long b,long long &x,long long &y)
{
    if(!b)
    {
        x=1,y=0;
        return a;
    }
    long long g=exc_gcd(b,a%b,x,y);
    long long t=x;
    x=y;
    y=t-(a/b)*y;
    return g;
}
int main()
{
    int t,n;
    long long b;
    while(~scanf("%d",&t))
    {
        while(t--)
        {
            scanf("%d%I64d",&n,&b);
            long long x=0,y=0;
            long long b1=exc_gcd(b,mod,x,y);
            if(x<0) x+=mod;
            long long ans=(n*(x%mod)+mod)%mod;
            printf("%I64d\n",ans);
        }
    }
    return 0;
}