spfa预处理出最短路数组dis，然后反向建边bfs出ok[u]表示u能到n点
然后发现有0环的话时候有inf解的，先dfs找0环判断即可
然后dfs，设状态f[u][v]为到u点，还可以跑最短路+v的长度，记忆化，然后直接搜即可
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<vector>
using namespace std;
const int N=100005,inf=1e9;
int T,n,m,k,mod,h[N],cnt,ans[N][55],dis[N],len,v[N],ti,fl;
bool f[N][55],ok[N];
vector<int>a[N];
struct qwe
{
int ne,to,va;
}e[N<<1];
int read()
{
int r=0,f=1;
char p=getchar();
while(p>'9'||p<'0')
{
if(p=='-')
f=-1;
p=getchar();
}
while(p>='0'&&p<='9')
{
r=r*10+p-48;
p=getchar();
}
return r*f;
}
inline void jia(int &x,int y)
{
x+=y;
if(x>=mod)
x-=mod;
}
void add(int u,int v,int w)
{
cnt++;
e[cnt].ne=h[u];
e[cnt].to=v;
e[cnt].va=w;
h[u]=cnt;
}
void zer(int u)
{
v[u]=ti;
for(int i=h[u];i;i=e[i].ne)
if(e[i].va==0)
{
if(v[e[i].to]==ti)
{
fl=1;
return;
}
else if(!v[e[i].to])
zer(e[i].to);
}
}
int dfs(int u,int v)
{
if(v<0)
return 0;
if(f[u][v])
return ans[u][v];
int sm=(u==n);
for(int i=h[u];i;i=e[i].ne)
if(ok[e[i].to])
sm=(sm+dfs(e[i].to,v-(e[i].va-(dis[e[i].to]-dis[u]))))%mod;
f[u][v]=1,ans[u][v]=sm;
return sm;
}
int main()
{
T=read();
while(T--)
{
memset(h,0,sizeof(h));
memset(ok,0,sizeof(ok));
memset(f,0,sizeof(f));
memset(ans,-1,sizeof(ans));
cnt=0;ti=0,fl=0;
n=read(),m=read(),k=read(),mod=read();
for(int i=1;i<=n;i++)
a[i].clear();
for(int i=1;i<=m;i++)
{
int x=read(),y=read(),z=read();
add(x,y,z);
a[y].push_back(x);
}
for(int i=1;i<=n&&!fl;i++)
if(!v[i])
ti++,zer(i);
if(fl)
{
puts("-1");
continue;
}
queue<int>q;
ok[n]=1,q.push(n);
while(!q.empty())
{
int u=q.front();
q.pop();
for(int i=0;i<a[u].size();i++)
if(!ok[a[u][i]])
{
ok[a[u][i]]=1;
q.push(a[u][i]);
}
}
for(int i=1;i<=n;i++)
dis[i]=inf,v[i]=0;
dis[1]=0,v[1]=1;
q.push(1);
while(!q.empty())
{
int u=q.front();
q.pop();
v[u]=0;
for(int i=h[u];i;i=e[i].ne)
if(dis[e[i].to]>dis[u]+e[i].va)
{
dis[e[i].to]=dis[u]+e[i].va;
if(!v[e[i].to])
{
v[e[i].to]=1;
q.push(e[i].to);
}
}
}
// for(int i=1;i<=n;i++)
// cerr<<ok[i]<<endl;
printf("%d\n",dfs(1,k));
}
return 0;
}