中国剩余定理2.0

问题概述：一个由很多人组成的队伍，可以进行n中排队方式，若这些人a1个人一排会剩下b1个人，a2个人一排会

剩下b2个人……求这个队伍至少多少人

输入样例：                               对应输出：

1                                               1033349243908123

4

3181 285

3203 12

7907 2952

16183 2905

设这个队伍总人数为n，且满足：

n % a[1] = b[1];

n % a[2] = b[2];

n % a[3] = b[3];

……;

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

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

则最中答案n（最小可能）= ∑(mx[i]*my[i]*b[i])%∏(a[i])

#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;
}