NOIP2017 Day1 T3 逛公园（最短路+拓扑排序+DP）

　　神tm比赛时多清个零就有60了T T

　　首先跑出1起点和n起点的最短路，因为k只有50，所以可以DP。设f[i][j]表示比最短路多走i的长度，到j的方案数。

　　我们发现如果在最短路上的和零边会有后向性，怎么办呢？拓扑排序。

　　把最短路上的点和零边的点拉出来跑拓扑排序，如果有零环的话必定度数不为0，而且要注意零环必须在<=最短路+k的路径上才输出-1，这个就用刚刚跑出来的1起点到n起点的最短路来判断就好了。

　　然后先按拓扑序DP出i相同的，然后再DP不在最短路上或者零边的。

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<queue>
#define MOD(x) ((x)>=mod?(x)-mod:(x))
using namespace std;
const int maxn=;
struct tjm{int too, dis, pre;}e[][maxn];
struct poi{int x, dis;};
priority_queue<poi>q;
bool operator<(poi a, poi b){return a.dis>b.dis;}
int T, x, y, z, n, m, K, mod, top, tot[], ans;
int f[][maxn], dist[][maxn], last[][maxn], d[maxn], st[maxn];
void read(int &k)
{
    int f=; k=; char c=getchar();
    while(c<'' || c>'') c=='-' && (f=-), c=getchar();
    while(c<='' && c>='') k=k*+c-'', c=getchar();
    k*=f;
}
inline void add(int x, int y, int z, int ty){e[ty][++tot[ty]]=(tjm){y, z, last[ty][x]}; last[ty][x]=tot[ty];}
inline void dij(int x, int ty)
{
    memset(dist[ty], , sizeof(dist[ty]));
    dist[ty][x]=; q.push((poi){x, });
    while(!q.empty())
    {
        poi now=q.top(); q.pop();
        if(dist[ty][now.x]!=now.dis) continue;
        for(int i=last[ty][now.x], too;i;i=e[ty][i].pre)
        if(dist[ty][too=e[ty][i].too]>dist[ty][now.x]+e[ty][i].dis)
        {
            dist[ty][too]=dist[ty][now.x]+e[ty][i].dis;
            q.push((poi){too, dist[ty][too]});
        }
    }
}
inline bool topo()
{
    memset(d, , sizeof(d));
    for(int i=;i<=n;i++)
    for(int j=last[][i], too;j;j=e[][j].pre)
    if(dist[][i]+e[][j].dis==dist[][too=e[][j].too]) d[too]++;
    top=; for(int i=;i<=n;i++) if(!d[i]) st[++top]=i;
    for(int i=;i<=top;i++)
    for(int j=last[][st[i]], too;j;j=e[][j].pre)
    if(dist[][st[i]]+e[][j].dis==dist[][too=e[][j].too])
    {
        d[too]--;
        if(!d[too]) st[++top]=too;
    }
    for(int i=;i<=n;i++) if(d[i] && dist[][i]+dist[i][n]<=dist[][n]+K) return ;
    return ;
}
int main()
{
    read(T);
    while(T--)
    {
        memset(last, , sizeof(last)); tot[]=tot[]=;
        read(n); read(m); read(K); read(mod);
        for(int i=;i<=m;i++) read(x), read(y), read(z), add(x, y, z, ), add(y, x, z, );
        dij(, ); dij(n, );
        if(!topo()) {puts("-1"); continue;}
        memset(f, , sizeof(f)); f[][]=; ans=;
        for(int i=;i<=K;i++)
        {
            for(int j=;j<=top;j++)
            for(int k=last[][st[j]], too;k;k=e[][k].pre)
            if(e[][k].dis+dist[][st[j]]==dist[][too=e[][k].too])
            f[i][too]+=f[i][st[j]], f[i][too]=MOD(f[i][too]);
            for(int j=;j<=n;j++)
            for(int k=last[][j], too, tmp;k;k=e[][k].pre)
            if((tmp=i+e[][k].dis+dist[][j]-dist[][too=e[][k].too])<=K && i!=tmp)
            f[tmp][too]+=f[i][j], f[tmp][too]=MOD(f[tmp][too]);
            ans+=f[i][n]; ans=MOD(ans);
        }
        printf("%d\n", ans);
    }
}