洛谷3953【NOIP2017DAY1T3】 逛公园（最短路+DP）

氵谷传送门

【题目描述】

很仙的一道最短路，想出状态定义后还是比较好做的。

首先跑出1到所有点的最短路，记作dis[i]。

定义状态方程dp[x][i]表示到点x比最短路大i的路径条数，根据这个数据规模明显直接O(nk)枚举就行。

DP转移方程：，这里因为一个点可以被多次搜到，所以记忆化可以做到剪枝效果。而后面这个式子看着长其实移项后就是，其中表示如果走这条边比最短路长多少，然后用当前枚举长的长度减去这么一坨就是剩下的长度，就是a。

注意初始化。。。

【代码~】

#include<bits/stdc++.h>
using namespace std;
const int MAXN=1e5+10;
const int MAXM=2e5+10;
const int INF=0x3f3f3f3f;
#define pii pair<int,int>

int T;
int n,m,p,k,cnt,flag,ans;
int head[MAXN],viss[MAXN],dis[MAXN],dp[MAXN][51],vis[MAXN][51];
int headd[MAXN];
int nxt[MAXM],to[MAXM],w[MAXM];
int nxtt[MAXM],too[MAXM],ww[MAXM];

int Read(){
	int i=0,f=1;
	char c=getchar();
	for(;(c>'9'||c<'0')&&c!='-';c=getchar());
	if(c=='-')  f=-1,c=getchar();
	for(;c>='0'&&c<='9';c=getchar())  i=(i<<3)+(i<<1)+c-'0';
	return i*f;
}

void add(int x,int y,int z){
	++cnt;
	nxt[cnt]=head[x];
	head[x]=cnt;
	to[cnt]=y;
	w[cnt]=z;
	nxtt[cnt]=headd[y];
	headd[y]=cnt;
	too[cnt]=x;
	ww[cnt]=z;
}

void dij(int s){
	memset(dis,INF,sizeof(dis));
	dis[s]=0;
	priority_queue<pii > q;
	q.push(make_pair(0,s));
	while(!q.empty()){
		int u=q.top().second;
		q.pop();
//		viss[u]=0;
		for(int i=head[u];i!=-1;i=nxt[i]){
			int v=to[i];
			if(dis[u]+w[i]<dis[v]){
				dis[v]=dis[u]+w[i];
//				if(!viss[v])
				q.push(make_pair(-dis[v],v));
			}
		}
	}
}

int dfs(int u,int k){
//	cout<<u<<' '<<k<<'\n';
	if(~dp[u][k])  return dp[u][k];
	dp[u][k]=0,vis[u][k]=1;
	for(int i=headd[u];i!=-1;i=nxtt[i]){
		int v=too[i];
		int delta=dis[u]+k-dis[v]-ww[i];
		if(delta<0)  continue;
		if(vis[v][delta])  flag=1;
		dp[u][k]=(dp[u][k]+dfs(v,delta))%p;
	}
	vis[u][k]=0;
	return dp[u][k];
}

void init(){
	memset(dp,-1,sizeof(dp));
	memset(head,-1,sizeof(head));
	memset(headd,-1,sizeof(headd));
	memset(vis,0,sizeof(vis));
	memset(viss,0,sizeof(viss));
	cnt=0,flag=0;
}

int main(){
	T=Read();
	while(T--){
		init();
		n=Read(),m=Read(),k=Read(),p=Read();
		for(int i=1;i<=m;++i){
			int x=Read(),y=Read(),z=Read();
			add(x,y,z);
		}
		dij(1);
//		for(int i=1;i<=n;++i)  cout<<dis[i]<<' ';
//		cout<<'\n';
		if(!dis[n]){puts("-1");return 0;}
		dp[1][0]=1;
		ans=0;
		for(int i=0;i<=k;++i)  ans=(ans+dfs(n,i))%p;
		if(!flag)  cout<<ans<<'\n';
		else  puts("-1");
	}
}