[BZOJ 2791]POI2012 Rendezvous

逗逗逗，一开始看错了题目，然后杯具了一个多小时。。。。。

看明白了题目还是思路还是很直观的。

首先若果两个点的在环之外就可以走到一起，很明显这就是最优解。

如果在环内的话，那么我们分找出两个点第一个进入环的祖先，很容易证明最优解一定是这两个祖先中的一个，然后。。。敲代码吧

#include <cstdio>
#include <cstring>
#include <algorithm>

using namespace std;
const int Maxn=500005;
int n,m,i,j,fa[Maxn][20],t,x,y,x1,y1,x2,y2,l,r,tot,N;
int v[Maxn],num[Maxn],nds[Maxn],cir[Maxn],lk[Maxn];
int q[Maxn],dep[Maxn],node[Maxn],next[Maxn],a[Maxn];

void add(int x,int y)
{ node[++tot]=y; next[tot]=a[x]; a[x]=tot; }

void init(){
  for (j=1;j<19;j++)
    for (i=1;i<=n;i++)
      fa[i][j]=fa[ fa[i][j-1] ][j-1];

  for (i=1;i<=n;i++)
  if (!v[i]){
  	for (j=i;!v[j];j=fa[j][0]) v[j]=i;
  	if (v[j]<i) continue;
  	for (nds[++N]=0;cir[j]==0;j=fa[j][0])
	  num[j] = (++nds[ cir[j]=N ]), lk[j]=j;
  }
  
  memset(dep,10,sizeof(dep));
  for (i=1;i<=n;i++)
    if (cir[i]>0) dep[ q[++r]=i ]=1;
  for (l=1;l<=r;l++){
  	for (i=a[q[l]];i;i=next[i])
  	if (dep[node[i]]>dep[q[l]]+1){
  	  dep[ q[++r]=node[i] ]=dep[q[l]]+1;
  	  cir[ q[r] ] = cir[ q[l] ];
  	  lk [ q[r] ] = lk [ q[l] ];
    }
  }
}

int LCA(int x,int y){
  if (dep[x]<dep[y]) swap(x,y);
  for (int i=18;i>=0;i--)
    if (dep[fa[x][i]]>dep[y]) x=fa[x][i];
  if (dep[x]>dep[y]) x=fa[x][0];
  for (int i=18;i>=0;i--)
   if (fa[x][i]!=fa[y][i])
    x=fa[x][i], y=fa[y][i];
  if (x==y) return x;
  return fa[x][0];
}

int way(int x,int y,int s){
  if (y>x) return y-x;
  return s+(y-x);
}

bool Judge(){
  if (max(x1,y1)<max(x2,y2)) return 0;
  if (max(x1,y1)>max(x2,y2)) return 1;
  if (min(x1,y1)<min(x2,y2)) return 0;
  if (min(x1,y1)>min(x2,y2)) return 1;
  return !(x1>=y1 && x2<y2);
}

int main(){
  freopen("ran.in","r",stdin);
  freopen("ran.out","w",stdout);
  scanf("%d%d",&n,&m);
  for (i=1;i<=n;i++){
    scanf("%d",&fa[i][0]);
    add(fa[i][0],i);
  }
  init();
  while (m--){
  	scanf("%d%d",&x,&y);
  	if (x==y) {printf("0 0\n"); continue;}
  	if (cir[x]!=cir[y])
  	  {printf("-1 -1\n");continue;}
  	if (lk[x]==lk[y]){
  	  t=LCA(x,y);
  	  printf("%d %d\n",dep[x]-dep[t],dep[y]-dep[t]);
  	  continue;
  	}
  	x1=dep[x]-1; y1=dep[y]-1+way(num[lk[y]],num[lk[x]],nds[cir[x]]);
  	x2=dep[x]-1+way(num[lk[x]],num[lk[y]],nds[cir[x]]); y2=dep[y]-1;
  	if (!Judge()) printf("%d %d\n",x1,y1);
  	 else printf("%d %d\n",x2,y2);
  }
  
  return 0;
}