[KMOI R1] 军事行动——最小生成树+bfs

[KMOI R1] 军事行动

题目背景

$$他们来了。$$

$$集结军队，干掉他们，一个不留。$$

$$是！$$

题目描述

喵星边境局势愈发紧张，导致发生边境冲突。喵星军队总司令小袁立即对 $Y$ 星采取军事行动。

整个宇宙可以看作一个平面直角坐标系，城市 $1,2,\dots ,n$ 的坐标分别为 $(x_1,y_1),(x_2,y_2),\dots (x_n,y_n)$。

现在小袁率领的若干支舰队（可以理解为小袁的军事力量是无限的）驻扎在边境要塞，城市 $1$ 处。他的舰队会进行以下移动：

• 如果舰队的坐标为 $(x,y)$，那么一天之后它可以移动到 $(x-1,y+2)$ 或 $(x+1,y+2)$ 或 $(x+2,y+1)$ 或 $(x-2,y+1)$ 或 $(x-1,y-2)$ 或 $(x+1,y-2)$ 或 $(x+2,y-1)$ 或 $(x-2,y-1)$ 处。

其中环绕在外的还有一条小行星带，当舰队的坐标 $(x,y)$ 且 $x\le 0$ 或 $y\le 0$ 或 $m<x$ 或 $m<y$ 时，舰队就会撞到小行星带。这是小袁所不想看到的。

现在小袁要攻打城市 $2,3,\dots ,n$，每一次他都会从一个已经占领的城市（城市 $1$ 也算），派出舰队前往城市 $i$ 并攻打它，舰队到达之后的第二天城市 $i$ 就被攻占了。

特别的，小袁在一个舰队前往攻打或攻打一个城市的时候不会派出另外一支舰队，在攻占一座城市后当天可以立即派出另外一支舰队。

小袁想问，最少要花多少时间才能攻占所有的城市。

攻打顺序可以不按照 $2,3\dots n$ 的顺序。

输入格式

第一行一个整数 $n,m$，表示城市个数和小行星带的范围。

接下来 $n$ 行，每一行两个正整数 $(x_i,y_i)$，表示城市 $i$ 的坐标。保证 $1\le x_i,y_i\le m$。

输出格式

一个整数 $ans$，表示最少要花的时间。

样例 #1

样例输入 #1

2 20
1 1
1 3

样例输出 #1

3

样例 #2

样例输入 #2

3 150
1 2
2 4
4 3

样例输出 #2

4

样例 #3

样例输入 #3

10 10
1 4
2 3
2 6
3 6
10 3
1 5
4 2
5 3
2 8
9 2

样例输出 #3

23

样例 #4

样例输入 #4

查看附件的 example4.in

样例输出 #4

查看附件的 example4.out

提示

样例一解释：

舰队在第一天来到了 $(3,2)$ 的位置，第二天到达了城市 $2$ 的位置，第三天占领了城市 $2$。总共花了 $3$ 天。

样例二解释：

舰队在第一天到达了城市 $2$ 的位置，第二天占领了城市 $2$。第三天到达了城市 $3$ 的位置，第四天占领了城市 $3$。总共花了 $4$ 天。

数据范围

本题采用 Subtask 捆绑测试。

子任务编号	测试点编号	$n$	$m$	特殊性质	分值
$1$	$1\sim 2$	$1\le n\le 7$	$4\le m\le 7$	无	$10$
$2$	$3\sim 7$	$1\le n\le 200$	$4\le m\le 70$	无	$25$
$3$	$8\sim 9$	$1\le n\le 150$	$4\le m\le 150$	有	$15$
$4$	$10\sim 20$	$1\le n\le 2000$	$4\le m\le 150$	无	$50$

特殊性质：对于每一个 $1\le i\le n-1$，都有 $x_i=x_{i+1}$。

数据严格保证不会有不同的城市拥有相同的坐标。

思路

这道题很明显bfs以及合并（直接用kruskal来做，每一次他都会从一个已经占领的城市（城市 1 也算），派出舰队前往城市i 并攻打它），但我们得注意一下最后的答案，因为我们根据样例可知，最后的答案要加上$n-1$。

- 特别注意： bool operator<(const E &a)const{ return step<a.step; }这句话多的大那么就是从小到大排序。

代码

#include<iostream>
#include<algorithm>
#include<cstring>
#include<queue>

using namespace std;

const int N = 4e3+10,M = N*N;

int n,m;
struct E{
    int x,y,step;
    bool operator<(const E &a)const{
        return step<a.step;
    }
}e[M];
int x[N],y[N];
int p[N];
bool st[N][N];
int cnt,ans;
int ma[N][N];
int dx[]={-1,1,2,-2,-1,1,2,-2},dy[]={2,2,1,1,-2,-2,-1,-1};

int find(int x){
    if(x!=p[x])p[x]=find(p[x]);
    return p[x];
}

void bfs(int x,int y){
    queue<E>q;
    q.push({x,y,0});
    
    st[x][y]=true;
    
    while(q.size()){
        E t=q.front();
        q.pop();
        
        int xx=t.x,yy=t.y,step=t.step;
        
        for(int i=0;i<8;i++){
            int a=xx+dx[i],b=yy+dy[i];
            if(a<=0||b<=0||a>m||b>m||st[a][b])continue;
            if(ma[a][b])e[++cnt]={ma[x][y],ma[a][b],step+1};
            st[a][b]=true;
            q.push({a,b,step+1});
        }
    }
}

int main(){
    cin>>n>>m;
    
    for(int i=1;i<=n;i++)p[i]=i;
    
    for(int i=1;i<=n;i++){
        scanf("%d%d",&x[i],&y[i]);
    }
    
    for(int i=1;i<=n;i++){
        ma[x[i]][y[i]]=i;
    }
    
    for(int i=1;i<=n;i++){
        memset(st,0,sizeof st);
        bfs(x[i],y[i]);
    }
    
    sort(e+1,e+1+cnt);
    for(int i=1;i<=cnt;i++){
        int pa=find(e[i].x),pb=find(e[i].y),step=e[i].step;
        if(pa!=pb){
            p[pa]=pb;
            ans+=step;
        }
    }
    
    cout<<ans+n-1;
    
    return 0;
    
}