BZOJ1875 HH去散步
这题 我承认 我过于NAIVE了 我把不能走回边 搞成了 走到一个点必须休息一下。。。
发现它的边数非常少 所以 按边建立矩阵 就可以十分科学的完成本题
矩乘我一定要封装进p,q
BZOJ 1875
1 #include<bits/stdc++.h> 2 #define mod 45989 3 #define ll long long 4 #define breaks printf("!\n"); 5 using namespace std; 6 7 inline int read() 8 { 9 int x=0,f=1;char ch=getchar(); 10 while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} 11 while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} 12 return x*f; 13 } 14 struct edge 15 { 16 int u,v,next; 17 }vs[300]; 18 int st[80],ee=1; 19 int n,m,t,A,B; 20 struct matrix 21 { 22 int mat[130][130]; 23 void clr() 24 { 25 for(int i=1;i<=ee;i++) 26 for(int j=1;j<=ee;j++) 27 mat[i][j]=0; 28 } 29 }M; 30 matrix operator*(matrix a,matrix b) 31 { 32 matrix tmp; tmp.clr(); 33 for(int i=1;i<=ee;i++) 34 for(int j=1;j<=ee;j++) 35 for(int k=1;k<=ee;k++) 36 (tmp.mat[i][j]+=1ll*a.mat[i][k]*b.mat[k][j]%mod)%=mod; 37 return tmp; 38 } 39 matrix operator^(matrix x,int p) 40 { 41 matrix tmp;tmp.clr(); 42 for(int i=1;i<=ee;i++) 43 tmp.mat[i][i]=1; 44 for(;p;p>>=1,x=x*x) 45 if(p&1)tmp=tmp*x; 46 return tmp; 47 } 48 inline void addedge(int u,int v) 49 { 50 vs[++ee].v=v;vs[ee].u=u; 51 vs[ee].next=st[u];st[u]=ee; 52 } 53 int main() 54 { 55 // freopen("a.out","w",stdout); 56 n=read();m=read();t=read(); 57 A=read();B=read(); 58 for(int i=1;i<=m;i++) 59 { 60 int a=read(),b=read(); 61 addedge(a,b);addedge(b,a); 62 } 63 matrix xx; xx.clr(); 64 for(int i=st[A];i;i=vs[i].next) 65 xx.mat[1][i]++; 66 for(int i=2;i<=ee;i++) 67 for(int j=2;j<=ee;j++) 68 if(vs[i].v==vs[j].u) 69 if(i!=(j^1))M.mat[i][j]++; 70 // xx.print(); M.print(); 71 xx=xx*(M^(t-1)); int ans=0; 72 for(int i=st[B];i;i=vs[i].next) 73 (ans+=xx.mat[1][i^1])%=mod; 74 cout<<ans<<endl; 75 return 0; 76 }
