なぜそれを並べ替えるTOPO？どのような判断スパイシーな複雑でなく、簡単にパンqwqに（私も料理でした）
直接tarjanのように、ああ！
我々が設定しなかった場合は、この質問は、しかし、DP式はいくつかのサイトで、良い非常に簡単ですそう（な空間制限の唯一の128メガバイトなど）カードとすることになっています
我々は、K値から開始することができますので、最大50があり、対象のK値が非常に小さいことがわかった、対象を観察しました。
私たちは、DP提供（i、j）は、私たちが来るのi番目の点を表し、最短「偏差値J」で、この時間は、プログラムの数です
DP面倒すぎるので、私は、検索を覚えていることを選択します
我々は、[i]はDISに、各点にダイクストラ最短1で最初に出て走っ
その後、我々は、検索を覚えておく必要があります。我々はポイント数に達したn個ならば、我々は今、プラスにお答えしますが、我々は最初の番に行ってNポイントして、n個のドットの波の後に走り回るために行ってきました、これも可能性があるため、直接ではなく返すように注意してくださいプログラムになります。
Vの一つは、我々は唯一の「ずれ」を計算する必要があることを命じた例譲渡、することができ、我々はxからVに到着した場合、その後、それはちょうど私達に依存して、この期間は、生産オフセット前回とは何の関係もないことがわかりましたX Vに、V 1から製造することができるどのくらい変化最短、すなわち：DIS [X] + W- DIS [V]
このように、単純な検索の波を実行することを忘れないでください、あなたは〜70ポイントを得ることができます
しかし、ガチョウ、我々は、ソリューションの無限の数があるかどうかを判断する必要があります。
どのようなオプションの無限の数を生成するのでしょうか？
リング0の実行可能な実施形態が存在する場合、明らかに、それは確実の解の対応する無限数を生成することができます
だから、我々は再びtarjan実行を描くからなるゼロの右側に端を開始することができます。
その後、0リングによって実現可能な解決策がある場合、我々は観察し、そうであれば、出力に-1。
コードは（一部起因してそのビット醜い、前に書かれていますが、ご容赦ください）：
#pragma GCC optimize(3,"inline","Ofast")//手动O3优化
#include<bits/stdc++.h>
using namespace std;
const int N=100001,P=5139222;
int n,m,k,p;
struct ver{
int u,v,w,next;
}t[N<<1];
priority_queue<pair<int,int> >s;
int las[N],resav[N];
int dis[N],d;
bool vis[N];
int tot;
inline void add(int u,int v,int w){
tot++;
t[tot].u=u,t[tot].v=v,t[tot].w=w;
t[tot].next=las[u],las[u]=tot;
return;
}
inline void dij(){
vis[1]=dis[1]=0;
s.push(make_pair(0,1));
while(!s.empty()){
int x=s.top().second,y=s.top().first;
s.pop();
if(vis[x]){
continue;
}
vis[x]=1;
for(int i=las[x];i;i=t[i].next){
int v=t[i].v,w=t[i].w;
if(dis[v]>dis[x]+w){
dis[v]=dis[x]+w;
s.push(make_pair(-dis[v],v));
}
}
}
}
int ans;int fa[N],sta[N],top,low[N],dfn[N];
int ji[N][51],cnt,tol[N];
bool vist[N][51],ip[N];
inline void Add(int x,int y){
q[++len]=(ver){x,y,0,pas[x]},pas[x]=len;
}
bool F;
inline int dfs(int x,int y){
if(y>k){
return 0;
}
if(vist[x][y]){//记忆化搜索
return ji[x][y];
}
vist[x][y]=1;
ji[x][y]=0;
if(x==n){
ji[x][y]++;
}
for(int i=las[x];i;i=t[i].next){
int nowp=dfs(t[i].v,y+(dis[x]+t[i].w-dis[t[i].v]));
if(nowp&&tol[fa[t[i].v]]>1){//判断是否有一种方案通过了零环
F=1;
return 0;
}
ji[x][y]+=nowp;
ji[x][y]%=p;
}
ji[x][y]%=p;
return ji[x][y];
}
inline void tarjan(int x){//tarjan
dfn[x]=low[x]=++cnt;sta[++top]=x;ip[x]=1;
for(int i=las[x];i;i=t[i].next){
if(t[i].w){
continue;
}
int v=t[i].v;
if(!dfn[v]){
tarjan(v);
low[x]=min(low[x],low[v]);
}
else if(ip[v]){
low[x]=min(low[x],dfn[v]);
}
}
if(low[x]==dfn[x]){
while(true){
int v=sta[top--];
ip[v]=0;fa[v]=x;++tol[x];
if(v==x){
break;
}
}
}
}
inline void check(){
cnt=0;
for(int i=1;i<=n;++i){
dfn[i]=low[i]=0,fa[i]=i,tol[i]=0;
}
for(int i=1;i<=n;++i){
if(!dfn[i]){
tarjan(i);
}
}
}
inline int read(){//fast read for int
int X=0,w=0; char ch=0;
while(!isdigit(ch)) {w|=ch=='-';ch=getchar();}
while(isdigit(ch)) X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
return w?-X:X;
}
int main(){
int T;
T=read();
bool flag=1;
while(T--){
ans=0;
n=read(),m=read(),k=read(),p=read();
las[1]=0;F=0;memset(pas,0,sizeof(pas)),len=0;
for(int i=2;i<=n;++i){
las[i]=0,dis[i]=P;
vis[i]=0;
}
tot=0;
for(int i=1;i<=m;++i){
int x,y,w;
x=read(),y=read(),w=read();
add(x,y,w);
}
check();
dij();
int rep=dfs(1,0)%p;
if(F){
puts("-1");
}else{
printf("%d\n",rep);
}
if(!T){
break;
}
memset(vist,0,sizeof(vist));
}
return 0;
}