解同余方程组

已知: *
N=a1x+b1;
N=a2y+b2;
。。。。。
N=anz+bn;
求N；

解：
a1x+b1=a2y+b2; ==> a1x - a2y=b2-b1; (因为求x 所以a2前正负不影响结果）；
令c=b2-b1;
通过扩展欧几里得算法对此式求解 求出x, 并得到gcd(a1,a2)=d;
如果c%d!=0,则此方程组无解；
如果有解则：
带入到 N=a1x+b1 中 ；
得到一个N的特解 ： N0=a1x+b1;
因为要同时满足式子： N=a1x+b1 N=a2y+b2;
所以 得到N的暂时的通解： N = N0 + (a1a2)/d * k; (k为正整数）；
此通解变形得到和已知相同形式的式子： N=( (a1a2)/d ) k + N0; 令a1= (a1a2)/d , b1=N0;
此时用此式和已知中下一个式子 联立 即可 求出 k 以及gcd,如此循环即可求解；

//例题   HDU 1573   最后求在已知范围内 此方程组有几个解 
#include<stdio.h>
#include<string.h>
#include<iostream>
using namespace std;


void look(int a,int b,int &d,int &x,int &y)
{
	if(!b)
	{
		x=1;y=0;d=a;
	}
	else
	{
		look(b,a%b,d,y,x);
		y-=(a/b)*x;
	}
}
int main()
{
	int n,m,a1,b1,a2,b2,flag=0,a[11],b[11],T,k,d,x,y,c;;
	scanf("%d",&T);
	while(T--)
	{
		cin>>n>>m;
	
		for(int i=0;i<m;i++)
			scanf("%d",&a[i]);
		for(int i=0;i<m;i++)
			scanf("%d",&b[i]);
		
		flag=0;
		
		a1=a[0];
		b1=b[0];
		
		for(int i=1;i<m;i++)
		{
			a2=a[i];
			b2=b[i];
			
			look(a1,a2,d,x,y);
			c=b2-b1;
			if(c%d)
			{
				flag=1;
				break;
			}
			int t=a2/d; 
			x=(x*(c/d)%t+t)%t;
			b1=a1*x+b1;
			a1=a1*a2/d;
		}
		if(flag||n<b1)
			printf("0\n");
		else
		{
			int ans=(n-b1)/a1+1;//+1是因为本身(b1)也为一个解 
			if(b1==0)
				ans--;//如果本身为0不满足题意正整数解  需要减掉 
			printf("%d\n",ans);
		}	
	}
	return 0;
}