P2151-[SDOI2009]HH去散步【矩阵乘法】
生活随笔
收集整理的這篇文章主要介紹了
P2151-[SDOI2009]HH去散步【矩阵乘法】
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
正題
題目鏈接:https://www.luogu.com.cn/problem/P2151
題目大意
nnn個點mmm條邊的一張無向圖,走過一條邊后下一次就不能再走這條,求AAA到BBB之間有多少長度為ttt的路徑。
解題思路
因為有不能再走同一條的限制,所以不能直接裸上矩乘,換一個狀態(tài)考慮。設(shè)fif_ifi?表示上次走到第iii條邊的結(jié)尾處,那么這樣就可以對于每條邊特化轉(zhuǎn)移了。
時間復(fù)雜度O((2m)3log?t)O(\ (2m)^3\log t\ )O(?(2m)3logt?)
codecodecode
#include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; const ll N=121,XJQ=45989; ll n,m,t,A,B; struct node{ll x,y; }a[N*2]; struct matrix{ll a[N][N]; }f,ans,c; matrix operator*(matrix &a,matrix &b){memset(c.a,0,sizeof(c.a));for(ll i=0;i<m;i++)for(ll j=0;j<m;j++)for(ll k=0;k<m;k++)(c.a[i][j]+=a.a[i][k]*b.a[k][j]%XJQ)%=XJQ;return c; } int main() {scanf("%lld%lld%lld%lld%lld",&n,&m,&t,&A,&B);for(ll i=0;i<m;i++){scanf("%lld%lld",&a[i].x,&a[i].y);a[i+m].y=a[i].x;a[i+m].x=a[i].y;}m<<=1;for(ll i=0;i<m;i++)for(ll j=0;j<m;j++)if(a[i].y==a[j].x&&i%(m/2)!=j%(m/2))f.a[i][j]++;for(ll i=0;i<m;i++)if(a[i].x==A)ans.a[0][i]=1;t--;while(t){if(t&1)ans=ans*f;f=f*f;t>>=1;}ll answer=0;for(ll i=0;i<m;i++)if(a[i].y==B)(answer+=ans.a[0][i])%=XJQ;printf("%lld\n",answer); }總結(jié)
以上是生活随笔為你收集整理的P2151-[SDOI2009]HH去散步【矩阵乘法】的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 1953年属相 在1953出生的属相
- 下一篇: 葛优个人资料 演员葛优介绍