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

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

[題解]bzoj1875(SDOI2009)HH去散步

2019-11-06 08:20:30
字體:
來源:轉載
供稿:網友

Description

HH有個一成不變的習慣,喜歡飯后百步走。所謂百步走,就是散步,就是在一定的時間 內,走過一定的距離。 但是同時HH又是個喜歡變化的人,所以他不會立刻沿著剛剛走來的路走回。 又因為HH是個喜歡變化的人,所以他每天走過的路徑都不完全一樣,他想知道他究竟有多 少種散步的方法。 現在給你學校的地圖(假設每條路的長度都是一樣的都是1),問長度為t,從給定地 點A走到給定地點B共有多少條符合條件的路徑

Input

第一行:五個整數N,M,t,A,B。其中N表示學校里的路口的個數,M表示學校里的 路的條數,t表示HH想要散步的距離,A表示散步的出發點,而B則表示散步的終點。 接下來M行,每行一組Ai,Bi,表示從路口Ai到路口Bi有一條路。數據保證Ai = Bi,但 不保證任意兩個路口之間至多只有一條路相連接。 路口編號從0到N ? 1。 同一行內所有數據均由一個空格隔開,行首行尾沒有多余空格。沒有多余空行。 答案模45989。

Output

一行,表示答案。

Sample Input

4 5 3 0 0 0 1 0 2 0 3 2 1 3 2

Sample Output

4

HINT

對于30%的數據,N ≤ 4,M ≤ 10,t ≤ 10。 對于100%的數據,N ≤ 20,M ≤ 60,t ≤ 2^30,0 ≤ A,B

Solution

這道題就是矩陣乘法,但是我們要把矩陣變一下。題目中要求不能立刻走回頭路,但是又有重邊,所以常規的鄰接矩陣是不行的。我們要用邊構矩陣。先把無向邊拆成兩條有向邊,然后得出矩陣A,其中Ai,j=1表示i號邊可以到達j號邊。然后求出At?1,此時A′i,j就表示從i號邊到j號邊走了t-1條邊的方案數。 然后我們枚舉每條起點的出邊i,對于每條出邊再枚舉每條終點的入邊j,∑A′i,j就是答案。

代碼:

#include<cstdio>#include<cstring>#include<algorithm>using namespace std;const int mod=45989,maxn=121;struct matrix{ int a[maxn][maxn]; matrix(){memset(a,0,sizeof a);}}temp,I;struct edge{ int same,to,next;}e[maxn];int n,m,t,A,B,num,head[maxn],in[maxn],ans;void add(int u,int v){ e[++num].to=v;e[num].next=head[u];head[u]=num; if(v==B)in[++in[0]]=num;}matrix cheng(matrix x,matrix y){ matrix re; for(int i=1;i<=num;i++){ for(int j=1;j<=num;j++){ for(int k=1;k<=num;k++){ (re.a[i][j]+=x.a[i][k]*y.a[k][j])%=mod; } } } return re;}matrix pow(matrix x,int y){ matrix re=I; while(y){ if(y&1)re=cheng(re,x); x=cheng(x,x);y>>=1; } return re;}int main(){ scanf("%d%d%d%d%d",&n,&m,&t,&A,&B); A++;B++; for(int i=1;i<=(m<<1);i++){ I.a[i][i]=1; } for(int i=1,u,v;i<=m;i++){ scanf("%d%d",&u,&v);u++;v++; add(u,v);e[num].same=num+1; add(v,u);e[num].same=num-1; } for(int i=1;i<=num;i++){ for(int j=head[e[i].to];j;j=e[j].next){ if(j!=e[i].same){ temp.a[i][j]=1; } } } temp=pow(temp,t-1); for(int i=head[A];i;i=e[i].next){ for(int j=1;j<=in[0];j++){ (ans+=(temp.a[i][in[j]]))%=mod; } }
發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
主站蜘蛛池模板: 青浦区| 公安县| 泉州市| 合水县| 治多县| 尼勒克县| 秦皇岛市| 邹平县| 都昌县| 喜德县| 台前县| 望城县| 南充市| 金溪县| 黑龙江省| 泰宁县| 连州市| 巴彦县| 房山区| 黎平县| 连州市| 道孚县| 东港市| 南雄市| 兴义市| 榆中县| 南和县| 佛冈县| 德惠市| 鄱阳县| 迁西县| 高雄市| 时尚| 肃北| 南充市| 阿勒泰市| 左云县| 萍乡市| 资阳市| 乌兰浩特市| 安国市|