Tarjan 學習筆記

coding_goat_qwq發表於2024-10-16
相關連結

相關概念:https://oi-wiki.org/graph/concept/

題單:https://www.luogu.com.cn/training/576466#problems

強連通分量

Tarjan 可以用於求強連通分量,然後就可以縮點了。

具體的,我們建出一個 dfs 樹,設 \(dfn_u\)\(u\) 被訪問的編號,\(low_u\) 表示 \(u\) 能連到最靠祖先的 \(dfn\)
那麼在從 \(u\) 搜尋 \(v\) 時,我們將 \(u\) 先入棧,可以知道:

  • 如果一個點未被訪問,那麼 \(low_u = low_v\),因為 \(v\) 可以到的地方,\(u\) 一定可以到。
  • 如果一個點已被訪問,且在棧中,那麼 \(low_u=dfn_v\),原因是這時 \(v\)\(u\) 的祖先。
  • 其他情況不用管。

注意到一個強連通分量有且僅有一個點使得 \(low_u=dfn_u\),那就是該分量的起始節點,那麼這個強連通分量的所有的點在 Tarjan 演算法的棧中在 \(u\) 之後。

相關題目:B3609 [圖論與代數結構 701] 強連通分量

程式碼
void dfs(int u) {
	vis[u]=1;
	low[u]=dfn[u]=++dfncnt;
	stk[++stktop]=u;
	for(auto v:G[u]) {
		if(!dfn[v]) dfs(v),low[u]=min(low[u],low[v]);
		else if(vis[v]) low[u]=min(low[u],dfn[v]);
	}
	if(dfn[u]==low[u]) {
		anscnt++;
		while(1) {
			int v=stk[stktop--];
			vis[v]=0;
			ans[anscnt].push_back(v);
			if(v==u) break;
		}
		sort(ans[anscnt].begin(),ans[anscnt].end());
	}
}

縮點

我們學會了如何去搜強連通分量,那麼把這些強聯通分量縮成一個點,就是縮點。
一張有向圖,縮點之後會變成一個 DAG。

相關題目:【模板】縮點

程式碼
int n,m;
int val[maxn],totval[maxn];
vector<int> G[maxn];
int low[maxn],dfn[maxn],dfncnt;
int stk[maxn],top,anscnt,bel[maxn];
bool vis[maxn];
int a[maxn],b[maxn];
int deg[maxn],dep[maxn],dp[maxn];
void dfs(int u) {
	vis[u]=1;
	dfn[u]=low[u]=++dfncnt;
	stk[++top]=u;
	for(int v:G[u]) {
		if(!dfn[v]) dfs(v),low[u]=min(low[u],low[v]);
		else if(vis[v]) low[u]=min(low[u],dfn[v]);
	}
	if(dfn[u]==low[u]) {
		anscnt++;
		while(1) {
			int v=stk[top--];
			vis[v]=0;
			bel[v]=anscnt;
			if(u==v) break;
		}
	}
}
signed main() {
	in2(n,m);
	For(i,1,n) in1(val[i]);
	For(i,1,m) {
		in2(a[i],b[i]);
		G[a[i]].push_back(b[i]);
	}
	For(i,1,n) if(!dfn[i]) dfs(i);
	For(i,1,n) G[i].clear(),totval[bel[i]]+=val[i];
	For(i,1,m) a[i]=bel[a[i]],b[i]=bel[b[i]];
	For(i,1,m) if(a[i]!=b[i]) G[a[i]].push_back(b[i]),deg[b[i]]++;
	queue<int> q;
	For(i,1,n) if(!deg[i]) q.push(i),dp[i]=totval[i];
	int ans=0;
	while(!q.empty()) {
		int u=q.front(); q.pop();
		for(int v:G[u]) {
			dep[v]=dep[u]+1;
			dp[v]=max(dp[v],dp[u]+totval[v]);
			deg[v]--;
			if(!deg[v]) q.push(v);
		}
		ans=max(ans,dp[u]);
	}
	cout<<ans<<'\n';
}

相關文章