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

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; // 找到了就一直return
    if(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++;} //如果目标已经使用并且相应位置是1
    if(count == 1) dfs(num + 1, in, out | cnt[num] ^ bit[j]); //则接着搜索时排除除了该1位置以外的1
    else 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 LOCAL
    freopen("ride.in","r",stdin);
    freopen("ride.out","w",stdout);
#endif
    int 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;
}