HDU 1573 多元同余线性方程

X问题

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)

Problem Description

求在小于等于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中的元素。

 

Output

对应每一组输入，在独立一行中输出一个正整数，表示满足条件的X的个数。

 

Sample Input

  
  

   
   3
10 3
1 2 3
0 1 2
100 7
3 4 5 6 7 8 9
1 2 3 4 5 6 7
10000 10
1 2 3 4 5 6 7 8 9 10
0 1 2 3 4 5 6 7 8 9
  
  

 

Sample Output

  
  

   
   1
0
3
  
  
  
  


  
  
  
  


  
  
  
  

   
   题解：解线性同余方程，然后找1-n之间有多少个数符合min(ans)+lcm(a1,a2...an).
  
  
  
  


  
  
  
  

   
   
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
typedef long long ll;
ll a[15],b[15];
void exgcd(ll a,ll b,ll &d,ll &x,ll &y){
	if(b==0){
		x=1;
		y=0;
		d=a;
		return ;
	}
	else{
		exgcd(b,a%b,d,x,y);
		ll temp=x;
		x=y;
		y=temp-(a/b)*y;
	}
}
int main(){
	int t;
	scanf("%d",&t);
	while(t--){
		ll n,m,i,j;
		ll lcm=1;
		scanf("%lld%lld",&n,&m);
		for(i=1;i<=m;i++){
			scanf("%lld",&a[i]);
			lcm=lcm*a[i]/__gcd(lcm,a[i]);
		}
		ll ifhave=1;
		for(i=1;i<=m;i++)scanf("%lld",&b[i]);
		ll d,x,y;
		for(i=2;i<=m;i++){
			ll al=a[1],bl=a[i],c=b[i]-b[1];
			exgcd(al,bl,d,x,y);
			if(c%d){
				ifhave=0;
				break;
			}
			ll t=bl/d;
			x=(x*(c/d)%t+t)%t;
			b[1]=a[1]*x+b[1];
			a[1]=a[1]*(a[i]/d);
		}
		if(!ifhave){
			printf("0\n");
		}
		else{
			ll ans=n-b[1];
			ans=ans/lcm+1;
			if(!b[1])ans--;
			if(b[1]>n)ans=0;//这里注意 如果最小的都不在n之内 ans就是0
			printf("%lld\n",ans);
		}
	}
	return 0;
}