P3758 [TJOI2017]可乐

P3758 [TJOI2017]可乐

---

这个题是利用了邻接矩阵的性质

即：对于一个邻接矩阵\(E\)，表示从\(E[i][j]\)到\(i\)到\(j\)路径长度为1的方案数是多少。那么\(E^k[i][j]\)表示从\(i\)到\(j\)路径长度为\(K\)

对这个题

对于走到其他城市，就可以直接使用邻接矩阵想乘

而对于原地不动，我们可以理解成向当前的城市添加一个自环就行了

然后对于自爆，我们就是前一秒的方案数都存下来，我们就可以设一个零号城市，所有城市都可以到他，然后这个城市无法到达其他城市

然后就可以愉快的快速幂了

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
const int maxn=110;
const int RQY=2017;
struct node
{
    int n,m;
    int map[maxn][maxn];
    node (int n=0,int m=0)
    {
        memset(map,0,sizeof(map));
    }
    node operator * (const node &a)const 
    {
        node res;
        res.n=n,res.m=a.m;
        for(int i=0;i<=n;i++)
            for(int j=0;j<=res.m;j++)
                for(int k=0;k<=m;k++)
                    res.map[i][j]=(res.map[i][j]+map[i][k]*a.map[k][j])%RQY;
        return res;
    }
};
node A,C;
node kasumi(node b,int t)
{
    node Res=b;t--;
    while(t)
    {
        if(t&1)
            Res=Res*b;
        t>>=1;
        b=b*b;
    }
    return Res;
}
int main()
{
    int n,m,T;
    scanf("%d%d",&n,&m);
    A.n=n;A.m=n;
    for(int i=0;i<=n;i++)   A.map[i][i]=A.map[i][0]=1;
    int a,b;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&a,&b);
        A.map[a][b]++;
        A.map[b][a]++;
    }
    scanf("%d",&T);
    C=kasumi(A,T);
    int tot=0;
    for(int i=0;i<=C.m;i++)
        tot=(tot+C.map[1][i])%RQY;
    printf("%d",tot);
}