bzoj1898: [Zjoi2005]Swamp 沼泽鳄鱼

一眼矩乘

把图分成12个，然后直接搞。

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const LL mod=;
 
int n;
struct Matrix
{
    LL a[][];
    Matrix(){}
    void clean()
    {
        memset(a,,sizeof(a));
    }
    void beone()
    {
        clean();
        for(int i=;i<=n;i++)a[i][i]=;
    }
    friend Matrix operator *(Matrix A,Matrix B)
    {
        Matrix C;C.clean();
        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.a[i][k]*B.a[k][j])%mod)%mod;
        return C;
    }
}G[];
 
LL mp[][];
int main()
{
    freopen("swamp.in","r",stdin);
    freopen("swamp.out","w",stdout);
    int m,st,ed,x,y;LL K;
    scanf("%d%d%d%d%lld",&n,&m,&st,&ed,&K);st++,ed++;
    memset(mp,,sizeof(mp));
    for(int i=;i<=m;i++)
    {
        scanf("%d%d",&x,&y);x++,y++;
        mp[x][y]=;mp[y][x]=;
    }
    for(int k=;k<=;k++)
        for(int i=;i<=n;i++)
            for(int j=;j<=n;j++)
                G[k].a[i][j]=mp[i][j];
 
    int fish;LL p;
    scanf("%d",&fish);
    for(int i=;i<=fish;i++)
    {
        scanf("%lld",&p);
        for(int j=;j<=p;j++)
        {
            scanf("%d",&x);x++;
            for(int k=j-;k<=;k+=p)
            {
                for(int u=;u<=n;u++)G[k].a[u][x]=;
            }
        }
    }
    Matrix YZH;YZH.beone();
    for(int i=;i<=;i++)
    {
        YZH=YZH*G[i%];
        if(K==i){printf("%lld\n",YZH.a[st][ed]);return ;}
    }
    p=K/;K%=;
    Matrix ans;ans.beone();
    while(p!=)
    {
        if(p%==)ans=ans*YZH;
        YZH=YZH*YZH;p/=;
    }
    for(int i=;i<=K;i++)ans=ans*G[i];
    printf("%lld\n",ans.a[st][ed]);
    return ;
}