LGOJ P3627 [APIO2009]抢掠计划 解题报告

最新推荐文章于 2020-05-02 17:15:04 发布
最新推荐文章于 2020-05-02 17:15:04 发布
阅读量170 收藏
点赞数
分类专栏： OI的那些事 文章标签： OI 题解 强连通分量
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_41996523/article/details/99482593
版权

文章目录

题目链接

解题方法

声明：部分参考@javalyc的题解：洛谷博客、博客园

A.解题思路1

首先，我们发现：如果我们抢掠了点 i i i，我们一定也可以抢掠掉 i i i所在的强连通分量，这样能够拿到更多的钱，结果肯定更优。故我们采用 T a r j a n Tarjan Tarjan算法缩点，求出图中每个强连通分量，计算每个强连通分量的总共的钱。

其次，缩完点后，我们可以在新图中，化点权为边权，然后使用SPFA算法跑最长路。然后判断酒吧位置的 d i s t dist dist的最大值。

注意：跑最长路时不要不要不要用Dijkstra，它的贪心在最长路中是 不！ 正！ 确！ 的！

为什么是不正确的？跑最长路，其实就相当于在边权全部变号，跑最短路。Dijkstra跑最短路有个前提：边权不能为负。

补充:我们不必从每个点开始跑个Tarjan，只要从起点开始跑。（反正其他点也到不了）。但这样就需要加些特殊判断。(见代码)

详细代码

#define rg register
#define inl inline
#define DEBUG printf("[Passing [%s] in line %d.]\n", __func__, __LINE__)
#define putline putchar('\n')
#define putsp putchar(' ')
#define Rep(a, s, t) for(rg int a = s; a <= t; a++)
#define Repdown(a, t, s) for(rg int a = t; a >= s; a--)
typedef long long ll;
#include<cstdio>
#define rs freopen("test.in", "r", stdin), freopen("test.out", "w", stdout)
struct IO
{
	IO(int set = 0) {if(set) rs;} 
	void RS() {rs;} 
	template<typename T> inline IO r(T& x)const
	{
	    x = 0; T f = 1; char ch = getchar();
	    for(; ch < '0' || ch > '9'; ch = getchar()) if(ch == '-') f = -1;
	    for(; ch >= '0' && ch <= '9'; ch = getchar()) x = x * 10 + int(ch - '0');
	    x *= f; return *this;
	}
	template<typename T> inline IO w(T x)const
	{
	    if(x < 0) {putchar('-'); x = -x;}
	    if(x >= 10) w(x / 10);
	    putchar(x % 10 + '0'); return *this;
	}
	template<typename T> inline IO wl(const T& x)const {w(x), putline; return *this;}
	template<typename T> inline IO ws(const T& x)const {w(x), putsp; return *this;}
	inline IO l() {putline; return *this;}
	inline IO s() {putline; return *this;}
}io;
template<typename T> inline T Max(const T& x, const T& y) {return y < x ? x : y;}
template<typename T> inline T Min(const T& x, const T& y) {return y < x ? y : x;}
template<typename T> inline void Swap(T& x, T& y) {T tmp = x; x = y; y = tmp;}
#include<stack>
#include<queue>
#include<cstring>
using namespace std;
const int MAXN = 500005, MAXM = 500005;
int N, M;
struct Edge
{
	int u, v, nxt, w;
}e[MAXM];
int head[MAXN], cnt;
void addedge(int u, int v, int w)
{
	e[++cnt].u = u;
	e[cnt].v = v;
	e[cnt].w = w;
	e[cnt].nxt = head[u];
	head[u] = cnt;
}
int mon[MAXN];
int dfn[MAXN], low[MAXN];
int ins[MAXN], tim;
int gg, g[MAXN];
int val[MAXN];
stack<int> stk;
void Tarjan(int u)
{
	dfn[u] = low[u] = ++tim;
	stk.push(u); ins[u] = 1;
	for(rg int i = head[u]; i; i = e[i].nxt)
	{
		int v = e[i].v;
		if(!dfn[v])
		{
			Tarjan(v);
			low[u] = Min(low[u], low[v]); 
		}
		else if(ins[v]) low[u] = Min(low[u], dfn[v]);
	}
	if(dfn[u] == low[u])
	{
		++gg;
		int now;
		do
		{
			now = stk.top(); stk.pop();
			ins[now] = 0;
			g[now] = gg;
			val[gg] += mon[now];
		}while(now != u);
	}
}
queue<int> q; 
int dist[MAXN], vist[MAXN];
void Spfa(int S)
{
	S = g[S];
	memset(dist, -127, sizeof dist);
	memset(vist, 0x00, sizeof vist);
	dist[S] = val[S];
	q.push(S);
	vist[S] = 1;
	while(!q.empty())
	{
		int u = q.front(); q.pop();
		vist[u] = 0;
		for(rg int i = head[u]; i; i = e[i].nxt)
		{
			int v = e[i].v, w = e[i].w;
			if(dist[v] < dist[u] + w)
			{
				dist[v] = dist[u] + w;
				if(!vist[v])
				{
					vist[v] = 1;
					q.push(v);
				}
			}
		} 
	}
}
int main()
{
    //io.RS();
	io.r(N).r(M);
	for(rg int i = 1; i <= M; i++)
	{
		int u, v;
		io.r(u).r(v);
		addedge(u, v, 0);
	}
	for(rg int i = 1; i <= N; i++) io.r(mon[i]);
	int S; io.r(S);
	Tarjan(S);
	cnt = 0;
	memset(head, 0x00, sizeof head);
	for(rg int i = 1; i <= M; i++)
	{
		int u = e[i].u, v = e[i].v;
		if(!g[u] || ! g[v]) continue;
		if(g[u] != g[v])
			addedge(g[u], g[v], val[g[v]]);
	}
	/*
	看到了吗，这种骚操作：
	只清空cnt和head数组，然后在老的e数组中枚举每条边
	（前向星存图的时候存个u是个好习惯）。
	这样能保证每条老边都能被考虑，并且
	不用存原来的边
	*/
	Spfa(S);
	int P; io.r(P);
	int ans = 0;
	for(rg int i = 1; i <= P; i++)
	{
		int k; io.r(k);
		k = g[k];
		ans = Max(ans, dist[k]);
	}
	io.wl(ans);
    return 0;
}

B.解题思路2

DP。

javalyc:DAG上搞DP是基本操作啊

嗯，所以这题我们可以先Tarjan缩点把图变DAG，再在新图上做DP。

由于状态 f [ i ] f[i] f[i]定义的是从 i i i点到 S S S点能得到的最大点权和，所以同上，我们建图只能建以 S S S为根的子树

一般来说，这种DP都有两种写法：

1）以BFS为基础的。 按拓扑序进行DP

2）以DFS为基础的。（其实本质上也是一种拓扑排序）。自然大部分人都喜欢这种咯，多好写。。。（本题解便以这种思路呈现）

同样，这道题的DP也有两种形式。（好像每个DP都有这两种形式）

1）顺推大法好： 由 f [ u ] f[u] f[u]的值，可以通过边 ( u − &gt; v ) (u-&gt;v) (u−>v)，更新 f [ v ] f[v] f[v]的值。

2）记忆化搜索： 考虑 f [ u ] f[u] f[u]的值，可以由所有边 ( k − &gt; u ) (k-&gt;u) (k−>u)得到，我们就可以先记搜出所有 f [ k ] f[k] f[k]的值,用所有 f [ k ] f[k] f[k]更新 f [ u ] f[u] f[u]的值。实际操作时考虑反向建边（本题解便采用这种方法）

详细代码

#define rg register
#define inl inline
#define DEBUG printf("[Passing [%s] in line %d.]\n", __func__, __LINE__)
#define putline putchar('\n')
#define putsp putchar(' ')
#define Rep(a, s, t) for(rg int a = s; a <= t; a++)
#define Repdown(a, t, s) for(rg int a = t; a >= s; a--)
typedef long long ll;
#include<cstdio>
#define rs freopen("test.in", "r", stdin), freopen("test.out", "w", stdout)
struct IO
{
	IO(int set = 0) {if(set) rs;} 
	void RS() {rs;} 
	template<typename T> inline IO r(T& x)const
	{
	    x = 0; T f = 1; char ch = getchar();
	    for(; ch < '0' || ch > '9'; ch = getchar()) if(ch == '-') f = -1;
	    for(; ch >= '0' && ch <= '9'; ch = getchar()) x = x * 10 + int(ch - '0');
	    x *= f; return *this;
	}
	template<typename T> inline IO w(T x)const
	{
	    if(x < 0) {putchar('-'); x = -x;}
	    if(x >= 10) w(x / 10);
	    putchar(x % 10 + '0'); return *this;
	}
	template<typename T> inline IO wl(const T& x)const {w(x), putline; return *this;}
	template<typename T> inline IO ws(const T& x)const {w(x), putsp; return *this;}
	inline IO l() {putline; return *this;}
	inline IO s() {putline; return *this;}
}io;
template<typename T> inline T Max(const T& x, const T& y) {return y < x ? x : y;}
template<typename T> inline T Min(const T& x, const T& y) {return y < x ? y : x;}
template<typename T> inline void Swap(T& x, T& y) {T tmp = x; x = y; y = tmp;}
#include<stack>
#include<queue>
#include<cstring>
using namespace std;
const int MAXN = 500005, MAXM = 500005;
int N, M;
struct Edge
{
	int u, v, nxt, w;
}e[MAXM];
int head[MAXN], cnt;
void addedge(int u, int v, int w)
{
	e[++cnt].u = u;
	e[cnt].v = v;
	e[cnt].w = w;
	e[cnt].nxt = head[u];
	head[u] = cnt;
}
int mon[MAXN];
int dfn[MAXN], low[MAXN];
int ins[MAXN], tim;
int gg, g[MAXN];
int val[MAXN];
stack<int> stk;
void Tarjan(int u)
{
	dfn[u] = low[u] = ++tim;
	stk.push(u); ins[u] = 1;
	for(rg int i = head[u]; i; i = e[i].nxt)
	{
		int v = e[i].v;
		if(!dfn[v])
		{
			Tarjan(v);
			low[u] = Min(low[u], low[v]); 
		}
		else if(ins[v]) low[u] = Min(low[u], dfn[v]);
	}
	if(dfn[u] == low[u])
	{
		++gg;
		int now;
		do
		{
			now = stk.top(); stk.pop();
			ins[now] = 0;
			g[now] = gg;
			val[gg] += mon[now];
		}while(now != u);
	}
}
int f[MAXN];
void dp(int u) //这里u是新图中点的编号 
{
	if(f[u]) return;
	f[u] = val[u];
	for(rg int i = head[u]; i; i = e[i].nxt)
	{
		int v = e[i].v;
		dp(v);
		f[u] = Max(f[u], f[v] + val[u]);
	}
}
int main()
{
    //io.RS();
	io.r(N).r(M);
	for(rg int i = 1; i <= M; i++)
	{
		int u, v;
		io.r(u).r(v);
		addedge(u, v, 0);
	}
	for(rg int i = 1; i <= N; i++) io.r(mon[i]);
	int S; io.r(S);
	Tarjan(S);
	cnt = 0;
	memset(head, 0x00, sizeof head);
	for(rg int i = 1; i <= M; i++)
	{
		int u = e[i].u, v = e[i].v;
		if(!g[u] || ! g[v]) continue;
		if(g[u] != g[v])
			addedge(g[v], g[u], val[g[v]]);
		//反向建边
	}
	for(rg int i = 1; i <= gg; i++)
		dp(i);
	int P; io.r(P);
	int ans = 0;
	for(rg int i = 1; i <= P; i++)
	{
		int k; io.r(k);
		k = g[k];
		ans = Max(ans, f[k]);
	}
	io.wl(ans);
    return 0;
}
非递归！APIO2009atm[抢掠计划]题解
t14t41t的专栏
10-13 1695
题目描述 输入描述 Input Description 第一行包含两个整数N、M。N 表示路口的个数，M 表示道路条数。接下来 M 行，每行两个整数，这两个整数都在1 到N 之间，第i+1 行的两个整数表示第 i 条道路的起点和终点的路口编号。接下来N 行，每行一个整数，按顺序表示每 个路口处的ATM 机中的钱数。接下来一行包含两个整数S、P，S 表示市中心的 编号，也就是出发的路口。P
P3627 [APIO2009]抢掠计划
starlords的博客
10-15 450
题目描述Siruseri 城中的道路都是单向的。不同的道路由路口连接。按照法律的规定， 在每个路口都设立了一个 Siruseri 银行的 ATM 取款机。令人奇怪的是，Siruseri 的酒吧也都设在路口，虽然并不是每个路口都设有酒吧。Banditji 计划实施 Siruseri 有史以来最惊天动地的 ATM 抢劫。他将从市中心 出发，沿着单向道路行驶，抢劫所有他途径的 ATM 机，最终他将在一个酒
【luogu P3627 [APIO2009]抢掠计划】 题解
07-27 117
题目链接：https://www.luogu.org/problemnew/show/P3627 把点权转化到边权上去。 #include <stack> #include <queue> #include <cstdio> #include <cstring> #include <iostream> #include <...
【题解】洛谷P3627 [APIO2009]抢掠计划（缩点+SPFA）
09-26 158
洛谷P3627：https://www.luogu.org/problemnew/show/P3627 思路 由于有强连通分量 所以我们可以想到先把整个图缩点 缩点完之后再建一次图 把点权改为边权 并把边权转为负数 即可用SPFA求最短路间接求最长路了 最后我们查询所有的酒吧 跳出最大的ans即可 思路简单 但是代码有些冗长 代码 #include<i...
洛谷P3627 [APIO2009]抢掠计划——题解
CleverLarry的博客
05-16 459
题目传送门 题目大意： 看题面吧。 思考过程&amp;具体做法： 先缩点，然后从市中心所在的新点开始做dp，最后对有酒吧的新点统计答案。 代码： #include &lt;bits/stdc++.h&gt; using namespace std; const int maxn=5e5+1000; struct stu { int to,next; }road[m...
[APIO2009]抢掠计划 解题报告
weixin_30906671的博客
09-18 71
这道题码的十分痛苦。 题目链接 思路：tarjan缩点，重新建图，点权转边权，边权变为负值，跑一遍spfa求最长路即可。 思路很简单，但是码量有点痛苦，打了一百行，打两个板子就发现了许多问题，例如：程序出错的时候不知道是哪里出了问题，找错十分麻烦。我自己出现的错误还是挺多的，模板不是很熟悉。 代码： #include<cstdio> #include<ios...
APIO2009 抢掠计划
AChunter的专栏
04-15 1840
【题意】 给一张图，每个点有一个ATM机，而其中的部分点是酒吧（可作为终点）。 抢掠计划就是从起点（1号点）出发，到终点的途中，将ATM机中的钱抢走。 问最多能抢到多少钱。 【题解】 很显然，由于图可能有环，所以自然地想到了SCC+DP。 题目不难，但数据很大，堆栈会溢出。表示懒得模拟堆栈了。 【代码】 {$M 100000000} uses math; type edge=r
C++解题报告：抢掠计划 [ APIO 2009 ] 解析
qq_44013342的博客
06-25 510
题目描述 Siruseri 城中的道路都是单向的。不同的道路由路口连接。按照法律的规定， 在每个路口都设立了一个 Siruseri 银行的 ATM 取款机。令人奇怪的是，Siruseri 的酒吧也都设在路口，虽然并不是每个路口都设有酒吧。 Banditji 计划实施 Siruseri 有史以来最惊天动地的 ATM 抢劫。他将从市中心 出发，沿着单向道路行驶，抢劫所有他途径的 ATM 机，最终...
【Bzoj1179】[Apio2009]抢掠计划atm
剑锋OI SharpSwordOI
11-08 286
时间限制：1S / 空间限制：256MB 【在线测试提交传送门】【问题描述】 ![这里写图片描述](http://www.lydsy.com/JudgeOnline/images/1179.jpg) 【输入格式】第一行包含两个整数N、M。N表示路口的个数，M表示道路条数。 接下来M行，每行两个整数，这两个整数都在1到N之间，第i+1行的两个整数表示第i条道路的起点和终点的路口编号。 接下来
【APIO2009】抢掠计划（scc+spfa）
Patrickpwq
10-14 258
&#13; 拿到题目首先想到：一个强连通分量里的点都可以互相走到&#13; 所以我们先缩点&#13; 缩完点后只需用spfa跑一遍最长路 最后再统计一下有哪些含有酒吧的强连通分量 取最大值&#13; 就是代码有点冗长&#13; &#13; #include&lt;bits/stdc++.h&gt;&#13; #define N 500005&#13; #define INF ...
洛谷 P3627 [APIO2009]抢掠计划 题解
weixin_30598225的博客
07-26 78
Analysis 建图+强连通分量+SPFA求最长路 但要保证最后到达的点中包含酒馆 虽然思路并不难想，但要求的代码能力很高。 1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<algorithm> 5 #i...
题解 P3627 [APIO2009]抢掠计划
ZHK's blog
05-02 598
Siruseri 城中的道路都是单向的。不同的道路由路口连接。按照法律的规定，在每个路口都设立了一个 Siruseri 银行的 ATM 取款机。令人奇怪的是，Siruseri 的酒吧也都设在路口，虽然并不是每个路口都设有酒吧。 Banditji 计划实施 Siruseri 有史以来最惊天动地的 ATM 抢劫。他将从市中心出发，沿着单向道路行驶，抢劫所有他途径的 ATM 机，最终他将在一个酒吧庆祝...
题解 P3627 【[APIO2009]抢掠计划】
weixin_30808575的博客
07-23 103
这道题前面的思路和其他题解一样，都是tarjan缩点。然后是要求金币的最大值，我用的是记忆化搜索，用jiyi[i]来表示从编号为i的联通块出发，最多可以抢到多少钱，搜索时如果遇到已经搜过的点（知道从那个点开始最多能抢多少钱）就直接return。因为记忆化搜索每个点只会遍历一次，所以记忆化搜索的复杂度应该是O(n)。 另外我的码风有点奇怪，先解释一下：stk[tk]分别是栈和栈的指针（tarjan...
【日常学习】【强连通分量tarjan缩点】codevs1611 抢掠计划题解
西家家
10-16 1770
题目描述 Description Siruseri 城中的道路都是单向的。不同的道路由路口连接。按照法律的规定， 在每个路口都设立了一个Siruseri 银行的ATM 取款机。令人奇怪的是，Siruseri 的酒吧也都设在路口，虽然并不是每个路口都设有酒吧。 Banditji 计划实施Siruseri 有史以来最惊天动地的ATM 抢劫。他将从市中心 出发，沿着单向道路行驶，抢劫所有他
抢掠计划
10-26 219
Siruseri 城中的道路都是单向的。不同的道路由路口连接。按照法律的规定，在每个路口都设立了一个 Siruseri 银行的 ATM 取款机。令人奇怪的是，Siruseri 的酒吧也都设在路口，虽然并不是每个路口都设有酒吧。Banditji 计划实施 Siruseri 有史以来最惊天动地的 ATM 抢劫。他将从市中心出发，沿着单向道路行驶，抢劫所有他途径的 ATM 机，最终他将在一个酒...
MongoDB分片集群搭建教程：副本集创建与数据分片
最新发布
11-12
内容概要：本文提供了详细的MongoDB分片集群的搭建指导，涵盖了从环境准备、配置文件编写、副本集的建立、主节点的选择、配置服务器和数据分片服务器的配置到最后的路由节点的搭建与操作整个流程，以及对数据库的哈希与范围两种分片策略的应用介绍和具体命令执行。 适合人群：熟悉NoSQL数据库概念并对MongoDB有一定了解的技术人员，尤其是在大型数据管理和分布式数据库架构设计中有需求的开发者。 使用场景及目标：帮助技术人员掌握构建高效能、高可用性的MongoDB分片集群的方法，适用于处理大规模、实时性强的数据存储与读取场景。 其他说明：文中通过实例演示了每个步骤的具体操作方法，便于跟随文档实操，同时也介绍了可能遇到的问题及其解决方案，如在没有正确配置的情况下试图写入数据时出现错误等情况的处理。
CPPC++_嵌入式硬件的物联网解决方案blinker库与Arduino ESP8266 ESP32一起工作.zip
11-12
CPPC++_嵌入式硬件的物联网解决方案blinker库与Arduino ESP8266 ESP32一起工作
CPPC++_逆向调用QQ Mojo IPC与WeChat XPlugin.zip
11-12
CPPC++_逆向调用QQ Mojo IPC与WeChat XPlugin
APIO2009竞赛试题解析-采油区域问题
"APIO2009试题-中文版，2014年的亚太地区信息学奥林匹克竞赛试题，包括三个题目：采油区域、会议中心、抢掠计划。" 在2009年的亚太地区信息学奥林匹克竞赛（APIO）中，参赛者面临的挑战之一是名为"采油区域"的问题...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

日居月诸Rijuyuezhu

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值