送外卖2（Floyd+三进制状压dp）

题目链接
链接：https://ac.nowcoder.com/acm/problem/13252
来源：牛客网

题目描述
美团外卖日订单数已经超过1200万，实时调度系统是背后的重要技术支撑，其中涉及很多复杂的算法。下面的题目是某类场景的抽象。

一张 nn 个点 mm 条有向边的图上，有 qq 个配送需求，需求的描述形式为( s_i , t_i , l_i , r_is
i
​
,t
i
​
,l
i
​
,r
i
​
)，即需要从点 s_is
i
​
送到 t_it
i
​
， 在时刻 l_il
i
​
之后（包括 l_il
i
​
）可以在 s_is
i
​
领取货物，需要在时刻 r_ir
i
​
之前（包括 r_ir
i
​
）送达 t_it
i
​
，每个任务只需完成一次。 图上的每一条边均有边权，权值代表外卖配送员通过这条边消耗的时间。在时刻 00 有一个配送员在 点 11 上，求他最多能完成多少个配送任务。
在整个过程中，我们忽略了取餐与最后给用户递餐的时间（实际场景中这两个时间是无法省略的），只考虑花费在路程上的时间。另外，允许在一个点逗留。

输入描述:
第一行，三个正整数 n , m , q (2 ≤ n ≤ 20, 1 ≤ m ≤ 400, 1 ≤ q ≤ 10)。
接下来 m 行，每行三个正整数 u_i , v_i , c_i (1 ≤ u_i,v_i ≤ n, 1 ≤ c_i ≤ 20000)，表示有一条从 u_i 到 v_i 耗时为 c_i 的有向边。
接下来 q 行，每行四个正整数 s_i , t_i , l_i , r_i (1 ≤ s_i,t_i ≤ n, 1 ≤ l_i ≤ r_i ≤ 10^6)，描述一个配送任务。
输出描述:
一个整数，表示最多能完成的任务数量。
示例1
输入
复制
5 4 3
1 2 1
2 3 1
3 4 1
4 5 1
1 2 3 4
2 3 1 2
3 4 3 4
输出
复制
2

分析：

因为n比较小，可以想到用状压做。先用floud计算点和点之间的最短时间。
dp[i][j]:状态为i，当前停留在j点，存储的是此时的时间
sta[i][j]:当前状态为i时，第j个任务的状态，如果sta[i][j]=1，表示第j个任务已经取餐，但没有送达终点。如果sta[i][j]=2,表示第j个任务已经完成。如果sta[i][j]=0,表示第j个任务还没有取餐。
plu[i]:表示第i个任务从没有取餐到完成取餐，或已取餐还没有送达到送达需要加上的值，即time+plu[i]即改变了第i个任务的状态

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 2e5+10 ;
const ll INF = 3e3+10,M = 3e4+10;
int n,m,q;
ll g[30][30],l[30],r[30],dp[N][30];
int s[30],t[30],sta[N][30],plu[30];
void init()
{
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++) g[i][j]=1e12;
	for(int i=1;i<=n;i++) g[i][i]=0;
	plu[0]=1;
	for(int i=1;i<=q;i++) plu[i]=plu[i-1]*3;
	for(int i=0;i<=plu[q];i++)//遍历每个状态 
	{
		int tmp=i;
		for(int j=0;j<q;j++)//对每个状态，找到当前状态各任务的完成进度 
		{
			sta[i][j]=tmp%3;
			tmp/=3;
		}
	}
	for(int i=0;i<=plu[q];i++)
	{
		for(int j=1;j<=n;j++)
		dp[i][j]=1e12;
	}
}
void floyd()
{
	for(int k=1;k<=n;k++)
	{
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=n;j++)
			{
				if(g[i][j]>g[i][k]+g[k][j])
				g[i][j]=g[i][k]+g[k][j];
			}
		}
	}	
}
int main()
{
	scanf("%d%d%d",&n,&m,&q);
	init();
	int a,b;
	ll c,d;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%lld",&a,&b,&c);
		g[a][b]=min(g[a][b],c);
	} 
	floyd();
	for(int i=0;i<q;i++)
	{
		scanf("%d%d%lld%lld",&s[i],&t[i],&l[i],&r[i]);
	}
	dp[0][1]=0;//初始在1这个点 
	for(int i=0;i<plu[q];i++)//遍历状态 
	{
		for(int j=1;j<=n;j++)//枚举当前状态下停留的点 
		{
			if(dp[i][j]==1e12) continue; 
			for(int k=0;k<q;k++)//枚举所有的任务 
			{
				if(sta[i][k]==0)//第k个任务没有取餐，则去取餐的点 
				{
					int tmp=i+plu[k];//完成取餐后的状态 
					ll time;
					if(dp[i][j]+g[j][s[k]]<l[k])  time=l[k];//到达s[k]时还没有到取餐的时间，则要等到取餐的时间才能走 
					else time=dp[i][j]+g[j][s[k]];// 达s[k]时已经到了取餐的时间 
					dp[tmp][s[k]]=min(dp[tmp][s[k]],time);
				}
				else if(sta[i][k]==1)//已取餐但是还没有送达 
				{
					int tmp=i+plu[k];
					ll time;
					if(dp[i][j]+g[j][t[k]]<=r[k])//判断能否在规定时间前到达t[k] 
					{
						time=dp[i][j]+g[j][t[k]];
						dp[tmp][t[k]]=min(dp[tmp][t[k]],time);
					}			
				}
				else continue;
			}
		}
	}
	int ans=0;
	for(int i=0;i<plu[q];i++)
	{
		for(int j=1;j<=n;j++)
		{
			int cnt=0;
			if(dp[i][j]==1e12) continue;
			for(int k=0;k<q;k++)
			{
				if(sta[i][k]==2) cnt++;
			}
			if(ans<cnt) ans=cnt;
		}
	}
	printf("%d",ans);
	return 0;
}