POJ 2676 3074 3076 数独3题 Dancing Links

最新推荐文章于 2019-08-15 15:24:54 发布
最新推荐文章于 2019-08-15 15:24:54 发布
阅读量112 收藏
点赞数
文章标签： c/c++ 数据结构与算法

今天就一直在做DLX的题目了，目前也只是入门而已，就做了几个数独的问题

将数独转化为DLX的模型是很繁琐的一部分。

首先，必须明白，有多少种选择就要建立多少行，有多少个条件约束就要建立多少个列

对于一个9*9的数独，我们可以在81个位置中选择一个地方，放1~9的数，所以选择数目就是9*9*9了，就要建立这么些个行

而约束呢，首先是81个格子，每个格子有且只能有1个数，那么就是9*9个列，然后对于某个位置的某个数，同一行的不能再有一样的数，也是9*9个列，同一列，同一宫都是一个原理，所以总共就是4*9*9个列了

然后问题就转化为了，从9*9*9行中选择某些行，使得每一列都只有1个1，这里的1代表的就是唯一性，比如每一宫中只能有1个1，如果你选的两行致使1个宫内有两个1，那么显然这是不符合约束的，必然在某一列中这个冲突会得到体现。同样的是同一列和同一行。

对于每个格子有且只能有1个数这个条件，也是不可少的。虽说插入的时候只要放置的位置不同，一定不会产生冲突，但是加入此条件就能保证最后输出的结果，每一格都充满了数字。

然后建图的时候就是按着行和列用十字循环链表建图了。


怎么说呢，今天这三道题做的，感觉学了不少东西，但是写起来，却又写挫了，16*16那个数独也不知他们怎么写的，我花了300+ms，前边一堆100ms以内的

首先是2676 

/*
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 9*9*9*9
#define INF 1000000000
#define N 9
#define M 9*9*9+5
using namespace std;
int L[MAXN], R[MAXN], C[MAXN], S[M], U[MAXN], D[MAXN], H[M], O[M], X[MAXN]; //X用来存储行，C用来存储列，H用来存储每行中的第一个结点，O用来存储结果
int cnt, head;
char mp[N + 5][N + 5], ans[N * N + 5];
bool vis[N * N * 4 + 5];
void link(int r, int c)
{
    S[c]++;
    C[cnt] = c;
    X[cnt] = r;
    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 <= N * N * 4; i++)
    {
        S[i] = 0;
        vis[i] = 0;
        D[i] = U[i] = i;
        R[i] = (i + 1) % (N * N * 4 + 1);
        L[i] = (i + N * N * 4) % (N * N * 4 + 1);
        cnt++;
    }
    memset(H, -1, sizeof(H));
}
void cal(int &r, int &cx, int &cy, int &ck, int &cg, int i, int j, int k)
{
    r = (i * N + j) * N + k - 1; //代表所属的行 
    cg = i * N + j + 1;   //代表的是数独中i,j位置所属的列
    cx = N * N + i * N + k;  //代表数独中同一行所属的列
    cy = N * N * 2 + j * N + k; //代表数独中同一列所属的列
    int LN = (int)sqrt(N * 1.0); 
    ck = N * N * 3 + (i / LN * LN + j / LN) * N + k;//代表数独中同一宫所属的列
}
void readdata()
{
    int r, cx, cy, ck, cg;
    for(int i = 0; i < N; i++)
        scanf("%s", mp[i]);
    for(int i = 0; i < N; i++)
        for(int j = 0; j < N; j++)
            if(mp[i][j] != '0')
            {
                cal(r, cx, cy, ck, cg, i, j, mp[i][j] - '0');
                link(r, cx);
                link(r, cy);
                link(r, ck);
                link(r, cg);
                vis[cx] = vis[cy] = vis[ck] = vis[cg] = 1;
            }
    for(int i = 0; i < N; i++)
        for(int j = 0; j < N; j++)
            for(int k = 1; k <= N; k++)
            {
                cal(r, cx, cy, ck, cg, i, j, k);
                if(vis[cx] || vis[cy] || vis[ck] || vis[cg]) continue;
                link(r, cx);
                link(r, cy);
                link(r, ck);
                link(r, cg);
            }
}
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)
    {
        for(int i = 0; i < k; i++) //根据所选择的行 能推出选择的位置和数字
            ans[X[O[i]] / N] = X[O[i]] % 9 + '1';
        for(int i = 0; i < N * N; i++)
        {
            printf("%c", ans[i]);
            if((i + 1) % N == 0) printf("\n");
        }
        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])
    {
        O[k] = 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()
{
    int T;
    scanf("%d", &T);
    while(T--)
    {
        init();
        readdata();
        dfs(0);
    }
    return 0;
}

然后是3074 


/*
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 9*9*9*9
#define INF 1000000000
#define N 9
#define M 9*9*9+5
using namespace std;
int L[MAXN], R[MAXN], C[MAXN], S[M], U[MAXN], D[MAXN], H[M], O[M], X[MAXN];
int cnt, head;
char mp[N * N + 5], ans[N * N + 5];
bool vis[N * N * 4 + 5];
void link(int r, int c)
{
    S[c]++;
    C[cnt] = c;
    X[cnt] = r;
    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 <= N * N * 4; i++)
    {
        S[i] = 0;
        vis[i] = 0;
        D[i] = U[i] = i;
        R[i] = (i + 1) % (N * N * 4 + 1);
        L[i] = (i + N * N * 4) % (N * N * 4 + 1);
        cnt++;
    }
    memset(H, -1, sizeof(H));
}
void cal(int &r, int &cx, int &cy, int &ck, int &cg, int i, int j, int k)
{
    r = (i * N + j) * N + k - 1;
    cg = i * N + j + 1;
    cx = N * N + i * N + k;
    cy = N * N * 2 + j * N + k;
    int LN = (int)sqrt(N * 1.0);
    ck = N * N * 3 + (i / LN * LN + j / LN) * N + k;
}
void readdata()
{
    int r, cx, cy, ck, cg;
    for(int i = 0; i < N; i++)
        for(int j = 0; j < N; j++)
            if(mp[i * N + j] != '.')
            {
                cal(r, cx, cy, ck, cg, i, j, mp[i * N + j] - '0');
                link(r, cx);
                link(r, cy);
                link(r, ck);
                link(r, cg);
                vis[cx] = vis[cy] = vis[ck] = vis[cg] = 1;
            }
    for(int i = 0; i < N; i++)
        for(int j = 0; j < N; j++)
            if(mp[i * N + j] == '.')
            for(int k = 1; k <= N; k++)
            {
                cal(r, cx, cy, ck, cg, i, j, k);
                if(vis[cx] || vis[cy] || vis[ck] || vis[cg]) continue;
                link(r, cx);
                link(r, cy);
                link(r, ck);
                link(r, cg);
            }
}
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)
    {
        for(int i = 0; i < k; i++)
            ans[X[O[i]] / N] = X[O[i]] % N + '1';
        for(int i = 0; i < N * N; i++)
        {
            printf("%c", ans[i]);
        }
        printf("\n");
        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])
    {
        O[k] = 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("%s", mp) != EOF && mp[0] != 'e')
    {
        init();
        readdata();
        dfs(0);
    }
    return 0;
}

最后是3076

这是16*16的数独

/*
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 16*16*16*16+16*16*4
#define INF 1000000000
#define N 16
#define M 16*16*16+5
using namespace std;
int L[MAXN], R[MAXN], C[MAXN], S[M], U[MAXN], D[MAXN], H[M], O[M], X[MAXN];
int cnt, head;
char mp[N + 5][N + 5], ans[N * N + 5];
bool vis[N * N * 4 + 5];
void link(int r, int c)
{
    S[c]++;
    C[cnt] = c;
    X[cnt] = r;
    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 <= N * N * 4; i++)
    {
        S[i] = 0;
        vis[i] = 0;
        D[i] = U[i] = i;
        R[i] = (i + 1) % (N * N * 4 + 1);
        L[i] = (i + N * N * 4) % (N * N * 4 + 1);
        cnt++;
    }
    memset(H, -1, sizeof(H));
}
void cal(int &r, int &cx, int &cy, int &ck, int &cg, int i, int j, int k)
{
    r = (i * N + j) * N + k - 1;
    cg = i * N + j + 1;
    cx = N * N + i * N + k;
    cy = N * N * 2 + j * N + k;
    int LN = (int)sqrt(N * 1.0);
    ck = N * N * 3 + (i / LN * LN + j / LN) * N + k;
}
void readdata()
{
    int r, cx, cy, ck, cg;
    for(int i = 1; i < N; i++)
        scanf("%s", mp[i]);
    for(int i = 0; i < N; i++)
        for(int j = 0; j < N; j++)
            if(mp[i][j] != '-')
            {
                cal(r, cx, cy, ck, cg, i, j, mp[i][j] - 'A' + 1);
                link(r, cx);
                link(r, cy);
                link(r, ck);
                link(r, cg);
                vis[cx] = vis[cy] = vis[ck] = vis[cg] = 1;
            }
    for(int i = 0; i < N; i++)
        for(int j = 0; j < N; j++)
            if(mp[i][j] == '-')
            for(int k = 1; k <= N; k++)
            {
                cal(r, cx, cy, ck, cg, i, j, k);
                if(vis[cx] || vis[cy] || vis[ck] || vis[cg]) continue;
                link(r, cx);
                link(r, cy);
                link(r, ck);
                link(r, cg);
            }
}
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)
    {
        for(int i = 0; i < k; i++)
            ans[X[O[i]] / N] = X[O[i]] % N + 'A';
        for(int i = 0; i < N * N; i++)
        {
            printf("%c", ans[i]);
            if((i + 1) % N == 0) printf("\n");
        }
        printf("\n");
        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])
    {
        O[k] = 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("%s", mp[0]) != EOF)
    {
        init();
        readdata();
        dfs(0);
    }
    return 0;
}


iteye_6233
关注
Dancing-Link+数独版
Luckfort
07-28 447
Dancing-Link ＆＆ Sodoku
【POJ 2676】Sudoku数独（DFS）
Chen_yuazzy的博客
07-28 1302
思路：DFS遍历解数独，很有意思的题目，由于行、列、每一个3*3小方格都不能重复，故需要三个数组row,col,dig来标记某1~9数字是否出现。DFS遍历，具体思路见注释。
java数独
ddige1793的博客
11-15 97
1.比较清晰，就是递归的求下一个元素是否能成功，这种肯定不是最优的。 http://www.cnblogs.com/lgfeng/p/3255566.html 2.Dancing Link解数独： Dancing Link是用来解精准覆盖的： http://www.cnblogs.com/grenet/p/3145800.html 这个是介绍的比较好的，还在研读中…… http...
(转)数独之Dancing Link解法
dreamhunter11770的专栏
11-07 1220
(转)数独之Dancing Link解法 默认分类 2009-08-24 19:04:14 阅读575 评论0 字号：大中小 订阅 Dancing Link算法(以下简称DLX)是解NPC难题中的精确覆盖(Exact Cover)的高效算法，一个问题，如果能转化成Exact Cover模型，则都能用DLX解。数独的解法也不列外。 对于一个N*N的(N=K*K)数独，我们可以用一个3位的N进制数rck(0 using namespace std; #define RR 7
HDU 1426 dancing links解决数独问题
weixin_30699463的博客
05-18 101
题目大意： 这是一个最简单的数独填充题目，题目保证只能产生一种数独，所以这里的初始9宫格较为稠密，可以直接dfs也没有问题 但最近练习dancinglinks，这类数据结构解决数独无疑效率会高很多 dancinglinks的数独限制条件是： 1.每行有9个元素，共9行对应dlx81列 2.每列有9个元素，共9行对应dlx81列 3.每个九宫格有9个元素，共9行对应...
舞蹈链 （dancing links X）
rgndao的博客
08-15 302
舞蹈链算法，用于求解精确覆盖或者重复覆盖问题，目的是保证每一列的值最后为1或大于1，搜索中选用某一行后就将这一行以及它所解决的所有列删去，最后如果发现有列无法删去则失败，回溯并将之前删去的都恢复；如果最后发现所有列都被删去了，则说明问题解决。 它首先可以用来解决数独问题：转化为每行，每列，每个宫都有1~9，每格只有一个数字这个条件。是精确覆盖问题。 HDU3111题目链接 #include...
POJ 2676 代码 搜索数独答案
09-15
简单搜索题 数独 答案 POJ 2676 也可以没事玩玩数独。
POJ2676-Sudoku
07-31
在提供的文件“POJ2676-Sudoku.cpp”中，应该包含了上述算法的C++实现，而“POJ2676-Sudoku.doc”可能包含了解题思路、代码解释或者测试数据的分析。详细阅读这两个文件，可以帮助我们更深入地理解这个问题的解决...
poj (3074：数独)Dancing Link算法神模版
Weiguang_123的专栏
09-02 4718
http://www.ocf.berkeley.edu/ jchu/publicportal/sudoku/0011047.pdf 题意：玩家需要根据9×9盘面上的已知数字，推理出所有剩余空格的数字，并满足每一行、每一列、每一个粗线宫内的数字均含1-9，不重复。 /* DLX解决9*9的数独问题，转化为729*324的精确覆盖问题 行： 一共9 * 9 * 9 == 729行。一共9 *
poj3074/3076 Dancing Links 解数独
10-06 2758
Dancing Links经典应用
一、Dancing Links：完全覆盖算法
李勇的专栏
06-08 1082
 1 问题描述 考虑一张棋盘，被分成6×6（6行6列）的36个正方形。设有形状一样的多米诺牌，每张牌恰好覆盖棋盘上相邻的两个方格。若把18张多米诺牌摆放在棋盘上，任何两张多米诺牌均不重叠，每张多米诺牌覆盖两个方格，并且棋盘上的所有的方格均被覆盖住，这样一种排列称为多米诺牌的完美覆盖。请问：6×6的棋盘有多少种不同的完美覆盖？ 2 功能要求 1) 良好的用户界面。要求能够显示所有的完美覆
POJ2676:Sudoku(DFS)
weixin_30835923的博客
11-29 119
Sudoku Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 10849 Accepted: 5414 Special Judge Description Sudoku is a very simple task. A square table with 9...
POJ 2676 解题报告
thestoryofsnow的专栏
10-02 463
这道题是DFS。看了discuss，说是倒着搜。从8,8到0,0，应该是根据测试数据来的。 thestoryofsnow 2676 Accepted 132K 16MS C++ 2172B /* ID: thestor1 LANG: C++ TASK: poj2676 */ #include #include #include #inclu
poj2676 - Sudoku
贰圣
04-14 2134
想看更多的解题报告: http://blog.csdn.net/wangjian8006/article/details/7870410                                      转载请注明出处：http://blog.csdn.net/wangjian8006 题目大意：                将数独补充完整 代码： #include usi
POJ2676
uyret的博客
05-24 541
数独题，9行9列，要求每个小九宫格中的数字都不一样 并且，每行每列数字都不一样，0代表没有填入数字。  #include int num[9][9],f;  int check(int n, int key) {      int  x,y,i,j;               for (j = 0; j     {           i =
poj3076(16*16数独)
lch的博客
09-01 1720
转载自：http://blog.csdn.net/lyhypacm/article/details/5923287 题意： 一个16*16的数独问题。 思路： 和poj3024一样，就改点东西。 代码： #include #include #include #include using namespace std; const int oo=0x3f
初次学习舞蹈链Dancing Links
qq_36386435的博客
10-01 704
介绍几个讲的比较好的博客： 这篇博客容易理解，可以先看： http://www.cnblogs.com/grenet/p/3145800.html 这篇博客阐述的也好，而且有实例分析（也是转载自此文）： https://blog.csdn.net/WhereIsHeroFrom/article/details/79220897?utm_source=copy Dancing Links正是十字交叉...
DFS-Sudoku（poj2676）
劫劫劫争
12-08 127
#include &lt;iostream&gt; #include &lt;cstring&gt; #include &lt;cstdio&gt; using namespace std; int map[10][10], hang[10][10], lie[10][10], kuai[10][10], cnt; int cal(int i, int j)//输入行列号，判断他在哪一块里 { ...
POJ编程挑战：精选50题解析与实战
"这是一个关于POJ在线判题系统推荐的50道编程题目，旨在帮助学习者提升编程和算法能力。" 这些题目涵盖了多种算法和编程技巧，是锻炼和提高编程能力的良好实践。以下是对每道题目特点和所涉及知识点的简要介绍： 1...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值