bzoj1787 紧急集合 lca

打完欧拉序列转rmq发现超了空间。。再打离线发现不能处理距离。。最后改了改欧拉序列rmq不MLE了。。结果发现跑的还没有倍增快。。。/(ㄒoㄒ)/~~以后一定用倍增了。

对于x,y,z，可以发现他们的lca只有两个，即有一个lca与其他两个不同（三个都相同除外），那么这个lca就是集合点。。没了。抑或一下可能会方便一点。

AC代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 1000005
using namespace std;

int n,m,cnt,tot,fst
,pnt
,nxt
,d
,f
[21],mi2[21],lg2
,pos
;
int read(){
int x=0; char ch=getchar();
while (ch<'0' || ch>'9') ch=getchar();
while (ch>='0' && ch<='9'){ x=x*10+ch-'0'; ch=getchar(); }
return x;
}
void add(int aa,int bb){
pnt[++tot]=bb; nxt[tot]=fst[aa]; fst[aa]=tot;
}
void dfs(int x,int last){
f[pos[x]=++cnt][0]=x; d[x]=d[last]+1; int p;
for (p=fst[x]; p; p=nxt[p])
if (pnt[p]!=last){ dfs(pnt[p],x); f[++cnt][0]=x; }
}
int lca(int x,int y){
x=pos[x]; y=pos[y]; if (x>y) swap(x,y);
int k=lg2[y-x+1]; x=f[x][k]; y=f[y-mi2[k]+1][k];
return (d[x]<d[y])?x:y;
}
int dist(int x,int y){ return d[x]+d[y]-(d[lca(x,y)]<<1); }
int main(){
n=read(); m=read(); int i,j;
for (i=1; i<n; i++){
int x=read(),y=read(); add(x,y); add(y,x);
}
dfs(1,0); for (i=2; i<=cnt; i++) lg2[i]=lg2[i>>1]+1;
mi2[0]=1; for (i=1; i<=lg2[cnt]; i++) mi2[i]=mi2[i-1]<<1;
for (i=1; i<=lg2[cnt]; i++)
for (j=1; j<=cnt; j++){
f[j][i]=f[j][i-1]; int k=j+mi2[i-1]; if (k>cnt) continue;
if (d[f[k][i-1]]<d[f[j][i]]) f[j][i]=f[k][i-1];
}
while (m--){
int x=read(),y=read(),z=read(),tmp=lca(x,y)^lca(x,z)^lca(y,z);
printf("%d %d\n",tmp,dist(x,tmp)+dist(y,tmp)+dist(z,tmp));
}
return 0;
}

by lych
2016.2.10