hdu 1573 X问题 两两可能不互质的中国剩余定理

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

 

Author

lwg

思路：不互质的中国剩余定理（和题目中的a，b数组反了，不要在乎这些细节）

　　　　x≡ a1（mod b1）

　　　　x≡ a2（mod b2)

　　　　a1+b1*m1=a2+b2*m2扩展欧几里德搞下

代码：

#include<iostream>
#include<cstdio>
#include<cmath>
#include<string>
#include<queue>
#include<algorithm>
#include<stack>
#include<cstring>
#include<vector>
#include<list>
#include<set>
#include<map>
using namespace std;
#define ll long long
#define mod 1000000007
#define inf 999999999
//#pragma comment(linker, "/STACK:102400000,102400000")
int scan()
{
    int res =  , ch ;
    while( !( ( ch = getchar() ) >= '' && ch <= '' ) )
    {
        if( ch == EOF ) return  <<  ;
    }
    res = ch - '' ;
    while( ( ch = getchar() ) >= '' && ch <= '' )
        res = res *  + ( ch - '' ) ;
    return res ;
}
int a[];
int b[];
int gcd(int x,int y)
{
    if(x%y==)
    return y;
    else
    return gcd(y,x%y);
}
void exgcd(int a, int b, int &x, int &y)
{
    if(b == )
    {
        x = ;
        y = ;
        return;
    }
    exgcd(b, a % b, x, y);
    int tmp = x;
    x = y;
    y = tmp - (a / b) * y;
}
int main()
{
    int x,y,z,i,t;
    scanf("%d",&x);
    while(x--)
    {
        scanf("%d%d",&y,&z);
        for(i=;i<z;i++)
        scanf("%d",&b[i]);
        for(i=;i<z;i++)
        scanf("%d",&a[i]);
        int a1=a[],b1=b[];
        int jie=;
        for(i=;i<z;i++)
        {
            int a2=a[i],b2=b[i];
            int xx,yy;
            int gys=gcd(b1,b2);
            if((a2-a1)%gys)
            {
                jie=;
                break;
            }
            exgcd(b1,b2,xx,yy);
            xx=(xx*(a2-a1))/gys;
            int gbs=b1*b2/gys;
            a1=(((xx*b1+a1)%gbs)+gbs)%gbs;
            b1=gbs;
        }
        if(!jie||y<a1)
        printf("0\n");
        else
        printf("%d\n",(y-a1)/b1+-((a1==)?:));
    }
    return ;
}