HDOJ---1573 X问题[中国剩余定理]

最新推荐文章于 2019-02-10 15:37:11 发布

weixin_30343157

最新推荐文章于 2019-02-10 15:37:11 发布

阅读量67

点赞数

原文链接：http://www.cnblogs.com/XBWer/archive/2012/08/10/2632122.html

版权

X问题

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

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 
   

Author

lwg

Source

HDU 2007-1 Programming Contest

Recommend

linle

转自：http://hi.baidu.com/nicker2010/item/da4754d5f5097c8c6cce3fed

解题思路：
1.因为(a1,a2,a3,a4,….,ak)不一定互质，所以不能够直接用中国剩余定理。
2.x=r1+a1*k1,x=r2+a2*k2,所以有r1+a1*k1=r2+a2*k2，化简后得到 a1*k1=(r2-r1) mod(a2);
用扩展欧几里得可以得到最小的k1,所以x=r1+a1*k1+a1*a2/gcd(a1,a2)，
就这样一直替换最后剩余一个同余方程。r1就是最后的解。

对于x=a1 mod b1,x= a2 mod b2，设x=a1+m*b1
所以b1*m=a2-a1 mod b2，利用欧几里德扩展定理求出最小的非负m，
那么x=a1+m*b1就已知,且x最小，如果无解，整个同余式组无解
同时，x+k*b1是所有满足x=a1 mod b1的解，而x+k'*b2又是所有满足x=a2 mod b2的解
那么，将x+k*b1与x+k'*b2合并，得到的式子就是x+k*lcm(b1,b2)
于是，上面两个式子可以用x'=x mod lcm(b1,b2)来替代
最后，就只剩下一个式子了，求得的最小的x就是答案
对于一次同余式ax=b mod n，设d=gcd(a,n)，则同余式有解的充要条件为d|b
假设d=a*x'+n*y',则x0=b/d*x'一定为方程组的一个解,且共有d个解，
（证明可以参考算法导论）最小正整数解为(x0%n/d+n/d)%(n/d)

code:

 1 #include <iostream>
 2 #include <vector>
 3 using namespace std;
 4 long long exgcd(long long a,long long b,long long &x,long long &y)
 5 {
 6     if(a==0)
 7     {x=0;y=1;return b;}
 8     long long g = exgcd(b%a,a,x,y);
 9     long long tem = y;
10     y=x;
11     x=tem-(b/a)*y;
12     return g;
13 }
14 long long CRT(const vector<long long>& m,const vector<long long>& b,long long& lcm)
15 {
16     bool flag = false;
17     long long x, y, i,d,result,a1,m1,a2,m2,Size=m.size();
18     m1 = m[0]; a1 = b[0];
19     for(i=1;i<Size;++i)
20     {
21         m2 = m[i]; a2 = b[i];
22         d = exgcd( m1, m2, x, y );
23         if((a2-a1)%d != 0) flag = true;
24         result = (x * ((a2-a1) / d ) % m2 + m2 ) % m2;
25         a1 = a1 + m1 * result; //对于求多个方程
26         m1 = (m1 * m2) / d;    //lcm(m1,m2)最小公倍数
27         a1 = (a1 % m1 + m1) % m1;
28     }
29     lcm = m1;
30     if (flag) return -1;
31     else return a1;
32 }
33 int main()
34 {
35     long long T,N,M,i,num,ans,lcm,result;
36     vector<long long> miVec,biVev;
37     cin>>T;
38     while(T--)
39     {
40         cin>>N>>M;
41         miVec.clear();biVev.clear();
42         for(i=0;i<M;++i)
43         {cin>>num;miVec.push_back(num);}
44         for(i=0;i<M;++i)
45         {cin>>num;biVev.push_back(num);}
46         ans = CRT(miVec,biVev,lcm);
47         if(ans<0){cout<<0<<endl;continue;}
48         if(N<ans) result = 0;
49         else result=(N-ans)/lcm + 1;
50         if(ans==0)--result;
51         cout<<result<<endl;
52     }
53     return 0;
54 }