【题解】洛谷P3953 逛公园（最短路 动态规划）

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;
}