hdoj 1576

最新推荐文章于 2019-01-27 03:07:20 发布

yyzs

最新推荐文章于 2019-01-27 03:07:20 发布

阅读量205

点赞数

分类专栏：数论

本文链接：https://blog.csdn.net/cdywy/article/details/53483393

版权

数论专栏收录该内容

2 篇文章 0 订阅

订阅专栏

//1Y真是爽啊

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

分析：根据题意a=b*x a=m*y+n

得b*x-m*y=n;

根据exgcd求得一组解x，y

#include "stdio.h"
#include "string.h"
#include "algorithm"
#define m 9973
using namespace std;
int exgcd(int a,int b,int &d,int &x,int &y){
    if(!b) {d=a;x=1;y=0;}
    else  {exgcd(b,a%b,d,y,x);y-=x*(a/b);}
}
int gcd(int a,int b)
{
    if(a%b==0)
           return b;
    return gcd(b,a%b);
}
int main()
{
    int x0,y0,x,y;
    int g,t;
    int n,a,b,m1,b1;
    scanf("%d",&t);
    while(t--){
        scanf("%d%d",&n,&b);
        g=gcd(b,m);
        m1=m/g;
        b1=b/g;
        exgcd(b,m,g,x0,y0);
        x0=x0*(n/g);
        x=x0+m1;
        while(x<0)
             x+=m;
        printf("%d\n",x%m);
    }
    return 0;
}