当前位置: 代码迷 >> 综合 >> POJ 3740 Easy Finding 位运算压缩+DFS or Dancing Links
  详细解决方案

POJ 3740 Easy Finding 位运算压缩+DFS or Dancing Links

热度:14   发布时间:2024-01-13 18:22:57.0

Easy Finding

这道题运用二进制DFS是比较好的方法,当然还可以用dancing links,这个比较高级,先把我看别人的二进制思路然后写的代码贴上来。

当然本题是体现不出DLX的高级之处的。用位运算和DLX的速度应该是相当

/*
ID: sdj22251
PROG: calfflac
LANG: C++
*/
#include <iostream>
#include <vector>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <queue>
#include <stack>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>
#include <sstream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <string>
#include <cstring>
#include <cmath>
#include <ctime>
#define MAX 2000000000
#define LOCA
#define PI acos(-1.0)
using namespace std;
int n, m, bit[18], a[18][303], cnt[303];
bool flag;
void dfs(int num, int in, int out) //in是已经使用的,out是排除了的,num是深度
{int i, j;if(flag) return; // 找到了就一直returnif(num >= n) {flag = true; return;} //找到了if(cnt[num] == (cnt[num] & out)) return; //如果深搜的目标都在排除的序列中,就return int count = 0;for(i = 0; i < m && count < 2; i++)if(cnt[num] & in & bit[i]) {j = i; count++;} //如果目标已经使用并且相应位置是1if(count == 1) dfs(num + 1, in, out | cnt[num] ^ bit[j]); //则接着搜索时排除除了该1位置以外的1else if (count == 0) //还没被使用过{for(i = 0; i < m && !flag; i++)if((cnt[num] & bit[i]) && !(out & bit[i])) // 相应位置是1 并且还没被排除过dfs(num + 1, in | bit[i], out | cnt[num] ^ bit[i]); }
}
int main()
{
#ifdef LOCALfreopen("ride.in","r",stdin);freopen("ride.out","w",stdout);
#endifint i, j;bit[0] = 1;for(i = 1; i < 18; i++) //将不同位置的1预存起来bit[i] = bit[i - 1] << 1; while(scanf("%d%d", &m, &n) != EOF){memset(cnt, 0, sizeof(cnt));for(i = 0; i < m; i++){for(j = 0; j < n; j++){scanf("%d", &a[i][j]);if(a[i][j])cnt[j] |= bit[i]; //如果这一个点是1,就将1及其位置存进二进制数组}}flag = false;for(i = 0; i < n; i++) //筛掉全是0这个情况if(!cnt[i]) break;if(i >= n) dfs(0, 0, 0);if(flag) puts("Yes, I found it");else puts("It is impossible");}return 0;
}


2012年1月

然后就是传说中的DLX了。

看了网上很多的资料,只是建图的时候比较麻烦吧

每一列先建立一个列头,然后往某一列中插入元素时只需要插到列头下面就行,因为是循环链表插到哪里都一样的

而行呢,不一定要行头,但是最好能把每一行的第一个1的编号记录下来,每次插入到它之后就行了。


/*
ID: sdj22251
PROG: inflate
LANG: C++
*/
#include <iostream>
#include <vector>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <queue>
#include <stack>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>
#include <sstream>
#include <iomanip>
#include <cstdio>
#include <cmath>
#include <cstdlib>
#include <cctype>
#include <string>
#include <cstring>
#include <cmath>
#include <ctime>
#define MAXN 5000
#define INF 1000000000
using namespace std;
int L[MAXN], R[MAXN], C[MAXN], S[MAXN], U[MAXN], D[MAXN], H[MAXN];//H数组用来存储每一行的第一个1的编号,C数组用来存列号,S数组用来存每一列中1的个数 
int n, m, cnt, head;//head表示列头,一般取0,cnt用来记录编号
void link(int r, int c)
{S[c]++;C[cnt] = c;U[cnt] = c;D[cnt] = D[c];U[D[c]] = cnt;D[c] = cnt;if(H[r] == -1){H[r] = cnt;L[cnt] = R[cnt] = cnt;}else{L[cnt] = H[r];R[cnt] = R[H[r]];L[R[H[r]]] = cnt;R[H[r]] = cnt;}cnt++;
}
void init()
{cnt = 0;head = 0;for(int i = 0; i <= m; i++){S[i] = 0;D[i] = U[i] = i;R[i] = (i + 1) % (m + 1);L[i] = (i + m) % (m + 1);cnt++;}memset(H, -1, sizeof(H));
}
void readdata()
{char t;getchar();for(int i = 1; i <= n; i++)for(int j = 1; j <= m; j++){t = getchar();getchar();if(t == '1') link(i, j);}
}
void removes(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]){U[D[j]] = U[j];D[U[j]] = D[j];S[C[j]]--;}
}
void resumes(int c)
{for(int i = U[c]; i != c; i = U[i])for(int j = L[i]; j != i; j = L[j]){U[D[j]] = j;D[U[j]] = j;S[C[j]]++;}L[R[c]] = c;R[L[c]] = c;
}
bool dfs(int k)
{if(R[head] == head) return true;int s = INF, c;for(int i = R[head]; i != head; i = R[i])if(s > S[i]){s = S[i];c = i;}removes(c);for(int i = U[c]; i != c; i = U[i]){for(int j = R[i]; j != i; j = R[j])removes(C[j]);if(dfs(k + 1)) return true;for(int j = L[i]; j != i; j = L[j])resumes(C[j]);}resumes(c);return false;
}
int main()
{while(scanf("%d%d", &n, &m) != EOF){init();readdata();if(dfs(0)) puts("Yes, I found it");else puts("It is impossible");}return 0;
}


  相关解决方案