ZJOI2005沼泽鳄鱼

Posted artlover

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了ZJOI2005沼泽鳄鱼相关的知识,希望对你有一定的参考价值。

矩阵优化dp

** 注意:矩阵乘法没有交换律 **

思路:类比P2151hh去散步

代码特点在一维的答案矩阵

1.矩阵优化两点间方案数不必赘述

2.注意2,3,4可以办到以他们的lcm为周期,正是因为如此我们可以矩阵加速(这样我们就可以化动为静,矩阵乘法了)

3.一维初始矩阵(一维邻接矩阵+第二个矩阵取交集)注意当前鳄鱼的位置与我们下一次走并无关,他们正好搓了一位,要小心

4.再次强调:矩阵乘法没有交换律:因此,我们在构造转移矩阵时要从2开始到12,再乘1(这里指的是矩阵)

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,m,S,E,k,t,mod=10000;
int record[20][60],Map[60][60];
struct M{
  int m[60][60];
  M(){
    memset(m,0,sizeof(m));
  }
  M operator *(const M &a)const {//重载定义矩阵乘法
    M ret;
    for(int i=1;i<=n;++i)
      for(int j=1;j<=n;++j)
        for(int k=1;k<=n;++k)
          ret.m[i][j]=(ret.m[i][j]+m[i][k]%mod*a.m[k][j]%mod)%mod;
    return ret;
  }
}c[20],e,A;
inline void build(){//构造转移矩阵
  for(int i=1;i<=n;++i)e.m[i][i]=1;
   for(int i=1;i<=12;++i){
      for(int j=1;j<=n;++j){
          for(int k=1;k<=n;++k){           
            if(Map[j][k]&&(!record[i][j]))c[i].m[k][j]=1;
          }
      }
   }
   for(int i=2;i<=12;i++)e=e*c[i];//矩乘没有交换律
   e=e*c[1];//因为我们的第二个矩阵跟第一步有关,所以先乘2矩阵
}
inline void power(M &a,M &b,int k){
  while(k){
    if(k&1)a=a*b;
    b=b*b;
    k>>=1;
  }
}//矩阵快速幂
inline void solve(){//极其繁琐需要头脑冷静
  if(k==1){printf("%d",Map[S][E]);return ;}
  for(int i=1;i<=n;++i)if(Map[S][i]&&c[2].m[S][i])A.m[S][i]=1;//先跟第二矩阵取交集,表示走第一步--此处表示构造初始矩阵
  if(k<=12)for(int i=3;i<=k+1;++i)A=A*c[i%12];//此处分类讨论
  else {//因为初始矩阵的存在,我们已经走了一步
    for(int i=3;i<=12;++i)A=A*c[i];//所以在这里我们先暴力走一遍
    A=A*c[1];//注意不能交换律
    power(A,e,(k/12)-1);//先乘了一遍所以-1
    for(int i=2;i<=(k-12)%12+1;++i)A=A*c[i];//最后再o再
  }
  printf("%d",A.m[S][E]);
}

int main(){
   scanf("%d%d%d%d%d",&n,&m,&S,&E,&k);
   S++,E++;//加1比较顺我的思路,下面的加1同
   for(int i=1;i<=m;++i){
    int x,y;
    scanf("%d%d",&x,&y);
    ++x,++y;
    Map[x][y]=1;
    Map[y][x]=1;
   }
   scanf("%d",&t);
   for(int i=1;i<=t;++i){
     int num;
     scanf("%d",&num);
     for(int j=1;j<=num;++j){
      int x;
      scanf("%d",&x);
      int k=j;
      x++;
      while(k<=12){
        record[k][x]=1;
            k+=num;
      }//12一循环
     }
   }
   build();
   solve();
   return 0;
}

以上是关于ZJOI2005沼泽鳄鱼的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ 1898: [Zjoi2005]Swamp 沼泽鳄鱼

BZOJ 1898: [Zjoi2005]Swamp 沼泽鳄鱼 [矩阵乘法]

BZOJ1898[Zjoi2005]Swamp 沼泽鳄鱼 矩阵乘法

bzoj1898: [Zjoi2005]Swamp 沼泽鳄鱼

bzoj1898 [Zjoi2005]沼泽鳄鱼

P2579 [ZJOI2005]沼泽鳄鱼