【中国剩余定理】洛谷P3868 [TJOI2009]猜数字

题目描述

现有两组数字，每组k个，第一组中的数字分别为：a1，a2，...，ak表示，第二组中的数字分别用b1，b2，...，bk表示。其中第二组中的数字是两两互素的。求最小的非负整数n，满足对于任意的i，n - ai能被bi整除。

格式输入

输入数据的第一行是一个整数k，（1 ≤ k ≤ 10）。接下来有两行，第一行是：a1，a2，...，ak，第二行是b1，b2，...，bk 

格式输出

输出所求的整数n。

输入样例

3
1 2 3
2 3 5

 输出样例

23

说明/提示

所有数据中，第一组数字的绝对值不超过10^9（可能为负数），第二组数字均为不超过6000的正整数，且第二组里所有数的乘积不超过10^18

每个测试点时限1秒

注意：对于C/C++语言，对64位整型数应声明为long long，如使用scanf, printf函数（以及fscanf, fprintf等），应采用%lld标识符。

题目分析：

其实是中国剩余定理，这里niiick解释的非常好，我是照他的思路写的，不会的去看看！

题目链接：https://www.luogu.org/problem/P3868

AC代码

#include<bits/stdc++.h>
using namespace std;
const int maxn=1e6+5;
typedef long long ll;
int a[20],b[20];
int n;
ll quick(ll a,ll b,ll mod)
{
	ll ans=0;
	while(b)
	{
		if(b%2)
			ans=(ans+a)%mod;
		a=(a+a)%mod;
		b/=2;
	}
	return ans;
}
ll x,y;
void exgcd(ll a1,ll &x,ll b1,ll &y)
{
	if(b1==0){
		x=1;y=0;return ;
	}
	exgcd(b1,y,a1%b1,x);
	y-=(a1/b1)*x;
}
ll calc()
{
	ll ans=0,tot=1,x,y;
	for(int i=1;i<=n;++i) 
		tot*=b[i];
    for(int i=1;i<=n;++i)
    {
        ll tmp=tot/b[i];
        exgcd(tmp,x,b[i],y);
        x=(x%b[i]+b[i])%b[i];
        ans=(ans+quick(quick(tmp,x,tot),a[i],tot))%tot;
    }
    return (ans+tot)%tot;
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%d",&a[i]);
	for(int i=1;i<=n;i++)
		scanf("%d",&b[i]);
	for(int i=1;i<=n;i++)
		a[i]=(a[i]%b[i]+b[i])%b[i];
	printf("%lld\n",calc());
	return 0;
}