codevs 3990 CRT中国剩余定理

题意:

给定n个质数，以及k mod 这些质数的余数。问：在闭区间[a,b]中，有多少个k？最小的k是多少？

中国剩余定理（CRT）的板子题

完整代码如下

#include<iostream>
#include<cstdio>
#include<algorithm>
#define ll long long
using namespace std;
ll l,r;
ll n,m[11],p[11];
void exgcd(ll a,ll b,ll &x,ll &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return;
	}
	else
	{
		exgcd(b,a%b,y,x);
		y-=a/b*x;
	}
}
int main()
{
	ll M=1,ans=0,cnt=0;//cnt 区间内解的个数 
	scanf("%lld%lld%lld",&n,&l,&r);
	for(int i=1;i<=n;++i)
	{
		scanf("%lld%lld",&p[i],&m[i]);
		M*=p[i];
	}
	for(int i=1;i<=n;++i)
	{
		ll x,y;
		ll Mi=M/p[i];
		exgcd(Mi,p[i],x,y);
		x=(x%p[i]+p[i])%p[i];
		if(!x)
			x+=p[i];
		ans=(ans+Mi*x*m[i])%M;
	}
	if(r>=ans)
		cnt=(r-ans)/M+1;
	if(l-1>=ans)
		cnt-=(l-ans)/M+1;
	if(ans<l)
		ans+=((l-ans)/M+1)*M;
	cout<<cnt<<endl;
	if(cnt==0)
		cout<<"0";
	else
		cout<<ans;
	return 0;
}

分析问题，我们发现：即是对这样一个方程组求解

我们定义

   

为 mod 意义下的逆元

根据百度百科可证

所以方程组所有的解的集合就是：

 

对于逆元，我们使用扩展欧几里得定理来求

下面背板

void exgcd(ll a,ll b,ll &x,ll &y)
{
	if(b==0)
	{
		x=1;
		y=0;
		return;
	}
	else
	{
		exgcd(b,a%b,y,x);
		y-=a/b*x;
	}
}

	exgcd(Mi,p[i],x,y);
	x=(x%p[i]+p[i])%p[i];
	if(!x)
		x+=p[i];

以下六行代码是本题的关键

	if(r>=ans)
		cnt=(r-ans)/M+1;
	if(l-1>=ans)
		cnt-=(l-ans)/M+1;
	if(ans<l)
		ans+=((l-ans)/M+1)*M;

如果我们得出的最小解ans在比区间的右端点r大

那么右端点r与ans的差值中M的个数再加上ans这个解(是+1不是+ans)就是[1,r]解的总数

同理

(l-ans)/M+1;

是[1,l-1]中解的个数，所以最终cnt就是[l,r]解的个数。

如果当前的ans不在区间[l,r]内，根据CRT通解，我们要加上k倍的M使其在区间内

而k的值就是l与ans的差值中有几倍的M再+1使其在*M+ans后在区间内。