HDU5046 Airport dancing links 重复覆盖+二分

这一道题和HDU2295是一样

是一个dancing links重复覆盖解决最小支配集的问题

在给定长度下求一个最小支配集，只要小于k就行

然后就是二分答案，每次求最小支配集

只不过HDU2295是浮点，这里是整数

我写的一个比较暴力

#include<cstdio>
#include<cstring>
#include<queue>
#include<cstdlib>
#include<algorithm>
#include<vector>
#include<cmath>
using namespace std;
typedef long long LL;
const int N=4e3;
int n,m,sz,k;
int u[N],l[N],r[N],d[N];
int h[],s[],col[N];
void init()
{
    for(int i=; i<=m; ++i)
    {
        s[i]=;
        u[i]=d[i]=i;
        l[i]=i-;
        r[i]=i+;
    }
    r[m]=;
    l[]=m;
    sz=m;
    for(int i=; i<=n; ++i)
        h[i]=-;
}
void link(int x,int y)
{
    ++sz;
    ++s[y],col[sz]=y;
    u[sz]=u[y],d[u[y]]=sz;
    d[sz]=y,u[y]=sz;
    if(h[x]==-)h[x]=l[sz]=r[sz]=sz;
    {
        l[sz]=l[h[x]];
        r[l[h[x]]]=sz;
        r[sz]=h[x];
        l[h[x]]=sz;
    }
}
void del(int y)
{
    for(int i=d[y]; i!=y; i=d[i])
        r[l[i]]=r[i],l[r[i]]=l[i];
}
void resume(int y)
{
    for(int i=d[y]; i!=y; i=d[i])
        r[l[i]]=l[r[i]]=i;
}
bool vis[];
int f()
{
    int ret=;
    for(int i=r[]; i; i=r[i])
        vis[i]=;
    for(int i=r[]; i; i=r[i])
    {
        if(vis[i])continue;
        vis[i]=;
        ++ret;
        for(int j=d[i]; j!=i; j=d[j])
            for(int k=r[j]; k!=j; k=r[k])
                vis[col[k]]=;
    }
    return ret;
}
bool dance(int pos)
{
    if(pos+f()>k)return ;
    if(!r[])
    {
        if(pos<=k) return ;
        return ;
    }
    int t=r[];
    for(int i=r[]; i!=; i=r[i])
        if(s[i]<s[t])t=i;
    for(int i=d[t]; i!=t; i=d[i])
    {
        del(i);
        for(int j=r[i]; j!=i; j=r[j])
            del(j);
        if(dance(pos+))return ;
        for(int j=l[i]; j!=i; j=l[j])
            resume(j);
        resume(i);
    }
    return ;
}
struct point
{
    LL x,y;
}o[];
LL ABS(LL x,LL y)
{
    if(x>=y)return x-y;
    return y-x;
}
LL dis(point a,point b)
{
    return ABS(a.x,b.x)+ABS(a.y,b.y);
}
LL D[N];
int main()
{
    int T,cas=;
    scanf("%d",&T);
    while(T--)
    {
      scanf("%d%d",&n,&k),m=n;
      for(int i=;i<=n;++i)
       scanf("%I64d%I64d",&o[i].x,&o[i].y);
      int cnt=;
      for(int i=;i<=n;++i)
           for(int j=i;j<=n;++j)
              D[++cnt]=dis(o[i],o[j]);
      sort(D+,D++cnt);
      cnt=unique(D+,D++cnt)-D-;
      int high=cnt,low=,mid;
      while(low<high)
      {
          mid=(low+high)>>;
          init();
          for(int i=;i<=n;++i)
           for(int j=;j<=n;++j)
              if(dis(o[i],o[j])<=D[mid])
                link(i,j);
          if(dance())high=mid;
          else low=mid+;
      }
      printf("Case #%d: %I64d\n",++cas,D[low]);
    }
    return ;
}

然后另一个是离散化的

#include<cstdio>
#include<cstring>
#include<queue>
#include<cstdlib>
#include<algorithm>
#include<vector>
#include<cmath>
using namespace std;
typedef long long LL;
const int N=4e3;
int n,m,sz,k;
int u[N],l[N],r[N],d[N];
int h[],s[],col[N];
void init()
{
    for(int i=; i<=m; ++i)
    {
        s[i]=;
        u[i]=d[i]=i;
        l[i]=i-;
        r[i]=i+;
    }
    r[m]=;
    l[]=m;
    sz=m;
    for(int i=; i<=n; ++i)
        h[i]=-;
}
void link(int x,int y)
{
    ++sz;
    ++s[y],col[sz]=y;
    u[sz]=u[y],d[u[y]]=sz;
    d[sz]=y,u[y]=sz;
    if(h[x]==-)h[x]=l[sz]=r[sz]=sz;
    {
        l[sz]=l[h[x]];
        r[l[h[x]]]=sz;
        r[sz]=h[x];
        l[h[x]]=sz;
    }
}
void del(int y)
{
    for(int i=d[y]; i!=y; i=d[i])
        r[l[i]]=r[i],l[r[i]]=l[i];
}
void resume(int y)
{
    for(int i=d[y]; i!=y; i=d[i])
        r[l[i]]=l[r[i]]=i;
}
bool vis[];
int f()
{
    int ret=;
    for(int i=r[]; i; i=r[i])
        vis[i]=;
    for(int i=r[]; i; i=r[i])
    {
        if(vis[i])continue;
        vis[i]=;
        ++ret;
        for(int j=d[i]; j!=i; j=d[j])
            for(int k=r[j]; k!=j; k=r[k])
                vis[col[k]]=;
    }
    return ret;
}
bool dance(int pos)
{
    if(pos+f()>k)return ;
    if(!r[])
    {
        if(pos<=k) return ;
        return ;
    }
    int t=r[];
    for(int i=r[]; i!=; i=r[i])
        if(s[i]<s[t])t=i;
    for(int i=d[t]; i!=t; i=d[i])
    {
        del(i);
        for(int j=r[i]; j!=i; j=r[j])
            del(j);
        if(dance(pos+))return ;
        for(int j=l[i]; j!=i; j=l[j])
            resume(j);
        resume(i);
    }
    return ;
}
struct point
{
    LL x,y;
}o[];
LL ABS(LL x,LL y)
{
    if(x>=y)return x-y;
    return y-x;
}
LL dis(point a,point b)
{
    return ABS(a.x,b.x)+ABS(a.y,b.y);
}
LL D[N];
int main()
{
    int T,cas=;
    scanf("%d",&T);
    while(T--)
    {
      scanf("%d%d",&n,&k),m=n;
      for(int i=;i<=n;++i)
       scanf("%I64d%I64d",&o[i].x,&o[i].y);
      int cnt=;
      for(int i=;i<=n;++i)
           for(int j=i;j<=n;++j)
              D[++cnt]=dis(o[i],o[j]);
      sort(D+,D++cnt);
      cnt=unique(D+,D++cnt)-D-;
      int high=cnt,low=,mid;
      while(low<high)
      {
          mid=(low+high)>>;
          init();
          for(int i=;i<=n;++i)
           for(int j=;j<=n;++j)
              if(dis(o[i],o[j])<=D[mid])
                link(i,j);
          if(dance())high=mid;
          else low=mid+;
      }
      printf("Case #%d: %I64d\n",++cas,D[low]);
    }
    return ;
}