正題

題目鏈接:https://www.luogu.org/problemnew/show/P2579

---

題目大意

一張無向圖，一個(gè)起點(diǎn)一個(gè)終點(diǎn)。

有食人魚，在若干個(gè)點(diǎn)之間有周期的移動(dòng)，周期為$2$或$3$或$4$個(gè)點(diǎn)為循環(huán)。

然后要求從起點(diǎn)到終點(diǎn)走$k$步且不碰到食人魚的方案數(shù)。

---

解題思路

因?yàn)槭橙唆~的周期為$2$或$3$或$4$，所以顯然所有食人魚的周期為$lcm(2,3,4)=12$一個(gè)循環(huán)。
然后$k$很大顯然矩陣乘法。

轉(zhuǎn)移矩陣$G(k)$中$[i,j]$表示在第$k$個(gè)時(shí)刻是否可以從$i$移動(dòng)到$j$。需要考慮連邊和食人魚具體分析。

構(gòu)建好轉(zhuǎn)移矩陣后$G(k)$我們開始轉(zhuǎn)移，我們發(fā)現(xiàn)如果不考慮快速冪轉(zhuǎn)移那么
$$Ans=G(1)?G(2)?G(3)...?G(12)?G(1)?G(2)...?G(k\%12)$$
中間有$?\frac{k}{12}?$個(gè)$G(1)?G(2)?G(3)?...G(12)$

那么我們讓$Gsum=G(1)?G(2)?G(3)...?G(12)$

然后十分顯然
$$Ans=Gsu{m}^{?\frac{k}{12}?}?\prod _{i=1}^{k\%12}G(i)$$
然后用$Ans$矩陣計(jì)算答案即可。

---

$code$

#include<cstdio> #include<algorithm> #include<cstring> using namespace std; const int Size=60,T=12,XJQ=10000; struct matrix{int a[Size][Size]; }ans,f[T]; int n,m,s,e,k,Nfish; matrix operator *(matrix a,matrix b) {matrix c;memset(c.a,0,sizeof(c.a));for(int i=0;i<Size;i++)for(int j=0;j<Size;j++)for(int k=0;k<Size;k++)(c.a[i][j]+=a.a[i][k]*b.a[k][j]%XJQ)%=XJQ;return c; } void power(int b) {matrix F=f[0];for(int i=1;i<T;i++)F=F*f[i];while(b){if(b&1) ans=ans*F;F=F*F;b>>=1;} } int main() {scanf("%d%d%d%d%d",&n,&m,&s,&e,&k);for(int i=1;i<=m;i++){int x,y;scanf("%d%d",&x,&y);for(int i=0;i<T;i++)f[i].a[x][y]=f[i].a[y][x]=1;}scanf("%d",&Nfish);for(int i=1;i<=Nfish;i++){int op,x,last;scanf("%d",&op);for(int j=0;j<op;j++){scanf("%d",&x);for(int k=j;k<T;k+=op){int now=k-1;if(now<0) now=T-1; for(int i=0;i<n;i++)f[now].a[i][x]=0;}}}ans.a[0][s]=1;power(k/T);for(int i=0;i<k%T;i++)ans=ans*f[i];printf("%d",ans.a[0][e]); }

總結(jié)

以上是生活随笔為你收集整理的P2579,jzoj2288-[ZJOI2005]沼泽鳄鱼【矩阵乘法】的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。