BZOJ 1008 越狱-阿里云开发者社区

BZOJ 1008 越狱

2017-03-04 988

版权

本文内容由阿里云实名注册用户自发贡献，版权归原作者所有，阿里云开发者社区不拥有其著作权，亦不承担相应法律责任。具体规则请查看《阿里云开发者社区用户服务协议》和《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容，填写侵权投诉表单进行举报，一经查实，本社区将立刻删除涉嫌侵权内容。

简介： 1008: [HNOI2008]越狱 Time Limit: 1 Sec Memory Limit: 162 MBSubmit: 8681 Solved: 3746[Submit][Status][Discuss] Description 　　监狱有连续编号为1...N的N个房间，每个房间关押一个犯人，有M种宗教，每个犯人可能信仰其中一种。

1008: [HNOI2008]越狱

Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 8681 Solved: 3746
[Submit][Status][Discuss]

Description

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

Input

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

Output

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

Sample Input

2 3

Sample Output

HINT

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

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1008

分析：快速幂取模运算，暂时先不解释，日后解释！

下面给出AC代码：

 1 #include <bits/stdc++.h>
 2 const int mod =100003;
 3 using namespace std;
 4 long long n,m;
 5 inline long long pow(long long a,long long b)
 6 {
 7     long long res=1;
 8     while(b)
 9     {
10         if(b&1) res=(res*a)%mod;
11         a=(a*a)%mod;
12         b>>=1;
13     }
14     return res;
15 }
16 int main()
17 {
18     cin>>m>>n;
19     long long a=pow(m,n);
20     long long b=(pow(m-1,n-1)*(m%mod))%mod;
21     a+=mod;
22     cout<<(a-b)%mod<<endl;
23    return 0;
24 }

更为简便的写法【之后补题的】

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 const int mod=100003;
 4 typedef long long ll;
 5 inline ll read()
 6 {
 7     ll x=0,f=1;
 8     char ch=getchar();
 9     while(ch<'0'||ch>'9')
10     {
11         if(ch=='-')
12             f=-1;
13         ch=getchar();
14     }
15     while(ch>='0'&&ch<='9')
16     {
17         x=x*10+ch-'0';
18         ch=getchar();
19     }
20     return x*f;
21 }
22 inline ll qpow(ll x,ll p)
23 {
24     ll ret=1;
25     for(;p;p>>=1,x=x*x%mod)
26     {
27         if(p&1)
28             ret=ret*x%mod;
29     }
30     return ret;
31 }
32 ll m,n;
33 int main()
34 {
35     m=read();
36     n=read();
37     ll a=qpow(m,n);
38     ll b=b=((m%mod)*qpow(m-1,n-1))%mod;
39     ll ans=(a+mod-b)%mod;
40     cout<<ans<<endl;
41     return 0;
42 }