POJ——1061（数论之线性同余方程求解）

题目地址：http://poj.org/problem?id=1061

参考地址：http://www.cnblogs.com/comeon4mydream/archive/2011/07/18/2109060.html（写的很详细，还有证明。）

哈工程地址：http://acm.hrbeu.edu.cn/index.php?act=status&cid=53（哈工程的OJ，就是那本基础题解）。

//不过这个借鉴的程序在哈工程的OJ上过不了，下面的我的程序过了。可能Poj上的数据弱吧。

#include <iostream>    //不过这个借鉴的程序在哈工程的OJ上过不了，下面的我的程序过了。
using namespace std;
long long extgcd(long long a, long long b, long long &x, long long &y)
{
    long long d, t;
    if (b == 0) { x = 1; y = 0; return a; }
    d = extgcd(b, a % b, x, y);
    t = x - a / b * y; x = y; y = t;
    return d;
}
int main()
{
    long long x, y, m, n, L, X, Y, d, r;
    while (cin >> x >> y >> m >> n >> L)
    {
        d = extgcd(n - m, L, X, Y); r = L / d;
        if ((x - y) % d) cout << "Impossible" << endl;
        else cout << ((x - y) / d * X % r + r) % r << endl;
    }
    return 0;
}

自己修改后的自己的代码：

#include <iostream>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <cstdio>
using namespace std;
typedef long long ll;
#define INF 0xfffffff
#define MAX(a,b) a>b?a:b
#define MIN(a,b) a>b?b:a
#define N 101
void extend_gcd(ll a,ll b,ll &d,ll &x,ll &y)
{
	if(!b) d=a,x=1,y=0;
	else{
		extend_gcd(b,a%b,d,y,x);
		y-=x*(a/b);
	}
}
int main()
{
	int i,j,k,t;
	ll x,y,m,n,l;
	while(cin>>x>>y>>m>>n>>l)
	 {
	    ll a=m-n;
	    ll c;
	 	if(a<0)   //保证a，b同时为正。 
		 {
		  a=-a;
	          c=(x-y+l)%l;
		 }
		 else 
		 {
 		 c=(y-x+l)%l;
 		 }
	 	ll b=l;	 
	 	ll xx,yy,dd;
	 	extend_gcd(a,b,dd,xx,yy);
	 	if(c%dd) cout<<"Impossible"<<endl;
	 	else 
	 	{
	 		 xx=xx*(c/dd);
	 	         b/=dd;
			 ll ans=(xx%b+b)%b;	
			 printf("%ld\n",ans);
	 	}
 	 }	
	return 0;
}