洛谷 题解 P3627 【[APIO2009]抢掠计划】
图论 tarjan缩点+最短路 的一道题
- tarjan求强连通分量(为以后缩点打下良好的基础)
(如果不会tarjan的请点击这儿)
你需要的东西:
(1)、dfn[],表示这个点在dfs时是第几个被搜到的。
(2)、low[],表示这个点以及其子孙节点连的所有点中dfn最小的值
(3)、stack[],表示当前所有可能能构成是强连通分量的点。
(4)、vis[],表示一个点是否在stack中。
(5)、color[],记录每一个点强连通分量的编号。
(6)、deep,记录dfs树的深度
inline void tarjan(int now)
{
dfn[now]=++deep;
low[now]=deep;
vis[now]=1;
st.push(now);
for(int i=0;i<ver[now].size();i++)
{
int x=ver[now][i];
if(!dfn[x])
{
tarjan(x);
low[now]=min(low[now],low[x]);
}
else
{
if(vis[x])
low[now]=min(low[now],low[x]);
}
}
if(dfn[now]==low[now])
{
color[now]=++sum;
vis[now]=0;
while(st.top()!=now)
{
color[st.top()]=sum;
vis[st.top()]=0;
st.pop();
}
st.pop();
}
}
- 缩点(去除图中的环)
//重点:这里建新图是依托强连通分量的编号来建的
for(int i=1;i<=n;i++)
{
for(int j=0;j<ver[i].size();j++)
{
int x=ver[i][j];
if(color[i]!=color[x])//如果不是属于同一个强连通分量中,那么就合并
{
g[color[i]].push_back(color[x]);
//千万不能写成g[i].push_back(x);坑死我了
}
}
}
//这一部分代码还可以适当优化...(想一想)
- 对于点权与酒馆的一些处理(方便求最短路)
for(int i=1;i<=n;i++)
{
ww[color[i]]+=w[i];//将这个强连通分量中所有的点权全部加起来
if(tf[i])tf[color[i]]=1;//只要这个强连通分量中有一个结点有酒馆,那么就设定为有酒馆
}
- 求最短路模板(然而实际是最长路)
//模板不做解释
inline void spfa()
{
d[color[s]]=ww[color[s]];
queue<int>q;
q.push(color[s]);
while(q.size())
{
int now=q.front();
q.pop();
for(int i=0;i<g[now].size();i++)
{
int x=g[now][i];
if(d[now]+ww[x]>d[x])
{
d[x]=d[now]+ww[x];
q.push(x);
}
}
}
}
所以...
在所有有酒馆的节点中选一个最大值输出就好了
for(int i=1;i<=sum;i++)
{
//cout<<d[i]<<" ";
if(tf[i])ans=max(ans,d[i]);
}
完整代码
#include<bits/stdc++.h>
using namespace std;
const int MAXN=500000+10;
int n,m;
vector<int>ver[MAXN];
vector<int>g[MAXN];
int w[MAXN],ww[MAXN];
bool tf[MAXN];
int s,p,ans=0;
int dfn[MAXN],color[MAXN],low[MAXN];
int deep,sum;
bool vis[MAXN];
int d[MAXN];
stack<int>st;
inline int read()
{
int tot=0;
char c=getchar();
while(c<'0'||c>'9')
c=getchar();
while(c>='0'&&c<='9')
{
tot=tot*10+c-'0';
c=getchar();
}
return tot;
}
inline void tarjan(int now)
{
dfn[now]=++deep;
low[now]=deep;
vis[now]=1;
st.push(now);
for(int i=0;i<ver[now].size();i++)
{
int x=ver[now][i];
if(!dfn[x])
{
tarjan(x);
low[now]=min(low[now],low[x]);
}
else
{
if(vis[x])
low[now]=min(low[now],low[x]);
}
}
if(dfn[now]==low[now])
{
color[now]=++sum;
vis[now]=0;
while(st.top()!=now)
{
color[st.top()]=sum;
vis[st.top()]=0;
st.pop();
}
st.pop();
}
}
inline void spfa()
{
d[color[s]]=ww[color[s]];
queue<int>q;
q.push(color[s]);
while(q.size())
{
int now=q.front();
q.pop();
for(int i=0;i<g[now].size();i++)
{
int x=g[now][i];
if(d[now]+ww[x]>d[x])
{
d[x]=d[now]+ww[x];
q.push(x);
}
}
}
}
int main()
{
//freopen("testdata.in","r",stdin);
n=read();m=read();
for(int i=1;i<=m;i++)
{
int x=read(),y=read();
ver[x].push_back(y);
}
for(int i=1;i<=n;i++)
w[i]=read();
s=read();p=read();
for(int i=1;i<=p;i++)
{
int x=read();
tf[x]=1;
}
for(int i=1;i<=n;i++)
{
if(!dfn[i])
tarjan(i);
}
/*cout<<endl;
for(int i=1;i<=n;i++)
cout<<color[i]<<" ";
cout<<endl;*/
for(int i=1;i<=n;i++)
{
for(int j=0;j<ver[i].size();j++)
{
int x=ver[i][j];
if(color[i]!=color[x])
{
//cout<<i<<" "<<x<<" "<<color[i]<<" "<<color[x]<<endl;
g[color[i]].push_back(color[x]);
}
}
}
for(int i=1;i<=n;i++)
{
ww[color[i]]+=w[i];
if(tf[i])tf[color[i]]=1;
}
/*cout<<color[s]<<endl;
cout<<endl;
for(int i=1;i<=sum;i++)
{
for(int j=0;j<g[i].size();j++)cout<<i<<" "<<g[i][j]<<"\n";
}
cout<<endl;
for(int i=1;i<=sum;i++)
cout<<ww[i]<<" ";cout<<endl;
for(int i=1;i<=sum;i++)
cout<<tf[i]<<" ";cout<<endl;
cout<<ww[color[s]]<<endl;*/
spfa();
for(int i=1;i<=sum;i++)
{
//cout<<d[i]<<" ";
if(tf[i])ans=max(ans,d[i]);
}
//cout<<endl;
cout<<ans<<endl;
return 0;
}
洛谷 题解 P3627 【[APIO2009]抢掠计划】的更多相关文章
- P3627 [APIO2009]抢掠计划
P3627 [APIO2009]抢掠计划 Tarjan缩点+最短(最长)路 显然的缩点...... 在缩点时,顺便维护每个强连通分量的总权值 缩完点按照惯例建个新图 然后跑一遍spfa最长路,枚举每个 ...
- 【题解】洛谷P3627 [APIO2009]抢掠计划(缩点+SPFA)
洛谷P3627:https://www.luogu.org/problemnew/show/P3627 思路 由于有强连通分量 所以我们可以想到先把整个图缩点 缩点完之后再建一次图 把点权改为边权 并 ...
- 洛谷 P3627 [APIO2009]抢掠计划 题解
Analysis 建图+强连通分量+SPFA求最长路 但要保证最后到达的点中包含酒馆 虽然思路并不难想,但要求的代码能力很高. #include<iostream> #include< ...
- 洛谷 P3627 [APIO2009]抢掠计划 Tarjan缩点+Spfa求最长路
题目地址:https://www.luogu.com.cn/problem/P3627 第一次寒假训练的结测题,思路本身不难,但对于我这个码力蒟蒻来说实现难度不小-考试时肛了将近两个半小时才刚肛出来. ...
- [洛谷P3627][APIO2009]抢掠计划
题目大意:给你一张$n(n\leqslant5\times10^5)$个点$m(m\leqslant5\times10^5)$条边的有向图,有点权,给你起点和一些可能的终点.问从起点开始,到任意一个终 ...
- 洛谷 P3627 [APIO2009]抢掠计划
这题一看就是缩点,但是缩完点怎么办呢?首先我们把所有的包含酒吧的缩点找出来,打上标记,然后建立一张新图, 每个缩点上的点权就是他所包含的所有点的点权和.但是建图的时候要注意,每一对缩点之间可能有多条边 ...
- 【luogu P3627 [APIO2009]抢掠计划】 题解
题目链接:https://www.luogu.org/problemnew/show/P3627 把点权转化到边权上去. #include <stack> #include <que ...
- 题解 P3627 【[APIO2009]抢掠计划】
咕了四个小时整整一晚上 P3627 [APIO2009] 抢掠计划(https://www.luogu.org/problemnew/show/P3627) 不难看出答案即为该有向图的最长链长度(允许 ...
- APIO2009 抢掠计划 Tarjan DAG-DP
APIO2009 抢掠计划 Tarjan spfa/DAG-DP 题面 一道\(Tarjan\)缩点水题.因为可以反复经过节点,所以把一个联通快中的所有路口看做一个整体,缩点后直接跑\(spfa\)或 ...
随机推荐
- Jquery检验输入值
1.检验邮件 function chkEmail(strEmail) { if (!/^\w+([-+.]\w+)*@\w+([-.]\w+)*\.\w+([-.]\w ...
- NetworkX系列教程(9)-线性代数相关
小书匠 Graph 图论 学过线性代数的都了解矩阵,在矩阵上的文章可做的很多,什么特征矩阵,单位矩阵等.grpah存储可以使用矩阵,比如graph的邻接矩阵,权重矩阵等,这节主要是在等到graph后 ...
- (5)打鸡儿教你Vue.js
条件与循环 条件判断使用 v-if 指令 <p v-if="seen"> <template v-if="ok"> <script ...
- java23种设计模式等等。。
23种设计模式http://www.cnblogs.com/maowang1991/archive/2013/04/15/3023236.html 提升Java代码性能和安全性https://blog ...
- Asp.Net.Mvc+MEF+EF 项目 ,源码在GitHub ..希望对大家有所帮助
自己开源的 一个 Asp.Net + MEF+ EF 的 项目 . 供大家学习和使用, 点击进入GitHub
- 小程序 image跟view标签上下会有间隙
图片文字等inline元素默许是跟父级元素的baseline对齐,而baseline又和父级底边有必定间距 我是使用: 加上这个消除了间隙,如果没有解决,你可以分别用 vertical-align:t ...
- 2018-2019-2 20165234 《网络对抗技术》 Exp9 Web安全基础
Exp9 Web安全基础 一. 实践内容 1. 安装JDK.Webgoat 2. SQL注入攻击 数字型注入(Numeric SQL Injection) 日志欺骗(Log Spoofing) 字符串 ...
- docker启动elasticsearch失败--jvm内存不足解决方案
centos下载完elasticsearch并修改完配置后运行docker命令: docker run –name es1 -p : -p : -d -v /docker/es/esmaster/es ...
- hibernate的各种查询
Hibernate Query Language(HQL)Criteria QueryNative SQL下面对其分别进行解释select子句:有时并不需要取得对象的所有属性,这时可以使用select ...
- 临界区代码 critical section Locks and critical sections in multiple threads
临界区 在同步的程序设计中,临界区段(Critical section)指的是一个访问共享资源(例如:共享设备或是共享存储器)的程序片段,而这些共享资源有无法同时被多个线程访问的特性. 当有线程进入临 ...