hdu1573 X问题中国剩余定理改进

最新推荐文章于 2023-07-21 20:41:03 发布

popcjz

最新推荐文章于 2023-07-21 20:41:03 发布

阅读量1.5k

点赞数

分类专栏：数学

数学专栏收录该内容

17 篇文章 0 订阅

订阅专栏

X问题

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

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

   310 31 2 30 1 2100 73 4 5 6 7 8 91 2 3 4 5 6 710000 101 2 3 4 5 6 7 8 9 100 1 2 3 4 5 6 7 8 9 
 

Sample Output

103

/*
正整数解=x0+lcm*k（x0最小正整数解，k任意非负整数）
坑点：
  1.中国剩余定理是有0解的
  2.当解为负数解的时候，你要通过加最小公倍数lcm将它调整为最小正整数解。
  3.0不是正整数。
  4.lcm要求，我的最开始写的模板没有lcm的求解。lcm(a,b)=a*b/gcd(a,b);
  5.尽量使用long long 数据类型，n不用。
*/

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

const int maxn=12;  
ll mod,lcm;  
int n,flag;  
ll a[maxn],m[maxn];//m为mod数组，a为余数数组  
  
ll exgcd(ll a,ll b,ll &x,ll &y)  //扩展欧几里得算法 求出a,b的最大公约数，且使x，y满足ax+by=GCD(a,b) 
{  
    if(b==0)
	{  
        x=1,y=0;  
        return a;  
    }  
    ll ans=exgcd(b,a%b,x,y);  
    ll temp=x;  
    x=y;  
    y=temp-a/b*y;  
    return ans;  
} 

ll GCD(ll a,ll b)  //最大公约数 
{  
    if(b==0) 
	return a;  
    return GCD(b,a%b);  
}  

int china2(int num)  //不互质的中国剩余定理  
{
    ll m1=m[0],a1=a[0],m2,a2,k1,k2,x0,gcd,c;  
    lcm=m[0];  
    for(int i=1;i<num;i++)
	{  
        m2=m[i],a2=a[i];  
        c=a2-a1;  
        gcd=exgcd(m1,m2,k1,k2);//解得：n1*k1+n2*k2=gcd(n1,n2)  
        lcm=lcm*m[i]/GCD(lcm,m[i]);  
        if(c%gcd)
		{  
            flag=1;  
            return 0;//无解  
        }  
        x0=c/gcd*k1;//n1*x0+n2*(c/gcd*k2)=c  PS:k1/gcd*c错误！  
        ll t=m2/gcd;  
        x0=(x0%t+t)%t;//求n1*x0+n2*y=c的x0的最小解  
        a1+=m1*x0;  
        m1=m2/gcd*m1;  
    }  
    return a1;  
}  
  
int main()  
{  
    int T;scanf("%d",&T);while(T--)
	{  
        scanf("%lld%d",&mod,&n);  
        for(int i=0;i<n;i++) scanf("%lld",&m[i]);  
        for(int i=0;i<n;i++) scanf("%lld",&a[i]);  
        flag=0;  
        ll ans=china2(n);  
        
        int sum=0;
        if(flag==1||ans>mod) 
		printf("0\n");//中国剩余定理解出来的值可能为0可能为负值可能为正，所以要加flag判断无解  
        else
		{  
            while(ans<=0) 
			ans+=lcm;//题目要求是正整数，所以要通过最小公倍数将负数不断变为正整数  
            while(ans<=mod)//当解小于范围mod，则不断往上加 
			{//计算正解的个数  
                sum++;  
                ans+=lcm;  
            }  
            printf("%d\n",sum);  
        }  
    }
    return 0;
}