hdoj 1576

时间:2024-09-26 14:38:17

  //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;
}