BZOJ1965洗牌(快速幂+快速乘+逆元)

思路:找规律题,可以发现,牌x移动后就会到2*x%(n+1)的位置,移动m次后就是2^{m}x\equiv L( mod n+1),则x=\tfrac{L}{2^m}(modn+1)

求2^m关于n+1的逆元就是2关于n+1的逆元。

#include<bits/stdc++.h>
using namespace std;
const int N = 2e5 + 7;
typedef long long ll;
ll n, m, l, mod;
ll a, b, x, y;
ll exgcd(ll a, ll b) {
    if(!b) {
        x = 1, y = 0;
        return a;
    }
    ll d = exgcd(b, a % b);
    ll t = x;
    x = y;
    y = t - (a / b) * y;
    return d;
}
ll mul(ll a, ll b, ll p) {
    ll ans = 0;
    for(ll i = b; i; i >>= 1, a = (a + a) % p)
        if(i & 1)
            ans = (ans + a) % p;
    return ans;
}
ll poww(ll a, ll b) {
    ll ans = 1;
    while(b) {
        if(b & 1) {
            ans = mul(ans, a, mod);//快速乘防止爆ll
        }
        a = mul(a, a, mod);
        b >>= 1;
    }
    return ans;
}
int main() {
    cin >> n >> m >> l;
    mod = n + 1;
    exgcd(2, n + 1);//计算逆元
    while(x < 0)
        x += (n + 1);
    ll niyuan = x;
    niyuan = poww(niyuan, m);
    x = mul(niyuan, l, mod);
    cout << x;
    return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值