[loj3525]喷泉公园

先将整张图$x$和$y$​都缩小一半，即"道路"长度变为1，"长椅"变为放在格子中心

如果在没有长椅的限制下也无解（直接dfs即可判定），显然原问题也无解

否则，将所有格子（注意不是点）黑白染色，并强制横向道路（即从$(x,y)$到$(x+1,y)$的道路，纵向同理）对应的长椅在黑色格子的中心、纵向道路对应的长椅在白色格子的中心

在这个限制（构造）下，长椅的限制即变为黑（白）色格子的上下（左右）两边不同时为横（纵）向道路

从上到下、从左到右依次贪心选边（能选即选），简单画图不难发现合法

为了方便，可以使用map维护，复杂度为$o(nlog n)$

#include<bits/stdc++.h>
#include"parks.h"
using namespace std;
#define N 200005
#define vi vector<int>
#define pii pair<int,int>
#define mp make_pair
#define fi first
#define se second
int n,id[N],f[N];
pii a[N];
vi ansu,ansv,ansa,ansb;
map<int,int>mat[N],vis[N];
bool cmp(int x,int y){
    return a[x]<a[y];
}
int find(int k){
    if (k==f[k])return k;
    return f[k]=find(f[k]);
}
void add(int x,int y,pii z){
    vis[z.fi][z.se]=1;
    f[find(x)]=find(y);
    ansu.push_back(x-1);
    ansv.push_back(y-1);
    ansa.push_back(z.fi*2+1);
    ansb.push_back(z.se*2+1);
}
int construct_roads(vi x,vi y){
    n=x.size();
    for(int i=1;i<=n;i++){
        a[i]=mp(x[i-1]/2,y[i-1]/2);
        id[i]=f[i]=mat[a[i].fi][a[i].se]=i;
    }
    sort(id+1,id+n+1,cmp);
    for(int j=1;j<=n;j++){
        int i=id[j],k=mat[a[i].fi+1][a[i].se];
        if ((k)&&(find(i)!=find(k))){
            if ((a[i].fi+a[i].se)&1)add(i,k,a[i]);
            else{
                if (!vis[a[i].fi][a[i].se-1])add(i,k,mp(a[i].fi,a[i].se-1));
            }
        }
        k=mat[a[i].fi][a[i].se+1];
        if ((k)&&(find(i)!=find(k))){
            if ((a[i].fi+a[i].se)%2==0)add(i,k,a[i]);
            else{
                if (!vis[a[i].fi-1][a[i].se])add(i,k,mp(a[i].fi-1,a[i].se));
            }
        }
    }
    for(int i=1;i<=n;i++)
        if (find(i)!=find(1))return 0;
    build(ansu,ansv,ansa,ansb);
    return 1;
}