poj 1084 舞蹈链（纠结题）

此题反正我自己是认为poj给的数据范围是有错的，不知道是不是自己太弱了，有大神在的话，欢迎来呸！

其实目的就在于建图，搞的我后来建了一个无比纠结的图，先建立了火柴棍和正方形的一个全图，然后再删除一些火柴和正方形

其实舞蹈链就是一个剪枝的深搜，好好理解就是的。数据n应该大于5的。

#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
int const N = 200;
int n;
int des,dnum;
int mark[N][N],map[N][N];
int rvis[N],cvis[N];

int L[N],R[N],U[N],D[N];
int S[N];
int nCol[N];
int nRow[N];
bool answer[N+1];
int best;
bool hash[N];

void init(int row,int col)
{
    int r=row;
    int w=col;
    for (int i=0;i<=w;i++)
    {
        L[i]=i-1; R[i]=i+1;
        U[i]=D[i]=i;
    }
    L[0]=w;
    R[w]=0;
    int cnt=w+1;
    for (int i=0;i<r;i++)
    {
        int head=cnt,tail=cnt;
        for (int j=0;j<w;j++)
        {
            int c = j+1;
            if(map[i][j]==1)
            {
                S[c]++;
                nCol[cnt]=c;
                nRow[cnt]=i;
                U[D[c]]=cnt;
                D[cnt]=D[c];
                U[cnt]=c;
                D[c]=cnt;
                L[cnt]=tail; R[tail]=cnt;
                R[cnt]=head; L[head]=cnt;
                tail=cnt;
                cnt++;
            }
        }
    }
}
void Remove(int x)
{
    for (int i=D[x];i!=x;i=D[i])
    {
        L[R[i]]=L[i];
        R[L[i]]=R[i];
        S[nCol[i]]--;
    }
}
void Resume(int x)
{
    for (int i=D[x];i!=x;i=D[i])
    {
        L[R[i]]=R[L[i]]=i;
        S[nCol[i]]++;
    }
}
int Hash()
{
    int ans=0;
    memset(hash,false,sizeof(hash));
    for (int c=R[0];c!=0;c=R[c])
    {
        if (!hash[c])
        {
            hash[c]=1;
            ans++;
            for (int i=D[c];i!=c;i=D[i])
                for (int j=R[i];j!=i;j=R[j])
                    hash[nCol[j]]=1;
        }
    }
    return ans;
}
void dfs(int ans)
{
    if (R[0]==0)
    {
        if(ans<best)
        best=ans;
        return;
    }
    int best2 = ans + Hash();
    if (best2 >= best)
        return;

    int c,minnum=100000;
    for (int i=R[0];i!=0;i=R[i])
    {
        if (S[i]==0) return;
        if (S[i]<minnum)
        {
            minnum=S[i];
            c=i;
        }
    }
    for (int i=U[c];i!=c;i=U[i])
    {
        answer[nRow[i]]=true;
        Remove(i);
        for (int j=R[i];j!=i;j=R[j])
            {Remove(j);}
        dfs(ans+1);
        for (int j=L[i];j!=i;j=L[j])
            {Resume(j);}
        Resume(i);
        answer[nRow[i]]=false;
    }
}
int main()
{
    int T;
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d",&n);
        memset(mark,0,sizeof(mark));
        int i,j,k,si,num=0,ct=0;
        for(si=1; si<=n; si++)
        {
            for(i=1; i<=n-si+1; i++)
            {
                for(j=1; j<=n-si+1; j++)
                {
                    for(k=0; k<si; k++)
                    {
                        mark[(i-1)*(2*n+1)+j+k-1][ct]=1;
                        mark[(i-1+si)*(2*n+1)+j+k-1][ct]=1;
                        mark[i*n+(i-1)*(n+1)+j+k*(2*n+1)-1][ct]=1;
                        mark[i*n+(i-1)*(n+1)+j+k*(2*n+1)+si-1][ct]=1;
                    }
                    ct++;
                }
            }
        }
        scanf("%d",&dnum);
        memset(rvis,0,sizeof(rvis));
        memset(cvis,0,sizeof(cvis));
        memset(map,0,sizeof(map));
        for(i=0;i<dnum;i++)
        {
            scanf("%d",&des);
            rvis[des-1]=1;
            for(j=0;j<ct;j++)
            {
                if(mark[des-1][j]==1)
                {
                   if(!cvis[j])
                   cvis[j]=1;
                }
            }
        }
        int row=0,col=0,maxr=2*n*n+2*n;
        for(i=0;i<maxr;i++)
        {
            if(rvis[i])continue;
            col=0;
            int f=1;
            for(j=0;j<ct;j++)
            {
                if(cvis[j])continue;
                if(mark[i][j])f=0;
                map[row][col++]=mark[i][j];
            }
            row++;
            if(f)row--;
        }
        init(row,col);
        best=10000;
        dfs(0);
        printf("%d
",best);
    }
    return 0;
}