HDU 1576 扩展欧几里得&&暴力

HDU 1576

A - A

Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u

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

这个题，题目很机智的给了一个很小的MOD，所以我们可以暴力

姿势如下：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <string>
#include <cstring>
#include <cmath>

using namespace std;
const int maxn=1e5+10;
const int INF=0x3f3f3f3f;
int main()
{
    int t;
    long long n,b;
    cin>>t;
    while(t--){
        cin>>n>>b;
        for(int i=1;i<9973;i++)
            if((i*b-n)%9973==0){
                cout<<i<<endl;
                break;
            }
    }

}

换一种正常的姿势，即扩展欧几里得：

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <string>
#include <cstdlib>
using namespace std;

const long long mod=9973;
void exgcd(long long a,long long &x,long long b,long long  &y){
    if(!b){
        x=1;
        y=0;
    }
    else{
        exgcd(b,y,a%b,x);
        y-=x*(a/b);
    }
}
int main(){
    long long t,N,B,x,y;
    cin>>t;
    while(t--){
        cin>>N>>B;
        exgcd(B,x,mod,y);
        long long ans=((x%mod)*(N%mod))%mod;
        while(ans<0) ans+=mod;
        cout<<ans<<endl;

    }
    return 0;
}