codeforces 100959B Airports 曼哈顿距离最大生成树

写了两种做法。

1. 对于每个点，向八个象限建最长边，注意最长边没有最短边的对称性，故每个点不能只枚举4个方向，要8个方向都枚举。这个做法速度较快，用的树状数组：

#include 
using namespace std;
typedef pair pii;
const int maxn=100005;
const int inf=2e9+1;
int n,pos[maxn],f[maxn],bit[maxn],res,x,y;
struct edge {
    int u,v,w;
    edge(){}
    edge(int u,int v,int w):u(u),v(v),w(w){}
    bool operator<(const edge &rhs)const{
        return w>rhs.w;
    }
};
vector v;
vector ID;
inline void init() {
    for (int i=0;ival) {
            val=bit[x];
            ret=pos[x];
        }
        x+=x&-x;
    }
    return ret;
}
struct P {
    int x,y,id;
    P(){}
    P(int x,int y,int id):x(x),y(y),id(id){}
    bool operator<(const P &rhs)const{
        return x==rhs.x?y=0;--i) {
        int o=id(p[i].y-p[i].x);
        int t=query(o);
        if (t!=-1)
            v.push_back(edge(p[i].id,p[t].id,manh(i,t)));
        update(o,p[i].x+p[i].y,i);
    }
}
inline void build() {
    for (int i=0;i<8;++i) {
        if (i&1)
            for (int j=0;j

2. Boruvka 算法

速度比上一个慢，用的 set 维护最远曼哈顿距离，去 begin() 或者 end()-1 均是最远点。

(1) begin():

#include 
using namespace std;
typedef pair pii;
typedef pair pip;
const int dx[4]={-1,-1,1,1};
const int dy[4]={1,-1,-1,1};
const int maxn=100005;
int n,x[maxn],y[maxn],f[maxn],val[maxn][4],res=2e9+1;
bool connected;
vector V[maxn];
vector A;
set st[4];
inline void init() {
    for (int i=0;iV[y].size())
        swap(x,y);
    for (int i=0;i<(int)V[x].size();++i)
        V[y].push_back(V[x][i]);
    V[x].clear();
    f[x]=y;
    return true;
}
int main()
{
    scanf("%d",&n);
    for (int i=0;i::iterator it;
            int v=-1,mx=-1;
            for (int j=0;j<(int)V[i].size();++j)
                for (int k=0;k<4;++k) {
                    it=st[(k+2)%4].begin();
                    if (v==-1||mx

(2) end()-1:

#include 
using namespace std;
typedef pair pii;
typedef pair pip;
const int dx[4]={-1,-1,1,1};
const int dy[4]={1,-1,-1,1};
const int maxn=100005;
int n,x[maxn],y[maxn],f[maxn],val[maxn][4],res=2e9+1;
bool connected;
vector V[maxn];
vector A;
set st[4];
inline void init() {
    for (int i=0;iV[y].size())
        swap(x,y);
    for (int i=0;i<(int)V[x].size();++i)
        V[y].push_back(V[x][i]);
    V[x].clear();
    f[x]=y;
    return true;
}
int main()
{
    scanf("%d",&n);
    for (int i=0;i::iterator it;
            int v=-1,mx=-1;
            for (int j=0;j<(int)V[i].size();++j)
                for (int k=0;k<4;++k) {
                    it=st[(k+2)%4].end();
                    --it;
                    if (v==-1||mx