Fast power, turtle speed multiplication, light speed multiplication template

fast power

abmodp ( 1 ≤ a , b , p ≤ 1 0 18 ) a^b mod p (1 \and a, b, p \and 10^{18})ab modp(1a,b,p1018)

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})ab m o d p ( 1 a,b,p1018)

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})ab m o d p ( 1 a,b,p1018)

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

Guess you like

Origin blog.csdn.net/m0_51156601/article/details/124181855