这几天终于鼓起勇气学了舞蹈链这种数据结构,在联系的过程中发现一般分为两种情况
精确覆盖和重复覆盖
首先将一下精确覆盖,顾名思义就是覆盖全部的元素且只能覆盖一次
在代码实现的方面讲即给出的0_1矩阵中找到对应的行号使这些行组成组成的新矩阵每列只有一个1
模板:
#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=1010;
int L[maxn*100],R[maxn*100],U[maxn*100],D[maxn*100];//节点的上下左右四个方向的链表
int C[maxn*100],H[maxn],cnt[maxn],vis[maxn],ans[maxn],Row[maxn*100];//C列H行cnt列链表中元素个数
int n,m,id,len;
void init(){
for(int i=0;i<=m;i++){
cnt[i]=0;U[i]=D[i]=i;
L[i+1]=i;R[i]=i+1;
}
R[m]=0;id=m+1;
memset(H,-1,sizeof(H));
}
void Link(int r,int c){
cnt[c]++;C[id]=c;Row[id]=r;
U[id]=U[c];D[U[c]]=id;
D[id]=c;U[c]=id;Row[id]=r;
if(H[r]==-1) H[r]=L[id]=R[id]=id;
else{
L[id]=L[H[r]];R[L[H[r]]]=id;
R[id]=H[r];L[H[r]]=id;
}
id++;
}
void Remove(int Size){
L[R[Size]]=L[Size];
R[L[Size]]=R[Size];
for(int i=D[Size];i!=Size;i=D[i]){
for(int j=R[i];j!=i;j=R[j]){
U[D[j]]=U[j];D[U[j]]=D[j];
cnt[C[j]]--;
}
}
}
void Resume(int Size){
for(int i=D[Size];i!=Size;i=D[i]){
for(int j=R[i];j!=i;j=R[j]){
U[D[j]]=j;D[U[j]]=j;
cnt[C[j]]++;
}
}
L[R[Size]]=Size;R[L[Size]]=Size;
}
int Dance(int k){
int pos,mm=maxn;
if(R[0]==0){
len=k;
return 1;
}
for(int i=R[0];i;i=R[i]){
if(mm>cnt[i]){
mm=cnt[i];pos=i;
}
}
Remove(pos);
for(int i=D[pos];i!=pos;i=D[i]){
ans[k]=Row[i];
for(int j=R[i];j!=i;j=R[j]) Remove(C[j]);
if(Dance(k+1)) return 1;
for(int j=L[i];j!=i;j=L[j]) Resume(C[j]);
}
Resume(pos);
return 0;
}
int main()
{
//len代表个数,具体存在ans数组中
}
2 就是重复覆盖,意思也是覆盖全部元素但可以重复
同上,只不过每列最少1个1
模板:
#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3fll;
const int maxn=60;
int L[maxn*maxn],R[maxn*maxn],U[maxn*maxn],D[maxn*maxn];//节点的上下左右四个方向的链表
int C[maxn*maxn],H[maxn],cnt[maxn],vis[maxn];//C列H行cnt列链表中元素个数
int n,m,id,fans;
vector<int>G[maxn];
void init(){
for(int i=0;i<=n;i++){
G[i].clear();
G[i].push_back(i);
cnt[i]=0;U[i]=D[i]=i;
L[i+1]=i;R[i]=i+1;
}
R[n]=0;id=n+1;
memset(H,-1,sizeof(H));
}
void Link(int r,int c){
cnt[c]++;C[id]=c;
U[id]=U[c];D[U[c]]=id;
D[id]=c;U[c]=id;
if(H[r]==-1) H[r]=L[id]=R[id]=id;
else{
L[id]=L[H[r]];R[L[H[r]]]=id;
R[id]=H[r];L[H[r]]=id;
}
id++;
}
void Remove(int Size){
for(int j=D[Size];j!=Size;j=D[j])
L[R[j]]=L[j],R[L[j]]=R[j];
}
void Resume(int Size){
for(int j=D[Size];j!=Size;j=D[j])
L[R[j]]=R[L[j]]=j;
}
int h()
{
int sum=0;
memset(vis,0,sizeof(vis));
for(int i=R[0];i;i=R[i]){
if(vis[i]) continue;
sum++;
for(int j=D[i];j!=i;j=D[j]){
for(int k=R[j];k!=j;k=R[k])
vis[C[k]]=1;
}
}
return sum;
}
void Dance(int k){
int mm=maxn,pos;
if(k+h()>=fans) return;
if(!R[0]){
if(k<fans) fans=k;
return;
}
for(int i=R[0];i;i=R[i])
if(mm>cnt[i]) mm=cnt[i],pos=i;
for(int i=D[pos];i!=pos;i=D[i])
{
Remove(i);
for(int j=R[i];j!=i;j=R[j]) Remove(j);
Dance(k+1);
for(int j=R[i];j!=i;j=R[j]) Resume(j);
Resume(i);
}
}
int main(){
int u,v;
while(scanf("%d%d",&n,&m)!=-1){
init();
for(int i=0;i<m;i++){
scanf("%d%d",&u,&v);
G[u].push_back(v);
G[v].push_back(u);
}
for(int i=1;i<=n;i++){
for(int j=0;j<G[i].size();j++){
int t=G[i][j];
Link(i,t);
printf("%d %d\n",i,t);
}
}
fans=n;
Dance(0);
printf("%d\n",fans);
}
return 0;
}