解决该题的关键是:
1、了解扩展欧几里德算法,可以运用其解出gcd(a,b)=ax1+by1中的x1、y1的值
2、由题可得以下内容:
n=A%9973,则n=A-A/9973*9973。又A/B=x,则A=Bx。所以Bx-A/9973*9973=n。即Bx-9973y=n。
到这里我们可以发现:只要求出x的值,即可算出x%9973,也就是(A/B)%9973了。顺利解决了!
3、题目关键转到如何求出x了。题目的输入是n和B,利用扩展欧几里德算法可求出gcd(B,9973)=Bx1+9973y1=1的x1。
等式两边同乘以n,得B(nx1)-9973(-ny1)=n。可知nx1就是Bx-9973y=n的解了!!!即x=nx1。
4、对于第三部得到的x可能是负数,由题这显然是不正确的。
可以做这样的转化:(x%9973+9973)%9973
(最后一点也不太懂,不懂转化后为啥任然正确!期待大神赐教)
AC代码:
#include<iostream>
using namespace std;
const int MOD=9973;
void extendGcd(int a,int b,int &x,int &y) //扩展gcd,可以求出gcd(a,b)以及ax+by=gcd(a,b)中x,y的值
{
if(b==0)
{
x=1;
y=0;
return;
}
else
{
extendGcd(b,a%b,x,y);
int tmp=x;
x=y;
y=tmp-a/b*y;
}
}
int main()
{
int t,n,b,x,y,tmp;
cin>>t;
while(t--)
{
cin>>n>>b;
extendGcd(b,MOD,x,y); //解出bx+9973y=1的x
x*=n; //此时的x为bx1-9973y1=n中的解x1了
tmp=(x%MOD+MOD)%MOD; //防止x为负,有题意x必为正数
cout<<tmp<<endl;
}
return 0;
}