SCC统计
Kosoraju
SCC总数及记录SCC所需要的最少边情况
#include<cstdio>
const int N = ;
int scc_cnt = ;
int Case, n, m, i, x, y, cnt, use[N], g[][N], nxt[][N], v[][N], ed, q[N], t, vis[N], e[N][];
inline void add(int x, int y)
{
v[][++ed] = y;
nxt[][ed] = g[][x];
g[][x] = ed;
v[][ed] = x;
nxt[][ed] = g[][y];
g[][y] = ed;
}
void dfs1(int x)
{
vis[x] = ;
for (int i = g[][x]; i; i = nxt[][i])
if (!vis[v[][i]])
{
use[i] = , dfs1(v[][i]);
}
q[++t] = x;
}
void dfs2(int x)
{
vis[x] = ;
for (int i = g[][x]; i; i = nxt[][i])
if (vis[v[][i]])
{
use[i] = , dfs2(v[][i]);
}
}
int main()
{
scanf("%d", &Case);
while (Case--)
{
scc_cnt = ;
scanf("%d%d", &n, &m);
for (ed = , i = ; i <= n; i++)
{
vis[i] = g[][i] = g[][i] = ;
}
for (i = ; i <= m; i++)
{
scanf("%d%d", &x, &y);
e[i][] = x, e[i][] = y;
add(x, y);
use[i] = ;
}
for (t = , i = ; i <= n; i++)
if (!vis[i])
{
dfs1(i);
}
for (i = n; i; i--)
if (vis[q[i]])
{
scc_cnt++;
dfs2(q[i]);
}
printf("%d\n",scc_cnt);
// cnt = n * 2;
// for (i = 1; i <= m; i++)if (use[i])
// {
// cnt--;
// }
// for (i = 1; i <= m; i++)if (!use[i] && cnt > 0)
// {
// use[i] = 1, cnt--;
// }
// for (i = 1; i <= m; i++)if (!use[i])
// {
// printf("%d %d\n", e[i][0], e[i][1]);
// }
}
}
Tarjan
转自：
https://www.cnblogs.com/stxy-ferryman/p/7779347.html
https://blog.csdn.net/sentimental_dog/article/details/53790582
用处：
1.有向图强连通分量 (一个有向图是强连通的当且仅当G中有一个回路，它至少包含每个节点一次)
2.无向图双连通分量(割点,桥)
3.离线LCA
割顶：若去掉一个点和与这个点相连的边后，图不再连通，则这个点是割顶。
求法：若节点uu存在一棵子树vv满足vv中所有节点的回边都指向uu及以下的节点（即low[v]≥pre[u]low[v]≥pre[u]），则uu是割顶；所以一次DFS即可求出所有割顶。但注意一个特殊情况！根节点是割顶当且仅当它的子节点数严格大于1。
桥：去掉这条边，图不再连通。
求法：若节点uu存在一棵子树vv满足vv中所有节点的回边都指向uu以下（不含uu）的节点（即low[v]>pre[u]low[v]>pre[u]），则边(u,v)(u,v)是桥；所以一次DFS即可求出所有桥。
BCC（点-双连通分量）：连通分量中任意两点间都至少有两条“点不重复路径”。
求法：点双等价于不含割顶的极大子图。但一个问题是割顶本身可能属于多个点双。解决方法是把边分开，即保存一个边的栈。
eBCC（边-双连通分量）：连通分量中任意两点间都至少有两条“边不重复路径”。
求法：边双等价于不含桥的极大子图。一条点不可能分属多个边双（否则可以将这些个边双合并）。所以一次DFS求出所有桥，再来一次DFS把边双分开即可（DFS到桥就不走）。
前向星版SCC
void tarjan(int x)
{
dfn[x] = ++deep;
low[x] = deep;
visit[x] = ;
sta[++top] = x;
for (int i = Head[x]; i; i = nxt[i])
{
int v = to[i];
if (!dfn[v])
{
tarjan(v);
low[x] = min(low[x], low[v]);
}
else
{
if (visit[v])
{
low[x] = min(low[x], low[v]);
}
}
}
if (dfn[x] == low[x])
{
color[x] = ++colorsum;
visit[x] = ;
while (sta[top] != x)
{
color[sta[top]] = colorsum;
visit[sta[top--]] = ;
}
top--;
}
}
[USACO06JAN]The Cow Prom
算强连通分量中点数不小于1的个数
/*#include<cstring>#include<algorithm>#include<queue>#include<vector>#include<cstdio>#include<cmath>#include<iostream>*/
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = ;
vector<int> gra[maxn];
int dfn[maxn];//表示这个点在dfs的时候是第几个搜到的；
int low[maxn];//表示这个点及其子节点连的所有点里dfn的最小值
int sta[maxn];//存着当前所有可能能构成强连通分量的点
int visit[maxn];//表示一个点目前是否在sta中
int cnt[maxn];//各个强连通分量中含点的数目
int color[maxn];//表示一个点属于哪个强连通分量
int deep;/*表示从前有多少个点被搜到*/
int top;/*sta目前的大小*/
int colorsum = ;/*目前强连通分量的数目*/
int ans = ;
int n, m;
void tarjan(int x)
{
dfn[x] = ++deep;
low[x] = deep;
visit[x] = ;
sta[++top] = x;
int sz = gra[x].size();
for (int i = ; i < sz; i++)
{
int to = gra[x][i];
if (!dfn[to])
{
tarjan(to);
low[x] = min(low[x], low[to]);
}
else
{
if (visit[to])
{
low[x] = min(low[x], low[to]);
}
}
}
if (dfn[x] == low[x])
{
color[x] = ++colorsum;
visit[x] = ;
while (sta[top] != x)
{
color[sta[top]] = colorsum;
visit[sta[top--]] = ;
}
top--;
}
}
int main()
{
cin >> n >> m;
int from, to;
for (int i = ; i <= m; i++)
{
cin >> from >> to;
gra[from].push_back(to);
}
for (int i = ; i <= n; i++)
{
if (!dfn[i])
{
tarjan(i);
}
}
for (int i = ; i <= n; i++)
{
cnt[color[i]]++;
}
for (int i = ; i <= colorsum; i++)
{
if (cnt[i] > )
{
ans++;
}
}
cout << ans << endl;
}
poj2186 Popular Cows
缩点后求出度为0的点的个数(如果求得的点为缩点后的超级点需要特判)
/*#include<cstring>#include<algorithm>#include<queue>#include<vector>#include<cstdio>#include<cmath>#include<iostream>*/
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = ;
vector<int> gra[maxn];
int du[maxn];
int dfn[maxn];//表示这个点在dfs的时候是第几个搜到的；
int low[maxn];//表示这个点及其子节点连的所有点里dfn的最小值
int sta[maxn];//存着当前所有可能能构成强连通分量的点
int visit[maxn];//表示一个点目前是否在sta中
int cnt[maxn];//各个强连通分量中含点的数目
int color[maxn];//表示一个点属于哪个强连通分量
int deep;/*表示从前有多少个点被搜到*/
int top;/*sta目前的大小*/
int colorsum = ;/*目前强连通分量的数目*/
int ans = ;
int n, m;
void tarjan(int x)
{
dfn[x] = ++deep;
low[x] = deep;
visit[x] = ;
sta[++top] = x;
int sz = gra[x].size();
for (int i = ; i < sz; i++)
{
int to = gra[x][i];
if (!dfn[to])
{
tarjan(to);
low[x] = min(low[x], low[to]);
}
else
{
if (visit[to])
{
low[x] = min(low[x], low[to]);
}
}
}
if (dfn[x] == low[x])
{
color[x] = ++colorsum;
visit[x] = ;
while (sta[top] != x)
{
color[sta[top]] = colorsum;
visit[sta[top--]] = ;
}
top--;
}
}
int main()
{
int temp = ;
cin >> n >> m;
int from, to;
for (int i = ; i <= m; i++)
{
cin >> from >> to;
gra[from].push_back(to);
}
for (int i = ; i <= n; i++)
{
if (!dfn[i])
{
tarjan(i);
}
}
for (int i = ; i <= n; i++)
{
int len = gra[i].size();
for (int j = ; j < len; j++)
{
int to = gra[i][j];
if (color[to] != color[i]) //如果缩点后两个点不在一个强连通分量内
{
du[color[i]]++;
}
}
cnt[color[i]]++;
}
for (int i = ; i <= colorsum; i++)
{
if (temp > )
{
cout << << endl;
return ;
}
if (du[i] == )
{
temp++;
ans = cnt[i];
}
}
cout << ans << endl;
}
洛谷3388 求割点并输出模板题（附带桥输出）
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = ;
vector<int> gra[maxn];
int dfn[maxn];//表示这个点在dfs的时候是第几个搜到的；
int low[maxn];//表示这个点及其子节点连的所有点里dfn的最小值
int iscut[maxn];//表示这个点是否是割点
int deep;/*表示从前有多少个点被搜到*/
int ans = ;
int n, m;
int tarjan(int x, int pre)
{
int child = ;
int lowx;
lowx = dfn[x] = ++deep;
int len = gra[x].size();
for (int i = ; i < len; i++)
{
int to = gra[x][i];
if (!dfn[to])
{
child++;
int lowto = tarjan(to, x);
lowx = min(lowx, lowto);
if (lowto >= dfn[x])
{
iscut[x] = ;
}
// if (lowto > dfn[x])
// {
// cout << "bridge " << x << " " << to << endl;
// }
}
else
{
if (to != pre && dfn[to] < dfn[x])
{
lowx = min(lowx, dfn[to]);
}
}
}
if (pre < && child == )
{
iscut[x] = ;
}
low[x] = lowx;
return lowx;
}
int main()
{
int temp = ;
cin >> n >> m;
int from, to;
for (int i = ; i <= m; i++)
{
cin >> from >> to;
gra[from].push_back(to);
gra[to].push_back(from);
}
for (int i = ; i <= n; i++)
{
if (!dfn[i])
{
tarjan(i, -);
}
}
for (int i = ; i <= n; i++)
{
if (iscut[i])
{
ans++;
}
}
cout << ans << endl;
for (int i = ; i <= n; i++)
{
if (iscut[i])
{
cout << i << " ";
}
}
cout << endl;
}
Codeforces 962 F
求点双连通分量中的边
/*Huyyt*/
#include<bits/stdc++.h>
#define mem(a,b) memset(a,b,sizeof(a))
#define pb push_back
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int dir[][] = {{, }, {, }, {, -}, { -, }, {, }, {, -}, { -, -}, { -, }};
const int mod = 1e9 + ;
const int gakki = + + + + 1e9;
const int MAXN = 2e5 + ;
const int MAXM = 2e5 + ;
int to[MAXM << ], nxt[MAXM << ], Head[MAXN], tot = ;
inline void addedge(int u, int v)
{
to[++tot] = v;
nxt[tot] = Head[u];
Head[u] = tot;
}
int n, m;
int dfn[MAXN], low[MAXN], dfs_clock = ;
int BCCcnt = , blong[MAXN], inque[MAXM << ];
int st[MAXN], l = , ans[MAXN], ansnum = ;
bool vis[MAXM << ];
void tarjanBCC(int x, int fa)
{
dfn[x] = low[x] = ++dfs_clock;
for (int i = Head[x]; i; i = nxt[i])
{
int v = to[i];
if (v == fa || vis[i])
{
continue;
}
vis[i] = vis[i ^ ] = ;
st[l++] = i;
if (!dfn[v])
{
tarjanBCC(v, x);
low[x] = min(low[v], low[x]);
if (dfn[x] <= low[v])
{
int now, vnumber = , enumber = ;
BCCcnt++;
while()
{
now = st[--l];
if (blong[to[now]] != BCCcnt)
{
blong[to[now]] = BCCcnt, ++vnumber;
}
if (blong[to[now ^ ]] != BCCcnt)
{
blong[to[now ^ ]] = BCCcnt, ++vnumber;
}
inque[++enumber] = now;
if(now==i)
break;
}
if (vnumber == enumber)
{
for (int i = ; i <= enumber; i++)
{
ans[++ansnum] = inque[i] / ;
}
}
}
}
else
{
low[x] = min(low[x], dfn[v]);
}
}
}
int main()
{
ios_base::sync_with_stdio();
cin.tie();
cin >> n >> m;
int u, v;
for (int i = ; i <= m; i++)
{
cin >> u >> v;
addedge(u, v), addedge(v, u);
}
for (int i = ; i <= n; i++)
{
if (!dfn[i])
{
tarjanBCC(i, -);
}
}
sort(ans + , ans + + ansnum);
cout << ansnum << endl;
for (int i = ; i <= ansnum; i++)
{
cout << ans[i] << " ";
}
return ;
}
附:求边双连通分量中的边
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1e5+;
int dfn[N],low[N],H[N],nxt[N<<],to[N<<];
int n,m,x,y,cnt,tot=;
bool ib[N],is[N];
void add(int x,int y){
to[++tot]=y;nxt[tot]=H[x];H[x]=tot;
}
void dfs(int u,int fa){
low[u]=dfn[u]=++cnt;
int chi=;
for(int i=H[u];i;i=nxt[i]){
int v=to[i];
if(v==fa) continue;
if(!dfn[v]) {
++chi;
dfs(v,u);
low[u]=min(low[u],low[v]);
if(low[v]>dfn[u]) ib[i]=ib[i^]=;
if(low[v]>=dfn[u]) is[u]=;
}
else low[u]=min(low[u],dfn[v]);
}
if(chi==&&fa==-) is[u]=;
}
int num,bnum,a[N],A[N],ans;
void dfs2(int u,int fa){
++num;for(int i=H[u];i;i=nxt[i]) {
int v=to[i];
if(ib[i]||ib[i^]||v==fa) continue;
ib[i]=ib[i^]=;
a[++bnum]=i>>;
if(!dfn[v]) dfn[v]=,dfs2(v,u);
}
}
int main(){
scanf("%d%d",&n,&m);
for(int i=;i<=m;++i) {
scanf("%d%d",&x,&y);
add(x,y);
add(y,x);
}
for(int i=;i<=n;++i) if(!dfn[i]) dfs(i,-);
for(int i=;i<=n;++i) dfn[i]=;
for(int i=;i<=n;++i) if(!dfn[i]) {
num=bnum=;
dfn[i]=;
dfs2(i,-);
if(num==bnum) for(int j=;j<=num;++j) A[++ans]=a[j];
}
printf("%d\n",ans);
sort(A+,A+ans+);
for(int i=;i<=ans;++i) printf("%d ",A[i]);
}
SCC统计的更多相关文章
- CCF CSP 201709-4 通信网络
CCF计算机职业资格认证考试题解系列文章为meelo原创,请务必以链接形式注明本文地址 CCF CSP 201709-4 通信网络 问题描述 某国的军队由N个部门组成,为了提高安全性,部门之间建立了M ...
- POJ 2186 Popular cows(Kosaraju+强联通分量模板）
题目链接:http://poj.org/problem?id=2186 题目大意:给定N头牛和M个有序对(A,B),(A,B)表示A牛认为B牛是红人,该关系具有传递性,如果牛A认为牛B是红人,牛B认为 ...
- Python-基于向量机SVM的文本分类
项目代码见 Github: 1.算法介绍 2.代码所用数据 详情参见http://qwone.com/~jason/20Newsgroups/ 文件结构 ├─doc_classification.py ...
- BZOJ 5450 轰炸 (强连通缩点+DAG最长路)
<题目链接> 题目大意: 有n座城市,城市之间建立了m条有向的地下通道.你需要发起若干轮轰炸,每轮可以轰炸任意多个城市.但每次轰炸的城市中,不能存在两个不同的城市i,j满足可以通过地道从城 ...
- atitit.提升备份文件复制速度(1) -----分析统计问题and解决方案
atitit.) -----分析统计问题and解决方案 1. 现在的情形 1 2. 硬盘信息 大概50mb/s, 50iops 1 3. 统计小的文件比率 2 4. 复制速度估计.. 2 5. 小文件 ...
- bzoj 2707 [SDOI2012]走迷宫（SCC+高斯消元）
Description Morenan被困在了一个迷宫里.迷宫可以视为N个点M条边的有向图,其中Morenan处于起点S,迷宫的终点设为T.可惜的是,Morenan非常的脑小,他只会从一个点出发随机沿 ...
- HDU-2767-tarjan/Kosaraju求scc
http://acm.hdu.edu.cn/showproblem.php?pid=2767 问最少添加几条边使得图为强连通. tarjan跑一下,然后对强连通分量缩点,找下此时出度为零和入度为零的点 ...
- poj 3592 Instantaneous Transference 【SCC +缩点 + SPFA】
Instantaneous Transference Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 6204 Accep ...
- HDU 3639 SCC Hawk-and-Chicken
求SCC缩点,统计出每个SCC中的点的个数. 然后统计能到达u的最多的点的个数,可以反向建图,再dfs一遍统计出来. 最后说一下,有必要开一个标记数组,因为测试数据中有重边,结果无限WA. #incl ...
随机推荐
- 【flask】表单-上传文件
依赖: flask-wtf upload_case.html <!DOCTYPE html> <html lang="en"> <head> & ...
- "Developer tools access" 需控制另一个进程才能继续调试 解决方案
解决方案: 打开终端输入下边命令: DevToolsSecurity --status 查看状态 DevToolsSecurity --enable 输入密码,修改为enable,即可用 DevToo ...
- jdbc简单连接oracle数据库
package com.shangsheng; import java.sql.*; public class UserOracle { public static void main(String[ ...
- Centos7 修改系统时间和硬件时间不一致的问题
查看系统时间 [root@localhost ~]# dateSat Feb 24 14:41:22 CST 2018 查看硬件时间 [root@localhost ~]# hwclock --sho ...
- SSM的开发步骤分析
完整开发步骤 导包 spring的jar包 mybatis的jar包 mybatis-Spring的jar包 aop的依赖jar包 oracle等数据库连接的jar包 DataSource的jar包 ...
- 【计算机视觉】阶编码本模型(Multi phase codebook model)
转自:http://www.cnblogs.com/xrwang/archive/2012/04/24/MPCBBGM.html 多阶编码本模型(Multi phase codebook model) ...
- [Vuejs] 在vue各个组件中应用全局scss变量
需要安装一个插件:sass-resources-loader 1.执行安装命令: npm i sass-resources-loader --save-dev 2.修改vue-cli环境下build文 ...
- NFC读写器调试总结20191203
以下为NFC读写器调试经验总结: 1.读写器部分,从TX1/TX2输出的13.56MHZ信号主要由L0/C0构成低通滤波器,用于滤除13.56MHZ的高次谐波,取值L0=1UH,C0=47PF时候,谐 ...
- Go语言中byte类型和rune类型（五）
本篇内容本来准备在上一篇写的,想了想还是拆开写. go语言中字符串需要使用用双引号,而单引号用来表示单个的字符,字符也是组成字符串的元素.go语言的字符有两种: uint8类型,或者叫 byte 型, ...
- map的常见用法
map的常见用法 map 是什么? map是一组键值对的组合,通俗理解类似一种特殊的数组,a[key]=val,只不过数组元素的下标是任意一种类型,而且数组的元素的值也是任意一种类型.有点类似pyth ...