国产探花免费观看_亚洲丰满少妇自慰呻吟_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);
發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 新乡市| 开阳县| 元江| 岑巩县| 临城县| 格尔木市| 枣庄市| 喀喇沁旗| 罗平县| 鄂托克前旗| 本溪| 长白| 铜梁县| 旺苍县| 新兴县| 南宁市| 临安市| 乐东| 鹤岗市| 江孜县| 太康县| 昌图县| 永善县| 陈巴尔虎旗| 定襄县| 博爱县| 呼伦贝尔市| 十堰市| 崇州市| 汶川县| 叶城县| 江城| 麻城市| 呼玛县| 武川县| 武安市| 临沭县| 柳州市| 临澧县| 志丹县| 北票市|