【NOIP2013模拟11.7A组】迷宫花坛

【NOIP2013模拟11.7A组】迷宫花坛

题目

圣玛格丽特学园的一角有一个巨大、如迷宫般的花坛。大约有一个人这么高的大型花坛，做成迷宫的形状，深受中世纪贵族的喜爱。维多利加的小屋就坐落在这迷宫花坛的深处。某一天早晨，久城同学要穿过这巨大的迷宫花坛，去探望感冒的维多利加。
整个迷宫可以用N个路口与M条连接两个不同路口的无向通道来描述。路口被标号为1到N，每条通道有各自的长度。整个迷宫一定是连通的，迷宫中可能存在若干个环路，但是，出于美观考虑，每个路口最多只会属于一个简单环路。例如，图1所示的迷宫是非常美观的，但图2则不符合我们的描述，因为3号路口同属于两个简单环。
你需要回答多个这样的询问：假如久城处在路口x，维多利加的小屋处在路口y，久城最短需要走多少距离才能到达小屋？

Input

第一行2个整数N,M，表示迷宫花坛的路口数和通道数；
接下来N行，每行3个整数x,y,z，描述一条连接路口x与路口y，长度为z的通道；
再接下来1行包含一个整数Q，表示询问数量；
之后Q行，每行2个整数x,y，描述一个询问。

Output

对于每个询问输出一行一个整数，表示最短距离。

Sample Input

4 4
1 2 1
2 3 2
1 3 2
3 4 1
2
2 4
1 3

Sample Output

3
2

Data Constraint

对于30%的数据，N≤100；
另有30%的数据，保证N=M；
对于100%的数据，1≤N≤100,000，Q≤200,000，1≤x,y≤N，1≤z≤1000.

解题思路

首先规定以DFS到的第一个在此环上的点视为父亲。
设$va{l}_i$表示第$i$个点向左或向右到父亲的距离（注意：同一环中，点到父亲的距离的方向是相同的，例如：
这里的$val$依次为$[0,1,3,0]$或$[0,4,2,0]$）
然后我们将每个环变成一个菊花图，其余点连向父亲，边权为这个点到父亲的最短路，将其命名为$Dis$
最后，在已经变好的图树上做lca，最后我们向上跳到的两个点（点$tx$和点$ty$），设$ans$为跳到点$tx$和点$ty$的路径和，两点就会有两种情况：1.两个点在原图中不在同一个环，处理情况和普通的lca一样；2.两个点在原图中在同一个环，那么之前定义的$Dis$和$val$就排上用场了。
怎么用呢？？？
看一看答案有什么情况：
一共有三种情况：1.答案路径不经过父亲且点$ty$到父亲的最短路（即$Dis$）不经过点$tx$，很容易发现答案路径的长度为$ans+\left|va{l}_{tx}-va{l}_{ty}\right|$；2.答案路径经过父亲且点$tx$和$ty$到父亲的最短路都不经过对方，那么答案路径的长度为$ans+Di{s}_{tx}+Di{s}_{ty}$；3.答案路径不经过父亲且点$ty$到父亲的最短路经过点$tx$，则答案路径的长度为$ans+\left|va{l}_{tx}-va{l}_{ty}\right|$。之后我们发现答案路径的长度总是优于另外一条路（这个的原因自己推导一下就能发现了）
综上所述，我们的两点$tx$和$ty$在情况2中答案为$ans+Min(\left|va{l}_{tx}-va{l}_{ty}\right|,Di{s}_{tx}+Di{s}_{ty})$。
就此，此题就轻松切掉啦！！！

Code

写得太丑了……

#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
using namespace std;
queue<int> q;
struct way
{
	int last,to,l;
}w[400005],eg[200005];
int f[100005][20],v[100005][20],deep[100005],fx,fy; 
int dfn[100005],low[100005],zh[100005],top,ti;
int c[100005],p[100005],dis[100005],val[100005];
int n,m,que,cnt,num,ans,head[100005],heade[100005];
bool bz[100005];
inline void adde(int from,int to,int l)
{
	eg[++cnt].l=l;
	eg[cnt].last=heade[from];
	eg[cnt].to=to;
	heade[from]=cnt;
}
inline void addw(int from,int to,int l)
{
	w[++cnt].last=head[from];
	w[cnt].to=to;
	w[cnt].l=l;
	head[from]=cnt;
}
inline void lca(int x,int y)
{
	if (deep[x]<deep[y])swap(x,y);
	for (register int i=19;i>=0;i--)if (deep[f[x][i]]>=deep[y])ans+=v[x][i],x=f[x][i];
	if (x==y)return;
	for (register int i=19;i>=0;i--)if (f[x][i]!=f[y][i])ans+=v[x][i]+v[y][i],x=f[x][i],y=f[y][i];
	fx=x;fy=y;
}
inline void st()
{
	for (register int j=1;j<20;j++)
	{
		for (register int i=1;i<=n;i++)
		{
			f[i][j]=f[f[i][j-1]][j-1];
			v[i][j]=v[i][j-1]+v[f[i][j-1]][j-1];
		}
	}
}
inline void dep(int u)
{
	for (register int i=heade[u];i;i=eg[i].last)
	{
		int to=eg[i].to;
		deep[to]=deep[u]+1;
		v[to][0]=eg[i].l;
		f[to][0]=u;
		dep(to);
	}
}
inline void ge(int u,int fa)
{
	int to;
	bz[u]=1;
	if (p[c[u]]!=u)adde(p[c[u]],u,dis[u]);
	for (register int i=head[u];i;i=w[i].last)
	{
		to=w[i].to;
		if (((i^1)^fa)&&(!bz[to]))
		{
			if (p[c[u]]^p[c[to]])adde(u,to,w[i].l);
			ge(to,i);
		}
	}
}
inline void tarjan(int u,int fa)
{
	int to;
	dfn[u]=low[u]=++ti;
	zh[++top]=u;
	bz[u]=1;
	for (register int i=head[u];i;i=w[i].last)
	{
		to=w[i].to;
		if (((i^1)^fa)&&(!dfn[to]))
		{
			tarjan(to,i);
			low[u]=min(low[u],low[to]);
		}
		else if ((bz[to])&&((i^1)^fa))low[u]=min(low[u],dfn[to]);
	}
	if (dfn[u]==low[u])
	{
		num++;
		int l;
		do
		{
			l=zh[top--];
			c[l]=num;
			bz[l]=0;
		}
		while(u!=l);
		p[num]=zh[top+1];
	}
}
inline void gv(int u,int fa)
{
	int to;bz[u]=1;
	for (register int i=head[u];i;i=w[i].last)
	{
		to=w[i].to;
		if (((i^1)^fa)&&(c[to]==c[u])&&(!bz[to]))
		{
			val[to]=val[u]+w[i].l;
			gv(to,i);
			break;
		}
	}
}
inline void spfa()
{
	val[1]=0;bz[1]=1;
	q.push(1);
	while (q.size())
	{
		int h=q.front();
		q.pop();
		bz[h]=0;
		for (register int i=head[h];i;i=w[i].last)
		{
			if (val[w[i].to]>val[h]+w[i].l)
			{
				val[w[i].to]=val[h]+w[i].l;
				if (!bz[w[i].to])
				{
					q.push(w[i].to);
					bz[w[i].to]=1;
				}
			}
		}
	}
}
inline int read()
{
	int an=0,ww=1;
	char ch=getchar();
	while ((ch<'0')||(ch>'9'))
	{
		if (ch=='-')ww=-1;
		ch=getchar();
	}
	while ((ch>='0')&&(ch<='9'))
	{
		an=(an<<1)+(an<<3)+(ch^48);
		ch=getchar();
	}
	return an*ww;
}
int main()
{
	freopen("garden.in","r",stdin);
	freopen("garden.out","w",stdout);
	int x,y,z;
	scanf("%d%d",&n,&m);
	cnt=1;ti=0;num=0;
	for (register int i=1;i<=n;i++)val[i]=2147483647;
	for (register int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&x,&y,&z);
		addw(x,y,z);
		addw(y,x,z);
	}
	spfa();
	memset(bz,0,sizeof(bz));
	memset(zh,0,sizeof(zh));
	for (register int i=1;i<=n;i++)
	{
		top=0;
		if (!dfn[i])tarjan(i,0);
	}
	memset(bz,0,sizeof(bz));
	for (register int i=1;i<=n;i++)if (p[c[i]]!=i)dis[i]=val[i]-val[p[c[i]]];
	for (register int i=1;i<=n;i++)
	{
		if (p[c[i]]==i)dis[i]=0,val[i]=0,gv(i,0);
	}
	memset(bz,0,sizeof(bz));
	cnt=1;ge(1,0);
	dep(1);st();
	scanf("%d",&que);
	for (register int i=1;i<=que;i++)
	{
		scanf("%d%d",&x,&y);
		ans=fx=fy=0;lca(x,y);
		if (!fx)printf("%d\n",ans);
		else if (c[fx]!=c[fy])printf("%d\n",ans+v[fx][0]+v[fy][0]);
		else printf("%d\n",ans+min(dis[fx]+dis[fy],abs(val[fx]-val[fy])));
	}
	return 0;
}