Strange Way to Express Integers——解题报告

一、题目链接：

http://poj.org/problem?id=2891

二、题目大意：

给定$2n$个正整数$a_1，a_2，\dots \dots ，a_n$和$m_1，m_2，...，m_n$，求一个最小的正整数$x$，满足任意$i\in [1,n]$，都存在$x\equiv a_i(mod{m}_i)$，若无解则输出$-1$。

三、题目分析：

1.因为此时没有保证$m_i$相互之间是互质关系，所以没有办法使用中国剩余定理，只能考虑其他的方法。
2.我们假设此时已经求出了前$k-1$个方程的解$x$。设$M=\begin{array}{c}{\prod }_{i=1}^{k-1}{x}_i\end{array}$，而此时显然对于前面的所有$k-1$个方程来说，$x+t\times M$都是方程的解，因为$t\times M$对任意$m_i(i\le k-1)$取余都为$0$。
3.而对于第k个方程来说，要是存在x+t*M是方程的解，那么此解就是前k个方程的通解。所以我们的操作如下：

• 因为此时要求$x+t\times M\equiv a_k(mod{m}_k)$，而方程等价于$t\times M\equiv a_k-x(mod{m}_k)$。
• 因为$a_k-x$已知，所以转化为了一个全新的线性同余方程，此时使用扩展欧几里得算法求出其中一个特解，若无解则直接输出$-1$。
• 而求出的特解即$t$，可以得到$x+t\times M\equiv a_k(mod{m}_k)$成立，同时$x^{\prime }=x+t\times M$同时也是前面$k$个方程的通解。
• 重复上述过程直到$k=n$运行结束。

4.由于$M$很容易超$longlong$，所以我们一定要注意取余的问题，详细的细节可以看代码。

四、正解代码：

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <cstdlib>

using namespace std;
typedef long long ll;
ll exgcd(ll a,ll b,ll &x,ll &y)
{
	if(!b)
	{
		x=1;
		y=0;
		return a;
	}
	ll d=exgcd(b,a%b,x,y);
	ll z=x;
	x=y;
	y=z-y*(a/b);
	return d;
}
ll n;
int main()
{
	while(~scanf("%lld",&n))
	{
		ll a,m;
		scanf("%lld%lld",&m,&a);
		ll x=a,M=m;
		bool flag=false;
		for(ll i=2;i<=n;i++)
		{
			ll t1,ty,New;
			scanf("%lld%lld",&m,&a);
			New=((a-x)%m+m)%m;
			ll d=exgcd(M,m,t1,ty);
			if(New%d==0)
				t1=t1*(New/d)%m;
			else
				flag=true;
			x+=t1*M;
			M=M/d*m;
			x=(x%M+M)%M;
		}
		if(!flag)
			printf("%lld\n",x);
		else
			printf("-1\n");
	}
	
	return 0;
}