HDU-1573 X问题(中国剩余定理)

原题链接

X问题

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 4923    Accepted Submission(s): 1641

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
  
  
  
  

   
   

  
  

    

这种情况就采用两两合并的思想，假设要合并如下两个方程

 

      

 

那么得到

 

       

 

在利用扩展欧几里得算法解出的最小正整数解，再带入

 

       

 

得到后合并为一个方程的结果为

 

       

 

这样一直合并下去，最终可以求得同余方程组的解。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#define INF 1e9
using namespace std;
typedef long long ll;

int a[15], b[15];
void gcd(int k1, int k2, int &x, int &y, int &d){
	
	if(k2 == 0){
		x = 1;
		y = 0;
		d = k1;
		return ;
	}
	int x1, y1;
	gcd(k2, k1%k2, x1, y1, d);
	x = y1;
	y = x1 - k1 / k2 * y1;
}
int main(){
	
//	freopen("in.txt", "r", stdin);
	int t;
	
	scanf("%d", &t);
	while(t--){
		int n, m;
		scanf("%d%d", &n, &m);
		for(int i = 0; i < m; i++)
		  scanf("%d", a+i);
		for(int i = 0; i < m; i++)
		  scanf("%d", b+i);
		int k1 = b[0], k2 = a[0];
		int sign = 0;
		for(int i = 1; i < m; i++){
		
			int h = b[i] - k1;
			int x, y, d;
			gcd(k2, a[i], x, y, d);
			if(h % d){
				sign = 1;
				break;
			}
			int v = a[i] / d;
			x *= h / d;
			x = (x % v + v) % v;
			x = k1 + k2 * x;
			k1 = x;
			k2 = k2 * a[i] / d;
			k1 %= k2;
		}
		if(sign){
			puts("0");
			continue;
		}
		int ans = n / k2;
		int h = ans * k2;
		if(h != n && (h+1)%k2 <= k1 && k1 <= n%k2)
		  ans++;
		printf("%d\n", ans);
	}
	return 0;
}