题单链接:
https://vjudge.net/article/752
推荐:1(强连通分量缩点),2(求割点),3(求桥边)
Tarjan 时间复杂度:O(n+e) n:点数,e:边数
割点学习 https://blog.csdn.net/weixin_44959460/article/details/89508621
1. POJ-1236 Network of Schools
Tarjan,缩点
缩完点后,
第一问就是 “入度为零的点” 的个数,
第二问就是 “入度为零的点的个数” 与 "出度为零的点的个数"的最大值,因为每次找一个出度为零的点 x,和一个入度为零的点 y,添加一条边 x->y,就各减少了一个点。总可以构造,使得最后变成一个强连通图。
参考这篇博文 https://xiaoxiaoh.blog.csdn.net/article/details/104284558
//POJ不支持 for (o: graph[X]) ? 所以遍历那里写的好长
#include <cstdio>
#include <vector>
#include <iostream>
using namespace std;
#define MAXN 105
int n;
vector<int> graph[MAXN];
int dfn[MAXN];
int low[MAXN];
int cnt; //dfn
int grp[MAXN];
int idx;
bool instack[MAXN];
int st[MAXN], top;
int in[MAXN], out[MAXN];
//求强连通分量
void Tarjan(int X)
{
low[X]=dfn[X]=++cnt;
st[++top]=X;
instack[X]=1;
for (int id=0,mxid=graph[X].size(),o; id<mxid && (o=graph[X][id]); id++) {
if (!dfn[o]) {
Tarjan(o);
low[X]=min(low[X],low[o]);
}
else {
if (instack[o])
low[X]=min(low[X],dfn[o]);
}
}
if (low[X]==dfn[X]) {
idx++;
int Y;
do {
Y=st[top--];
instack[Y]=0;
grp[Y]=idx;
} while (X!=Y);
}
}
void solve()
{
scanf("%d",&n);
for (int u=1,v; u<=n; u++) {
scanf("%d",&v);
while (v!=0) {
graph[u].push_back(v);
scanf("%d",&v);
}
}
for (int i=1; i<=n; i++)
if (!dfn[i])
Tarjan(i);
if (idx==1) { //原本就是强连通图
printf("1\n0\n");
return;
}
for (int u=0; u<=n; u++) {
for (int id=0,mxid=graph[u].size(),v; id<mxid && (v=graph[u][id]); id++)
if (grp[u]!=grp[v]) {
in[grp[v]]++;
out[grp[u]]++;
}
}
int cnt1=0,cnt2=0,ans1,ans2;
for (int i=1; i<=idx; i++) {
cnt1+= (in[i]==0);
cnt2+= (out[i]==0);
}
ans1=cnt1;
ans2=max(cnt1,cnt2);
printf("%d\n%d\n",ans1,ans2);
}
int main()
{
solve();
return 0;
}
2. UVA-315 Network
无向图 Tarjan 求割点。(和求强连通分量的Tarjan算法差不多,都利用了 dfn 和 low)。割点有两种情况:
1.是根节点,且儿子数量超过 1;
2.不是根节点,且dfn[X] >= low[o]
,其中 X 是割点, o 是其儿子节点。这样一旦删除 X ,其 o 子树就被锁住,不能与 X 上面的节点相连了。
#include <bits/stdc++.h>
using namespace std;
#define MAXN (105)
int dfn[MAXN], low[MAXN], cnt;
vector<int> graph[MAXN];
int n;
bool is_gedian[MAXN];
//求割点、桥边
void Tarjan(int X, int Fa)
{
int son=0; //子节点个数(对根节点才有用)
dfn[X]=low[X]=++cnt;
for (int o: graph[X]) {
if (!dfn[o]) { //没走过的点才是儿子
son++;
Tarjan(o,X);
low[X]=min(low[X],low[o]);
if (X!=Fa && low[o]>=dfn[X]) { //X是非根节点的割点
is_gedian[X]=1;
}
if (low[o]>dfn[X]) { //桥边
}
}
else if (o!=Fa) {
low[X]=min(low[X],dfn[o]);
}
}
if (X==Fa && son>1) { //X是根节点,且是割点
is_gedian[X]=1;
}
}
void init()
{
for (int i=0; i<MAXN; i++) graph[i].clear();
memset(dfn,0,sizeof(dfn));
memset(low,0,sizeof(low));
memset(is_gedian,0,sizeof(is_gedian));
cnt=0;
}
bool solve()
{
init();
scanf("%d",&n);
if (n==0) return false;
int uu,vv;
while (1) {
scanf("%d",&uu);
if (uu==0) break;
while (1) {
scanf("%d",&vv);
graph[uu].push_back(vv);
graph[vv].push_back(uu);
if (getchar()=='\n') break;
}
}
for (int i=1; i<=n; i++)
if (!dfn[i])
Tarjan(i,i);
int ans=0;
for (int i=1; i<=n; i++) ans+=(is_gedian[i]);
printf("%d\n",ans);
return true;
}
int main()
{
while (solve())
continue;
}
3. UVA-796 Critical Links
无向图 Tarjan 求桥边。基本上和求割点一样,就是那个大于等于变成了大于号,不能等于。
#include <bits/stdc++.h>
using namespace std;
#define MAXN (1005)
int n;
int dfn[MAXN], cnt;
int low[MAXN];
vector<int> graph[MAXN];
vector<pair<int,int> > ans;
bool vis[MAXN][MAXN];
int ans_num;
//求桥边
void Tarjan(int X, int Fa)
{
dfn[X]=low[X]=++cnt;
for (int o: graph[X]) {
if (!dfn[o]) { //没走过的点才是儿子
Tarjan(o,X);
low[X]=min(low[X],low[o]);
if (low[o]>dfn[X] && !vis[X][o]) { //桥边
vis[X][o]=vis[o][X]=true;
ans.push_back(make_pair(min(X,o),max(X,o)));
ans_num++;
}
}
else if (o!=Fa) {
low[X]=min(low[X],dfn[o]);
}
}
}
void init()
{
ans_num=0;
cnt=0;
memset(dfn,0,sizeof(dfn));
memset(low,0,sizeof(low));
for (int i=0; i<MAXN; i++) graph[i].clear();
ans.clear();
memset(vis,0,sizeof(vis));
}
bool solve()
{
init();
if (scanf("%d",&n)!=1)
return 0;
int uu,vv,tt;
for (int i=0; i<n; i++) {
scanf("%d (%d)", &uu, &tt);
while (tt--) {
scanf("%d",&vv);
graph[uu].push_back(vv);
}
}
for (int i=0; i<n; i++)
if (!dfn[i])
Tarjan(i,i);
sort(ans.begin(),ans.end());
printf("%d critical links\n",ans_num);
for (auto i: ans)
printf("%d - %d\n", i.first, i.second);
puts("");
return 1;
}
int main()
{
while (solve())
continue;
return 0;
}