社区讨论

这份边双为什么T了

学术版参与者 1已保存回复 0

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
0 条
当前快照
1 份
快照标识符
@loxkj9vl
此快照首次捕获于
2023/11/14 08:01
2 年前
此快照最后确认于
2023/11/14 13:29
2 年前
查看原帖
经下载数据,n>50000的点都T了
CPP
#include<iostream>
#include<cstring>
#include<vector>
#include<algorithm>
#include<stack>
using namespace std;
struct edge {
	int u;int v;
};
stack<edge> s;
int pre[500000],bccno[500000],dfs_clock,vis[500000];
vector<int> g[500000],dfn[500000];
int bcc (int u,int fa) {
	int lowu=pre[u]=++dfs_clock;
	int child=0;
	for (int i=0;i<g[u].size();i++) {
		int v=g[u][i];
		if (!pre[v]) {
			child++;
			int lowv=bcc(v,u);
			lowu=min(lowv,lowu);
			if (lowv>pre[u]) {
				s.push((edge){u,v});
			}
		}
		else if (pre[v]<pre[u]&&v!=fa) lowu=min(lowu,pre[v]);
	}	
	return lowu;
}
int fb (int u,int dc) {
	bccno[u]=dc;
	vis[u]=1;
	for (int i=0;i<g[u].size();i++) {
		int v=g[u][i];
		if (v==-1) continue;
		if (vis[v]) continue;
		fb(v,dc);
	}
}
int find_bcc (int n) {
	vector<int> del[500000];
	while(!s.empty()) {
		edge e=s.top();
		s.pop();
		int u=e.u,v=e.v;
		//cout<<u<<' '<<v<<endl;
		del[u].push_back(v);
		del[v].push_back(u);
		for (int i=0;i<n;i++) {
			sort(g[i].begin(),g[i].end());
			sort(del[i].begin(),del[i].end());
		}
	}
	for (int i=1;i<=n;i++) {
		int p=0;
		for (int j=0;j<del[i].size();j++) {
			while (g[i][p]!=del[i][j]) {
				p++;
				if (p>=g[i].size()) break;
			}
			if (p<g[i].size()) g[i][p]=-1;
		}
	}
	dfs_clock=0;
	for (int i=1;i<=n;i++) if (!vis[i]) {
		//cout<<i<<' '<<dfs_clock<<endl;
		fb(i,++dfs_clock);
	}
}
int main() {
	int n,m;
	cin>>n>>m;
	for (int i=0;i<m;i++) {
		int f,t;
		cin>>f>>t;
		g[f].push_back(t);
		g[t].push_back(f);
	}
	for (int i=1;i<=n;i++) if (!pre[i]) bcc(i,-1);
	find_bcc(n);
	for (int i=1;i<=n;i++) dfn[bccno[i]].push_back(i);
	//for (int i=1;i<=n;i++) cout<<pre[i]<<' ';
	//cout<<endl;
	cout<<dfs_clock<<endl;
	for (int i=1;i<=dfs_clock;i++) {
		cout<<dfn[i].size()<<' ';
		for (int j=0;j<dfn[i].size();j++) cout<<dfn[i][j]<<' ';
		cout<<endl;
	}
	return 0;
}
	

回复

0 条回复,欢迎继续交流。

正在加载回复...