HH有個一成不變的習慣,喜歡飯后百步走。所謂百步走,就是散步,就是在一定的時間 內,走過一定的距離。 但是同時HH又是個喜歡變化的人,所以他不會立刻沿著剛剛走來的路走回。 又因為HH是個喜歡變化的人,所以他每天走過的路徑都不完全一樣,他想知道他究竟有多 少種散步的方法。 現在給你學校的地圖(假設每條路的長度都是一樣的都是1),問長度為t,從給定地 點A走到給定地點B共有多少條符合條件的路徑
第一行:五個整數N,M,t,A,B。其中N表示學校里的路口的個數,M表示學校里的 路的條數,t表示HH想要散步的距離,A表示散步的出發點,而B則表示散步的終點。 接下來M行,每行一組Ai,Bi,表示從路口Ai到路口Bi有一條路。數據保證Ai = Bi,但 不保證任意兩個路口之間至多只有一條路相連接。 路口編號從0到N ? 1。 同一行內所有數據均由一個空格隔開,行首行尾沒有多余空格。沒有多余空行。 答案模45989。
一行,表示答案。
4 5 3 0 0 0 1 0 2 0 3 2 1 3 2
4
對于30%的數據,N ≤ 4,M ≤ 10,t ≤ 10。 對于100%的數據,N ≤ 20,M ≤ 60,t ≤ 2^30,0 ≤ A,B
這道題就是矩陣乘法,但是我們要把矩陣變一下。題目中要求不能立刻走回頭路,但是又有重邊,所以常規的鄰接矩陣是不行的。我們要用邊構矩陣。先把無向邊拆成兩條有向邊,然后得出矩陣A,其中
代碼:
#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; } }新聞熱點
疑難解答