洛谷P1074 靶形数独（跳舞链）

传送门

坑着，等联赛之后再填（联赛挂了就不填了233）

//minamoto
#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
char buf[1<<21],*p1=buf,*p2=buf;
template<class T>inline bool cmax(T&a,const T&b){return a<b?a=b,1:0;}
inline int read(){
    #define num ch-'0'
    char ch;bool flag=0;int res;
    while(!isdigit(ch=getc()))
    (ch=='-')&&(flag=true);
    for(res=num;isdigit(ch=getc());res=res*10+num);
    (flag)&&(res=-res);
    #undef num
    return res;
}
int score[100]={0,6,6,6,6,6,6,6,6,6,

                6,7,7,7,7,7,7,7,6,
                
                6,7,8,8,8,8,8,7,6,
                
                6,7,8,9,9,9,8,7,6,
                
                6,7,8,9,10,9,8,7,6,
                
                6,7,8,9,9,9,8,7,6,
                
                6,7,8,8,8,8,8,7,6,
                
                6,7,7,7,7,7,7,7,6,

                6,6,6,6,6,6,6,6,6};
const int N=9,mm=N*N*N*N*N*4+N,mn=N*N*N+N;
int mp[10][10];
int ans=-1,sz;
int U[mm],D[mm],L[mm],R[mm],C[mm],X[mm];
int S[mn],Q[mn],H[mn];bool v[mn];
void init(int r,int c){
    //建好虚拟节点 
    for(int i=0;i<=c;++i){
        S[i]=0,U[i]=D[i]=i,
        L[i+1]=i,R[i]=i+1;
    }
    R[sz=c]=0,L[0]=sz;
    while(r) H[r--]=-1;//判断每列是否有节点 
}
void place(int &r,int &c1,int &c2,int &c3,int &c4,int i,int j,int k){
    //看不懂 
    r=((i-1)*N+j-1)*N+k;
    c1=(i-1)*N+j;
    c2=N*N+(i-1)*N+k;
    c3=N*N*2+(j-1)*N+k;
    c4=N*N*3+(((i-1)/3)*3+(j-1)/3)*N+k;
}
void link(int r,int c){
    //S记录每列的元素个数，C是个队列，记录总的节点个数（大概） 
    //好像看不太懂这跳舞链怎么连的…… 
    ++S[C[++sz]=c];
    X[sz]=r,D[sz]=D[c],U[D[c]]=sz,
    U[sz]=c,D[c]=sz;
    if(H[r]==-1) H[r]=L[sz]=R[sz]=sz;//这行没有的话就都先连起来 
    else{
        R[sz]=R[H[r]],L[R[H[r]]]=sz,
        L[sz]=H[r],R[H[r]]=sz;
    }
}
void remove(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])
    D[U[j]]=D[j],U[D[j]]=U[j],--S[C[j]];
}
void resume(int c){
    for(int i=U[c];i!=c;i=U[i])
    for(int j=L[i];j!=i;j=L[j])
    ++S[C[D[U[j]]=U[D[j]]=j]];
    L[R[c]]=R[L[c]]=c;
}
void dance(int k){
    if(!R[0]){
        int res=0;
        for(int i=0;i<k;++i)
        res+=score[(X[Q[i]]-1)/N+1]*((X[Q[i]]-1)%N+1);
        cmax(ans,res);
        return;
    }
    int tmp=mm,c;
    for(int i=R[0];i;i=R[i])
    if(S[i]<tmp) tmp=S[c=i];
    remove(c);
    for(int i=D[c];i!=c;i=D[i]){
        Q[k]=i;
        for(int j=R[i];j!=i;j=R[j]) remove(C[j]);
        dance(k+1);
        for(int j=L[i];j!=i;j=L[j]) resume(C[j]);
    }
    resume(c);
}
int main(){
//    freopen("testdata.in","r",stdin);
    int r,c1,c2,c3,c4;
    init(mn,N*N*4);
    for(int i=1;i<=N;++i)
    for(int j=1;j<=N;++j){
        mp[i][j]=read();
        if(mp[i][j]){
            place(r,c1,c2,c3,c4,i,j,mp[i][j]);
            link(r,c1),link(r,c2),link(r,c3),link(r,c4);
            v[c1]=v[c2]=v[c3]=v[c4]=1;
        }
    }
    for(int i=1;i<=N;++i)
    for(int j=1;j<=N;++j)
    for(int k=1;k<=N;++k){
        place(r,c1,c2,c3,c4,i,j,k);
        if(v[c1]||v[c2]||v[c3]||v[c4])continue;
        link(r,c1),link(r,c2),link(r,c3),link(r,c4);
    }
    dance(0);
    printf("%d\n",ans);
    return 0;
}