bzoj 1898 [Zjoi2005]Swamp 沼泽鳄鱼——矩阵快速幂

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=1898

当然是邻接矩阵做转移矩阵来快速幂。

对于鳄鱼，好在它们周期的lcm是12，也就是每12次就又一样了。

所以把12个转移矩阵合成一下，就可以每次乘一样的，进而快速幂。%12剩下的次数暴力一下。

学到了一些方便的东西，比如struct的构造函数没有参数之类的。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=,mod=1e4;
int n,m,st,ed,k,fs;
struct Matrix{
    int a[N][N];
    Matrix(){memset(a,,sizeof a);}
    Matrix operator*(const Matrix &b)const
    {
        Matrix c;
        memset(c.a,,sizeof c.a);//
        for(int i=;i<=n;i++)
            for(int k=;k<=n;k++)
                for(int j=;j<=n;j++)
                    (c.a[i][j]+=a[i][k]*b.a[k][j]%mod)%=mod;
        return c;
    }
    void init()
    {
        for(int i=;i<=n;i++)a[i][i]=;
    }
}r[],ans;
int main()
{
    scanf("%d%d%d%d%d",&n,&m,&st,&ed,&k);st++;ed++;
    int x,y;
    for(int i=;i<=m;i++)
    {
        scanf("%d%d",&x,&y);x++;y++;
        for(int j=;j<=;j++)
            r[j].a[x][y]=,r[j].a[y][x]=;
    }
    scanf("%d",&fs);
    for(int i=;i<=fs;i++)
    {
        scanf("%d",&x);
        for(int j=;j<=x;j++)
        {
            scanf("%d",&y);y++;
            for(int k=j;k<=;k+=x)
                memset(r[k].a[y],,sizeof r[k].a[y]);
        }
    }
    r[].init();
    for(int i=;i<=;i++)r[]=r[]*r[i];
    ans.a[][st]=;x=k/;
    while(x)
    {
        if(x&)ans=ans*r[];
        r[]=r[]*r[];x>>=;
    }
    x=k%;
    for(int i=;i<=x;i++)ans=ans*r[i];
    printf("%d",ans.a[][ed]);
    return ;
}