如题= =自用模板

精确覆盖问题 选取几行使每一列有且仅有1个1
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
import java.io.*;
import java.util.*;
public class Main {

public static void main(String[] args) {
FastScanner sc=new FastScanner();
PrintWriter pw=new PrintWriter(System.out);
DancingLinks dl=new DancingLinks();
while(sc.hasNext()){
int n=sc.nextInt();
int m=sc.nextInt();
dl.init(n, m);
for(int i=1;i<=n;i++){
int k=sc.nextInt();
for(int o=1;o<=k;o++){
int j=sc.nextInt();
dl.link(i, j);
}
}
boolean judge=dl.dance(0);
if(!judge){
pw.println("NO");
}else{
pw.print(dl.ansd);
for(int i=0;i<dl.ansd;i++){
pw.print(" "+dl.ans[i]);
}
pw.println();
}
pw.flush();
}

}

}
class DancingLinks{
/*图的下标为1~n 1~m*/
int n,m,size;
int maxn=100100;
int[]U,D,R,L,Row,Col;
int[]H;//行头节点
int[]S;//每列有多少个节点
int ansd;
int[]ans;//如果有答案,则选了ansd行,具体是哪几行放在ans[ ]数组里面,ans[0~ansd-1];

DancingLinks(){
U=new int[maxn];
D=new int[maxn];
R=new int[maxn];
L=new int[maxn];
Row=new int[maxn];
Col=new int[maxn];
H=new int[maxn];
S=new int[maxn];
ans=new int[maxn];
}
void init(int n,int m){
this.n=n;
this.m=m;
for(int i=0;i<=m;i++){
S[i]=0;
U[i]=D[i]=i;//初始状态下,上下自己指向自己
L[i]=i-1;
R[i]=i+1;
}
R[m]=0;
L[0]=m;
size=m;//编号,每列都有一个头节点,编号1-m
for(int i=1;i<=n;i++)
H[i]=-1;//每一行的头节点
}

void link(int r,int c){//第r行,第c列
++S[Col[++size]=c];//第size个节点所在的列为c,当前列的节点数++
Row[size]=r;//第size个节点行位置为r
D[size]=D[c];//下面这四句头插法(图是倒着的?)
U[D[c]]=size;
U[size]=c;
D[c]=size;
if(H[r]<0){
H[r]=L[size]=R[size]=size;
}else {
R[size]=R[H[r]];
L[R[H[r]]]=size;
L[size]=H[r];
R[H[r]]=size;
}
}

void remove(int c){//删除节点c,以及c上下节点所在的行,每次调用这个函数,都是从列头节点开始向下删除,这里c也可以理解为第c列
L[R[c]]=L[c]; //因为第c列的列头节点编号为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]){
U[D[j]]=U[j];
D[U[j]]=D[j];
--S[Col[j]];
}
}
}

void resume(int c){//恢复节点c,以及c上下节点所在的行(同上,也可以理解为从第c列的头节点开始恢复
for(int i=U[c];i!=c;i=U[i])
for(int j=L[i];j!=i;j=L[j])
++S[Col[U[D[j]]=D[U[j]]=j]]; //打这一行太纠结了 T T
L[R[c]]=R[L[c]]=c;
}

boolean dance(int d){//递归深度
if(R[0]==0){
ansd=d;
return true;
}
int c=R[0];
for(int i=R[0];i!=0;i=R[i])
if(S[i]<S[c])
c=i;
remove(c);//找到节点数最少的列,当前元素不是原图上0,1的节点,而是列头节点
for(int i=D[c];i!=c;i=D[i]){
ans[d]=Row[i];//列头节点下面的一个节点
for(int j=R[i];j!=i;j=R[j])
remove(Col[j]);
if(dance(d+1))//找到,返回
return true;
for(int j=L[i];j!=i;j=L[j])
resume(Col[j]);
}
resume(c);
return false;
}
}