poj2891（中国剩余定理模版：：不要求模数互质）

2016.8.23完全按照思路来写的中国剩余定理

#include<cstdio>  
#include<algorithm>  
#include<cstring>  
#define ll long long   
using namespace std;  
ll k;  
void exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
	if (b==0){d=a;x=1;y=0;return ;}
	exgcd(b,a%b,d,y,x);y-=x*(a/b);
}
int main()  
{  
    while (scanf("%lld",&k)!=EOF)  
    {  
        ll a1,a2,r1,r2;  
        bool o=true;  
       for (int i=1;i<=k;i++)
       {
       		if (i==1)
       		{
       			scanf("%lld%lld",&a1,&r1);
       			continue;
			}
			scanf("%lld%lld",&a2,&r2);
			
			ll x,y,d;
			exgcd(a1,a2,d,x,y);
			if ((r2-r1)%d) o=false;
			x=(r2-r1)/d*x;
			
			ll k=a2/d;
			x=((x%k)+k)%k;
			
			ll X=x*a1+r1;//这是每次求解后的X，就是那个解
			a1=a1/d*a2;//模数求出来
			r1=X%a1;//知道解，知道模数，直接mod就出来了
	   }
        if (o) printf("%lld\n",r1);  
        else printf("-1\n");  
    }  
    return 0;  
}  

写法一

#include<cstdio>
#include<algorithm>
#include<cstring>
#define ll long long 
using namespace std;
ll k;
void exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
	if (b==0){d=a;x=1;y=0;return ;}
	exgcd(b,a%b,d,y,x);y-=x*(a/b);
}
int main()
{
	while (scanf("%lld",&k)!=EOF)
	{
		ll a1,a2,r1,r2;
		bool o=true;
		for (int i=1;i<=k;i++)
		{
			if (i==1)
			{
				scanf("%lld%lld",&a1,&r1);
				continue;
			}
			scanf("%lld%lld",&a2,&r2);
			ll x,y,d;
			exgcd(a1,a2,d,x,y);//
			if ((r2-r1)%d) o=false;//如果不符合要求，标记
			x=x*((r2-r1)/d);
			ll m=a2/d;
			x=(x%m+m)%m;//1
			r1=r1+a1*x;//更新余数
			a1=a1/d*a2;//更新模数
			r1=(r1%a1+a1)%a1;//2，两次取模保证是最小正整数解
		}
		if (o) printf("%lld\n",r1);
		else printf("-1\n");
	}
	return 0;
}

写法二

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define ll long long
using namespace std;
void exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
	if (b==0){x=1;y=0;d=a;return;}
	exgcd(b,a%b,d,y,x);y-=x*(a/b);
}
int main()
{
	int k;
	while (scanf("%d",&k)!=EOF)
	{
		ll a1,r1,a2,r2;
		scanf("%lld%lld",&a1,&r1);
		bool o=true;
		for (int i=1;i<k;i++)
		{
			scanf("%lld%lld",&a2,&r2);
			ll x,y,d;
			exgcd(a1,a2,d,x,y);
			if ((r2-r1)%d) o=false;
			x*=(r2-r1)/d;
			ll h=a2/d;
			x=(x%h+h)%h;
			r1=x*a1+r1;
			a1=a1/d*a2;
			r1=(r1%a1+a1)%a1;
		}
		if (o) printf("%lld\n",r1);else printf("-1\n");
	}
	return 0;
}