bzoj 4242: 水壶

题意：

给出一些平面上的点，两两间的距离为最短路。求两点间的路径中最大值的最小值。

题解：

容易想到，就是最小生成树上的路径的最大值。
关键怎么用各种奇技淫巧求这个最小生成树呢？
orz题解
那么每个点最多有四个区域重合，即建四条边。
然后就是普通的MST+倍增。
就这么一道题做+调了一下午+一晚上是不是很菜
code：

#include<queue>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
struct node{
    int x,y,c,next;
}b[200010],e[20000010],a[800010];int Len=0,len=0,last[200010];
struct trnode{
    int dep,fa[25],max[25];
    bool vis;
    trnode(){vis=false;}
}tr[200010];
queue<pair<int,int> > q;
int yh[7]={-1,0,1,0};
int yl[7]={0,1,0,-1};
int n,m,p,Q;
int bel[2010][2010],dis[2010][2010];
int map[2010][2010],fa[200010];
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
void insert(int x,int y,int c){e[++Len].x=x;e[Len].y=y;e[Len].c=c;}
bool cmp(node a,node b){return a.c<b.c;}
void ins(int x,int y,int c)
{
    a[++len].y=y;a[len].c=c;
    a[len].next=last[x];last[x]=len;
}
int findfa(int x)
{
    if(fa[x]!=x) fa[x]=findfa(fa[x]);
    return fa[x];
}
void dfs(int x,int fa,int c)
{
    tr[x].dep=tr[fa].dep+1;tr[x].fa[0]=fa;tr[x].max[0]=c;
    tr[x].vis=true;
    for(int i=1;(1<<i)<=tr[x].dep;i++)
    {
        tr[x].fa[i]=tr[tr[x].fa[i-1]].fa[i-1];
        tr[x].max[i]=max(tr[tr[x].fa[i-1]].max[i-1],tr[x].max[i-1]);
    }
    for(int i=last[x];i;i=a[i].next)
    {
        int y=a[i].y;
        if(y==fa) continue;
        dfs(y,x,a[i].c);
    }
}
int solve(int x,int y)
{
    int tx=findfa(x),ty=findfa(y);
    if(tx!=ty) return -1;
    if(tr[x].dep<tr[y].dep) swap(x,y);
    int ans=0;
    for(int i=19;i>=0;i--)
        if((1<<i)<=tr[x].dep-tr[y].dep) ans=max(ans,tr[x].max[i]),x=tr[x].fa[i];
    if(x==y) return ans;
    for(int i=19;i>=0;i--)
        if((1<<i)<=tr[x].dep&&tr[x].fa[i]!=tr[y].fa[i])
        {
            ans=max(ans,tr[x].max[i]);
            ans=max(ans,tr[y].max[i]);
            x=tr[x].fa[i];y=tr[y].fa[i];
        }
    ans=max(ans,tr[x].max[0]);
    ans=max(ans,tr[y].max[0]);
    return ans;
}
int main()
{
    scanf("%d %d %d %d",&n,&m,&p,&Q);
    memset(map,-1,sizeof(map));
    for (int i=1;i<=n;i++)
    {
        char ch;
        while (ch!='.' && ch!='#') ch=getchar();
        for (int j=1;j<=m;j++,ch=getchar())
            map[i][j]=(ch=='.')?0:-1;
    }
    memset(dis,-1,sizeof(dis));
    for(int i=1;i<=p;i++)
    {
        b[i].x=read();b[i].y=read();
        int x=b[i].x,y=b[i].y;
        bel[x][y]=i;dis[x][y]=0;
        q.push(make_pair(x,y));
    }
    while(!q.empty())
    {
        int x=q.front().first,y=q.front().second;q.pop();
        int t=bel[x][y];
        for(int i=0;i<4;i++)
        {
            int tx=x+yh[i],ty=y+yl[i];
            if(map[tx][ty]==-1) continue;
            if(bel[tx][ty]==0)
            {
                bel[tx][ty]=t;
                dis[tx][ty]=dis[x][y]+1;
                q.push(make_pair(tx,ty));
            }
            else
                if(bel[tx][ty]!=bel[x][y]) insert(bel[tx][ty],t,dis[tx][ty]+dis[x][y]);
        }
    }
    sort(e+1,e+Len+1,cmp);
    for(int i=1;i<=p;i++) fa[i]=i;
    for(int i=1;i<=Len;i++)
    {
        int x=e[i].x,y=e[i].y,tx=findfa(x),ty=findfa(y);
        if(tx!=ty) fa[tx]=ty,ins(x,y,e[i].c),ins(y,x,e[i].c);
    }
    tr[0].dep=-1;
    for(int i=1;i<=p;i++) if(!tr[i].vis) dfs(i,0,0);
    while(Q--)
    {
        int S,T;S=read();T=read();
        printf("%d\n",solve(S,T));
    }
}