2016.7.13 刷题打卡

到22:29 才写完一道题：
bzoj 1898
主要是 搞了搞矩阵乘法在 信息学中的应用 没错就是 国家集训队论文
另外翻出来 lyd讲的 沼泽鳄鱼，写了些：

这是一道 dp题 但是 光看k就知道这个要 优化，怎么优化呢？

首先我们考虑 没有 鳄鱼的情况那问题就变成 了 从s到e路径为k 的 路径。
但是现在 我们有一群烦人的 小生物 鳄鱼。
没关系，不虚他！ 我们只要建立以下 Gi 即删掉 i时刻的不能经过的 点的相关的边。

那么此问题变成： π（i=1-k） Gi

那我如何快速的乘出这些东西呢？

这个时候题目里的重要信息 T=2，3,4 就要用上了！

这些的最小公倍数 是 12
所以 我们有 字丑勿怪。。。

那么上代码：

#include<cstdio>  #include<cstdlib>  #include<cstring>  // by mar_ch#define maxn 60  #define mod 10000  typedef int pp[maxn][maxn];  int n,m,start,end,k,Nfish,w[21][5],ans=0,len[maxn];  pp t,f[12],guy,all,a,apple;  //矩乘；   inline void mul(pp x,pp &y)  {      memset(t,0,sizeof(t));      for(int i=0;i<n;i++)       for(int j=0;j<n;j++)        for(int k=0;k<n;k++)         t[i][j]=(t[i][j]+x[i][k]*y[k][j]%mod)%mod;      memcpy(y,t,sizeof(y));  }  inline void fast(int y)  {       for(;y;y>>=1,mul(all,all))      if(y&1)mul(all,guy);  }  int main()  {      scanf("%d%d%d%d%d",&n,&m,&start,&end,&k);      for(int i=1,x,y;i<=m;i++){          scanf("%d%d",&x,&y);          a[x][y]=a[y][x]=1;      }      scanf("%d",&Nfish);      for(int i=1,t;i<=Nfish;i++){          scanf("%d",&len[i]);          for(int j=0;j<len[i];j++)          scanf("%d",&w[i][j]);      }      for(int i=0;i<n;i++)all[i][i]=1;      for(int i=0;i<12;i++){          memcpy(f[i],a,sizeof(f[i]));          for(int j=1;j<=Nfish;j++){              for(int k=0;k<n;k++)               f[i][w[j][(i+1)%len[j]]][k]=f[i][k][w[j][i%len[j]]]=0;          }          mul(f[i],all);      }      for(int i=0;i<n;i++)guy[i][i]=1;      fast(k/12);      for(int i=0;i<k%12;i++)mul(f[i],guy);      printf("%d\n",guy[end][start]);      return 0;  }  

嗯嗯。

然后 容我洗个澡，回来再大战 poj 3613

嗯哼，题解明天。。

#include<cstdio>  #include<cstring>  #include<climits>  #define min(a,b) ((a)<(b))?(a):(b)  const int N = 205;  const int INF = 1000000100;  int m,p[1001];  int ans[N][N];  void matrixmult(int a[N][N],int b[N][N]){      int i,j,k;      int c[N][N];      for(i=0;i<m;i++){          for(j=0;j<m;j++){              c[i][j]=INF;              for(k=0;k<m;k++)                  c[i][j]=min( a[i][k]+b[k][j] , c[i][j] );          }      }      for(i=0;i<m;i++){          for(j=0;j<m;j++)a[i][j]=c[i][j];      }  }  void mult(int s[N][N],int k){      int i,j;      for(i=0;i<m;i++){          for(j=0;j<m;j++)              ans[i][j]=s[i][j];      }      k--;      while(k>0){          if(k&1) matrixmult(ans,s);          k=k>>1;          matrixmult(s,s);      }  }  int main()  {      int mat[N][N];      int i,j,n,t,s,e,from,to,val;      while(~scanf("%d%d%d%d",&n,&t,&s,&e))      {          for(i=0;i<N;i++){              for(j=0;j<N;j++)mat[i][j]=INF;          }          memset(p,-1,sizeof(p));          m=0;          for(i=0;i<t;i++){              scanf("%d%d%d",&val,&from,&to);              if(p[from]<0) p[from] = m++;              if(p[to]<0) p[to] = m++;              if(val<mat[ p[from] ][ p[to] ])                  mat[ p[from] ][ p[to] ] = mat[ p[to] ][ p[from] ] = val;          }          s=p[s]; e=p[e];          mult(mat,n);          printf("%d\n",ans[s][e]);      }      return 0;  }