TJOI 2009 猜数字 题解

题目传送门

题目大意： 求下面这个方程组的最小正整数解。
$$\{\begin{array}{l}(n-a_1)|b_1\\ (n-a_2)|b_2\\ ⋮\\ (n-a_k)|b_k\end{array}$$

题解

我们把柿子变化一下：
$$\{\begin{array}{l}(n-a_1)|b_1\\ (n-a_2)|b_2\\ ⋮\\ (n-a_k)|b_k\end{array}\Rightarrow \{\begin{array}{l}n-a_1\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}{b}_1)\\ n-a_2\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}{b}_2)\\ ⋮\\ n-a_k\equiv 0(\mathrm{m}\mathrm{o}\mathrm{d}{b}_k)\end{array}\Rightarrow \{\begin{array}{l}n\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{b}_1)\\ n\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{b}_2)\\ ⋮\\ n\equiv a_k(\mathrm{m}\mathrm{o}\mathrm{d}{b}_k)\end{array}$$

然后就是个中国剩余定理裸题了。

但是题目里面说 $\prod b_i\le 10^{18}$，这意味着中间的运算依然可能爆 $longlong$，所以要用快速乘。

以及，$b_i$ 不一定是个质数，所以求逆元需要用 exgcd。

代码如下：

#include <cstdio>
#include <cstring>
#define ll long long

int n;
ll a[11],b[11],t[11],M=1,ans=0;
void exgcd(ll a,ll b,ll &x,ll &y)
{
	if(b==0){x=1,y=0;return;}
	exgcd(b,a%b,y,x);
	y-=a/b*x;
}
ll inv(ll a,ll mod)
{
	ll x,y;
	exgcd(a,mod,x,y);
	return (x+mod)%mod;
}
ll ksc(ll x,ll y)
{
	ll re=0,tot=x;
	while(y)
	{
		if(y&1ll)re=(re+tot)%M;
		tot=(tot+tot)%M;
		y>>=1ll;
	}
	return re;
}

int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	scanf("%lld",&a[i]);
	for(int i=1;i<=n;i++)
	scanf("%lld",&b[i]),a[i]=(a[i]%b[i]+b[i])%b[i],M*=b[i];
	for(int i=1;i<=n;i++)
	t[i]=inv(M/b[i],b[i]);
	for(int i=1;i<=n;i++)
	ans=(ans+ksc(ksc(a[i],t[i]),M/b[i]))%M;
	printf("%lld",ans);
}