Luogu P3387

强连通分量的定义如下:

有向图强连通分量:在有向图G中,如果两个顶点vi,vj间(vi>vj)有一条从vi到vj的有向路径,同时还有一条从vj到vi的有向路径,则称两个顶点强连通(strongly connected)。如果有向图G的每两个顶点都强连通,称G是一个强连通图。有向图的极大强连通子图,称为强连通分量(strongly connected components)。

来源于百度百科

我本人的理解:有向图内的一个不能再拓展得更大的强连通子图叫做这个有向图的一个强连通分量(也可以说是一个环)

注意:单独的一个孤立的点也会是一个强连通分量

求出强连通分量以后有什么用呢?

很显然,我们可以把整个强连通分量作为单独的一个点,权值按照题目要求取(在这题就是取所有点权的总和),这样就可以让这一个有向有环图转化成一个有向无环图。

Tarjan算法

Tarjan算法(在这里指Tarjan对于强连通分量提出的算法)就是用于求出一个有向图内的所有强连通分量的有效算法。

基本思想就是利用DFS往下搜索,标记顺序,如果找到返回祖先的一条边,则说明会构成一个环(强连通分量)。

这里要引入几个Tarjan算法必备的数组

数组名 作用
dfn[i] 用于记录节点i的dfs序
stk[i] 一个栈,用于记录当前搜索的这一条链上的节点
low[i] 用于记录节点i能访问到的节点中最小的dfs序 ,也就是最上层的祖先

关键点:如果dfn[i]==low[i],意味着在节点i的子树中没有任何的节点可以访问到节点i的祖先,说明节点i与仍然在栈内子节点(必要条件)构成了一个强连通分量

不在栈内的子节点无法与节点i构成强连通分量,原因是不在栈内则说明它本身已经作为一个强连通分量被弹出栈了。

void tarjan(int now)
{
dfn[now]=++tim;//记录dfs序
low[now]=tim;//当前能访问到dfs序最小的点就是自己
stk[++cnt]=now;
vis[now]=true;
for (int i=head[now];i;i=e[i].nxt)
{
int to=e[i].to;
if (!dfn[to])
{
tarjan(to);
low[now]=min(low[now],low[to]);
//如果该点没被遍历过,那么就进行遍历。
}
else
{
if (vis[to]) low[now]=min(low[now],dfn[to]);
//必须判断是否在栈中。只有在同时在栈内的点才有可能构成强连通分量。
}
}
if (low[now]==dfn[now])
{
tot++;//强连通分量的编号
while (stk[cnt]!=now)
{
scc[stk[cnt]]=tot;
val[tot]+=a[stk[cnt]];
vis[stk[cnt]]=false;
cnt--;
}
scc[stk[cnt]]=tot;
val[tot]+=a[stk[cnt]];
vis[stk[cnt]]=false;
cnt--;
//将栈中比u后进入的点和u本身出栈,这些点构成一个强联通分量,打上标记
}
}

结合代码进行理解。

拓扑排序和缩点操作

对一个有向无环图(Directed Acyclic Graph简称DAG)G进行拓扑排序,是将G中所有顶点排成一个线性序列,使得图中任意一对顶点u和v,若边<u,v>∈E(G),则u在线性序列中出现在v之前。通常,这样的线性序列称为满足拓扑次序(Topological Order)的序列,简称拓扑序列。简单的说,由某个集合上的一个偏序得到该集合上的一个全序,这个操作称之为拓扑排序。

来源于百度百科

我个人的理解:对有向无环图中所有节点排序形成一个合法的访问次序。

做一个比喻:吃饭之前要端盘子,端盘子之前要炒菜——那么对于这三件事的拓扑次序就是炒菜→端盘子→吃饭。

那么具体应该如何处理呢?

事实上有两种实现方法,但是我个人暂时只会一种。

利用队列的方式,把所有入度为0的点入队,然后把这几个点对其他点入度的贡献删除,然后再把入度为0的点入队,直到排序完成为止。

完成拓扑排序后就可以利用拓扑次序进行动态规划了。

完整代码如下:

#include<cstdio>
#include<queue>
using namespace std;
queue<int> que;
struct data
{
int sta,to,nxt;
}e[500005],ine[500005],oute[500005];
int dfn[100005],tim,low[100005],stk[100005],cnt,cnti,cnto,head[100005],inhead[100005],outhead[100005];
bool vis[10005];
int tot,scc[100005],val[100005],a[100005],n,in[100005],order[100005],m,u,v,f[100005],ans;
void tarjan(int now)
{
dfn[now]=++tim;//记录dfs序
low[now]=tim;//当前能访问到dfs序最小的点就是自己
stk[++cnt]=now;
vis[now]=true;
for (int i=head[now];i;i=e[i].nxt)
{
int to=e[i].to;
if (!dfn[to])
{
tarjan(to);
low[now]=min(low[now],low[to]);
//如果该点没被遍历过,那么就进行遍历。
}
else
{
if (vis[to]) low[now]=min(low[now],dfn[to]);
//必须判断是否在栈中。只有在同时在栈内的点才有可能构成强连通分量。
}
}
if (low[now]==dfn[now])
{
tot++;//强连通分量的编号
while (stk[cnt]!=now)
{
scc[stk[cnt]]=tot;
val[tot]+=a[stk[cnt]];
vis[stk[cnt]]=false;
cnt--;
}
scc[stk[cnt]]=tot;
val[tot]+=a[stk[cnt]];
vis[stk[cnt]]=false;
cnt--;
//将栈中比u后进入的点和u本身出栈,这些点构成一个强联通分量,打上标记
}
}
void topo()//拓扑排序
{
cnt=0,cnti=0,cnto=0;
for (int i=1;i<=n;i++)
{
for (int j=head[i];j;j=e[j].nxt)
{
if (scc[i]!=scc[e[j].to])
{
oute[++cnto].to=scc[e[j].to];
oute[cnto].nxt=outhead[scc[i]];
outhead[scc[i]]=cnto;
in[scc[e[j].to]]++;
ine[++cnti].sta=scc[i];
ine[cnti].nxt=inhead[scc[e[j].to]];
inhead[scc[e[j].to]]=cnti;
//out前缀的变量是出边的记录
//in前缀的变量是入边的记录,使用了一种另类的链式前向星
}
}
}
for (int i=1;i<=tot;i++)
if (in[i]==0) que.push(i);//入度为零则入队
cnt=0;
while (!que.empty())
{
int u=que.front();
que.pop();
order[++cnt]=u;//记录顺序
for (int i=outhead[u];i;i=oute[i].nxt)
{
int v=oute[i].to;
in[v]--;
if (in[v]==0) que.push(v);
}
}
}
int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++)
scanf("%d",&a[i]);
for (int i=1;i<=m;i++)
{
scanf("%d%d",&u,&v);
e[i].to=v;
e[i].nxt=head[u];
head[u]=i;
//链式前向星存原图
}
tim=0;
for (int i=1;i<=n;i++) if (!dfn[i]) tarjan(i);
//如果没有被遍历过的点要继续遍历。
topo();
for (int i=1;i<=tot;i++)
{
f[order[i]]=val[order[i]];
for (int j=inhead[order[i]];j;j=ine[j].nxt)
f[order[i]]=max(f[order[i]],f[ine[j].sta]+val[order[i]]);
//很容易的一个动态规划
}
for (int i=1;i<=tot;i++) ans=max(f[i],ans);//统计答案
printf("%d",ans);
return 0;
}

【Luogu P3387】缩点模板(强连通分量Tarjan&拓扑排序)的更多相关文章

  1. 【BZOJ2330】糖果(差分约束系统,强连通分量,拓扑排序)

    题意: 幼儿园里有N个小朋友,lxhgww老师现在想要给这些小朋友们分配糖果,要求每个小朋友都要分到糖果.但是小朋友们也有嫉妒心,总是会提出一些要求,比如小明不希望小红分到的糖果比他的多,于是在分配糖 ...

  2. 强连通分量(tarjan求强连通分量)

    双DFS方法就是正dfs扫一遍,然后将边反向dfs扫一遍.<挑战程序设计>上有说明. 双dfs代码: #include <iostream> #include <cstd ...

  3. 【bzoj1093】[ZJOI2007]最大半连通子图 Tarjan+拓扑排序+dp

    题目描述 一个有向图G=(V,E)称为半连通的(Semi-Connected),如果满足:对于u,v∈V,满足u→v或v→u,即对于图中任意两点u,v,存在一条u到v的有向路径或者从v到u的有向路径. ...

  4. 【bzoj5017】[Snoi2017]炸弹 线段树优化建图+Tarjan+拓扑排序

    题目描述 在一条直线上有 N 个炸弹,每个炸弹的坐标是 Xi,爆炸半径是 Ri,当一个炸弹爆炸时,如果另一个炸弹所在位置 Xj 满足:  Xi−Ri≤Xj≤Xi+Ri,那么,该炸弹也会被引爆.  现在 ...

  5. 强连通分量(Tarjan)模板

    贴模板,备忘. 模板1: #include<iostream> #include<cstring> #include<cmath> #include<cstd ...

  6. 强连通分量tarjan缩点——POJ2186 Popular Cows

    这里的Tarjan是基于DFS,用于求有向图的强联通分量. 运用了一个点dfn时间戳和low的关系巧妙地判断出一个强联通分量,从而实现一次DFS即可求出所有的强联通分量. §有向图中, u可达v不一定 ...

  7. POJ 1236 Network of Schools(强连通分量/Tarjan缩点)

    传送门 Description A number of schools are connected to a computer network. Agreements have been develo ...

  8. 强连通分量Tarjan模板

    #include<iostream> #include<stdio.h> #include<string.h> #include<stack> #inc ...

  9. POJ1236_A - Network of Schools _强连通分量::Tarjan算法

    Time Limit: 1000MS   Memory Limit: 10000K Description A number of schools are connected to a compute ...

随机推荐

  1. 一个html,3D 标签 鼓励自己

    效果如图: <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w ...

  2. 一般链表实现集合运算(C语言)

    最近在学习数据结构,遇到以下问题: 假设集合A = (c, b, e, g, f, d),B = (a, b, n, f),利用一般线性链表实现集合运算(A-B)∪(B-A). 分析: 上面的问题只要 ...

  3. ubuntu12.04 添加程序启动器

    方法: 1. 在/usr/share/applications目录下创建eclipse启动器配置文件 cd /usr/share/applications vim eclipse.desktop 2. ...

  4. (JavaScript) 时间转为几天前、几小时前、几分钟前

    // 时间戳转多少分钟之前 getDateDiff(dateTimeStamp) { // 时间字符串转时间戳 var timestamp = new Date(dateTimeStamp).getT ...

  5. 发布兼容TS的JS库到nexus和npmjs

    一. 前言 由于node以及绝大多数前端库都是用JavaScript(以下简称JS)语言实现,而Angular是用TypeScript(以下简称TS)实现,虽然TS是JS的超集,但是由于TS和JS对于 ...

  6. NOIP模拟 39

    考的嘛也不是. 伤心(怎么可能) T1稍想想组合数,然后牢记: 取模题随时取模,包括刚刚读入的数据  T2想到了基环树,然而不会打QAQ.. 非常简洁但非常大神的做法:随便断掉环上的一条边 利用“这条 ...

  7. 基于xposed逆向微信、支付宝、云闪付来实现个人免签支付功能

    我的个人网站如何实现支付功能? 想必很多程序员都有过想开发一个自己的网站来获得一些额外的收入,但做这件事会遇到支付这个问题.目前个人网站通过常规手法是无法实现支付管理的,所有支付渠道都需要以公司的身份 ...

  8. Maven多模块项目介绍和搭建

    http://www.open-open.com/lib/view/open1418263515855.html

  9. go xml 序列化

    /** 序列化 */ package main import ( "encoding/xml" "fmt" ) // Person 结构 type Person ...

  10. 更新linux时候提示“由于没有公钥,无法验证下列签名".

    本文链接:https://blog.csdn.net/loovejava/article/details/21837935 新安装的Ubuntu在使用sudo apt-get update更新源码的时 ...