P2579 [ZJOI2005]沼泽鳄鱼（矩乘）

P2579 [ZJOI2005]沼泽鳄鱼

没有食人鱼：直接矩乘优化

有食人鱼：食人鱼周期2、3、4，公倍数12，可以以12为一个周期矩乘，剩下的暴力

注意矩乘不满足乘法交换律，一定要注意乘的顺序

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int n,m,S,T,K,d[];
struct mat{
    int a[][];
    mat(){memset(a,,sizeof(a));}
    mat operator * (const mat &G) const{
        mat c;
        for(int i=;i<n;++i)
            for(int j=;j<n;++j)
                for(int k=;k<n;++k)
                    c.a[i][j]=(c.a[i][j]+a[i][k]*G.a[k][j]%)%;
        return c;
    }
}s,f[],w;
mat Pow(mat x,int y){
    mat re;
    for(int i=;i<n;++i) re.a[i][i]=;
    for(;y;y>>=,x=x*x) if(y&) re=re*x;
    return re;
}
int main(){
    scanf("%d%d%d%d%d",&n,&m,&S,&T,&K);
    for(int i=,u,v;i<=m;++i){
        scanf("%d%d",&u,&v);
        s.a[u][v]=s.a[v][u]=;
    }
    for(int i=;i<;++i) f[i]=s;
    scanf("%d",&m);
    for(int i=,t;i<=m;++i){
        scanf("%d",&t);
        for(int u=;u<t;++u) scanf("%d",&d[u]);
        for(int j=;j<;++j)
            for(int k=;k<n;++k)
                f[j].a[k][d[j%t]]=;
    }
    for(int i=;i<n;++i) w.a[i][i]=;
    for(int i=;i<;++i) w=w*f[i];
    w=w*f[]; w=Pow(w,K/);
    for(int i=;i<=K%;++i) w=w*f[i];
    printf("%d",w.a[S][T]);
    return ;
}