hdu 1573(中国剩余定理)

Posted AC菜鸟机

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了hdu 1573(中国剩余定理)相关的知识,希望对你有一定的参考价值。

X问题

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


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.没有整数解输出0
2.求出满足条件的最小非负整数,设该数为X
X+k*mod<=N
k<=(n-x)/mod 当X为0的时候,答案为k,如果不为0,那么将自己也算进去,答案是k+1
#include <stdio.h>
#include <string.h>
#include <algorithm>
#include <iostream>
#include <stdlib.h>
#include <math.h>
using namespace std;
typedef long long LL;
const int N = 10;
LL extend_gcd(LL a,LL b,LL &x,LL &y)
{
    if(b==0)
    {
        x=1,y=0;
        return a;
    }
    else
    {
        LL x1,y1;
        LL d = extend_gcd(b,a%b,x1,y1);
        x = y1;
        y = x1-a/b*y1;
        return d;
    }
}
LL m[N],a[N];///模数为m,余数为a, X % m = a
bool solve(LL &m0,LL &a0,LL m,LL a)
{
    long long y,x;
    LL g = extend_gcd(m0,m,x,y);
    LL t = a-a0>0?a-a0:a0-a;
    if( t%g )return false;
    x *= (a - a0)/g;
    x %= m/g;
    a0 = (x*m0 + a0);
    m0 *= m/g;
    a0 %= m0;
    if( a0 < 0 )a0 += m0;
    return true;
}
/**
* 无解返回false,有解返回true;
* 解的形式最后为 a0 + m0 * t (0<=a0<m0)
*/
bool MLES(LL &m0 ,LL &a0,LL n)///解为 X = a0 + m0 * k
{
    bool flag = true;
    m0 = 1;
    a0 = 0;
    for(int i = 0; i < n; i++)
        if( !solve(m0,a0,m[i],a[i]) )
        {
            flag = false;
            break;
        }
    return flag;
}
int main()
{
    int n;
    LL N;
    int tcase;
    scanf("%d",&tcase);
    while(tcase--)
    {
        scanf("%lld%d",&N,&n);
        for(int i=0; i<n; i++)
        {
            scanf("%lld",&m[i]);
        }
        for(int i=0; i<n; i++)
        {
            scanf("%lld",&a[i]);
        }
        LL m0,a0;
        bool flag = MLES(m0,a0,n);
        LL x = (a0%m0+m0)%m0;
        if(!flag) printf("0\n");
        else
        {
            if(x>N) printf("0\n");
            else{
                LL ans = (N-x)/m0;
                if(x==0) printf("%lld\n",ans);
                else printf("%lld\n",ans+1);
            }
        }
    }
    return 0;
}

 

以上是关于hdu 1573(中国剩余定理)的主要内容,如果未能解决你的问题,请参考以下文章

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

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

hdu1573-X问题-(扩展欧几里得定理+中国剩余定理)

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

hdu1573 X问题中国剩余定理

HDU1573X问题