#include #include #include #include #include #include #include using namespace std; #define N 1005*1005 #define RN 1005 #define M 1005 struct DLX { int n,m,C; int U[N],D[N],L[N],R[N],Row[N],Col[N]; int H[M],S[M],cnt,ans[M]; void init(int _n,int _m) { n=_n; m=_m; for(int i=0; i<=m; i++) { U[i]=D[i]=i; L[i]=(i==0?m:i-1); R[i]=(i==m?0:i+1); S[i]=0; } C=m; for(int i=1; i<=n; i++) H[i]=-1; } void link(int x,int y) { C++; Row[C]=x; Col[C]=y; S[y]++; U[C]=U[y]; D[C]=y; D[U[y]]=C; U[y]=C; if(H[x]==-1) H[x]=L[C]=R[C]=C; else { L[C]=L[H[x]]; R[C]=H[x]; R[L[H[x]]]=C; L[H[x]]=C; } } void del(int x) { R[L[x]]=R[x]; L[R[x]]=L[x]; for(int i=D[x]; i!=x; i=D[i]) { for(int j=R[i]; j!=i; j=R[j]) { U[D[j]]=U[j]; D[U[j]]=D[j]; S[Col[j]]--; } } } void rec(int x) { for(int i=U[x]; i!=x; i=U[i]) { for(int j=L[i]; j!=i; j=L[j]) { U[D[j]]=j; D[U[j]]=j; S[Col[j]]++; } } R[L[x]]=x; L[R[x]]=x; } int dance(int x) { if(R[0]==0 || R[0]>16) { cnt=x; return 1; } int now=R[0]; for(int i=R[0]; i!=0 && i<=16; i=R[i]) { if(S[i]