扩展欧几里得:
已知:设:K = (A/B)%9973 , 则A/B = k + 9973x (x未知), 因此A == kB + 9973xB
A%9973==n,所以 KB%9973==n;
故(K/n)B +(-y/n)*9973 = gcd(B,9973) ;
转化为 ax+by=gcd(a,b);
所以直接用扩展欧几里得求出x==K/n;
再让x*n即为所求;
#include<stdio.h> #include<iostream> #include<algorithm> #include<math.h> #include<string.h> using namespace std; int exgcd(int a,int b,int &x,int &y) { if(b==0) { x=1; y=0; return a; } int r=exgcd(b,a%b,x,y); int t=x; x=y; y=t-a/b*y; return r; } int main() { int n,b,ans,x,y; int T; scanf("%d",&T); while(T--) { scanf("%d%d",&n,&b); exgcd(b,9973,x,y); x=(x+9973)%9973; printf("%d\n",(x*n)%9973); } }
