[HNOI2008]越狱

题面

<center> 1917: [HNOI2008]越狱 <center>
<center>时间限制:1秒 内存限制:162MB<center>

题目描述
　　监狱有连续编号为1...N的N个房间，每个房间关押一个犯人，有M种宗教，每个犯人可能信仰其中一种。如果
相邻房间的犯人的宗教相同，就可能发生越狱，求有多少种状态可能发生越狱

输入
　　输入两个整数M，N.1<=M<=10^8^,1<=N<=10^12^

输出
　　可能越狱的状态数，模100003取余
　　
样例输入

2 3

样例输出

6

提示
　　6种状态为(000)(001)(011)(100)(110)(111)

题解

　　每个房间的一个人信仰M个宗教中的一种，则N个房间共有M^N^种状态，其中不可越狱的状态：第一个房间的人可从M个宗教中选一，剩余N-1个房间的人不能和前一个人的宗教相同(否则会越狱)，有(M-1)个宗教选择，共M*(M-1)^(N-1)^种状态。则最终可越狱的状态：M^N^- M*(M-1)^(N-1)^。最终答案：((M^N^- M*(M-1)^(N-1)^))%mod。运用费马小定理(a^b^%p=a^b%(p-1)^%p)进行降幂,并运用各种求余公式整理得：
　　 (M^N%(mod-1)^%mod- M%mod*((M-1)^(N-1)%(mod-1)^%mod)+mod)%mod
最后，对其进行计算即可。
## 源码

#include <iostream>
using namespace std;
typedef long long ll;
const ll mod=1e5+3;
ll fun(ll a,ll b){
    ll i;
    ll ans=1;
    for(i=0;i<b;i++){
        ans*=a;
        ans%=mod;
    }
    return ans;
}
int main(int argc, char** argv) {
ll n,m,a,b,c,d,e;
cin>>m>>n;
a=n%(mod-1);
b=fun(m,a);
c=(n-1)%(mod-1);
d=fun(m-1,c);
m%=mod;
d*=m;
d%=mod;
e=(b-d+mod)%mod;
cout<<e;
return 0;
}