ZOJ 3609 Modular Inverse

最新推荐文章于 2019-04-03 19:06:54 发布

LsFlyt

最新推荐文章于 2019-04-03 19:06:54 发布

阅读量454

点赞数

分类专栏：寒假练习题数学文章标签：扩展欧几里得 zoj

本文链接：https://blog.csdn.net/LsFlyt/article/details/50670129

版权

寒假练习题同时被 2 个专栏收录

30 篇文章 0 订阅

订阅专栏

数学

10 篇文章 0 订阅

订阅专栏

Modular Inverse

Time Limit:2000MS MemoryLimit:65536KB 64bitIO Format:%lld& %llu

Description

The modularmodular multiplicative inverse of aninteger a modulo m isan integer x suchthat a^-1≡x (mod m). Thisisequivalent toax≡1 (mod m).

Input

There aremultiple test cases. The first line of input isan integer T ≈2000 indicating the number of test cases.

Each test casecontains two integers 0 < a ≤1000 and 0 < m ≤1000.

Output

For each testcase, output the smallest positive x.If such x doesn'texist, output "Not Exist".

Sample Input

3 11

4 12

5 13

Sample Output

Not Exist

References

http://en.wikipedia.org/wiki/Modular

题目数据不大，可以直接暴力过

训练时候用的扩展欧几里得，在这里顺便回顾一下：

对于求最大公约数有gcd(a,b)=gcd(b,a%b),gcd(a,0)=a;

而扩展欧几里得可以求形如ax+by=gcd(a,b)的方程的解，所以对于ax+by=c;如果c!=gcd(a,b)则无解。

设

ax₁+by₁=gcd(a,b)

则

bx₂+(a%b)y₂=gcd(b,a%b)=ax₁+by₁

即

bx₂+(a-a/b)y₂=ax₁+by₁

ax₁+by₁=ay₂+b(x₂-a/b*y₂)

所以x₁₌y₂ y₁=x₂-a/b*y₂

而b=0时，x=1,y=0

代码实现如下：

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

对于本题求ax+bm=1

无解情况：m!=1 && gcd(a,m)!=1或者m==1 && a==1

代码如下：

#include<cstdio>

using namespace std;

	int T;
	int a,m;
	int x,y;
	
int gcd(int a,int b)
{
	if (!b)
	{
		return a;
	}
	return gcd(b,a%b);
}	

void extend_gcd(int a,int b)
{
    if(b==0)
	{  
      x=1;  
      y=0;  
    }  
    else
	{  
      extend_gcd(b,a%b);  
      int temp=x;  
      x=y;  
      y=temp-a/b*y;  
    }  
}  
int main()
{
	scanf("%d",&T);
	while (T--)
	{
		scanf("%d%d",&a,&m);
		if (m==1)
		{
			if (!(a%m))
			{
				printf("1\n");
			}
			else
			{
				if (!(m%a))
				{
					printf("%d\n",m/a);
				}
				else
				{
					printf("Not Exist\n");
				}
			}
		}
		else
		{
			if (gcd(a,m)==1)
			{
				extend_gcd(a,m);
				while (x<0)
				{
					x+=m;
				}
				printf("%d\n",x);
			}
			else
			{
				printf("Not Exist\n");
			}
		}
		
	}
	return 0;
}