fast power
求abmodp ( 1 ≤ a , b , p ≤ 1 0 18 ) a^b mod p (1 \and a, b, p \and 10^{18})ab modp(1≤a,b,p≤1018)
ll quick_power(ll a, ll b, ll p) {
// ll 是long long类型
ll ans = 1 % p;
for(; b; b >>= 1) {
if(b & 1) ans = ans * a % p;
a = a * a % p;
}
return ans;
}
Turtle Speed Multiplication (Solution 1)
Let a ∗ b modp ( 1 ≤ a , b , p ≤ 1 0 18 ) a*b \ mod p (1 \and a, b, p \ and 10^{18})a∗b m o d p ( 1 ≤a,b,p≤1018)
ll quick_mul(ll a, ll b, ll p) {
// ll 是long long类型
ll ans = 0;
for(; b; b >>= 1) {
if(b & 1) ans = (ans + a) % p;
a = a * 2 % p;
}
return ans;
}
Multiplication of the speed of light (solution 2)
Let a ∗ b modp ( 1 ≤ a , b , p ≤ 1 0 18 ) a*b \ mod p (1 \and a, b, p \ and 10^{18})a∗b m o d p ( 1 ≤a,b,p≤1018)
ull mul(ull a, ull b, ull p) {
// ull 是 unsigned long long 类型,ll 是long long类型
a %= p, b %= p;
ull c = (long double)a * b / p;
ull x = a * b, y = c * p;
ll ans = (long long)(x % p) - (long long)(y % p);
if(ans < 0) ans += p;
return ans;
}