HDOJ 1576 A/B(扩展欧几里得算法或者数学转化解法)



A/B

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 3436    Accepted Submission(s): 2608

Problem Description

要求(A/B)%9973，但由于A很大，我们只给出n(n=A%9973)(我们给定的A必能被B整除，且gcd(B,9973) = 1)。

 

Input

数据的第一行是一个T，表示有T组数据。
每组数据有两个数n(0 <= n < 9973)和B(1 <= B <= 10^9)。

 

Output

对应每组数据输出(A/B)%9973。

 

Sample Input

  
  

   
   2
1000 53
87 123456789
  
  

 

Sample Output

  
  

   
   7922
6060
  
  

 

第一种，扩展欧几里得算法:

运用扩展欧几里德算法能解出gcd(a,b)=a*x1+b*y1的x1和y1的值。

根据题目我们知道：

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的值。

利用扩展欧几里德算法可求出gcd(B,9973)=Bx1+9973y1=1的x1。题中说gcd(B,9973)=1；所以等式两边同乘以n，得B(n*x1)-9973(-n*y1)=n。可知n*x1就是B*x-9973*y=n的解了！！！即x=n*x1。

在扩展欧几里得算法中得到的x可能为负值，所以还需要x=(x%9973+9973)%9973。

代码如下：

<span style="font-size:18px;">#include<cstdio>
#include<cstring>

void extgcd(int a,int b,int &x,int &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return ;
	}
	else
	{
		extgcd(b,a%b,y,x);
		y-=(a/b)*x;
	} 
}

int main()
{
	int t,n,b,x,y;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d%d",&n,&b);
		extgcd(b,9973,x,y);
		x*=n;
		x=(x%9973+9973)%9973;
		printf("%d\n",x);
	}
	return 0;
} </span>

第二种，数学思维解法：

思路：设X=（A/B）%9973
          因为n=A%9973，所以A=k*9973+n (k为一常数)
          又因为X=（A/B）%9973，所以A/B=d*9973+X （d为一常数）       两边同乘以B，得：A=B*d*9973+B*X

即              B*d*9973+B*X=k*9973+n
移项得：   B*d*9973+B*X-n=k*9973
所以题意转为  只要满足（B*X-n）%9973=0，X即为要求的结果。n的值知道，B的值知道，又因为x的取值范围是0到9972，因此枚举x的值即可，满足条件的就是答案。

代码如下：

<span style="font-size:18px;">#include<cstdio>
int main()
{
	int n,a,i,t;
	long long b;
	scanf("%d",&t);
	while(t--)
	{
		scanf("%d%lld",&n,&b);
		for(i=0;i<9973;++i)
		{
			if((b*i-n)%9973==0)
				break;
		}
		printf("%d\n",i);
	}
	return 0;
}</span>