DLX模板

const int maxnode=100010;
const int maxn=1010;
const int maxm=1010;
struct DLX
{
    int L[maxnode],R[maxnode],U[maxnode],D[maxnode],Row[maxnode],Col[maxnode],C[maxm],H[maxn],cnt;
    int ans[maxn],ansd;
    void Init(int n,int m)
    {
        for(int i=0;i<=m;i++)
        {
            L[i]=i-1;R[i]=i+1;
            U[i]=D[i]=i;C[i]=0;
        }
        cnt=m;L[0]=m;R[m]=0;

        ansd=0;

        for(int i=1;i<=n;i++)
            H[i]=0;
    }
    void Link(int x,int y)
    {
        C[Col[++cnt]=y]++;
        Row[cnt]=x;

        U[cnt]=y;
        U[D[y]]=cnt;
        D[cnt]=D[y];
        D[y]=cnt;

        if(H[x])
            L[R[H[x]]]=cnt,R[cnt]=R[H[x]],R[H[x]]=cnt,L[cnt]=H[x];
        else
            H[x]=L[cnt]=R[cnt]=cnt;
    }

    void Delete(int c)
    {
        L[R[c]]=L[c];R[L[c]]=R[c];
        for(int i=D[c];i!=c;i=D[i])
            for(int j=R[i];j!=i;j=R[j])
                --C[Col[j]],U[D[j]]=U[j],D[U[j]]=D[j];
    }

    void Resume(int c)
    {
        L[R[c]]=c;R[L[c]]=c;
        for(int i=U[c];i!=c;i=U[i])
            for(int j=L[i];j!=i;j=L[j])
                ++C[Col[j]],U[D[j]]=j,D[U[j]]=j;
    }

    bool Solve()
    {
        if(!R[0])return true;
        int p=R[0];
        for(int i=R[p];i;i=R[i])
            if(C[p]>C[i])
                p=i;
        Delete(p);
        ansd++;
        for(int i=D[p];i!=p;i=D[i]){
            ans[ansd]=Row[i];
            for(int j=R[i];j!=i;j=R[j])
                Delete(Col[j]);
            if(Solve())
                return true;

            for(int j=L[i];j!=i;j=L[j])
                Resume(Col[j]);
        }
        Resume(p);
        --ansd;
        return false;
    }
}DLX;