雙倍經驗題

題意很簡單，求仙人掌兩點間距離

各種情況討論了我一天多...

在膜拜了一種新的建圖方式後終於思路清晰了

思路：首先樹上的很好搞，dis[x]+dis[y]-2*dis[lca(x,y)]

先用SPFA求出1號點到任意點的距離dis[x]，這個待會求答案時要用

仙人掌也是一種植物，所以我們要把他變成一棵樹

dfs時把每個簡單環取出，每個環總有一個dfn最小的點，把它記做這個環的根

把所有環上的點向環根連邊，記錄它屬於哪個環bel[x]（根不記錄，因為它還有可能是另一個環的節點），順便求出環的長度

環邊去掉，其他邊不變，這就形成了一棵樹

這棵樹有什麼用呢？

首先這棵樹上的點沒有變，只是邊較原仙人掌有所變化

那麼我們就可以求任意兩點x,y,的lca了

接下來就是愉悅的分類討論了

設dep[x]>dep[y]

分別記錄x和y的祖先且是lca的兒子的點fx和fy

1.x和y中一個是lca

畫圖可知這個lca是不是原圖的環都沒有關係

直接return dis[x]-dis[y]

2.fx和fy都在同一個環上

即bel[x]!=0&&bel[x]==bel[y]

那麼ans=dis[x]-dis[fx]+dis[y]-dis[fy]+d，d是fx,和fy在環上的距離

這個畫圖也能知道

d怎麼求，可以在dfs時得到每個點的另一個距離rdis，這個rdis不是到1的最小距離，而是每個環都按一個方向走的距離

或者你可以在取環的時候順便求出環上點到環根的一個方向的距離rdis，這種更好理解一些，雖然我腦子一抽用了上一種方法

這樣就滿足可減了，可以直接用abs(rdis[fx]-rdis[fy])得出fx和fy在環上的一段距離，用len[bel[fx]]-d可得另一種距離，取min即可

3.fx和fy有一個不在環上或x和y不在同一個環上

這個簡單，畫圖即可 ans=dis[x]+dis[y]-dis[lca]*2

#include<cstdio>
#include<vector>
#include<cstring>
#include<iostream>
#include<algorithm>
#define abs(a) (a>0?a:(-(a)))
const int maxn=50010,maxm=100010,maxk=22;
using namespace std;
int n,m,Q,dis[maxn],head,tail,q[maxm+10],sta[maxn],top,rdis[maxn],rcnt,fa[maxn][maxk],dfn[maxn],last[maxn],bel[maxn],rlen[maxn],dep[maxn];
int pre[maxm],now[maxn],son[maxm],val[maxn],tot,tim;bool bo[maxn],del[maxm],vis[maxn];
void add(int a,int b,int c){pre[++tot]=now[a],now[a]=tot,son[tot]=b,val[tot]=c;}
void ins(int a,int b,int c){add(a,b,c),add(b,a,c);}
void spfa(){
	memset(dis,63,sizeof(dis));
	head=dis[1]=0,q[tail=1]=1,bo[1]=1;
	while (head<=tail){
		if (++head>maxm) head=1;
		int x=q[head];
		for (int y=now[x];y;y=pre[y]){
			if (dis[son[y]]>dis[x]+val[y]){
				dis[son[y]]=dis[x]+val[y];
				if (!bo[son[y]]){
					if (++tail>maxm) tail=1;
					q[tail]=son[y],bo[son[y]]=1;
				}
			}
		}
		bo[x]=0;
	}
}

void getring(int st,int ed,int id){
	del[id]=del[id^1]=1,rlen[++rcnt]+=val[id];
	for (int x=ed;x!=st;x=son[last[x]^1]){
		bel[x]=rcnt,del[last[x]]=del[last[x]^1]=1;
		ins(st,x,0),rlen[rcnt]+=val[last[x]];
	}
}

void dfs(int x){
	dfn[x]=++tim;
	for (int y=now[x];y;y=pre[y]) if (y!=(last[x]^1)&&(y<=(m*2+1))){
		if (!dfn[son[y]]) last[son[y]]=y,rdis[son[y]]=rdis[x]+val[y],dfs(son[y]);
		else if (dfn[son[y]]<dfn[x]) getring(son[y],x,y);
	}
}

void dfs2(int x){
	vis[x]=1;
	for (int i=1;i<=18;i++) fa[x][i]=fa[fa[x][i-1]][i-1];
	for (int y=now[x];y;y=pre[y]) if (!del[y]&&son[y]!=fa[x][0])
		dep[son[y]]=dep[x]+1,fa[son[y]][0]=x,dfs2(son[y]);
}

int query(int x,int y){
	if (dep[x]<dep[y]) swap(x,y);
	//printf("step1: %d %d\n",x,y);
	int a=x,b=y;
	for (int h=dep[x]-dep[y],i=18;h&&i>=0;i--) if (h&(1<<i)) h-=(1<<i),x=fa[x][i];
	//printf("step2: %d %d\n",x,y);
	if (x==y) return dis[a]-dis[b];//puts("cas1:鏈 "),
	for (int i=18;i>=0;i--)
		if (fa[x][i]!=fa[y][i])
			x=fa[x][i],y=fa[y][i];
	int lca=fa[x][0];
	//if (!x||!y||!lca) return 0;
	if (bel[x]&&bel[x]==bel[y]){
		//puts("cas2:相同的環  ");
		int d=abs(rdis[x]-rdis[y]);d=min(d,rlen[bel[x]]-d);
		return dis[a]+dis[b]-dis[x]-dis[y]+d;
	}
	//puts("cas3：不同的環或不是環  ");
	return dis[a]+dis[b]-2*dis[lca];
}

int main(){
	scanf("%d%d%d",&n,&m,&Q),tot=1;
	for (int i=1,x,y,z;i<=m;i++) scanf("%d%d%d",&x,&y,&z),ins(x,y,z);
	spfa(),last[1]=-1,dfs(1),dfs2(1);
	for (int i=1,x,y;i<=Q;i++) scanf("%d%d",&x,&y),printf("%d\n",query(x,y));
	return 0;
}

/*
13 18 3
1 2 5

2 14 3
2 9 5
14 9 2

9 15 4
9 10 6
15 10 3

10 11 7
11 12 6
12 13 5
13 10 4

2 3 2
3 4 3
4 5 2
5 6 3
6 7 2
7 8 3
8 3 2

15 12
13 2
15 12
*/