洛谷P6512 【[QkOI#R1] Quark and Flying Pigs】题解

最新推荐文章于 2020-05-10 18:27:33 发布

原创最新推荐文章于 2020-05-10 18:27:33 发布 · 325 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

本文详细解析了一种用于抓取飞猪的算法，通过排序、最短路径计算及动态规划，实现对飞猪出现时间及位置的有效追踪，旨在最大化捕获数量。文章介绍了算法的四个部分：排序飞猪出现时间、计算最短路径、动态规划求解最长路径及调整起始条件。

蛮综合的一道好题！

Solution

Part 1

首先，我们把所有猪猪出现的时间从小到大排序。

考虑: 如果要抓第 $i$ 只猪猪，那么第 $j (j > i)$ 只猪猪能否也被抓到。注意这里眼里只有第i只和第j只猪猪，没有别的猪猪。

假定第 $i$ 只猪猪在 $t 1$ 时刻， $x 1$ 号节点上出现；第 $j$ 只猪猪在第 $t 2$ 时刻， $x 2$ 号节点上出现 $(j > i)$ 。如果从 $x 1$ 到 $x 2$ 的最短耗时不大于 $t 2$ 与 $t 1$ 的差(时间差)，那么它们就可能一起被抓到；否则不可能一起被抓到（因为抓了一只就来不及赶过去抓另一只了呀）。

注意我们这里从 $x 1$ 到 $x 2$ 的最短耗时就是从 $x 1$ 到 $x 2$ 的最短路径；由于 $n$ 不是很大，没必要用 $P 5905$ Johnson全源最短路的板子，因为我们用简单的 $F l o y d$ 就可以轻松跑过。然后两重循环枚举飞猪，再用上面的方法看看它们可不可能一起被抓到即可。

同时发现，如果第 $i$ 只飞猪与第 $j$ 只飞猪能够同时被抓到，就相当于从 $i$ 到 $j$ 连了一条单向边。

Part 2

仔细研究上一个Part的最后一句话。

恍然大悟的您发现: 我们就是要找到从 $1$ 到 $k$ 的最长路径，并且每个节点只能经过一次（一只可爱的猪猪只能被抓一次啊），答案就是上述路径的长度(每条边的长度均为 $1$ )。而且，因为总是从编号小的点连到编号大的点，所以说不可能有环。

下面就是一个显而易见的 $d p$ 了。拓扑排序没有必要，毕竟数据不大。

状态设计，即 $dp_i$ 表示从 $i$ 到 $k$ 号节点的最长路径( $k$ 为飞猪的数量)。

状态转移需要想一想。每一个节点都继承了它指向的点，那么我们就可以找到它所有指向的点中 $d p$ 值最大的一个，并让这个节点跟在它后面。

for (int i=head[now];i;i=e[i].next)
{
	if (e[i].to>now)  maxv=max(maxv,dp[e[i].to]);
}//找到最大值
dp[now]=maxv+1;//跟在它后面

Part 3

哎呀哎样例过不了了。怎么回事呢？

容易发现，从 $1$ 号节点到 $1$ 号猪猪(最早出现的猪猪)也需要时间的。

那么怎么办呢？很简单啊。

只需要加一个 $0$ 号猪猪，它在 $0$ 时刻， $1$ 号节点出现(即开始就能抓到它)。然后其他的正常操作就可以啦，记得带上可怜的 $0$ 号猪猪哦。

而由于多了一只猪猪并且它肯定能够被抓到，所以答案还要减一。综上所述，答案为 $dp_0-1$ 。

Part 4

①时间复杂度: $O(m+n^3+k^2)$

②做题时间: $27 m i n$

③综合难度: 绿+/蓝-

④个人思路标签: 图论+最短路+Floyd+快速排序+动规dp(递推)。

注意不要开long long，搞不好这样可能会导致MLE甚至TLE。

Code

#include <bits/stdc++.h>
#define inf 5000005
using namespace std;

int n,m,k,cnt=0;
int GA[205][205],head[5005],dp[5005];

struct edge
{
	int next;
	int to;
}e[30000000];

struct node
{
	int t;
	int rt;
}pig[5005];

bool cmp(node a,node b)
{
	return a.t<b.t;
}

inline void add_edge(int u,int v)
{
	cnt++;
	e[cnt].to=v;
	e[cnt].next=head[u];
	head[u]=cnt;
}

signed main()
{
	cin>>n>>m>>k;
	for (int i=1;i<=n;i++)
	{
		for (int j=1;j<=n;j++)
		{
			if (i!=j)  GA[i][j]=inf;
			else GA[i][j]=0;
		}
	}
	for (int i=1;i<=m;i++)
	{
		int u,v,w;
		cin>>u>>v>>w;
		GA[u][v]=w;
		GA[v][u]=w;
	}
	for (int k=1;k<=n;k++)
	{
		for (int i=1;i<=n;i++)
		{
			for (int j=1;j<=n;j++)
			{
				if (GA[i][k]+GA[k][j]<GA[i][j])
				{
					GA[i][j]=GA[i][k]+GA[k][j];
					GA[j][i]=GA[i][k]+GA[k][j];
				}
			}
		}
	}
	for (int i=1;i<=k;i++)  cin>>pig[i].t>>pig[i].rt;
	
	sort(pig+1,pig+k+1,cmp);
	pig[0].t=0,pig[0].rt=1;
	
	for (int i=0;i<=k;i++)
	{
		for (int j=i+1;j<=k;j++)
		{
			if (GA[pig[i].rt][pig[j].rt]<=pig[j].t-pig[i].t)  add_edge(i,j);
		}
	}
	for (int now=k;now>=0;now--)
	{
		int maxv=0;
		for (int i=head[now];i;i=e[i].next)
		{
			if (e[i].to>now)  maxv=max(maxv,dp[e[i].to]);
		}
		dp[now]=maxv+1;
	}
	cout<<dp[0]-1<<endl;
	
	return 0;
}

撒花✿✿ヽ(°▽°)ノ✿撒花