hdu 3498 whosyourdaddy 重复覆盖

#include <bits/stdc++.h>
using namespace std;
struct DLX
{
    const static int maxn=64,maxm=64,maxnode=maxn*maxn;
    int n,m,sz,anssz;
    int u[maxnode],d[maxnode],l[maxnode],r[maxnode],row[maxnode],col[maxnode];
    //row col 记录该节点i坐标
    int h[maxn],colcnt[maxm];
    void init(int _n,int _m) //数据位于1,1...n,m
    {
        n=_n,m=_m;
        for(int i=0;i<=m;i++)//1..m是每一列的辅助结点 0代表(0,0)即head
        {
            u[i]=d[i]=i;
            l[i]=i-1;r[i]=i+1;
            col[i]=i;row[i]=0;
            colcnt[i]=0;
        }
        l[0]=m;r[m]=0;
        sz=m;anssz=INT_MAX;
        for(int i=1;i<=n;i++)h[i]=-1;
    }
    void push(int x,int y)
    {
        colcnt[y]++;
        col[++sz]=y;
        row[sz]=x;
        d[sz]=d[y]; //加入列首元素后面
        u[sz]=y;
        u[d[y]]=sz;
        d[y]=sz;
        if(h[x]<0) h[x]=r[sz]=l[sz]=sz; //成为行首元素
        else{
            int hd=h[x]; //加入行首节点右边
            l[sz]=hd;
            r[sz]=r[hd];
            l[r[hd]]=sz;
            r[hd]=sz;
        }
    }
    void repeat_remove(int x){
        for(int i=d[x];i!=x;i=d[i])
            l[r[i]]=l[i],r[l[i]]=r[i];
    }
    void repeat_resume(int x){
        for(int i=u[x];i!=x;i=u[i])
            l[r[i]]=r[l[i]]=i;
    }
    int H(){
        static bool vis[maxm];
        memset(vis,0,sizeof vis);
        int res=0;
        for(int i=r[0];i!=0;i=r[i]){ //如果该列未被标标记,则标记该列中每行的列
            if(vis[i]) continue;
            res++;vis[i]=1;
            for(int j=d[i];j!=i;j=d[j])
                for(int k=r[j];k!=j;k=r[k]) vis[col[k]]=1;
         }
        return res;//易证res <= ans-dep
    }
    void repeat_dance(int dep)
    {
        if(r[0]==0){
            anssz=min(anssz,dep); return ;
        }
        if(dep+H()>=anssz) return ;
        int cur=r[0];
        for(int i=r[0];i!=0;i=r[i]) if(colcnt[i]<colcnt[cur])cur=i;
//        repeat_remove(cur);
        for(int i=d[cur];i!=cur;i=d[i]){
            repeat_remove(i);
            for(int j=r[i];j!=i;j=r[j]) repeat_remove(j);
            repeat_dance(dep+1);
            for(int j=l[i];j!=i;j=l[j]) repeat_resume(j);
            repeat_resume(i);
        }
//        repeat_resume(cur);
    }
};
DLX dlx;
int main()
{
#ifdef shuaishuai
    freopen("in.txt","r",stdin);
#endif // shuaishuai
//把这个图存为矩阵，会发现最少的次数就等于找最少的行使每列都至少有1
    int n,m;
    while(~scanf("%d%d",&n,&m))
    {
        dlx.init(n,n);
        for(int x,y;m--;)
        {
            scanf("%d%d",&x,&y);
            dlx.push(x,y);dlx.push(y,x);
        }
        for(int i=1;i<=n;i++)dlx.push(i,i);
        dlx.repeat_dance(0);
        printf("%d\n",dlx.anssz);
    }
    return 0;
}