GCD
LLI gcd(LLI a, LLI b) {
return b == 0 ? a : gcd(b, a % b);
}
LCM
关于最小公倍数的结论:如果lcm(a,b) == n
其中
n=p1e1* p2e2 * p3e3…… * pnex
a=p1a1* p2a2 * p3a3…… * pnax
b=p1b1 * p2b2 * p3b3…… * pnbx
则有e1=max(a1,b1);……ex=max(ax,bx);
扩展欧几里得
LLI Ex_gcd(LLI a,LLI b,LLI &x,LLI &y) {
if(b == 0) {
x = 1;
y = 0;
return a;
}
LLI ans = Ex_gcd(b,a % b,x,y);
LLI temp = x;
x = y;
y = temp - a / b * y;
return ans;
}
求出a * x + b * y = gcd(a,b)的一个特解x0,y0,通解
x = x0 + b / gcd(a,b) * t, y = y0 – a / gcd(a,b) * t (t是任意整数);
而对于a * x + b * y = c的一个特解x1 = x0 * c / gcd(a,b), y1 = y0 * c / gcd(a,b);
通解:x = x1 + b / gcd(a,b) * t, y = y1 – a / gcd(a,b) * t
中国剩余定理
LLI Chinese_Remainder(LLI n) {
LLI M = m[1],R = r[1],x,y;//m[]和r[]数组都是从1到n的
for(int i = 2; i <= n; i ++) {
LLI d = Ex_gcd(M,m[i],x,y);
LLI c = r[i] - R;
if(c % d) return -1;
LLI temp = Ex_gcd(M / d,m[i] / d,x,y);//调用扩展欧几里得
x = (c / d * x) % (m[i] / d);
R = R + x * M;
M = M / d * m[i];
R = R % M;
}
if(R < 0) R += M;
return R;
}