信息学奥赛一本通 1618：越狱

早睡身体好呀

已于 2022-07-25 12:50:57 修改

阅读量236

点赞数

分类专栏：快速幂文章标签： c++ 快速幂

于 2019-05-05 22:14:51 首次发布

本文链接：https://blog.csdn.net/qq_42815188/article/details/89856571

版权

快速幂专栏收录该内容

6 篇文章 0 订阅

订阅专栏

题目链接：http://ybt.ssoier.cn:8088/problem_show.php?pid=1618

思路：

总方案数为 $m*m*…*m=m^n$

不发生越狱的方案数为 $m*(m-1)*...*(m-1) = m*(m-1)^{n-1}$

故发生越狱的方案数为 $m^n-m*(m-1)^{n-1}$ 。

对答案取余 $m^n-m*(m-1)^{n-1}) \ \% \ mod$ ，即 $m^n \ \% \ mod - m*(m-1)^{n-1} \ \% \ mod) \ \% \ mod$ ，

但是 $m^n \ \% \ mod$ 减去 $m*(m-1)^{n-1} \ \% \ mod$ 可能是负数，

所以还需要把负余数转换成正余数，即 $m^n \ \% \ mod - m*(m-1)^{n-1} \ \% \ mod + mod) \ \% \ mod$ 。

#include <iostream>

using namespace std;

typedef long long ll;

const int mod = 100003;

int qmi(int a, ll b, int m)
{
    int res = 1;
    while (b)
    {
        if (b & 1) res = (ll)res * a % m;
        a = (ll)a * a % m;
        b >>= 1;
    }
    return res;
}

int main()
{
    int m;
    ll n;
    cin >> m >> n;
    
    cout << (qmi(m, n, mod) - (ll)m * qmi(m - 1, n - 1, mod) % mod + mod) % mod << endl;
    
    return 0;
}