P2579 [ZJOI2005]沼泽鳄鱼（邻接矩阵，快速幂）

题目简洁明了（一点都不好伐）

照例，化简题目

给一张图，每一个时间点有一些点不能走，（有周期性），求从起点第k秒恰好在终点的方案数，可重复，不可停留。

额dp实锤

于是就被打脸了....

有一种东西叫做邻接矩阵，还有一种东西叫做矩阵乘法，bk201大仙曾经讲过，能用邻接矩阵的k次方求这个东西。

那，难度下降了很多了。

但是，对于那周期是3,4,6的鳄鱼怎么办呢？

答案就是：

观察3,4,6，lcm（3,4,6）=12，所以，当k%12==0，鳄鱼们都回到了原状态。

所以，我们就暴力地对于每一种情况建一个邻接矩阵，把不能走的点刨掉，然后把他们依次乘起来（因为不满足交换律）

然后把它给k次方就行了。。。

值得注意的地方：

1、因为点是从0开始，所以根据我的习惯，都加上了1；

2、不满足交换律不满足交换律不满足交换律

3、k不一定整除12，所以在剩余的那一部分要手动乘一下初始矩阵

代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=55;
const int mod=10000;
int m,n,st,en,K,fn;
int map[maxn][maxn];
int unable[maxn][maxn];
struct node
{
    int a[maxn][maxn];
    node()
    {
        memset(a,0,sizeof(a));
    }
}q,w[55];
inline node operator * (const node& x,const node& y)
{
    node z;
    for(int i=1;i<=n;i++)
    {
        for(int k=1;k<=n;k++)
        {
            for(int j=1;j<=n;j++)
            {
                z.a[i][j]=(z.a[i][j]+x.a[i][k]*y.a[k][j])%mod;
            }
        }
    }
    return z;
}
node ksm(node x,int y)
{
    node z;
    for(int i=1;i<=n;i++)
    z.a[i][i]=1;
    while(y)
    {
        if(y&1)
        z=z*x;
        x=x*x;
        y>>=1;
    }
    return z;
}

int main()
{
    scanf("%d%d%d%d%d",&n,&m,&st,&en,&K);
    st++;en++;
    for(int i=1;i<=m;i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        x++;y++;
        map[x][y]=map[y][x]=1;
    }
    scanf("%d",&fn);
    while(fn--)
    {
        int x;
        scanf("%d",&x);
        for(int j=0;j<x;j++)
        {
            int y;
            scanf("%d",&y);
            y++;
            for(int k=j;k<12;k+=x)
            {
                unable[k][y]=true;
            }
        }
    }

    for(int i=1;i<=n;i++)
    {
        q.a[i][i]=1;
    }
    for(int i=0;i<12;i++)
    {
        for(int j=1;j<=n;j++)
        {
            for(int k=1;k<=n;k++)
            {
                w[i].a[j][k]=(map[j][k]&&!unable[i][k]);
            }
        }
    }
    for(int i=1;i<12;i++)
    {
        q=q*w[i];
    }
    q=q*w[0];
    int k=K/12,t=K-k*12;
    node ans;
    for(int i=1;i<=n;i++)
    {
        ans.a[i][i]=1;
    }
    ans=ans*ksm(q,k);
    for(int i=1;i<=t;i++)
    ans=ans*w[i];
    printf("%d ",ans.a[st][en]);
    return 0;
}

（完）