AcWing 383. 观光(上课摸鱼(×)开摆(√))

题目

题意: 给定n个点，m条有向边。给定起点S和终点T。求最短路的方案数 + 仅比最短路长度大1的路径的方案数.

思路: Dij + dp.涉及到dp我就不会了。
最短路方案数易求,比最短路长度大1的路径方案数是关键。Dij预处理距离数组dist[]
考虑dp: f[u][k]: 表示到u点长度为dist[u] + k的方案数
答案为f[T][0] + f[T][1]
转移: f[j][k] = f[j][k] + f[u][k], if(dist[j] == dist[u] + w[i]). 到达u长度为dist[u] + k的方案数,会被到j长度为dist[j] + k所继承，因为j恰好可以通过u到达。
f[j][1] = f[j][1] + f[u][0], if(dist[j] + 1 == dist[u] + w[i]) 到达j长度为dist[j] + 1的方案数，可以继承到达u长度为dist[u]的方案数。

// Problem: 观光
// Contest: AcWing
// URL: https://www.acwing.com/problem/content/385/
// Memory Limit: 64 MB
// Time Limit: 1000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<complex>
#include<cstring>
#include<cmath>
#include<vector>
#include<map>
#include<unordered_map>
#include<list>
#include<set>
#include<queue>
#include<stack>
#define OldTomato ios::sync_with_stdio(false),cin.tie(nullptr),cout.tie(nullptr)
#define fir(i,a,b) for(int i=a;i<=b;++i) 
#define mem(a,x) memset(a,x,sizeof(a))
#define p_ priority_queue
// round() 四舍五入 ceil() 向上取整 floor() 向下取整
// lower_bound(a.begin(),a.end(),tmp,greater<ll>()) 第一个小于等于的
// #define int long long //QAQ
using namespace std;
typedef complex<double> CP;
typedef pair<int,int> PII;
typedef long long ll;
// typedef __int128 it;
const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const ll inf = 1e18;
const int N = 1e3+10;
const int M = 2e5+10;
const int mod = 1e9+7;
const double eps = 1e-6;
inline int lowbit(int x){ return x&(-x);}
template<typename T>void write(T x)
{
    if(x<0)
    {
        putchar('-');
        x=-x;
    }
    if(x>9)
    {
        write(x/10);
    }
    putchar(x%10+'0');
}
template<typename T> void read(T &x)
{
    x = 0;char ch = getchar();ll f = 1;
    while(!isdigit(ch)){if(ch == '-')f*=-1;ch=getchar();}
    while(isdigit(ch)){x = x*10+ch-48;ch=getchar();}x*=f;
}
int n,m,k,T;
int h[N],e[M],w[M],ne[M],idx = 0;
void add(int a,int b,int c)
{
	e[idx] = b,w[idx] = c,ne[idx] = h[a],h[a] = idx++;
}
int dist[N];
bool vis[N];
int cnt[N];
int st,ed;
void Dij(int S)
{
	mem(dist,0x3f);
	mem(vis,false);
	mem(cnt,0);
	priority_queue<PII,vector<PII>,greater<PII> > q;
	cnt[S] = 1;
	dist[S] = 0;
	q.push({0,S});
	while(q.size())
	{
		PII tmp = q.top(); q.pop();
		int u = tmp.second;
		if(vis[u]) continue;
		vis[u] = true;
		for(int i=h[u];~i;i=ne[i])
		{
			int j = e[i];
			if(dist[j] > dist[u] + w[i])
			{
				dist[j] = dist[u] + w[i];
				q.push({dist[j],j});
				cnt[j] = cnt[u];
			}
			else if(dist[j] == dist[u] + w[i])
			{
				cnt[j] += cnt[u];
			}
		}
	}
}
struct s2 //按照dist数组排序
{
	int id;
	int dis;
	bool operator<(const s2&rhs)
	{
		return dis < rhs.dis;
	}
}va[N];
int f[N][2]; //表示到达i点，且距离为 dist[i] + j的道路数量
void solve()
{
   mem(f,0);
   mem(h,-1); idx = 0;
   read(n); read(m);
   for(int i=0;i<m;++i)
   {
   	   int x,y,z; read(x),read(y),read(z);
   	   add(x,y,z);
   }
   read(st),read(ed);
   Dij(st);
   
   for(int i=1;i<=n;++i) 
   {
   	  va[i].id = i;
   	  va[i].dis = dist[i];
   }
   sort(va+1,va+n+1);
   
   f[st][0] = 1;
   for(int k=0;k<=1;++k)
   {
   	  for(int t=1;t<=n;++t)
   	  {
   	  	 int u = va[t].id;
   	  	 for(int i=h[u];~i;i=ne[i])
   	  	 {
   	  	 	int j = e[i];
   	  	 	if(dist[u] + w[i] == dist[j]) f[j][k] += f[u][k];
   	  	 	if(k == 0 && dist[u] + w[i] == dist[j]+1) f[j][1] += f[u][0];
   	  	 }
   	  }
   }
   int ans = f[ed][0] + f[ed][1];
   
   
   write(ans); puts("");
}
signed main(void)
{  
   // T = 1;
   // OldTomato; cin>>T;
   read(T);
   while(T--)
   {
   	 solve();
   }
   return 0;
}