题解-[ZJOI2005]沼泽鳄鱼
题解-[ZJOI2005]沼泽鳄鱼
前置知识:
邻接矩阵
矩阵乘法
矩阵快速幂
给一个有 \(N\) 个点,\(M\) 条双向边的图 \(G\),其中有 \(NFish\) 只鳄鱼以 \(T\) 个点 \(P_0\sim P_t\) 为周期运动。求从 \(Start\) 出发到 \(End\) 不停留走 \(K\) 步每步不碰到鳄鱼的方案数(节点下标从 \(0\) 开始编号)。
数据范围:\(1\le N\le 50\),\(1\le K\le 2\times 10^9\),\(1\le NFish\le 20\),\(2\le T\le 4\)。
看到这个 \(K\) 的范围就知道要带个 \(\log\),而在图上可以带 \(\log\) 的算法,唯有二分、倍增(树形图)和邻接矩阵快速幂(小图)。看到这题 \(N\) 很小,又无从二分,便知道可以邻接矩阵快速幂。
首先要知道什么是矩阵乘法。大小为 \(a\times b\) 的矩阵 \(A\) 和大小为 \(b\times c\) 的矩阵方可相乘,乘积矩阵 \(C\) 大小为 \(a\times c\),满足 \(C_{i,j}=\sum\limits_{k=1}^bA_{i,k}\cdot B_{k,j}(1\le i\le a,1\le j\le c)\)。
矩阵乘法可以解决集合的定向求和变换操作。其中原矩阵乘以变换矩阵变成目标矩阵。
矩阵快速幂:原矩阵多次定向变换,可以通过多次乘以变换矩阵解决。矩阵乘法满足结合律,所以可以先求出变换矩阵的幂。矩阵也可以快速幂,时间复杂度 \(\Theta(n^3\log k)\)。
邻接矩阵:用于表示图边,在无权图中,如果矩阵 \(M\) 的元素 \(M_{i,j}=1\),表示存在边 \((i,j)\)。
而邻接矩阵快速幂就是建立在邻接矩阵和矩阵快速幂上的。原矩阵表示起点状态,变换矩阵为邻接矩阵,每乘一次就表示不停留地走一步后的状态,乘 \(K\) 次就成了最终状态,取终点矩阵值则为答案。
状态:到每个节点的方案数。
Example
讲解时暂时初始下标为 \(1\) 吧,暂时不考虑鳄鱼。
例如 \(N=3\),\(Start=1\),\(End=3\),有双向边 \({(1,3),(3,2)}\),求走 \(K=3\) 步后到终点 \(End\) 的方案数。
原始状态:
\begin{bmatrix}
1 & 0 & 0\\
\end{bmatrix}
\]
变换矩阵为邻接矩阵:
\begin{bmatrix}
0 & 0 & 1\\
0 & 0 & 1\\
1 & 1 & 0\\
\end{bmatrix}
\]
目标状态:
C=&A\times B^K\\
=&
\begin{bmatrix}
1 & 0 & 0\\
\end{bmatrix}
\times
\begin{bmatrix}
0 & 0 & 1\\
0 & 0 & 1\\
1 & 1 & 0\\
\end{bmatrix}
^3\\
=&
\begin{bmatrix}
1 & 0 & 0\\
\end{bmatrix}
\times
\begin{bmatrix}
0 & 0 & 2\\
0 & 0 & 2\\
2 & 2 & 0\\
\end{bmatrix}\\
=&
\begin{bmatrix}
0 & 0 & 2\\
\end{bmatrix}
\end{split}
\]
所以到终点 \(End=3\) 的方案数为 \(2\)。
但是鳄鱼怎么办呢?
可以发现 \(2\le T\le 4\),所以所有鳄鱼运动的总周期是 \(12\)。可以计算 \(12\) 种变换邻接矩阵,表示到每个周期时可以走的边。然后顺次相乘,求乘积的 \(\lfloor\frac K{12}\rfloor\) 次幂(可以用矩阵快速幂)乘以 前 \(\left(K\bmod 12\right)\) 种变换矩阵的乘积,就是总变换矩阵。用表示起点的原矩阵乘以总变换矩阵,即可得答案。
时间复杂度 \(\Theta(NFish+N^3\log K)\)(注意,\(NFish\) 为一个完整变量名)。
真的难讲,还是放代码吧。
Code
#include <bits/stdc++.h>
using namespace std;
//&Start
#define inf 0x3f3f3f3f
#define re register
#define il inline
#define hash unorded_map
typedef long long lng;
typedef unsigned long long ulng;
typedef vector<int> veci;
#define fo(i,st,xb,y) for(re int i=st;i xb;i y)
//&Data
#define N 50
#define mod 10000
int n,m,s,t,k,fish,p[25][12];
//&Matrix
struct Matrix{ //矩阵
int arr[N+5][N+5];
Matrix(re int op=0){
memset(arr,0,sizeof arr);
if(op==1) fo(i,1,<=N,++) arr[i][i]=1; //构造单位矩阵 E 满足 E*A=A*E=A
}
il int* operator[](re int x){return arr[x];}
il friend Matrix operator*(re Matrix x,re Matrix y){
re Matrix res;
fo(k,1,<=N,++)fo(i,1,<=N,++)fo(j,1,<=N,++)
(res[i][j]+=x[i][k]*y[k][j]%mod)%=mod;
return res;
}
il void print(re char*s){
puts(s);
fo(i,1,<=n,++)fo(j,1,<=n,++)
printf("%d%c",arr[i][j],"\n "[j<n]);
}
}st,e,g[12],all,ans;
il Matrix Pow(re Matrix a,re int x){ //矩阵快速幂
re Matrix res(1);
for(;x;a=a*a,x>>=1)if(x&1) res=res*a;
return res;
}
//&Main
int main(){
scanf("%d%d%d%d%d",&n,&m,&s,&t,&k);
s++,t++,st[1][s]=1; //原矩阵
fo(i,1,<=m,++){
re int u,v;scanf("%d%d",&u,&v);
e[u+1][v+1]=e[v+1][u+1]=1; //邻接矩阵
}
scanf("%d",&fish);
fo(i,1,<=fish,++){
re int tmp; scanf("%d",&tmp);
fo(j,0,<tmp,++) scanf("%d",p[i]+j),p[i][j]++;
fo(j,tmp,<12,++) p[i][j]=p[i][j%tmp];
}
re int fb[N+5];
fo(i,0,<12,++){
fill(fb+1,fb+n+1,0);
fo(j,1,<=fish,++) fb[p[j][i]]=1;
fo(j,1,<=n,++)fo(k,1,<=n,++)
if(!fb[k]) g[i][j][k]=e[j][k]; //12种变换矩阵
}
all=Matrix(1);
fo(i,1,<12,++) all=all*g[i];
all=all*g[0]; //12个矩阵顺次相乘(1,2,...,11,0)
ans=st*Pow(all,k/12); //原矩阵乘以变换矩阵
fo(i,1,<=k%12,++) ans=ans*g[i]; //乘以剩余 K%12 个矩阵
printf("%d\n",ans[1][t]); //答案
return 0;
}
祝大家学习愉快!
题解-[ZJOI2005]沼泽鳄鱼的更多相关文章
- 【BZOJ1898】[ZJOI2005]沼泽鳄鱼(矩阵快速幂,动态规划)
[BZOJ1898][ZJOI2005]沼泽鳄鱼(矩阵快速幂,动态规划) 题面 BZOJ 洛谷 题解 先吐槽,说好了的鳄鱼呢,题面里面全是食人鱼 看到数据范围一眼想到矩乘. 先不考虑食人鱼的问题,直接 ...
- P2579 [ZJOI2005]沼泽鳄鱼(矩乘)
P2579 [ZJOI2005]沼泽鳄鱼 没有食人鱼:直接矩乘优化 有食人鱼:食人鱼周期2.3.4,公倍数12,可以以12为一个周期矩乘,剩下的暴力 注意矩乘不满足乘法交换律,一定要注意乘的顺序 #i ...
- [ZJOI2005]沼泽鳄鱼 矩阵乘法
---题面--- 题解: 乍一看还是挺懵逼的.和HH去散步很像,思路也是类似的. 复制一段我在HH去散步的题解里面写的一段话吧: 考虑f[i][j]表示i和j是否右边相连,有为1,否则为0,那么f同时 ...
- [ZJOI2005]沼泽鳄鱼
题目描述 潘塔纳尔沼泽地号称世界上最大的一块湿地,它地位于巴西中部马托格罗索州的南部地区.每当雨季来临,这里碧波荡漾.生机盎然,引来不少游客. 为了让游玩更有情趣,人们在池塘的中央建设了几座石墩和石桥 ...
- [luogu2579 ZJOI2005] 沼泽鳄鱼(矩阵快速幂)
传送门 题目描述 潘塔纳尔沼泽地号称世界上最大的一块湿地,它地位于巴西中部马托格罗索州的南部地区.每当雨季来临,这里碧波荡漾.生机盎然,引来不少游客. 为了让游玩更有情趣,人们在池塘的中央建设了几座石 ...
- bzoj1898 [Zjoi2005]沼泽鳄鱼
Description 潘塔纳尔沼泽地号称世界上最大的一块湿地,它地位于巴西中部马托格罗索州的南部地区.每当雨季来临,这里碧波荡漾.生机盎然,引来不少游客.为了让游玩更有情趣,人们在池塘的中央建设了几 ...
- P2579 [ZJOI2005]沼泽鳄鱼
传送门 话说邻接矩阵居然还能快速幂的么-- 把原图的邻接矩阵\(G\)打出来,那么\(G[u][v]\)表示一秒后\(u\)到\(v\)的方案数,\(G^2[u][v]\)表示\(2\)秒后的方案数- ...
- ZJOI2005沼泽鳄鱼
矩阵优化dp ** 注意:矩阵乘法没有交换律 ** 思路:类比P2151hh去散步 代码特点在一维的答案矩阵 1.矩阵优化两点间方案数不必赘述 2.注意2,3,4可以办到以他们的lcm为周期,正是因为 ...
- P2579 [ZJOI2005]沼泽鳄鱼(邻接矩阵,快速幂)
题目简洁明了(一点都不好伐) 照例,化简题目 给一张图,每一个时间点有一些点不能走,(有周期性),求从起点第k秒恰好在终点的方案数,可重复,不可停留. 额dp实锤 于是就被打脸了.... 有一种东西叫 ...
随机推荐
- yum安装Ceph指定Jewel版本
前言 通过yum安装指定的rpm包,这个一般是 yum --showduplicates list ceph | expand ,然后去通过yum安装指定的版本即可,这个在hammer下是没有问题的, ...
- Redis订阅
1.Redis订阅简介 进程间的一种消息通信模式:发送者(pub)发送消息,订阅者(sub)接收消息. 2.Redis订阅命令 3.Redis订阅的使用 先订阅后发布后才能收到消息, 1 可以一次性订 ...
- 学会网页制作,web app开发,必须先从语法基础开始学习
做软件开发,是从事编程开发工作,必须先从语法基础开始学习,通过语法组成产品效果. 前端开发的基础语法,由HTML+CSS+JavaScript组成,这是前端开发最基本的3个语言. ...
- 关于steam平台“wallpaper engine”软件出现界面黑屏,但壁纸能播放的问题
前阵子重装电脑后,在使用wallpaper engine这款软件时发现了以下令人疑惑的画面: 点击"设置"和"壁纸选择"界面全是黑的......这还没完,更气人 ...
- 来吧,展示!SpringBoot OSS 整合全过程,没见过比这更详细的了
前言 阿里云对象存储服务(Object Storage Service,简称 OSS),是阿里云提供的海量.安全.低成本.高可靠的云存储服务.其数据设计持久性不低于 99.9999999999%(12 ...
- 新鲜出炉!2020年最新java面试题大全,面试突击必备!
前言 发现网上很多Java面试题都没有答案,所以花了很长时间搜集整理出来了一套Java面试题,希望对大家有帮助哈~ 打算这几天每天更新15~20题.(这样有助于你们阅读和理解!)我们先从简单的开始 1 ...
- Model class apps.goods.models.GoodsType doesn't declare an explicit app_label and isn't in an application in INSTALLED_APPS
在admin.py注册这个model时,报了个错: RuntimeError: Model class apps.goods.models.GoodsType doesn't declare an e ...
- Appium上下文和H5测试(一)
坚持原创输出,点击蓝字关注我吧 作者:清菡 博客:oschina.云+社区.知乎等各大平台都有. 目录 一.混合应用-H5 1.混合应用是什么? 2.怎么样分辨一个 App 页面究竟是原生的还是 We ...
- Java 枚举 enum 详解
本文部分摘自 On Java 8 枚举类型 Java5 中添加了一个 enum 关键字,通过 enum 关键字,我们可以将一组拥有具名的值的有限集合创建为一种新的类型,这些具名的值可以作为常规的程序组 ...
- 在IDM上设置防止过度抓取网站信息
在使用Internet Download Manager(IDM)下载器时,有时会发现IDM自带的抓取功能过于强大,以至于有时会抓取一些无效的链接.那么,该如何避免IDM的过度抓取呢? 图1:IDM的 ...