我在有0环的图里跑了最短路计数
我可能已经是个废物啦
很早之前就想写这道题啦，但是太菜了发现自己不会，今天终于写啦
首先我们建图的时候建出一个正图还有一个反图，我们对着这两个图分别跑两次最短路，求出\(1\)到所有点的最短路，以及所有点到\(n\)的最短路
如果不考虑无解的情况，我们现在就可以大力记搜了
我们设\(dp[x][j]\)表示从\(x\)走到\(n\)相比\(x\)到\(n\)的最短路多走了\(j\)的路径条数
之后大力记搜即可，\(ans=\sum_{i=0}^Kdp[1][i]\)
我们去记忆化搜索就行了，比如说从\(dp[i][now]\)转移到\(j\)这个点，我们转移的条件就是从\(i\)到\(j\)有一条边相连，之后由于\(dp[i][now]\)这个状态要走的路程是\(dis[i]+now\)，之后走了一条\(dis(i,j)\)的路径，到了\(j\)点还要走的就是\(dis[i]+now-dis(i,j)\)，相比最短路多走的就是\(dis[i]+now-dis(i,j)-dis[j]\)，于是我们直接在记搜里这样转移就好了
之后边界条件就是\(dp[n][0]=1\)
再来考虑一下如何判断无解的情况
显然如果有\(0\)环，而且环里有一个点到\(1\)和\(n\)的最短路的和走的冤枉路小于\(K\)，那么就可以在这个环里一直绕下去，就会出现无解了
显然我们只需要考虑\(0\)环的话，我们只需要建一张只有\(0\)边的图就行了，对这张图跑一个\(topsort\)或者\(tarjan\)，之后判断环里是否有这种点就可以啦
其实不不是非常难，主要就是\(dp\)的思想还有最简单的图论
但是我非常sb，我对于存在\(0\)环的图跑了最短路计数，显然这种图的最短路条数是无限的，这样跑出来显然是错的
代码
#include<iostream>
#include<cstring>
#include<cstdio>
#include<queue>
#define re register
#define maxn 100005
#define LL long long
#define mp std::make_pair
int n,m,K,num[2],cnt;
int mod;
typedef std::pair<int,int> pii;
std::priority_queue<pii,std::vector<pii>,std::greater<pii> > q;
struct node
{
int v,nxt,w;
}e[2][maxn<<1];
struct ZE
{
int v,nxt;
}E[maxn<<1];
int Head[maxn];
inline void write(int x)
{
if(x>9) write(x/10);
putchar(x%10+'0');
}
inline void add(int x,int y)
{
E[++cnt].v=y;
E[cnt].nxt=Head[x];
Head[x]=cnt;
}
int dp[maxn][51],Q[maxn],r[maxn];
int head[2][maxn],dis[2][maxn];
int f[maxn];
inline int read()
{
char c=getchar();
int x=0;
while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9')
x=(x<<3)+(x<<1)+c-48,c=getchar();
return x;
}
inline void add_edge(int x,int y,int z,int opt)
{
e[opt][++num[opt]].v=y;
e[opt][num[opt]].nxt=head[opt][x];
e[opt][num[opt]].w=z;
head[opt][x]=num[opt];
}
inline void dij(int s,int o)
{
memset(f,0,sizeof(f));
while(!q.empty()) q.pop();
dis[o][s]=0;
if(o) dp[s][0]=1;
q.push(mp(dis[o][s],s));
while(!q.empty())
{
int k=q.top().second;
q.pop();
if(f[k]) continue;
f[k]=1;
for(re int i=head[o][k];i;i=e[o][i].nxt)
if(dis[o][e[o][i].v]>dis[o][k]+e[o][i].w)
{
dis[o][e[o][i].v]=dis[o][k]+e[o][i].w;
q.push(mp(dis[o][e[o][i].v],e[o][i].v));
}
}
}
int dfs(int x,int now)
{
if(dp[x][now]) return dp[x][now]%mod;
for(re int i=head[0][x];i;i=e[0][i].nxt)
if(dis[1][x]-dis[1][e[0][i].v]-e[0][i].w+now>=0)
dp[x][now]=(dp[x][now]+dfs(e[0][i].v,dis[1][x]-dis[1][e[0][i].v]-e[0][i].w+now))%mod;
return dp[x][now]%mod;
}
int main()
{
int T=read();
while(T--)
{
cnt=num[0]=num[1]=0;
memset(Head,0,sizeof(Head));
memset(head,0,sizeof(head));
memset(dis,0x7f,sizeof(dis));
memset(dp,0,sizeof(dp));
memset(Q,0,sizeof(Q));
memset(r,0,sizeof(r));
n=read(),m=read(),K=read(),mod=read();
int x,y,z;
for(re int i=1;i<=m;i++)
{
x=read(),y=read(),z=read();
if(!z) add(x,y),r[y]++;
add_edge(x,y,z,0),add_edge(y,x,z,1);
}
dij(1,0),dij(n,1);
int tot=0;
for(re int i=1;i<=n;i++)
if(!r[i]) Q[++tot]=i;
for(re int i=1;i<=tot;i++)
{
int t=Q[i];
for(re int j=Head[t];j;j=E[j].nxt)
{
r[E[j].v]--;
if(!r[E[j].v]) Q[++tot]=E[j].v;
}
}
int flag=0;
for(re int i=1;i<=n;i++)
if(r[i]&&dis[0][i]+dis[1][i]<=dis[0][n]+K)
{
puts("-1");
flag=1;
break;
}
if(flag) continue;
int ans=0;
for(re int i=0;i<=K;i++) ans=(ans+dfs(1,i))%mod;
write(ans),putchar(10);
}
return 0;
}