HDU 5046 Airport （DLX可重复覆盖+二分）

题意：

在N个城市选出K个城市，建飞机场(1 ≤ N ≤ 60,1 ≤ K ≤ N)，N个城市给出坐标，选择这K个机场，使得从城市到距离自己最近的机场的 最大的距离 最小。

输出这个最小值。

思路：

我居然想歪，歪到想到了模式识别的K均值算法，简直瞎搞2333333.（虽说和DLX一样都是暴力）

DLX可重复覆盖+二分距离制关系表即可

代码：

kuangbin模板，注意爆int

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=65;
const double eps=1e-7;
int L[maxn*maxn],R[maxn*maxn],U[maxn*maxn],D[maxn*maxn];//节点的上下左右四个方向的链表
int C[maxn*maxn],H[maxn],cnt[maxn],vis[maxn];//C列H行cnt列链表中元素个数
int n,m,id,fans,T,lim,cut;
long long le,ri,mid;
typedef struct Node{
    long long x;long long y;
    Node(long long xx=0,long long yy=0){x=xx;y=yy;}
}Node;
Node city[maxn];
void init(){
    for(int i=0;i<=cut;i++){
        cnt[i]=0;U[i]=D[i]=i;
        L[i+1]=i;R[i]=i+1;
    }
    R[cut]=0;id=cut+1;
    memset(H,-1,sizeof(H));
}
void Link(int r,int c){
    cnt[c]++;C[id]=c;
    U[id]=U[c];D[U[c]]=id;
    D[id]=c;U[c]=id;
    if(H[r]==-1) H[r]=L[id]=R[id]=id;
    else{
        L[id]=L[H[r]];R[L[H[r]]]=id;
        R[id]=H[r];L[H[r]]=id;
    }
    id++;
}
void Remove(int Size){
    for(int j=D[Size];j!=Size;j=D[j])
        L[R[j]]=L[j],R[L[j]]=R[j];
}
void Resume(int Size){
    for(int j=D[Size];j!=Size;j=D[j])
        L[R[j]]=R[L[j]]=j;
}
int h()
{
    int sum=0;
    memset(vis,0,sizeof(vis));
    for(int i=R[0];i;i=R[i]){
        if(vis[i]) continue;
        sum++;
        for(int j=D[i];j!=i;j=D[j]){
            for(int k=R[j];k!=j;k=R[k])
                vis[C[k]]=1;
        }
    }
    return sum;
}
bool Dance(int k){
    int mm=maxn,pos;
    if(k+h()>lim) return 0;
    if(!R[0]){
        return k<=lim;
    }
    for(int i=R[0];i;i=R[i])
        if(mm>cnt[i]) mm=cnt[i],pos=i;

    for(int i=D[pos];i!=pos;i=D[i])
    {
        Remove(i);
        for(int j=R[i];j!=i;j=R[j]) Remove(j);
        if(Dance(k+1)) return 1;
        for(int j=R[i];j!=i;j=R[j]) Resume(j);
        Resume(i);
    }
    return 0;
}
bool get(long long dis){
    init();
    for(int i=0;i<n;i++){
        for(int j=0;j<n;j++)
        if((llabs(city[i].x-city[j].x)+llabs(city[i].y-city[j].y))<=dis){
            Link(i+1,j+1);
        }
    }
    return Dance(0);
}
int main(){
    scanf("%d",&T);
    for(int Case=1;Case<=T;Case++){
        scanf("%d%d",&n,&lim);
        cut=n;
        for(int i=0;i<n;i++)
            scanf("%lld%lld",&city[i].x,&city[i].y);
        le=0;
        ri=0;
        for(int i=0;i<n;i++){
            for(int j=0;j<n;j++){
                ri=max(ri,(llabs(city[i].x-city[j].x)+llabs(city[i].y-city[j].y)));
            }
        }
        while(le<ri){
            mid=(le+ri)/2;
            if(get(mid)){
                ri=mid;
            }else{
                le=mid+1;
            }
        }
        printf("Case #%d: %lld\n",Case,le);
    }
}