hdu 2157 How many ways??

来源：互联网发布：流程图软件推荐编辑：程序博客网时间：2024/05/16 06:52

这是一道经典的矩阵快速幂的题（可参考acm题集方向中的矩阵快速幂第8条），只需套一下模板，对应求一下数据就可以了，矩阵快速幂有点像disjkstra或者FLOYD 懂模板会用解题很方便（不太像万恶的背包或者说dp hh）

AC：

#include<iostream>  #include<cstdio>  #include<cstdlib>  #include<cstring>  #include<string>    #include<algorithm>      using namespace std;    const int MAX=20+10;  const int mod=1000;  int f[MAX][MAX],sum[MAX][MAX],array[MAX][MAX];  int n,m;    void MatrixMult(int a[MAX][MAX],int b[MAX][MAX]){//矩阵相乘     int c[MAX][MAX]={0};      for(int i=0;i<n;++i){          for(int j=0;j<n;++j){              for(int k=0;k<n;++k){                  c[i][j]=(c[i][j]+a[i][k]*b[k][j])%mod;              }          }      }      for(int i=0;i<n;++i){          for(int j=0;j<n;++j)a[i][j]=c[i][j];      }  }    int Matrix(int a,int b,int k){//矩阵快速幂       for(int i=0;i<n;++i){//初始化矩阵sum,使sum*b=b(b是矩阵)           for(int j=0;j<n;++j){              if(i == j)sum[i][j]=1;              else sum[i][j]=0;          }      }      for(int i=0;i<n;++i){//将最初矩阵f复制给array           for(int j=0;j<n;++j)array[i][j]=f[i][j];      }      while(k){          if(k&1)MatrixMult(sum,array);          MatrixMult(array,array);          k>>=1;      }      return sum[a][b];  }int main(){int s,t,T,a,b,k;while(scanf("%d %d",&n,&m)==2&&(n||m)){memset(f,0,sizeof(f));for(int i=0;i<m;i++){scanf("%d %d",&s,&t);f[s][t]=1;}scanf("%d",&T);while(T--){scanf("%d %d %d",&a,&b,&k);cout<<Matrix(a,b,k)%1000<<endl;}}return 0;}

其中非main函数部分是模板可以收集一下嘿嘿

0 0