CF100015C
主要找到环上任意一条边,有比较dis(u,v),dis(u,a)+w+dis(b,v),dis(u,b)+w+dis(a,u)
然后,然后没了
lca求dis(u,v);dis(u,v)=dis[u]+dis[v]-2*dis[lca(u,v)];
我加油慢慢写
咳咳,我可真是写了半天才写过了样例,lca经测试无误。
可是啊
哭了...
还是附上代码
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;
struct node{
int u,v,w,nxt;
}e[1000101];
int fir[1000101],cnt=0,top[1000101];
void add(int x,int y,int z){
e[++cnt].nxt=fir[x];e[cnt].u=x;e[cnt].v=y;e[cnt].w=z;fir[x]=cnt;
}
int dis[1000101]={0};
int vis[1000101]={0},id[100101],tot=0,p=0,a=0,b=0,ww,siz[1000101],dep[1000101],son[1000101]={0},fa[1000101];
void dfs(int x,int f){
for(int i=fir[x];i;i=e[i].nxt){
int vv=e[i].v;//cout<<vv<<' '<<x<<' '<<i<<endl;
if(vv==f)continue;
if(vis[vv]==1&&a==0&&b==0){
p=i;a=e[i].u;b=e[i].v;ww=e[i].w;e[i].w=-1;return;
}
if(!vis[vv]){
vis[vv]=1; dfs(vv,x);
}
}
}
void dfs1(int x,int f){
siz[x]=1;dep[x]=dep[f]+1;fa[x]=f;
for(int i=fir[x];i;i=e[i].nxt){
int vv=e[i].v;//cout<<vv<<' '<<x<<' '<<i<<endl;
if((a==x&&vv==b)||(vv==a&&x==b))e[i].w=-1;
if(vv==f||(a==x&&vv==b)||(vv==a&&x==b))continue;
if(!vis[vv]&&e[i].w!=-1){
vis[vv]=1;dis[vv]=dis[x]+e[i].w;dfs1(vv,x);
}
siz[x]+=siz[vv];
if(siz[vv]>siz[son[x]])son[x]=vv;
}
}
void dfs2(int x,int topp){
top[x]=topp;
if(son[x])dfs2(son[x],topp);
for(int i=fir[x];i;i=e[i].nxt){
int vv=e[i].v;
if(vv==fa[x]||vv==son[x]||(a==x&&vv==b)||(vv==a&&x==b))continue;
dfs2(vv,vv);
}
}
int lca(int x,int y){
while(top[x]!=top[y]){
if(dep[top[x]]>=dep[top[y]])x=fa[top[x]];
else y=fa[top[y]];
}
return dep[x]<dep[y]?x:y;
}
int diss(int x,int y,int zz){
return dis[x]+dis[y]-2*dis[zz];
}
int n,m;
int main(){
scanf("%d",&n);
for(int i=1;i<=n;i++){
int x,y,z;
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);add(y,x,z);
}
dfs(0,0);
//cout<<a<<' '<<b<<endl;
memset(vis,0,sizeof(vis));siz[0]=0;dep[0]=0;
dfs1(0,0);
dfs2(0,0);
int q;
cin>>q;
//for(int i=1;i<=n;i++)cout<<dis[i]<<' ';
while(q--){
int x,y;
scanf("%d%d",&x,&y);
int zz=lca(x,y),tt=lca(x,a),kk=lca(b,y);
int ans=min(diss(x,y,zz),diss(x,a,tt)+ww+diss(b,y,kk));
tt=lca(x,b);kk=lca(a,y);
ans=min(ans,diss(x,b,tt)+ww+diss(a,y,kk));
printf("%d\n",ans);
}
int lla;
cin>>lla;
return 0;
}
我觉得,照理,树剖的lca是跑得草鸡快的,至少比神马倍增跑得快
可是啊...不想说了