bzoj1875 HH去散步

题目链接

一开始特别天真，

$Folyd$传递闭包写了一下过了样例就交上去了

然后$gg$

$qwq$

想了想$……$

还要写一点东西

最后统计答案的时候有细节要注意

/**************************************************************
    Problem: 1875
    User: zhangheran
    Language: C++
    Result: Accepted
    Time:17164 ms
    Memory:3876 kb
****************************************************************/

// luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
//#include"suqingnian.h"
#include<algorithm>
#include<cstring>
using namespace std;
const int mod=;
int n,m,t,a,b;
struct data{
    int v;int u;int next;
}edge[];int alist[];int cnt;
void add(int u,int v)
{
    edge[++cnt].u=u;
    edge[cnt].v=v;
    edge[cnt].next=alist[u];
    alist[u]=cnt;
    return ;
}
struct Martix{
    long long num[][];
    void hint(){memset(num,,sizeof(num));for(int i=;i<=m*;i++) num[i][i]=;}
    void clear(){memset(num,,sizeof(num));}
    friend Martix operator *(const Martix &a,const Martix &b)
    {
        Martix c;
        memset(c.num,,sizeof(c.num));
        for(int k=;k<=m*;k++)
          for(int i=;i<=m*;i++)
            for(int j=;j<=m*;j++)
              c.num[i][j]=(c.num[i][j]+a.num[i][k]*b.num[k][j]%mod)%mod;
        return c;
    }
    friend Martix operator +(const Martix &a,const Martix &b)
    {
        Martix c;memset(c.num,,sizeof(c.num));
        for(int i=;i<=m*;i++)
          for(int j=;j<=m*;j++)
            c.num[i][j]=(a.num[i][j]+b.num[i][j])%mod;
        return c;
    }
};
Martix
    _pow(Martix _a,int _b)
    {
        Martix _res;_res.hint();
        for(;_b;_b >>= , _a= _a * _a )
          if(_b & ) _res = _res * _a ;
        return _res ;
    }
int u,v;
Martix map,dis;
int ans=;
int main()
{
    memset(alist,-,sizeof(alist));
    scanf("%d%d%d%d%d",&n,&m,&t,&a,&b);
    a++,b++;
    for(int i=;i<=m;i++) scanf("%d%d",&u,&v),add(u+,v+),add(v+,u+);
    int next=alist[a];
    while(next!=-) map.num[][next]++,next=edge[next].next;
    for(int i=;i<=cnt;i++){
        v=edge[i].v,next=alist[v];
        while(next!=-){
            if(next!=(i+((i&)?:-)))
              dis.num[i][next]++;
            next=edge[next].next;
        }
    }
    dis=_pow(dis,t-);
    map=map*dis;
//  for(int i=0;i<=m*2;i++,puts(""))
//    for(int j=0;j<=m*2;j++)
//      printf("%d ",map.num[i][j]);
    for(int i=;i<=cnt;i++) if(edge[i].v==b) ans=(ans+map.num[][i])%mod;
    printf("%d",ans);
    return ;
}