日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問(wèn) 生活随笔!

生活随笔

當(dāng)前位置: 首頁(yè) > 编程资源 > 编程问答 >内容正文

编程问答

bzoj1297 [SCOI2009]迷路——拆点+矩阵快速幂

發(fā)布時(shí)間:2024/9/5 编程问答 26 豆豆
生活随笔 收集整理的這篇文章主要介紹了 bzoj1297 [SCOI2009]迷路——拆点+矩阵快速幂 小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.

題目:https://www.lydsy.com/JudgeOnline/problem.php?id=1297

一看感覺是矩陣快速冪之類的,但邊權(quán)不好處理啊;

普通的矩陣快速冪只能處理邊權(quán)為1的,所以想辦法把邊權(quán)處理成1;

仔細(xì)一看還有一個(gè)條件是邊權(quán)小于10;

所以拆點(diǎn)!把一個(gè)點(diǎn)拆成10個(gè)點(diǎn)表示到它不同的距離,那么和它相連的那些點(diǎn)就可以跟某個(gè)距離的點(diǎn)連邊權(quán)為1的邊;

雖然沒有自己想出來(lái),不過(guò)1A還是極好的!(因?yàn)樘?jiǎn)單了)

代碼如下:

#include<iostream> #include<cstdio> #include<cstring> using namespace std; int n,T,w[15][15],tot,id[15][15],mod=2009; struct Matrix{int a[105][105];Matrix(){memset(a,0,sizeof a);}Matrix operator * (const Matrix &y) const{Matrix x;for(int i=1;i<=tot;i++)for(int k=1;k<=tot;k++)for(int j=1;j<=tot;j++)(x.a[i][j]+=a[i][k]*y.a[k][j])%=mod;return x;}void init(){for(int i=1;i<=tot;i++)a[i][i]=1;} }f,ans; Matrix pw(Matrix x,int k) {Matrix ret; ret.init();for(;k;k>>=1,x=x*x)if(k&1)ret=ret*x;return ret; } int main() {scanf("%d%d",&n,&T);char ch[15];for(int i=0;i<n;i++){scanf("%s",&ch);for(int j=0;j<n;j++)w[i][j]=ch[j]-'0';}for(int i=0;i<n;i++)for(int j=0;j<=9;j++){id[i][j]=++tot;if(j)f.a[tot-1][tot]++;}for(int i=0;i<n;i++)for(int j=0;j<n;j++)if(w[i][j])f.a[id[i][w[i][j]-1]][id[j][0]]++;ans.a[1][id[0][0]]=1;Matrix fn=pw(f,T);ans=ans*fn;printf("%d",ans.a[1][id[n-1][0]]);return 0; }

?

轉(zhuǎn)載于:https://www.cnblogs.com/Zinn/p/9254628.html

總結(jié)

以上是生活随笔為你收集整理的bzoj1297 [SCOI2009]迷路——拆点+矩阵快速幂的全部?jī)?nèi)容,希望文章能夠幫你解決所遇到的問(wèn)題。

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