`
阅读更多

/*
无向图边的双连通分量,在同一个连通分量里的边之间的路径不会有必须经过的点
点和边的下标从1开始。
*/
#include <cstdio>
#include <vector>
using namespace std;
const int N=10005; //点数
const int M=100005; //边数,边的双连通分量的点数是O(M)
const int QN=10005;
const int NUM=N+M; //双连通分量的点数
typedef pair<int, int> pii;
#define mkpii make_pair<int, int>

struct e{
	int v;
	e* nxt;
}es[M<<1], *fir[NUM];
int n, m, en, scnt;//scnt是连通分量的个数
int tes[M][2];
vector<int> has[N]; //每个点相关的 连通分量
int stk[N], minc[N], dep[NUM], vid[N], eid[M];
bool vis[NUM];
int top;
pii qs[QN];
vector<pii> qlca[NUM];
e estk[NUM];  //用于模拟dfs的数据结构
int par[NUM], fa[NUM], anc[QN]; //par[i]为i的直接前驱, fa用于并查集, anc[i]为第i次查询的结果


void add_e(int u, int v){
	es[en].v=v;
	es[en].nxt=fir[u];
	fir[u]=&es[en++];
}
void insert(int u, int v){
	add_e(u, v);
	add_e(v, u);
}
void dfs(int u, int f, int d){
	vis[u]=true;
	minc[u]=dep[u]=d;
	stk[++top]=u;
	e* cur;
	int v;
	for(cur=fir[u]; cur; cur=cur->nxt){
		if(!vis[v=cur->v]){
			dfs(v, u, d+1);
			if(minc[u]>minc[v]){
				minc[u]=minc[v];
			}
			if(minc[v]>dep[u]){ //边u-->v作为一个连通分量
				has[v].push_back(++scnt);
				has[u].push_back(scnt);
				--top;
			}else if(minc[v]==dep[u]){
				++scnt;
				do{
					has[stk[top--]].push_back(scnt);
				}while(stk[top]!=u);
				has[u].push_back(scnt);
			}
		}else if(v!=f){
			if(minc[u]>minc[v]){
				minc[u]=minc[v];
			}
		}
	}
}
void edgeSCC(){
	int i;
	for(i=1; i<=n; i++){
		has[i].clear();
		vis[i]=false;
	}
	for(top=-1, scnt=0, i=1; i<=n; i++){
		if(vis[i]) continue;
		dfs(i, -1, 1);
	}
}
//构建新树
void reBuild(){
	int i, s, j;
	for(en=0, i=1; i<=scnt; i++){
		fir[i]=NULL;
	}
	for(i=1; i<=n; i++){
		s=has[i].size();
		if(s==1){
			vid[i]=has[i][0];
		}else if(s>1){
			vid[i]=++scnt;
			fir[scnt]=NULL;
			for(j=0; j<s; j++){
				insert(scnt, has[i][j]);
			}
		}else{
			vid[i]=0;
		}
	}
	int su, sv, u, v, iu, iv;
	bool flag;
	//给每一条边编号
	for(i=1; i<=m; i++){
		u=tes[i][0];
		v=tes[i][1];
		for(flag=true, su=has[u].size(), iu=0; iu<su&&flag; iu++){
			for(sv=has[v].size(), iv=0; iv<sv&&flag; iv++){
				if(has[u][iu]==has[v][iv]){
					flag=false;
					eid[i]=has[u][iu];
					break;
				}
			}
		}
	}
}

//向栈里添加节点u
void push(int u, int p, int& top){
	top++;
	estk[top].v = u;
	estk[top].nxt = fir[u];
	par[u] = p;
	fa[u] = u;
	vis[u] = true;
}
int findFa(int u){
	int k = u;
	while(k != fa[k]) k = fa[k];
	while(u != k){
		int tf = fa[u];
		fa[u] = k;
		u = tf;
	}
	return k;
}
//qn为查询lca的次数,qs记录查询lca的两个几点,anc记录每次查询的结果
void lca(pii* qs, int qn, int* anc, int n){
	int i, r, u, v, top;
	for(i = 1; i <= n; i++){
		qlca[i].clear();
		vis[i] = false;
	}
	for(i = 1; i <= qn; i++){
		//这里的pii的first表示第几个查询,second表示另一个点
		qlca[qs[i].first].push_back(mkpii(i, qs[i].second));
		qlca[qs[i].second].push_back(mkpii(i, qs[i].first));
		if(qs[i].first == qs[i].second){
			anc[i] = qs[i].first;
		}
	}
	for(i=1; i<= n; i++){
		if(vis[i]) continue;
		r = i;
		top = -1;
		push(r, -1, top);
		dep[r]=1;
		while(top >= 0){
			u = estk[top].v;
			if(!estk[top].nxt){ //u为根的子树已经访问完毕
				if(par[u] >= 0){
					fa[u] = par[u];
				}
				top--;
				continue;
			}
			v = estk[top].nxt->v;
			estk[top].nxt = estk[top].nxt->nxt;
			if(v != par[u]){
				push(v, u, top);
				dep[v]=dep[u]+1;
				//处理查询
				int size = qlca[v].size();
				for(i = 0; i < size; i++){
					int k = qlca[v][i].second;
					if(v != k && vis[k]){ //要求另一个点已经被访问了
						anc[qlca[v][i].first] = findFa(k);
					}
				}
			}
		}
	}
}

bool input(){
	scanf("%d%d", &n, &m);
	if(n==0) return false;
	int i;
	for(en=0, i=1; i<=n; i++){
		fir[i]=NULL;
	}
	for(i=1; i<=m; i++){
		scanf("%d%d", tes[i], tes[i]+1);
		insert(tes[i][0], tes[i][1]);
	}
	return true;
}
void solve(){
	edgeSCC();
	reBuild();
	int i, qn, d;
	scanf("%d", &qn);
	for(i=1; i<=qn; i++){
		scanf("%d%d", &qs[i].first, &qs[i].second);
		qs[i].first=eid[qs[i].first];
		qs[i].second=eid[qs[i].second];
	}
	lca(qs, qn, anc, scnt);
	for(i=1; i<= qn; i++){
		d=dep[qs[i].first]+dep[qs[i].second]-dep[anc[i]]*2;
		printf("%d\n", d>>1);
	}
}
int main(){
	while(input()) solve();
	return 0;
}
 
0
1
分享到:
评论

相关推荐

    hdu.rar_hdu

    HDU(杭州电子科技大学在线评测系统)是一个深受程序员喜爱的在线编程练习平台,它提供了丰富的算法题目供用户挑战,帮助他们提升编程技能和算法理解能力。"hdu.rar_hdu"这个压缩包文件很可能是某位程序员整理的他在...

    HDU_2010.rar_hdu 2010_hdu 20_hdu acm20

    【标题】"HDU_2010.rar"是一个压缩包文件,其中包含了与"HDU 2010"相关的资源,特别是针对"HDU ACM20"比赛的编程题目。"hdu 2010"和"hdu 20"可能是该比赛的不同简称或分类,而"hdu acm20"可能指的是该赛事的第20届...

    HDU题目java实现

    【标题】"HDU题目java实现"所涉及的知识点主要集中在使用Java编程语言解决杭州电子科技大学(HDU)在线评测系统中的算法问题。HDU是一个知名的在线编程竞赛平台,它提供了大量的算法题目供参赛者练习和提交解决方案...

    ACM HDU题目分类

    ACM HDU 题目分类 ACM HDU 题目分类是指对 HDU 在线判题系统中题目的分类,总结了大约十来个分类。这些分类将有助于编程选手更好地理解和解决问题。 DP 问题 DP(Dynamic Programming,动态规划)是一种非常重要...

    hdu1250高精度加法

    ### hdu1250高精度加法 #### 背景介绍 在计算机科学与编程竞赛中,处理大整数运算(特别是加法、减法、乘法等)是常见的需求之一。当数字的位数超过了标准数据类型(如`int`、`long`等)所能表示的最大值时,就需要...

    hdu.rar_HDU 1089.cpp_OJ题求和_hdu_horsekw5_杭电obj

    【标题】"hdu.rar_HDU 1089.cpp_OJ题求和_hdu_horsekw5_杭电obj" 提供的信息是关于一个压缩文件,其中包含了一个名为 "HDU 1089.cpp" 的源代码文件,这个文件是为了解决杭州电子科技大学(Hangzhou Dianzi ...

    HDU DP动态规划

    【标题】"HDU DP动态规划"涉及到的是在算法领域中的动态规划(Dynamic Programming,简称DP)技术,这是解决复杂问题的一种高效方法,尤其适用于有重叠子问题和最优子结构的问题。动态规划通常用于优化多阶段决策...

    HDU1059的代码

    HDU1059的代码

    hdu1001解题报告

    hdu1001解题报告

    hdu 1574 passed sorce

    hdu 1574 passed sorce

    ACM HDU

    【ACM HDU】指的是在ACM(国际大学生程序设计竞赛,International Collegiate Programming Contest)中,参赛者在杭州电子科技大学(Hangzhou Dianzi University,简称HDU)的在线评测系统上完成并已解决的题目集合...

    hdu2101解决方案

    hdu2101AC代码

    杭电ACMhdu1163

    【标题】:杭电ACMhdu1163 【描述】:这是一道源自杭州电子科技大学(Hangzhou Dianzi University,简称HDU)的ACM编程竞赛题目,编号为1163。这类问题通常需要参赛者利用计算机编程解决数学、逻辑或算法上的挑战,...

    HDU.rar_hdu_hdu07_com_shownv9b_www.563hdu.

    【标题】"HDU.rar_hdu_hdu07_com_shownv9b_www.563hdu." 暗示这是一个与HDU(杭州电子科技大学在线编程平台)相关的压缩包,其中可能包含了该平台上的编程竞赛题目或练习题目的源代码。"hdu07"可能是某个特定题目的...

    Hdu1000—2169部分代码

    HDU是杭州电子科技大学(Hangzhou Dianzi University)举办的一个在线编程竞赛平台,全称为HDU Online Judge。ACM是国际大学生程序设计竞赛(International Collegiate Programming Contest)的缩写,是一个全球性的...

    hdu 5007 Post Robot

    hdu 5007 Post Robot 字符串枚举。 暴力一下就可以了。

    HDU acm-PPT课件

    【ACM入门与提高:HDU ACM竞赛课程详解】 ACM(国际大学生程序设计竞赛,International Collegiate Programming Contest,简称ICPC或ACM/ICPC)是一项全球性的竞赛,旨在激发大学生对计算机科学的兴趣,提升他们的...

    hdu_acm_1084.rar_ACM_HDU10_acm10_hdu_hdu 1084

    【标题】"hdu_acm_1084.rar_ACM_HDU10_acm10_hdu_hdu 1084" 提供的是一个关于杭电(HDU)ACM竞赛第1084题的解决方案。该题目可能是在编程竞赛中常见的算法问题,而ACM(国际大学生程序设计竞赛)是全球知名的编程...

    HDU最全ac代码

    HDU(Hangzhou Dianzi University)是国内外知名的在线编程竞赛平台,主要服务于ACM/ICPC(国际大学生程序设计竞赛)以及相关的算法训练。"HDU最全ac代码"这个压缩包很可能是包含了在HDU平台上解题通过的完整源代码...

    hdu动态规划算法集锦

    根据提供的信息,我们可以总结出以下关于“hdu动态规划算法集锦”的知识点: ### 动态规划基础概念 动态规划是一种解决多阶段决策问题的方法,它通过将原问题分解为互相重叠的子问题,利用子问题的解来构建原问题...

Global site tag (gtag.js) - Google Analytics