BZOJ（1008）

bzoj1008: [HNOI2008]越狱

Description

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

Input

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

Output

可能越狱的状态数，模100003取余

Sample Input

2 3

Sample Output

6

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

虽然这题被称为bzoj上的水题，但是本弱鸡还是做了挺久的。不久前还写过一个快速幂，但还是不能独立完成，决定写个博客加深一下印象。

先算出所有的情况，因为一共有N个房间，M种宗教，

所有每个房间都有M种可能，所以总共有M^N种可能。

然后算出不可能越狱的情况，第一个房间有M种可能，因为相邻房间的人信仰的宗教不能相同，剩下的每个房间都有M-1种可能，所有不可能越狱的情况共有

M*(M-1)^(N-1)

可能越狱的情况是所有的情况减去不可能越狱的情况

M^N- M*(M-1)^(N-1)

然后写个快速幂就可以了。快速幂的时间复杂度是

O(logn)。下面讲一下快速幂是如何实现的。

比如说 a^10,因为10写成二进制形式是1010。

a^10= a^(1010)//(1010是二进制形式）

a^10= a^(1* 2 ^ 3) *a^(0 *2 ^ 2) *a^(1 * 2 ^ 1) *a^(0 *2 ^0) = a ^ 8 * a^(4 * 0) *a ^ (2) * a^(0 *1)

=a^8 * 1* a ^ 2 * 1

考虑到a ^ 2= a * a, a^4 = a^2 *a^2, a^8 =a ^ 4 * a ^ 4。

我们只要定义一个 t ,t=a,每次循环 t = t * t， 同时定义指数为b ,把b看成二进制形式，每次循环b都向右移一位，同时定义一个ans来保存结果，如果b&1==1时说明这时b的二进制表示形式的末尾是1，ans=ans * t，否则ans=ans。

下面粘贴一下AC代码

#include<bits/stdc++.h>
using namespace std;
const int MOD=100003;
long long powmod(long long x,long long y)
{
	long long ans=1,t=x;
	while(y)
	{
		if(y&1)
		ans=(ans*t)%MOD;
		t=(t*t)%MOD;
		y=y>>1;
	}
	return ans;
}


int main()
{
	long long m,n;
	cin>>m>>n;
	long long kiss1=powmod(m,n)%MOD;
	long long kiss2=((m%MOD)*(powmod(m-1,n-1))%MOD)%MOD;
	cout<<(MOD+kiss1-kiss2)%MOD<<endl;//加一个MOD防止kiss1-kiss2为负数
	return 0;
}

特别鸣谢https://blog.csdn.net/mosbest/article/details/69264953
和http://www.cnblogs.com/qscqesze/p/4348573.html