博客园 首页 私信博主 显示目录 隐藏目录 管理 动画

POJ.2891.Strange Way to Express Integers(扩展CRT)

题目链接
扩展中国剩余定理:1(直观的)2(详细证明)

[Upd:]https://www.luogu.org/problemnew/solution/P4774


#include <cstdio>
#include <cctype>
#define gc() getchar()
typedef long long LL;
const int N=1e6+5;

LL n,m[N],r[N];

inline LL read()
{
	LL now=0,f=1;register char c=gc();
	for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
	for(;isdigit(c);now=now*10+c-'0',c=gc());
	return now*f;
}
LL Exgcd(LL a,LL b,LL &g,LL &x,LL &y)
{
	if(!b) g=a, x=1ll, y=0ll;
	else Exgcd(b,a%b,g,y,x),y-=a/b*x;
}
LL Ex_CRT()
{
	LL M=m[1],R=r[1],x,y,g,t;
	for(int i=2; i<=n; ++i)
	{
		Exgcd(M,m[i],g,x,y);
		if((r[i]-R)%g) return -1ll;
		x*=(r[i]-R)/g, t=m[i]/g, x=(x%t+t)%t;//相当于M*(x/((ri-R)/g)) ≡ g(mod mi/g)?不管了就这么理解吧 
		R+=M*x, M*=t, R%=M;
	}
	return (R%M+M)%M;
}

int main()
{
	while(~scanf("%lld",&n))
	{
		for(int i=1; i<=n; ++i) m[i]=read(),r[i]=read();
		printf("%lld\n",Ex_CRT());
	}
	return 0;
}
posted @ 2018-02-19 17:31  SovietPower  阅读(282)  评论(0编辑  收藏  举报