51Nod-中国剩余定理（加优化版）

最新推荐文章于 2021-09-25 17:00:54 发布

信仰..

最新推荐文章于 2021-09-25 17:00:54 发布

阅读量568

点赞数

分类专栏：数论文章标签：数论

本文链接：https://blog.csdn.net/haut_ykc/article/details/52601538

版权

数论专栏收录该内容

131 篇文章 0 订阅

订阅专栏

1079 中国剩余定理

基准时间限制：1 秒空间限制：131072 KB 分值: 0 难度：基础题

关注

一个正整数K，给出K Mod 一些质数的结果，求符合条件的最小的K。例如，K % 2 = 1, K % 3 = 2, K % 5 = 3。符合条件的最小的K = 23。

Input

第1行：1个数N表示后面输入的质数及模的数量。（2 <= N <= 10)
第2 - N + 1行，每行2个数P和M，中间用空格分隔，P是质数，M是K % P的结果。（2 <= P <= 100, 0 <= K < P)

Output

输出符合条件的最小的K。数据中所有K均小于10^9。

Input示例

Output示例

李陶冶 (题目提供者)

先规定一个mx[i] = (a[1]*a[2]*…*a[end])/a[i];

和一个my[i] = 使得(mx[i]*my[i])%a[i]==1（my[i]尽可能的小，但它可能为负数）

则最中答案n（最小可能）= ((mx[1]*my[1]*b[1]+mx[2]*my[2]*b[2]+…+mx[end]*my[end]*b[end])%(a[1]*a[2]*…

*a[end])+(a[1]*a[2]*…*a[end]))%(a[1]*a[2]*…*a[end]);

至今不知道如何证明。233333

代码如下：

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<limits.h>
#include<algorithm>
#include<queue>
#include<stack>
#include<vector>
#include<math.h>
#include<map>
using namespace std;
#define maxn 105
long long n,m,sum=1;
long long a[maxn],b[maxn];
long long x,y;
void EX_gcd(long long t1,long long t2)//求逆元 
{
	long long temp;
	if(t2==0)
	{
		x=1;
		y=0;
		return ;
	}
	EX_gcd(t2,t1%t2);
	temp=x;
	x=y;
	y=temp-t1/t2*y;
}
long long CRT(long long t1,long long t2)
{
	long long i;
	long long ans=0,d,mi;
	for(i=1;i<=t2;i++)
	{
		mi=t1/a[i];
		EX_gcd(mi,a[i]);
		ans=(ans+x*mi*b[i])%t1;
	}
	return (ans+t1)%t1;
}																					
int  main()
{
	int i,j;
	scanf("%lld",&n);
	for(i=1;i<=n;i++)
	{
		scanf("%lld%lld",&a[i],&b[i]);
		sum*=a[i];
	}
	long long ans=CRT(sum,n);
	printf("%lld\n",ans);
}

优化版（防爆longlong（别人的））

#include<stdio.h>  
#define LL long long  
LL x, y, a[10005], p[10005];  
LL Powto(LL a, LL b, LL mod)            /*乘法优化，防爆long long*/  
{  
    LL aim;  
    aim = 0;  
    while(b>0)               /*原理：将a和b拆成a*(b1+b2+…+bn)，所以b一定要是正数*/  
    {  
        if(b&1)   
            aim = (aim+a)%mod;  
        a = (a+a)%mod;  
        b >>= 1;  
    }  
    return aim;  
}  
LL Ex_Gcd(LL a, LL b)  
{  
    LL d, temp;  
    if(b==0)  
    {  
        x = 1, y = 0;  
        return a;  
    }  
    d = Ex_Gcd(b, a%b);  
    temp = x;  
    x = y;  
    y = temp-a/b*y;  
    return d;  
}  
LL CRT(LL sum, LL k)  
{  
    int i;  
    LL ans, d, mi;  
    ans = 0;  
    for(i=1;i<=k;i++)  
    {  
        mi = sum/a[i];  
        d = Ex_Gcd(mi, a[i]);  
        ans = (ans+Powto(Powto(x, mi, sum), p[i], sum))%sum;            /*这里的x可能为负数，所以不能放后面乘*/  
    }  
    return (ans+sum)%sum;  
}  
int main(void)  
{  
    int T, n, i;  
    LL sum, ans;  
    scanf("%d", &T);  
    while(T--)  
    {  
        sum = 1;  
        scanf("%d", &n);  
        for(i=1;i<=n;i++)  
        {  
            scanf("%lld%lld", &a[i], &p[i]);  
            sum *= a[i];  
        }  
        ans = CRT(sum, n);  
        printf("%lld\n", ans);  
    }  
    return 0;  
}