BZOJ.4793.[CERC2016]Hangar Hurdles(Kruskal重构树 BFS)

题目链接

\(Description\)

有一个\(n\times n\)的正方形网格，上面有若干障碍点。\(q\)次询问，每次询问把一个正方形箱子从\((x1,y1)\)推到\((x2,y2)\) ，不能接触障碍点且不超出边界，箱子的边长最大能是多少。
\(n\leq1000,q\leq3\times10^5\)。

\(Solution\)

我们可以直接将非障碍点看做顶点，向周围四个非障碍点连边。这样是要求一棵最大生成树。
从每个障碍点做八连通BFS，可以得到通过每个非障碍点位置的最大直径。
然后按最大直径从大到小依次加入点并标为访问过，每次与周围之前被访问过的点连边（也就是枚举权值小的与大的连边，优先连长边，和Kruskal一样）。
其实就是Kruskal重构树。连边时因为权值在点上，所以不需要新建点来表示边了，直接将权值大的点连到权值小的点作为其子节点即可。
两点间路径询问就是直接求LCA的权值了。

复杂度\(O(n^2+q\log n^2)\)。

注意初始化dis为INF（否则无障碍点时所有dis都为0）。

//126688kb  9332ms
#include <queue>
#include <cstdio>
#include <cctype>
#include <vector>
#include <cstring>
#include <algorithm>
#define mp std::make_pair
#define pr std::pair<int,int>
//#define gc() getchar()
#define MAXIN 1000000
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
const int N=1003,M=1e6+5;
const int dd[]={1,0,-1,0,1};
const int dx[]={0,0,-1,-1,-1,1,1,1};
const int dy[]={-1,1,-1,0,1,-1,0,1};

int n,id[N][N],A[M],fa[M];
std::vector<int> v[N];
char IN[MAXIN],*SS=IN,*TT=IN;

namespace HLD
{
    int Enum,H[M],nxt[M],to[M],fa[M],sz[M],son[M],dep[M],top[M];

    inline void AE(int u,int v)
    {
        to[++Enum]=v, nxt[Enum]=H[u], H[u]=Enum;
    }
    inline int LCA(int u,int v)
    {
        while(top[u]!=top[v]) dep[top[u]]>dep[top[v]]?u=fa[top[u]]:v=fa[top[v]];
        return dep[u]<dep[v]?u:v;
    }
    void DFS1(int x)
    {
        int mx=0; sz[x]=1;
        for(int i=H[x],v; i; i=nxt[i])
            if((v=to[i])!=fa[x])
            {
                fa[v]=x, dep[v]=dep[x]+1, DFS1(v), sz[x]+=sz[v];
                if(sz[v]>mx) mx=sz[v], son[x]=v;
            }
    }
    void DFS2(int x,int tp)
    {
        top[x]=tp;
        if(son[x])
        {
            DFS2(son[x],tp);
            for(int i=H[x]; i; i=nxt[i])
                if(to[i]!=fa[x]&&to[i]!=son[x]) DFS2(to[i],to[i]);
        }
    }
}

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}
int Find(int x)
{
    return x==fa[x]?x:fa[x]=Find(fa[x]);
}
void Merge(int x,int y)//x<-y
{
    y=Find(y);
    if(x!=y) fa[y]=x, HLD::AE(x,y);
}
void BFS()
{
    static int dis[N][N];
    static bool vis[N][N];
    static std::queue<pr> q;

    memset(dis,0x3f,sizeof dis);
    for(int i=1; i<=n; ++i)
        for(int j=1; j<=n; ++j)
            if(!id[i][j]) q.push(mp(i,j)),vis[i][j]=1,dis[i][j]=0;
    while(!q.empty())
    {
        int x=q.front().first,y=q.front().second;
        q.pop();
        for(int i=0,xn,yn,d=dis[x][y]; i<8; ++i)
            if((xn=x+dx[i])>0&&xn<=n&&(yn=y+dy[i])>0&&yn<=n&&!vis[xn][yn])
                vis[xn][yn]=1, dis[xn][yn]=d+1, q.push(mp(xn,yn));
    }
    for(int i=1; i<=n; ++i)
        for(int j=1,ds,t; j<=n; ++j)
            if(t=id[i][j])
            {
                ds=std::min(dis[i][j],std::min(std::min(i,n-i+1),std::min(j,n-j+1)));
                ds=(ds<<1)-1;//maxd = 2r-1
                A[t]=ds, v[ds].push_back(t);
            }
}

int main()
{
    static pr ref[M];
    static bool vis[N][N];

    n=read(); int cnt=0;
    for(int i=1; i<=n; ++i)
    {
        register char c=gc(); while(c!='.'&&c!='#') c=gc();
        if(c=='.') id[i][1]=++cnt, ref[cnt]=mp(i,1);
        for(int j=2; j<=n; ++j) if(gc()=='.') id[i][j]=++cnt, ref[cnt]=mp(i,j);
    }

    BFS();
    for(int i=1; i<=cnt; ++i) fa[i]=i;
    for(int i=n; i; --i)
    {
        const std::vector<int> &vec=v[i];
        for(int j=0,l=vec.size(); j<l; ++j)
        {
            int t=vec[j],x=ref[t].first,y=ref[t].second;
            vis[x][y]=1;
            for(int i=0,xn,yn; i<4; ++i)
                if((xn=x+dd[i])>0&&xn<=n&&(yn=y+dd[i+1])>0&&yn<=n&&vis[xn][yn])
                    Merge(t,id[xn][yn]);
        }
    }
    for(int i=1; i<=cnt; ++i) if(fa[i]==i) Merge(0,i);
    HLD::DFS1(0), HLD::DFS2(0,0);

    for(int q=read(),p1,p2,x1,y1,x2,y2; q--; )
    {
        p1=id[read()][read()], p2=id[read()][read()];
        if(!p1||!p2) puts("0");
        else printf("%d\n",A[HLD::LCA(p1,p2)]);
    }
    return 0;
}

转载于:https://www.cnblogs.com/SovietPower/p/9764716.html