国产探花免费观看_亚洲丰满少妇自慰呻吟_97日韩有码在线_资源在线日韩欧美_一区二区精品毛片,辰东完美世界有声小说,欢乐颂第一季,yy玄幻小说排行榜完本

首頁 > 學院 > 開發設計 > 正文

[BZOJ1898][Zjoi2005]Swamp 沼澤鱷魚(dp+矩陣優化)

2019-11-10 17:19:21
字體:
來源:轉載
供稿:網友

題目描述

傳送門

題解

f(i,j)表示時間為i時到達j的方案數,答案即為f(k,t) condi(i,j)表示時間為i時j有沒有鱷魚0/1 那么f(i,j)=∑(v,j)∈Ef(i?1,v) 并且當condi(i,j)=1時f(i,j)=0 時間復雜度O(kn2),考慮怎么優化 因為鱷魚的周期只可能是234,所以f的轉移是以12為周期循環的 那么我們可以把這12次轉移用矩陣處理出來,由于矩陣滿足結合律,可以先將它們乘起來然后快速冪k/12次,最后在乘上剩余的 矩陣的構造方法就是一個n*n的矩陣,如果f(i-1,j)可以轉移到f(i,v),那么在a[j,v]置1

ljw學姐的測試題,完完全全想出來了正解,然而被埋在了一個奇怪的地方——long long不能做位運算! GG

代碼

#include<algorithm>#include<iostream>#include<cstring>#include<cstdio>#include<cmath>using namespace std;#define N 55#define Mod 10000int n,m,s,t,k,x,y,nfish;int T[N],p[N][5];bool condi[N][N],flag[N][N];struct data{int a[N][N];}unit,st,A,trans[N],ans;data cheng(data a,data b){ data ans;memset(ans.a,0,sizeof(ans.a)); for (int k=1;k<=n;++k) for (int i=1;i<=n;++i) for (int j=1;j<=n;++j) ans.a[i][j]=(ans.a[i][j]+a.a[i][k]*b.a[k][j])%Mod; return ans;}data fast_pow(data a,int p){ data ans=unit; for (;p;p>>=1,a=cheng(a,a)) if (p&1) ans=cheng(ans,a); return ans;}int main(){ scanf("%d%d%d%d%d",&n,&m,&s,&t,&k);++s,++t; for (int i=1;i<=m;++i) { scanf("%d%d",&x,&y);++x,++y; flag[x][y]=flag[y][x]=1; } scanf("%d",&nfish); for (int i=1;i<=nfish;++i) { scanf("%d",&T[i]); for (int j=1;j<=T[i];++j) scanf("%d",&p[i][j]),++p[i][j]; } for (int i=1;i<=min(k,12);++i) for (int j=1;j<=nfish;++j) { int now=p[j][i%T[j]+1]; condi[i][now]=1; } for (int i=1;i<=n;++i) unit.a[i][i]=1; st.a[1][s]=1; for (int i=1;i<=min(k,12);++i) { for (int j=1;j<=n;++j) if (!condi[i][j]) for (int l=1;l<=n;++l) if (flag[l][j]) trans[i].a[l][j]=1; } A=unit; for (int i=1;i<=min(k,12);++i) A=cheng(A,trans[i]); if (k>12) { A=fast_pow(A,k/12); for (int i=1;i<=k%12;++i) A=cheng(A,trans[i]); } ans=cheng(st,A);
發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 屏边| 珠海市| 松桃| 田林县| 定结县| 常山县| 伊吾县| 井陉县| 嘉善县| 九江市| 大同县| 荣昌县| 金川县| 枣阳市| 洪江市| 东至县| 佳木斯市| 高台县| 醴陵市| 乌拉特前旗| 淅川县| 淮安市| 阳江市| 宁河县| 昆明市| 贺州市| 瑞昌市| 南昌县| 安福县| 盐边县| 克什克腾旗| 东乡| 扎赉特旗| 固阳县| 高州市| 曲阳县| 中宁县| 昆山市| 新余市| 东港市| 神农架林区|