NOIP2017最难的题目。。。
这里给一种比较方便理解的做法。
拿到这个题目,啥也不用想,首先得把最短路求出来。然而求最短路时我们要反着建图,也就是求出n到其他所有点的最短路。为什么这样做呢?因为这样可以避免正向某个点无法到达n点的情况。求出最短路后,我们可以利用动态规划解决这个问题。首先考虑没有0边的情况。我们开一个二维数组f[u][know],代表在反向图中从u到n与从u到n的最短路径之差等于know的路径的条数。那么我们得到结果就是∑f[1][know] (0<=know<=k)。为了求出f,我们可以通过正向图dfs的方法。接下来我们就需要得到know与转移后得到新的know值(以下称为know')的关系。这里我们可以得到一个表达式:dis[u]+know=dis[y]+know'+length[i]。只有满足这个条件时,我们更新得到的新的路径才对原来的know有贡献。求和即可。注意初始化f[n][know]为1。对于有0边的情况,我们开w数组记录是否有0边。如果成功更新就没有,反之就有,flag=1,返回主函数输出-1.为了卡过这道题,我们应当用读入优化,而且还要在主函数里计算f从0到k的累加值时判断如果出现flag为true就输出-1;
#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<queue>
#define ll long long
#define INF 0x3f3f3f3f
using namespace std;
const int maxn=100010;
const int maxm=200010;
int n,m,T,k,mod;
int head[maxn],nnext[maxm*2],to[maxm*2],length[maxm*2],tot;
int head1[maxn],nnext1[maxm*2],to1[maxm*2],length1[maxm*2],tot1;
int dis[maxn];
bool b[maxn];
int wd[maxn][55],f[maxn][55];
bool flag;
inline int read()
{
char c=getchar();
while(c<'0'||c>'9') c=getchar();
int ans=0;
while(c>='0'&&c<='9')
{
ans=ans*10+c-'0';
c=getchar();
}
return ans;
}
void add(int x,int y,int l)
{
tot++;
nnext[tot]=head[x];
head[x]=tot;
to[tot]=y;
length[tot]=l;
}
void add1(int x,int y,int l)
{
tot1++;
nnext1[tot1]=head1[x];
head1[x]=tot1;
to1[tot1]=y;
length1[tot1]=l;
}
void csh()
{
memset(head,0,sizeof(head));
memset(head1,0,sizeof(head1));
tot=tot1=0;
flag=false;
memset(f,0,sizeof(f));
}
void spfa()
{
queue<int> q;
memset(b,false,sizeof(b));
memset(dis,INF,sizeof(dis));
dis[n]=0;
q.push(n);
b[n]=true;
while(!q.empty())
{
int now=q.front();
q.pop();
b[now]=false;
for(int i=head1[now];i;i=nnext1[i])
{
int y=to1[i];
if(dis[y]>dis[now]+length1[i])
{
dis[y]=dis[now]+length1[i];
if(!b[y])
{
b[y]=true;
q.push(y);
}
}
}
}
}
int dfs(int u,int know)
{
if(wd[u][know])
{
flag=true;
return 0;
}
if(f[u][know]>0) return f[u][know];
wd[u][know]=1;
int sum=0;
for(int i=head[u];i;i=nnext[i])
{
// int y=;
int tmp=know+dis[u]-dis[to[i]]-length[i];
if(tmp<0||tmp>k) continue;
sum=(sum+dfs(to[i],tmp))%mod;
if(flag) return 0;
}
if(u==n&&know==0) sum=1;
wd[u][know]=0;
f[u][know]=sum;
return sum;
}
int main()
{
// freopen("park.in","r",stdin);
// freopen("park.out","w",stdout);
scanf("%d",&T);
while(T--)
{
csh();
scanf("%d%d%d%d",&n,&m,&k,&mod);
for(int i=1;i<=m;i++)
{
int x,y,z;
x=read(),y=read(),z=read();
add(x,y,z);
add1(y,x,z);
}
spfa();
int ans=0;
memset(wd,0,sizeof(wd));
for(int i=0;i<=k;i++)
{
ans=(ans+dfs(1,i))%mod;
if(flag) break;
}
if(flag) cout<<"-1"<<endl;
else printf("%d\n",ans);
}
fclose(stdin);
fclose(stdout);
return 0;
}