P2783有机化学之神36pts求助
  • 板块学术版
  • 楼主PCCP
  • 当前回复3
  • 已保存回复3
  • 发布时间2022/11/16 20:33
  • 上次更新2023/10/27 02:44:36
查看原帖
P2783有机化学之神36pts求助
310773
PCCP楼主2022/11/16 20:33

蒟蒻 #5 #6TLETLE,#4 和 #8-11WAWA ,找到的数据和自建的数据都过了,实在不知道哪里有问题,于是求助谷内大佬,能否帮忙看看哪里有问题,或者给蒟蒻一些警示(

代码如下 (马蜂不行,敬请原谅)

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<stack>
using namespace std;
const int N=1e6+10;
int n,m,o,in[N];
int he[N<<1],ne[N<<1],to[N<<1],tot=1;
void addedge(int x,int y){
	to[++tot]=y;
	ne[tot]=he[x];
	he[x]=tot; 
}
int her[N<<1],ner[N<<1],tor[N<<1];
void add(int x,int y){
	tor[++tot]=y;
	ner[tot]=her[x];
	her[x]=tot;
}
int dfo[N],low[N],rootid[N],cnt=0,edcccnt=0;
stack <int> q;
vector <int> edcc[N];
void tarjan(int now,int inedge){
	dfo[now]=low[now]=++cnt;
	q.push(now);
	for(int i=he[now];i;i=ne[i]){
		int v=to[i];
		if(i==(inedge^1)){
			continue;
		}
		else if(!dfo[v]){
			tarjan(v,i);
			low[now]=min(low[now],low[v]);
			if(low[v]>dfo[now]){
				++edcccnt;
				int p;
				do{
					p=q.top();
					edcc[edcccnt].push_back(p);
					rootid[p]=edcccnt;
					q.pop();
				}
				while(p!=v);
			}
		}
		else{
			low[now]=min(low[now],dfo[v]);
		}
	}
}
int fa[N],dep[N],f[30][N];
void dfs(int now,int father,int deep){
	fa[now]=father;
	dep[now]=deep;
	for(int i=her[now];i;i=ner[i]){
		int v=tor[i];
		if(v==fa[now]){
			continue;
		}
		dfs(v,now,deep+1);
		f[0][v]=now;
	}
}
int lca(int x,int y){
	int l=0,ll=x,rr=y;
	while((1<<l)<=n){
		l++;
	}
	if(dep[x]<dep[y]){
		swap(x,y);
	}
	for(int i=l;i>=0;i--){
		if(dep[y]<=dep[x]-(1<<i)){
			x=f[i][x];
		}
	}
	if(x==y){
		return (dep[ll]-dep[f[0][x]])+(dep[rr]-dep[f[0][x]])-1;
	}
	for(int i=l;i>=0;i--){
		if(f[i][x]!=f[i][y]){
			x=f[i][x];
			y=f[i][y];
		}
	}
	return (dep[ll]-dep[f[0][x]])+(dep[rr]-dep[f[0][x]])-1;
}
void print(int x){
	int xx=0,t[50];
    for(;x;x>>=1){
    	xx++;
    	if(x&1){
    		t[xx]=1;
		}
    	else{
    		t[xx]=0;
		}
    }
    for(int i=xx;i>=1;i--){
    	printf("%d",t[i]);
	}
    printf("\n");
}
int main(){
	scanf("%d%d",&n,&m);
	int x,y,flag;
	for(int i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		flag=0;
		for(int j=he[x];j;j=ne[j]){
			if(to[j]==y){
				flag++;
				break;
			}
		}
		if(x==y||flag>0){
			continue; 
		}
		addedge(x,y);
		addedge(y,x);
	}
	for(int i=1;i<=n;i++){
		if(!dfo[i]){
			addedge(i+n,i);
			tarjan(i+n,0);
		}
	}
	tot=0;
	for(int i=1;i<=n;i++){
		for(int j=he[i];j;j=ne[j]){
			int v=to[j];
			if(rootid[i]!=rootid[v]){
				add(rootid[i],rootid[v]);
				add(rootid[v],rootid[i]);
			}
		}
	}
	dfs(1,0,1);
	scanf("%d",&o);
	int ans=0;
	for(int i=1;i<=o;i++){
		scanf("%d%d",&x,&y);
		ans=lca(rootid[x],rootid[y]);
		print(ans);
	}
}
2022/11/16 20:33
加载中...