【DLX】【vijos P1345】【codevs1174】舞蹈链(Dancing Links)——数独问题

最新推荐文章于 2020-08-19 21:25:11 发布
最新推荐文章于 2020-08-19 21:25:11 发布
阅读量756 收藏
点赞数
分类专栏： DLX 搜索 文章标签： vijos 数独问题 dlx
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u011327397/article/details/52760255
版权
搜索 同时被 2 个专栏收录
12 篇文章 0 订阅
订阅专栏
DLX
3 篇文章 0 订阅
订阅专栏

好久都没有写博客了，要不起来

在写之前，先安利另外一个作者写的
跳跃的舞者，舞蹈链（Dancing Links）算法——求解精确覆盖问题：http://www.cnblogs.com/grenet/p/3145800.html
算法实践——舞蹈链（Dancing Links）算法求解数独：http://www.cnblogs.com/grenet/p/3163550.html
然后vijos上有道题http://vijos.org/p/1345

我们将数独转换为精确覆盖，然后用舞蹈链来求解求精确覆盖，它虽然用了和DFS一样的方法，但却有着变态的效率。
vijos上那个题dfs也可以过，通过程序大概在3000ms~5500ms，但是用DLK来写却只要 0ms~60ms，足见效率的提升

下面是代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<set>
#include<queue>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<stack>
#define INF 2100000000
#define ll long long
#define clr(x)  memset(x,0,sizeof(x));
#define M 1000

using namespace std;

struct link2
{
    link2 *up,*down,*right,*left,*col;
    int row,count,lie;
}*head,*c[M],*link[5*M],dizhi[5*M];

int T,a[90],ans[M],sta,n;

void remove(link2 *C)
{
    C->left->right=C->right;
    C->right->left=C->left;

    link2 *i=C->down;
    while(i!=C)
    {
        link2 *j=i->right;
        while(j!=i)
        {
            j->up->down=j->down;
            j->down->up=j->up;
            j=j->right;
            j->col->count--;
        }
        i=i->down;
    }
}

void resume(link2 *C)
{
    C->left->right=C;
    C->right->left=C;

    link2 *i=C->down;
    while(i!=C)
    {
        link2 *j=i->right;
        while(j!=i)
        {
            j->up->down=j;
            j->down->up=j;
            j=j->right;
            j->col->count++;
        }
        i=i->down;
    }
}

int t;

void addline(int *index,int row)
{
    for(int i=1;i<=4;i++)
    {
        link[++t]=&dizhi[t+T];
        if(i==1)
        {
            link[t]->left=link[t];
            link[t]->right=link[t];
        }
        else
        {
            link[t]->left=link[t-1];
            link[t]->right=link[t-1]->right;
            link[t-1]->right->left=link[t];
            link[t-1]->right=link[t];
        }
        int j=index[i];
        c[j]->count++;
        link[t]->up=c[j]->up;
        c[j]->up->down=link[t];
        link[t]->down=c[j];
        c[j]->up=link[t];
        link[t]->col=c[j];
        link[t]->row=row;
        link[t]->lie=j;
    }
}

void bulid()
{
    int m=324;
    head=&dizhi[T++];
    for(int i=1;i<=m;i++)c[i]=&dizhi[T++];
    head->right=c[1];
    head->left=c[m];
    for(int i=1;i<=m;i++)
    {
        c[i]->left=i==1?head:c[i-1];
        c[i]->right=i==m?head:c[i+1];
        c[i]->up=c[i]->down=c[i];
        c[i]->count=0;
        c[i]->lie=i;
    }
    for(int i=1;i<=81;i++)
    {
        if(!a[i])continue;
        int y=(i-1)/9+1,x=(i-1)%9+1;
        int p1=i,p2=81+(y-1)*9+a[i];
        int p3=162+(x-1)*9+a[i],p4=243+((y-1)/3*3+(x-1)/3)*9+a[i];
        c[p1]->count=c[p2]->count=c[p3]->count=c[p4]->count=-1;
        c[p1]->left->right=c[p1]->right;c[p1]->right->left=c[p1]->left;
        c[p2]->left->right=c[p2]->right;c[p2]->right->left=c[p2]->left;
        c[p3]->left->right=c[p3]->right;c[p3]->right->left=c[p3]->left;
        c[p4]->left->right=c[p4]->right;c[p4]->right->left=c[p4]->left;
    }

    for(int i=1;i<=9;i++)
        for(int j=1;j<=9;j++)
            for(int k=1;k<=9;k++)
            {
                if(a[(i-1)*9+j])continue;
                int p1=(i-1)*9+j;
                int p2=81+(i-1)*9+k;
                int p3=162+(j-1)*9+k;
                int p4=243+((i-1)/3*3+(j-1)/3)*9+k;
                if(c[p1]->count==-1||c[p2]->count==-1||c[p3]->count==-1||c[p4]->count==-1)continue;
                int index[5]={0,p1,p2,p3,p4};
                addline(index,((i-1)*9+j)*10+k);
            }
}

bool dance(int x)
{
    if(head->right==head)
    {
        for(int i=1;i<=sta;i++)
            a[ans[i]/10]=ans[i]%10;
        //printf("%d:\n",n+1);
        for(int i=1;i<=81;i++)
        {
            printf("%d",a[i]);
        }

        cout<<'\n';
        return 1;
    }
    link2 *i=head->right->right;
    link2 *temp=head->right;
    while(i!=head)
    {
        if(i->count<temp->count)
            temp=i;
        i=i->right;
    }
    if(temp->down==temp)return 0;
    remove(temp);
    i=temp->down;
    while(i!=temp)
    {
        link2 *j=i->right;
        ans[++sta]=i->row;
        while(j!=i)
        {
            remove(j->col);
            j=j->right;
        }
        if(dance(x+1))return 1;
        j=i->left;
        while(j!=i)
        {
            resume(j->col);
            j=j->left;
        }
        ans[sta]=0;
        sta--;
        i=i->down;
    }
    resume(temp);
    return 0;
}

int main()
{
    cin>>n;
    while(n--)
    {
        char ch;
        getchar();
        for(int i=1;i<=81;i++)
        {
            ch=getchar();
            a[i]=ch-'0';
        }
        t=T=0;
        clr(ans);sta=0;
        clr(dizhi);clr(link);clr(c);
        bulid();
        dance(0);
    }

}

然后另外一道题
http://codevs.cn/problem/1174/
也是数独，效率我没研究，这里主要想说的是怎么讲我们刚才的求单解数独的DLK转换为多解（但是这样DLK强劲的速度就体现不出来了）
在上面，只要求出一个解，就会return 1，我们将这个1改成0就可以了

然后是这道题的代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<set>
#include<queue>
#include<algorithm>
#include<vector>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<stack>
#define INF 2100000000
#define ll long long
#define clr(x)  memset(x,0,sizeof(x));
#define M 325

using namespace std;

struct link2
{
    link2 *up,*down,*right,*left,*col;
    int row,count,lie;
}*head,*c[M],*lik[5*M],dizhi[5*M];

int T,a[90],ans[M],sta,n,Ans;

int get_num(int a)
{
    int x=((a-1)/9)+1,y=((a-1)%9)+1;
    int m1=min(x-1,9-x),m2=min(y-1,9-y),mi=min(m1,m2);
    return mi+6;
}

void remove(link2 *C)
{
    C->left->right=C->right;
    C->right->left=C->left;

    link2 *i=C->down;
    while(i!=C)
    {
        link2 *j=i->right;
        while(j!=i)
        {
            j->up->down=j->down;
            j->down->up=j->up;
            j=j->right;
            j->col->count--;
        }
        i=i->down;
    }
}

void resume(link2 *C)
{
    C->left->right=C;
    C->right->left=C;

    link2 *i=C->down;
    while(i!=C)
    {
        link2 *j=i->right;
        while(j!=i)
        {
            j->up->down=j;
            j->down->up=j;
            j=j->right;
            j->col->count++;
        }
        i=i->down;
    }
}

int t;

void addline(int *index,int row)
{
    for(int i=1;i<=4;i++)
    {
        t++;
        lik[t]=&dizhi[t+T];
        if(i==1)
        {
            lik[t]->left=lik[t];
            lik[t]->right=lik[t];
        }
        else
        {
            lik[t]->left=lik[t-1];
            lik[t]->right=lik[t-1]->right;
            lik[t-1]->right->left=lik[t];
            lik[t-1]->right=lik[t];
        }
        int j=index[i];
        c[j]->count++;
        lik[t]->up=c[j]->up;
        c[j]->up->down=lik[t];
        lik[t]->down=c[j];
        c[j]->up=lik[t];
        lik[t]->col=c[j];
        lik[t]->row=row;
        lik[t]->lie=j;
    }
}

void bulid()
{
    int m=324;
    head=&dizhi[T++];
    for(int i=1;i<=m;i++)c[i]=&dizhi[T++];
    head->right=c[1];
    head->left=c[m];
    for(int i=1;i<=m;i++)
    {
        c[i]->left=i==1?head:c[i-1];
        c[i]->right=i==m?head:c[i+1];
        c[i]->up=c[i]->down=c[i];
        c[i]->count=0;
        c[i]->lie=i;
    }
    for(int i=1;i<=81;i++)
    {
        if(!a[i])continue;
        int y=(i-1)/9+1,x=(i-1)%9+1;
        int p1=i,p2=81+(y-1)*9+a[i];
        int p3=162+(x-1)*9+a[i],p4=243+((y-1)/3*3+(x-1)/3)*9+a[i];
        c[p1]->count=c[p2]->count=c[p3]->count=c[p4]->count=-1;
        c[p1]->left->right=c[p1]->right;c[p1]->right->left=c[p1]->left;
        c[p2]->left->right=c[p2]->right;c[p2]->right->left=c[p2]->left;
        c[p3]->left->right=c[p3]->right;c[p3]->right->left=c[p3]->left;
        c[p4]->left->right=c[p4]->right;c[p4]->right->left=c[p4]->left;
    }

    for(int i=1;i<=9;i++)
        for(int j=1;j<=9;j++)
            for(int k=1;k<=9;k++)
            {
                if(a[(i-1)*9+j])continue;
                int p1=(i-1)*9+j;
                int p2=81+(i-1)*9+k;
                int p3=162+(j-1)*9+k;
                int p4=243+((i-1)/3*3+(j-1)/3)*9+k;
                if(c[p1]->count==-1||c[p2]->count==-1||c[p3]->count==-1||c[p4]->count==-1)continue;
                int index[5]={0,p1,p2,p3,p4};
                addline(index,((i-1)*9+j)*10+k);
            }
}

bool dance(int x)
{
    if(head->right==head)
    {
        for(int i=1;i<=sta;i++)
            a[ans[i]/10]=ans[i]%10;
        //printf("%d:\n",n+1);
        int temp=0;
        for(int i=1;i<=81;i++)
            temp+=get_num(i)*a[i];
        if(temp>Ans)
            Ans=temp;
        return 0;   //就是这里
    }
    link2 *i=head->right->right;
    link2 *temp=head->right;
    while(i!=head)
    {
        if(i->count<temp->count)
            temp=i;
        i=i->right;
    }
    if(temp->down==temp)return 0;
    remove(temp);
    i=temp->down;
    while(i!=temp)
    {
        link2 *j=i->right;
        ans[++sta]=i->row;
        while(j!=i)
        {
            remove(j->col);
            j=j->right;
        }
        if(dance(x+1))return 1;
        j=i->left;
        while(j!=i)
        {
            resume(j->col);
            j=j->left;
        }
        ans[sta]=0;
        sta--;
        i=i->down;
    }
    resume(temp);
    return 0;
}

int main()
{
    for(int i=1;i<=81;i++)
        scanf("%d",&a[i]);
    t=T=0;
    bulid();
    dance(0)；
    if(Ans==0)cout<<-1;
    else cout<<Ans;
}

大概就是这个样子，如果有什么问题，或错误，请在评论区提出，谢谢。

anantheparty
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
TCP Ports list (3498 ports in list)
zyb378747350的专栏
01-11 41万+
参考地址:https://www.gasmi.net/tcp.php TCP Ports list (3498 ports in list) 1 tcpmux TCP Port Service Multiplexer 2 compressnet Management Utility 3 compressnet Compression Process
洛谷P4929 【模板】舞蹈链（DLX）
悲伤无声地涌出来
08-20 448
本题是舞蹈链模板——精确覆盖问题
舞蹈链数独求解算法
04-10
舞蹈链算法求解数独
用舞蹈链解数独问题
零昔无浪
02-23 1513
转载自：http://fayaa.com/code/view/1738/ 把debug文件夹中的exe文件放在指定路径，这里放在F盘 然后打开cmd文件输入 F:\53.exe 030000000010200060000600400004000800007000000000001003000460005000800000050000001 输出结果 数独测试案例 0
洛谷 P1074 靶形数独 —— 舞蹈链（DLX）精确覆盖
天翼之城的博客
07-29 222
This way 题意： 给你一个9*9，并且至少有24个位置已经填上的数独矩阵，让你去补全这个数独。要求每一个格子只能有一个数，同一个数只能在一行，一列，一个九宫格中出现一次。并且它是带有计分的数独： 填在相应格子里的数对于答案的贡献是数值乘上相应的权值。问你最终答案最大是多少。 题解： 那么首先9*9的数独在上一篇博客里面已经说了： This way 9*9的格子中数独的解法有6.6E21种，这个已经有证明了。 但是这道题它有一个很重要的限制就是已经填了24个数，虽然我无法证明之后的复杂度，但是
dancing links 骑士遍历
07-10
dancing links 解决骑士遍历问题，实验报告，含源代码
DLX_Links_DLX_
09-28
标题中的"DLX_Links_DLX_"暗示了我们讨论的主题是关于一种名为"Dancing Links"（跳舞链）的算法，通常简称为DLX。这个算法是由Donald Knuth提出的一种高效解决完全背包问题、数独问题和其他一系列组合优化问题的方法...
sudoku:使用Dancing Links（DLX）的Sudoku求解器
05-09
使用Dancing Links（DLX）的Sudoku求解器 数独可以简化为精确的掩护问题，该问题被认为是NP完全的。 NP-complete的分类仅适用于广义nxn数独，而不适用于9x9数独，因为它是有限实例。 确切的掩护问题是一个决策问题...
基于舞蹈链的数独算法，以Python实现
最新发布
12-26
舞蹈链算法，又称Dancing Links（DLX）算法，是由日本数学家和计算机科学家Donald Knuth提出的一种用于解决完全多归类问题的有效方法。在数独求解中，我们可以将每个空格看作一个“列”，每个可能填入的数字看作一个...
dlx.js:Knuth的Dancing Links与Algorithm X JS实现
05-13
Dancing Links是一种解决精确覆盖问题的有效算法，它在计算机科学中主要用于求解组合优化问题，如回文分割、集合完全覆盖等。 **描述详解：** 描述中提到的"Donald Knuth的Dancing Links与算法X（DLX）JavaScript...
【NOIP2009】靶形数独 DLX（Dancing Links）
辗转山河弋流歌
10-30 2139
[NOIP2009]靶形数独 T4 Time Limit: 2 Sec  Memory Limit: 128 MB Description 　　小城和小华都是热爱数学的好学生，最近，他们不约而同地迷上了数独游戏，好胜的他们想用数独来一比高低。但普通的数独对他们来说都过于简单了，于是他们向Z 博士请教，Z 博士拿出了他最近发明的“靶形数独”，作为这两个孩子比试的题目。 　　靶形数
P4929 【模板】舞蹈链（DLX）
aishang6250的博客
04-26 148
题目背景 本题是舞蹈链模板——精确覆盖问题 题目描述 给定一个N行M列的矩阵，矩阵中每个元素要么是1，要么是0你需要在矩阵中挑选出若干行，使得对于矩阵的每一列j，在你挑选的这些行中，有且仅有一行的第j个元素为1 输入输出格式 输入格式： 第一行两个数N,M 接下来N行，每行M个数字0或1，表示这个矩阵 输出格式： 一行输出若干...
靶形数独
Uhen's Blog
03-24 497
靶形数独原题地址：洛谷（听说数独是一个NP完全问题） 本题使用搜索，且需要一定的优化，从网上百度到了以下优化：剪掉不可行解为了提高效率，可以像八皇后一样，开三个数组，记录行列以及九宫格的状态，只有在三个数组中都未出现的值才可以被填入。 可以参考以下实现：//定义数组，以及判断格子是否可行，设置格子已被占用及回溯的方法（使用宏定义） #define get_block(a, b) ((a/3*3)
DLX
wen_xp的专栏
12-23 456
// hdu 2828 // poj 3740 int mx[N][M]; int U[nm],D[nm],L[nm],R[nm],X[nm],Y[nm]; int S[M],H[N]; int sz,ans; void ver(int x,int y){///上 下 D[x]=y;U[y]=x; } void hor(int x,int y){////左 右 R[x]=y;L[y]
【Dancing Links舞蹈链】poj 3076 Sudoku
Hacker_vision
09-03 928
Source : poj 3076 Sudoku http://poj.org/problem?id=3076 Problem Description A Sudoku grid is a 16x16 grid of cells grouped in sixteen 4x4 squares, where some cells are filled with letters from A to
洛谷 P4929 【模板】舞蹈链（DLX）
wineandchord
08-19 183
精确覆盖问题 使用 dancing links 本质思想是十字链表 // P4929 模板 舞蹈链 DLX #include<iostream> #include<cstdio> #define MAXN 510 #define MAXM 510 #define MAX (MAXN*MAXM) using namespace std; int n,m; struct DLX{ int n,m,sz; // n为行数，m为列数，sz表示总节点数 int U[MAX].
C++求解数独（舞蹈链）算法实例「超详细步骤+学习文档」
12-23 1180
求解数独，包括了数独所有解和舞蹈链的源码文件，源码主要是数独求解（舞蹈链算法）由于回朔法易于理解，不做具体介绍（便于求出所有可行解），本次仅具体介绍舞蹈链算法。 　　程序设计采用了数组代替十字交叉链表。主要步骤为：问题转化-->舞蹈链算法求解-->转化输出可行解。 　　回溯法与舞蹈链算法求解同一问题时间对比，c++实现，存在误差（操作系统自身，输入输出），但...
poj (3074：数独)Dancing Link算法神模版
Weiguang_123的专栏
09-02 4602
http://www.ocf.berkeley.edu/ jchu/publicportal/sudoku/0011047.pdf 题意：玩家需要根据9×9盘面上的已知数字，推理出所有剩余空格的数字，并满足每一行、每一列、每一个粗线宫内的数字均含1-9，不重复。 /* DLX解决9*9的数独问题，转化为729*324的精确覆盖问题 行： 一共9 * 9 * 9 == 729行。一共9 *
ZOJ-3209 Treasure Map（dancing links，舞蹈链算法，DLX）
int_Brosea的博客
08-09 419
ZOJ - 3209 题目描述 Your boss once had got many copies of a treasure map. Unfortunately, all the copies are now broken to many rectangular pieces, and what make it worse, he has lost some of the pieces. L...
用舞蹈链算法（Dancing Links）求解俄罗斯方块覆盖问题matlab代码
06-07
dl = DancingLinks(max_nodes); % 构建 Dancing Links 矩阵 for i = 1:m for j = 1:n for k = 1:numel(blocks) block = blocks{k}; if i + size(block, 1) - 1 > m || j + size(block, 2) - 1 > n continue; ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值