设为首页 加入收藏

TOP

BZOJ 1008 HNOI2008 越狱 快速幂
2015-07-20 17:42:13 来源: 作者: 【 】 浏览:1
Tags:BZOJ 1008 HNOI2008 越狱 快速

题目大意:有三名越狱犯杀死狱警后穿警服逃狱,现已被警方全部抓获……啊呸呸呸呸呸,扯远了

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

m<=10^8,n<=10^12

一看这数据范围就是logn嘛!于是我们来分析一下

f[i]表示连续i个房间可能发生越狱的情况数

当状态从f[i]转移到f[i+1]时有两种情况:

1.原先状态就会发生越狱 那么我们再塞一个人进去也一定会发生越狱 情况数为f[i]*m

2.原先状态不会发生越狱 那么我们只有塞进去一个与旁边人信仰相同的人才会越狱 情况数为m^i-f[i]
故得到递推式 f[i+1]=f[i]*(m-1)+m^i

这一看就是矩阵乘法!麻烦!不会写!妈蛋!

我们换一个思路,求一下补集

f[i]表示连续i个房间无法发生越狱的情况数

当状态从f[i]转移到f[i+1]时有两种情况:

1.原先状态就会发生越狱 那么我们再塞一个人进去也一定会发生越狱 情况数为0

2.原先状态不会发生越狱 那么我们只有塞进去一个与旁边人信仰不同的人才不会发生越狱 情况数为f[i]*(m-1)

故f[i+1]=f[i]*(m-1)

初值f[1]=m,最后得到结果用m^n相减即可

#include
  
   
#include
   
     #include
    
      #include
     
       #define M 100003ll using namespace std; typedef long long ll; ll m,n,ans; ll ksm(ll x,ll y) { ll re=1; while(y) { if(y&1)re*=x,re%=M; x*=x,x%=M; y>>=1; } return re; } int main() { cin>>m>>n; m%=M; ans=ksm(m,n); ans-=m*ksm(m-1,n-1); ans%=M; ans+=M; ans%=M; cout<
      
       

】【打印繁体】【投稿】【收藏】 【推荐】【举报】【评论】 【关闭】 【返回顶部
分享到: 
上一篇HUST - 1010 The Minimum Length 下一篇线性表 - 模拟指针实现线性表

评论

帐  号: 密码: (新用户注册)
验 证 码:
表  情:
内  容:

·工业机器人TCP校准中 (2025-12-25 05:19:17)
·opc 通讯协议与 TCP (2025-12-25 05:19:15)
·labview中tcp/ip通信 (2025-12-25 05:19:13)
·新书介绍《Python数 (2025-12-25 04:49:47)
·怎么利用 Python 进 (2025-12-25 04:49:45)