题解-[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\) 的方案数。

原始状态:

\[A=
\begin{bmatrix}
1 & 0 & 0\\
\end{bmatrix}
\]

变换矩阵为邻接矩阵:

\[B=
\begin{bmatrix}
0 & 0 & 1\\
0 & 0 & 1\\
1 & 1 & 0\\
\end{bmatrix}
\]

目标状态:

\[\begin{split}
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]沼泽鳄鱼的更多相关文章

  1. 【BZOJ1898】[ZJOI2005]沼泽鳄鱼(矩阵快速幂,动态规划)

    [BZOJ1898][ZJOI2005]沼泽鳄鱼(矩阵快速幂,动态规划) 题面 BZOJ 洛谷 题解 先吐槽,说好了的鳄鱼呢,题面里面全是食人鱼 看到数据范围一眼想到矩乘. 先不考虑食人鱼的问题,直接 ...

  2. P2579 [ZJOI2005]沼泽鳄鱼(矩乘)

    P2579 [ZJOI2005]沼泽鳄鱼 没有食人鱼:直接矩乘优化 有食人鱼:食人鱼周期2.3.4,公倍数12,可以以12为一个周期矩乘,剩下的暴力 注意矩乘不满足乘法交换律,一定要注意乘的顺序 #i ...

  3. [ZJOI2005]沼泽鳄鱼 矩阵乘法

    ---题面--- 题解: 乍一看还是挺懵逼的.和HH去散步很像,思路也是类似的. 复制一段我在HH去散步的题解里面写的一段话吧: 考虑f[i][j]表示i和j是否右边相连,有为1,否则为0,那么f同时 ...

  4. [ZJOI2005]沼泽鳄鱼

    题目描述 潘塔纳尔沼泽地号称世界上最大的一块湿地,它地位于巴西中部马托格罗索州的南部地区.每当雨季来临,这里碧波荡漾.生机盎然,引来不少游客. 为了让游玩更有情趣,人们在池塘的中央建设了几座石墩和石桥 ...

  5. [luogu2579 ZJOI2005] 沼泽鳄鱼(矩阵快速幂)

    传送门 题目描述 潘塔纳尔沼泽地号称世界上最大的一块湿地,它地位于巴西中部马托格罗索州的南部地区.每当雨季来临,这里碧波荡漾.生机盎然,引来不少游客. 为了让游玩更有情趣,人们在池塘的中央建设了几座石 ...

  6. bzoj1898 [Zjoi2005]沼泽鳄鱼

    Description 潘塔纳尔沼泽地号称世界上最大的一块湿地,它地位于巴西中部马托格罗索州的南部地区.每当雨季来临,这里碧波荡漾.生机盎然,引来不少游客.为了让游玩更有情趣,人们在池塘的中央建设了几 ...

  7. P2579 [ZJOI2005]沼泽鳄鱼

    传送门 话说邻接矩阵居然还能快速幂的么-- 把原图的邻接矩阵\(G\)打出来,那么\(G[u][v]\)表示一秒后\(u\)到\(v\)的方案数,\(G^2[u][v]\)表示\(2\)秒后的方案数- ...

  8. ZJOI2005沼泽鳄鱼

    矩阵优化dp ** 注意:矩阵乘法没有交换律 ** 思路:类比P2151hh去散步 代码特点在一维的答案矩阵 1.矩阵优化两点间方案数不必赘述 2.注意2,3,4可以办到以他们的lcm为周期,正是因为 ...

  9. P2579 [ZJOI2005]沼泽鳄鱼(邻接矩阵,快速幂)

    题目简洁明了(一点都不好伐) 照例,化简题目 给一张图,每一个时间点有一些点不能走,(有周期性),求从起点第k秒恰好在终点的方案数,可重复,不可停留. 额dp实锤 于是就被打脸了.... 有一种东西叫 ...

随机推荐

  1. FreeSql接入CAP的实践

    CAP CAP 是一个基于 .NET Standard 的 C# 库,它是一种处理分布式事务的解决方案,同样具有 EventBus 的功能,它具有轻量级.易使用.高性能等特点. https://git ...

  2. python3中的os.path模块

    os.path模块主要用于获取文件的属性,这里对该模块中一些常用的函数做些记录. os.abspath(path):获取文件的绝对路径.这里path指的是路径,例如我这里输入"data.cs ...

  3. TCP/IP模型简介和/etc/hosts文件说明

    软件=协议的实现. IP决定了主机的位置.端口号决定了进程的位置. 两台主机上的通讯实际是两台主机上两个具体进程的通讯. TCP/IP模型分四层: TCP/IP模型:应用层---传输层----网络层- ...

  4. fio测试ceph的filestore

    前言 fio是一个适应性非常强的软件,基本上能够模拟所有的IO请求,是目前最全面的一款测试软件,之前在看德国电信的一篇分享的时候,里面就提到了,如果需要测试存储性能,尽量只用一款软件,这样从上层测试到 ...

  5. 重构rbd镜像的元数据

    这个已经很久之前已经实践成功了,现在正好有时间就来写一写,目前并没有在其他地方有类似的分享,虽然我们自己的业务并没有涉及到云计算的场景,之前还是对rbd镜像这一块做了一些基本的了解,因为一直比较关注故 ...

  6. Git-stash(暂存)

    修改某文件后,不想commit,使用stash保存在本地的某分支内 # 暂存 git stash ## 可暂存新增文件 git stash -u ## 为此次暂存添加标识 git stash save ...

  7. linux 身份鉴别口令复杂度整改

    口令复杂度: 1.首先安装apt install libpam-cracklib -y2.vim /etc/pam.d/common-password3.在第2步末尾添加password requis ...

  8. phpstorm 远程调式 php

    https://cloud.tencent.com/developer/article/1561767 超时设置 fastcgi: ``` 1. apache module的情况下: 修改配置文件 h ...

  9. Python自动玩"别再踩白块了"小游戏

    开发工具 Python版本:3.6.4 相关模块: numpy模块: CV2模块: pyautogui模块: pyscreenshot模块: 以及一些Python自带的模块. 环境搭建 安装Pytho ...

  10. Markdown的应知应会

    Markdown介绍 什么是Markdown Markdown是一种纯文本.轻量级的标记语言,常用作文本编辑器使用.和记事本.notepad++相比,Markdown可以进行排版:和Word相比,Ma ...