洛谷P3387 【模板】缩点 题解
背景
今天\(loj\)挂了,于是就有了闲情雅致来刷\(luogu\)
题面
洛谷P3387 【模板】缩点传送门
题意
给定一个\(n\)个点\(m\)条边有向图,每个点有一个权值,求一条路径,使路径经过的点权值之和最大。你只需要求出这个权值和。
允许多次经过一条边或者一个点,但是,重复经过的点,权值只计算一次。
思路
~~点击获得更好的阅读体验~~
\(Kosaraju\)缩点+记忆化搜索
什么是\(Kosaraju\)缩点?
与\(Tarjan\)缩点的时间复杂度差不多,都是\(O(n+m)\)。
其主要思想就是\(dfs\)两次。
核心代码如下:
void dfs_1(int x){
vis[x]=1;
for(int i=fir[x];i;i=nxt[i]){
if(vis[son[i]]==0) dfs_1(son[i]);
}
d[++t]=x;
}
void dfs_2(int x){
vis[x]=t;
s[t]++;
for(int i=fir2[x];i;i=nxt2[i]){
if(vis[son2[i]]==0) dfs_2(son2[i]);
}
}
void Kosaraju(){
t=0;
memset(vis,0,sizeof(vis));
for(int i=1;i<=n;i++){
if(vis[i]==0) dfs_1(i);
}
memset(vis,0,sizeof(vis));t=0;
for(int i=n;i>=1;i--){
if(vis[d[i]]==0) t++,dfs_2(d[i]);
}
}
当然,我认为\(Kosaraju\)有一定的缺点,比如说要建反边。
具体\(Kosaraju\)算法的介绍戳这里
什么是记忆化搜索
记忆化搜索(Memory search)心理学是指搜索信息的流程,但是搜索到的一些解用动态规划的那种思想和模式作一些保存。
——百度百科
其实说白了就是\(dp\)的进阶。多了一个搜索。。。
只不过在\(dfs\)时开了个数组存储了当前的状态,以便以后搜索需要。
关于这道题
那么这道题就愉快的解决了。。。
思路如下:
为什么有这丑陋无比的水印
CODE:(PS:请忽略边权)
#include<algorithm>
#include<bitset>
#include<complex>
#include<deque>
#include<exception>
#include<fstream>
#include<functional>
#include<iomanip>
#include<ios>
#include<iosfwd>
#include<iostream>
#include<istream>
#include<iterator>
#include<limits>
#include<list>
#include<locale>
#include<map>
#include<memory>
#include<new>
#include<numeric>
#include<ostream>
#include<queue>
#include<set>
#include<sstream>
#include<stack>
#include<stdexcept>
#include<streambuf>
#include<string>
#include<typeinfo>
#include<utility>
#include<valarray>
#include<vector>
#include<cctype>
#include<cerrno>
#include<cfloat>
#include<ciso646>
#include<climits>
#include<clocale>
#include<cmath>
#include<csetjmp>
#include<csignal>
#include<cstdarg>
#include<cstddef>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#define E 200010
using namespace std;//丑陋无比的头文件
inline int read(){
int res=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=getchar();}
while(ch>='0'&&ch<='9') res=res*10+ch-'0',ch=getchar();
return res*f;
}//丑陋无比的读优
inline void write(int x){
if(x<0) putchar('-'),x=-x;
if(x<10) putchar(x+'0');
else{
write(x/10);
putchar(x%10+'0');
}
}//丑陋无比的输优
int n,m,a[E];
int fir[E],nxt[E],w[E],son[E],tot;
int vis[E],d[E],top,f[E],ans;
int fir2[E],nxt2[E],w2[E],son2[E],tot2;
void add(int x,int y,int z){
++tot;
nxt[tot]=fir[x];
fir[x]=tot;
son[tot]=y;
w[tot]=z;
}//建正边
void add2(int x,int y,int z){
++tot2;
nxt2[tot2]=fir2[x];
fir2[x]=tot2;
son2[tot2]=y;
w2[tot2]=z;
}//建反边
void dfs1(int x){//Kosaraju算法第一次dfs
vis[x]=1;
for(int i=fir[x];i;i=nxt[i]){//正向
int to=son[i];
if(!vis[to]) dfs1(to);
}
d[++top]=x;
}
void dfs2(int x){//Kosaraju算法第二次dfs
vis[x]=top;//标记第几个强连通分量
for(int i=fir2[x];i;i=nxt2[i]){//反向
int to=son2[i];
if(!vis[to]) dfs2(to);
}
}
void dp(int x){//记忆化搜索
if(f[x]!=0) return ;//记忆化体现
int ss=0;
f[x]=d[x];
for(int i=fir2[x];i;i=nxt2[i]){//使用新的DAG
int to=son2[i];
dp(to);
ss=max(ss,f[to]);//取最大值
}
f[x]+=ss;//别忘记加回去
}
int main(){
n=read();m=read();//读入
for(int i=1;i<=n;i++) a[i]=read();//读入
for(int i=1;i<=m;i++){
int x=read(),y=read();
add(x,y,1);//建正边
add2(y,x,1);//建反边
}
//Kosaraju Start
memset(vis,0,sizeof(vis));
memset(d,0,sizeof(d));top=0;//清空
for(int i=1;i<=n;i++){
if(!vis[i]) dfs1(i);//dfs第一次(正向)
}
memset(vis,0,sizeof(vis));top=0;//清空
for(int i=n;i>=1;i--){
if(!vis[d[i]]) ++top,dfs2(d[i]);//dfs第二次(反向)
}
//Kosaraju End
memset(d,0,sizeof(d));//清空
for(int i=1;i<=n;i++){
d[vis[i]]+=a[i];//权值累计
}
memset(fir2,0,sizeof(fir2));tot2=0;//清空
for(int i=1;i<=n;i++){//重新建DAG
for(int j=fir[i];j;j=nxt[j]){
int to=son[j];
if(vis[i]!=vis[to]) add2(vis[i],vis[to],1);//需要两个不同的强连通分量
}
}
for(int i=1;i<=top;i++){//记忆化搜索
if(!f[i]){//记忆化
dp(i);
ans=max(ans,f[i]);//取最大值
}
}
write(ans);putchar('\n');//输出
return 0;//结束
}
洛谷P3387 【模板】缩点 题解的更多相关文章
- 洛谷 P3387 【模板】缩点 DAGdp学习记
我们以洛谷P3387 [模板]缩点 来学习DAGdp 1.这道题的流程 //伪代码 for i->n if(i未被遍历) tarjan(i) 缩点() DAGdp() 完成 首先tarjan这部 ...
- tarjan缩点练习 洛谷P3387 【模板】缩点+poj 2186 Popular Cows
缩点练习 洛谷 P3387 [模板]缩点 缩点 解题思路: 都说是模板了...先缩点把有环图转换成DAG 然后拓扑排序即可 #include <bits/stdc++.h> using n ...
- 洛谷P1783 海滩防御 分析+题解代码
洛谷P1783 海滩防御 分析+题解代码 题目描述: WLP同学最近迷上了一款网络联机对战游戏(终于知道为毛JOHNKRAM每天刷洛谷效率那么低了),但是他却为了这个游戏很苦恼,因为他在海边的造船厂和 ...
- 洛谷P3373 [模板]线段树 2(区间增减.乘 区间求和)
To 洛谷.3373 [模板]线段树2 题目描述 如题,已知一个数列,你需要进行下面两种操作: 1.将某区间每一个数加上x 2.将某区间每一个数乘上x 3.求出某区间每一个数的和 输入输出格式 输入格 ...
- 洛谷P4047 [JSOI2010]部落划分题解
洛谷P4047 [JSOI2010]部落划分题解 题目描述 聪聪研究发现,荒岛野人总是过着群居的生活,但是,并不是整个荒岛上的所有野人都属于同一个部落,野人们总是拉帮结派形成属于自己的部落,不同的部落 ...
- 洛谷P1155 双栈排序题解(图论模型转换+二分图染色+栈)
洛谷P1155 双栈排序题解(图论模型转换+二分图染色+栈) 标签:题解 阅读体验:https://zybuluo.com/Junlier/note/1311990 原题地址:洛谷P1155 双栈排序 ...
- 【AC自动机】洛谷三道模板题
[题目链接] https://www.luogu.org/problem/P3808 [题意] 给定n个模式串和1个文本串,求有多少个模式串在文本串里出现过. [题解] 不再介绍基础知识了,就是裸的模 ...
- 洛谷10月月赛II题解
[咻咻咻] (https://www.luogu.org/contestnew/show/11616) 令人窒息的洛谷月赛,即将参加NOIp的我竟然只会一道题(也可以说一道也不会),最终145的我只能 ...
- 洛谷P3375 [模板]KMP字符串匹配
To 洛谷.3375 KMP字符串匹配 题目描述 如题,给出两个字符串s1和s2,其中s2为s1的子串,求出s2在s1中所有出现的位置. 为了减少骗分的情况,接下来还要输出子串的前缀数组next.如果 ...
随机推荐
- TCP协议 连接三次握手
TCP(Transmission Control Protocol) 传输控制协议 TCP是主机对主机层的传输控制协议,提供可靠的连接服务,采用三次握手确认建立一个连接: 位码即tcp标志位,有6种标 ...
- Vasya and Multisets CodeForces - 1051C(英语限制了我的想象力)
题意: 对输入的每个数都进行分配到两个多重集中,问分配完以后 每个多重集中的个数为1的元素是否相等 相等则输出一个分配方式, 不想等为NO 解析: 三种情况 1.原数列中个数为1 的数的个数为偶数 ...
- Ryuji doesn't want to study 2018 徐州赛区网络预赛
题意: 1.区间求 a[l]×L+a[l+1]×(L−1)+⋯+a[r−1]×2+a[r](L is the length of [ l, r ] that equals to r - l + 1) ...
- emWin notes
emwin 学习小记 @2018-7-5 ## 在使用 graph 控件时,需要在坐标上显示波形,波形刷新函数 GRAPH_DATA_YT_AddValue() 放在一个 while(1) 中,这样便 ...
- 第一、介绍Canvas
canvas能做什么? canvas是HTML5中的新元素,你可以使用javascript用它来绘制图形.图标.以及其它任何视觉性图像.它也可用于创建图片特效和动画.如果你掌握了完整的命令,你可以用c ...
- AngularJs -- 模 块
在JavaScript中,将函数代码全部定义在全局命名空间中绝对不是什么好主意,这样做会导致冲突从而是调试变得非常困难,浪费宝贵的时间. 上一章介绍数据绑定时,就是写在全局命名空间中定义的函数. 在A ...
- CSS3 响应式布局: @media (min/max-width:***) @font-face
响应式布局 responsive design @media 属性 bootstrap css 分析: @media (min-width:768px){ body{***} } use @medi ...
- 使用js获取浏览器地址栏里的参数
用JS获取地址栏参数的方法(超级简单) 方法一:采用正则表达式获取地址栏参数:( 强烈推荐,既实用又方便!) function GetQueryString(name) { var reg = new ...
- 第11月第23天 markedTextRange 崩溃
1. 在对textView.textField限制文字长度时,如果不做特殊处理,当联想文字加上已输入文字超出设定长度时,iOS 7.0系统会崩溃(ios 8.0以上系统做了处理,不会崩溃). http ...
- 第11月第11天 avplayer循环播放
1. /* Setting actionAtItemEnd to None prevents the movie from getting paused at item end. A very sim ...