中国剩余定理（模数互质 究极模板 快速乘防爆long long 附赠拓展欧几里德求逆元 ）

#include<iostream>
#include<algorithm>
#include<iomanip>
#include<set>
#include<queue>
#include<stack>
#include<vector>
#include<map>
#include<string.h>
#include<string>
#include<math.h>
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define VNAME(value) (#value)
#define debug(x) {cout<<VNAME(x)<<"="<<(x)<<endl;} 
using namespace std;
const int N = 200000;
const int inf=0x3f3f3f3f;
typedef long long ll;
inline ll ksc(ll x, ll y, ll mod)//快速乘 
{
    return ( x * y - (ll) ( (long double) x / mod*y )*mod + mod ) % mod;     
}
void exgcd(ll a,ll b,ll& d,ll& x,ll& y)//拓展欧几里得
{
	if(!b){d=a;x=1;y=0;}
	else {exgcd(b,a%b,d,y,x);y-=x*(a/b);} 
}
ll inverse(ll a,ll b)//求a对于b的逆元 
{
	ll d,x,y;
	exgcd(a,b,d,x,y);
	return d==1?(x%b+b)%b:-1;//如果两个数不互质，那么不存在逆元 
}
//中国剩余定理,pri[]存放互质的模数，r[]存放余数，n为个数 
ll china(ll pri[],ll r[],ll n){
	ll M=1,m,d,res=0;
	for(int i=0;i<n;i++) 
		M*=pri[i]; //所有数的最小公倍数
	for(int i=0;i<n;i++){
		m=M/pri[i]; //除以当前模数。 
		ll x,y;
		exgcd(pri[i],m,d,x,y); // pri[i]*x+m*y=gcd(pri[i],m])
		 y=(y%pri[i]+pri[i])%pri[i];
		res=(res+ksc(r[i],ksc(y,m,M),M))%M;//快速乘，防爆long long.  
	} 
	return (res%M+M)%M;  //满足的最小解 
}
ll a[15];
ll b[15]; 
int main()
{
	int n;
	cin>>n;
	for(int i=0;i<n;i++)
	{
		cin>>a[i]>>b[i];//a[]模数 b[]余数 
	}
	cout<<china(a,b,n)<<endl;
}

稍加更改可AC
poj 1006
牛客oj 曹冲养猪
洛谷oj P3868 [TJOI2009]猜数字