hdu 1576

设A=B*k1=k2*99973+n
则B*k1-k2*9973=n
因为gcd(B,9973)=1;
用扩展欧几里得算法可以求出x1,x2使
x1*b+x2*9973=1;
n*x1*b+n*x2*9973=n;
求(n*x1)%9973;

因为求出的x1可能有负数,(n*k%9973+9973)%9973;






#include<stdio.h>
#include<math.h>
#define inf 9973
int q,x,y;
void gcd(int a,int b)
{
	if(b==0)
	{
		q=a;x=1;y=0;return;
	}
	else 
	{
		gcd(b,a%b);
		int temp=x;
		x=y;
		y=temp-a/b*x;
	}
}
int main()
{
	int i,j,b,n,t,k1;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d%d",&n,&b);
		if(inf<b)
		{
			gcd(b,inf);
			k1=x;
		}
		else
		{
		   gcd(inf,b);
		   k1=y;
		}
		//printf("%d=(%d)*(%d)+(%d)*(%d)\n",q,x+y-k1,inf,k1,b);
		k1=(n*k1)%inf;
		k1=(k1+inf)%inf;//求出的k1有负数,刚开始求绝对值发现不对,跟答案相加正好等于9973
		printf("%d\n",k1);		
	}
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值