社区讨论

30pts求助

P7924「EVOI-RD2」旅行家参与者 1已保存回复 1

讨论操作

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

当前回复
1 条
当前快照
1 份
快照标识符
@lozmvfjb
此快照首次捕获于
2023/11/15 18:42
2 年前
此快照最后确认于
2023/11/15 20:13
2 年前
查看原帖
CPP
#include<bits/stdc++.h>
using namespace std;
const int maxn=5*1e5+10,maxm=2*1e6+10;
struct node{
	int ne,bi;
};
int n,m,Q,a[maxn],dfn[maxn],low[maxn],num,sum1,dcc[maxn],dc,c[maxn],sum,d[maxn],f[maxn][22],sun[maxn][22],tot,vis[maxn];
vector<node> t[maxn];
vector<int> t1[maxn];
bool v[maxm];
queue<int> q;
void tarjan(int x,int fa){
	dfn[x]=low[x]=++num;
	for(int i=0;i<t[x].size();i++){
		int y=t[x][i].ne;
		if(!dfn[y]){
			tarjan(y,x);
			low[x]=min(low[x],low[y]);
			if(low[y]>dfn[x]) v[t[x][i].bi]=1;
		}else if(y!=fa) low[x]=min(low[x],dfn[y]);
	}
}
void dfs(int x){
	c[x]=dc;
	dcc[dc]+=a[x];
	for(int i=0;i<t[x].size();i++){
		int y=t[x][i].ne;
		if(c[y]||v[t[x][i].bi]) continue;
		dfs(y);
	}
}
void bfs(){
	q.push(1),d[1]=0;
	while(!q.empty()){
		int x=q.front();
		q.pop();
		for(int i=0;i<t1[x].size();i++){
			int y=t1[x][i];
			if(d[y]) continue;
			d[y]=d[x]+1;
			f[y][0]=x;
			sun[y][0]=dcc[x];
			for(int j=1;j<=tot;j++){
				f[y][j]=f[f[y][j-1]][j-1];
			}
			for(int j=1;j<=tot;j++){
				sun[y][j]=sun[y][j-1]+sun[f[y][j-1]][j-1];
			}
		}
	}
}
int lca(int x,int y){
	if(d[x]<d[y]) swap(x,y);
	for(int i=tot;i>=0;i--) if(d[x]-(1<<i)>=d[y]) x=f[x][i];//sum+=sun[x][i];
	if(x==y) return x;
	for(int i=tot;i>=0;i--) if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];//sum+=sun[x][i]+sun[y][i];
	//sum+=sun[x][0];
	return f[x][0];
}
void dfs1(int x){
	for(int i=0;i<t1[x].size();i++){
		int y=t1[x][i];
		dfs(y);
		vis[x]+=vis[y];
	}
	if(vis[x]) sum1+=dcc[x];
}
int main(){
	scanf("%d%d",&n,&m);	
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	for(int x,y,i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		node temp;
		temp.bi=i,temp.ne=y;
		t[x].push_back(temp);
		node temp1;
		temp1.bi=i,temp1.ne=x;
		t[y].push_back(temp1);
   }
	for(int i=1;i<=n;i++){
		if(!dfn[i]) tarjan(i,0);
	}
	for(int i=1;i<=n;i++){
		if(!c[i]){
			dc++;
			dfs(i);
		}
	}
	for(int i=1;i<=n;i++){
		for(int j=0;j<t[i].size();j++){
			int y=t[i][j].ne;
			if(c[i]==c[y]) continue;
			t1[c[i]].push_back(c[y]);
			t1[c[y]].push_back(c[i]);
		}
	}
	bfs();
	scanf("%d",&Q);
	tot=log(dc)/log(2)+1;
	for(int x,y,i=1;i<=Q;i++){
		scanf("%d%d",&x,&y);
		if(c[x]==c[y]){
			if(!vis[c[x]]) /*sum+=dcc[c[x]],*/vis[c[x]]++;
			continue;
		}else{
			int kkk=lca(c[x],c[y]);
			vis[c[x]]++,vis[c[y]]++,vis[kkk]--,vis[f[kkk][0]]--;
		}
	}
	dfs1(1);
	//for(int i=1;i<=dc;i++) if(vis[i]) sum+=dcc[i];
	cout<<sum1;
}
我发现题解中找联通分量都是用栈弹出后找的
但是为什么找出割边再删去的重建图方式不行

回复

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

正在加载回复...