hdu 1576 求解模线性方程

 

题意：
已知  gcd(B, mod)=1, A%B==0,A/B *B=n(mod mod);
给出B,n;求 A/B= ans (mod mod)中的ans

题解：
求解模线性方程
先求A/B，令A/B=x
则有B*x=n(mod mod),可想到扩展欧几里德
求出x,因为这里的x是B*x=1(mod 9973)的解，x之后还要乘n
代码如下：

#include<iostream>
#include<cstdio>

using namespace std;

const int mod=9973;

__int64 a,b,c,d,x,y;
__int64 exgcd(__int64 a,__int64 b)
{
	if(b==0)
	{
		x=1;y=0;
		return a;
	}
	d=exgcd(b,a%b);
	__int64 xx=y,yy=x-a/b*y;
	x=xx;y=yy;
	return d;
}

int main()
{
	int cas;
	scanf("%d",&cas);
	while(cas--)
	{
		__int64 B,n;
		scanf("%I64d%I64d",&n,&B);
		d=exgcd(B,mod);
		x*=n;
		printf("%d\n",(x%mod+mod)%mod);
	}
	return 0;
}