hdu1573（中国剩余定理，在区间内解的数量）

求在小于等于N的正整数中有多少个X满足：X mod a[0] = b[0], X mod a[1] = b[1], X mod a[2] = b[2], …, X mod a[i] = b[i], … (0 < a[i] <= 10)。

 

Input

输入数据的第一行为一个正整数T，表示有T组测试数据。每组测试数据的第一行为两个正整数N，M (0 < N <= 1000,000,000 , 0 < M <= 10)，

表示X小于等于N，数组a和b中各有M个元素。接下来两行，每行各有M个正整数，分别为a和b中的元素。

因为，解是以lcm（m1，m2。。。mn） 作为循环的

所以 用n/lcm，再看一看余数是否大于等于a1（就是最小正整数答案），如果这样，ans++

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#define ll long long 
using namespace std;

ll n,m,b[10009],a[10009];
ll gcd(ll a,ll b)
{
	return b? gcd(b,a%b):a;
}
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()
{
	int T;
	scanf("%d",&T);
	for (int kk=1;kk<=T;kk++)
	{
		scanf("%lld%lld",&n,&m);
		bool o=true;
		ll m1,m2,a1,a2;
		for (int i=1;i<=m;i++) scanf("%lld",&b[i]);
		for (int i=1;i<=m;i++) scanf("%lld",&a[i]);
		a1=a[1];m1=b[1];
		for (int i=2;i<=m;i++)
		{
			ll d,x,y;
			m2=b[i];a2=a[i];//还是按这样写中国剩余定理吧！
			exgcd(m1,m2,d,x,y);
			if ((a2-a1)%d) o=false;
			x*=(a2-a1)/d;
			ll k=m2/d;
			x=(x%k+k)%k;
			a1=m1*x+a1;
			m1=m1/d*m2;
			a1=(a1%m1+m1)%m1;
		}
		if (a1==0)
		{
			a1=1;
			for (int i=1;i<=m;i++) a1=a1/gcd(a1,b[i])*b[i];
		}//特殊情况，如果余数都是0
		ll ans=0;
		ans=n/m1;
		n=n%m1;
		if(n>=a1) ans++;//计算答案
		if (o) printf("%lld\n",ans);else printf("0\n");
	}
	return 0;
}