POJ 3074&&3076 舞蹈链

最新推荐文章于 2021-10-09 16:30:12 发布
最新推荐文章于 2021-10-09 16:30:12 发布
阅读量2.1k 收藏
点赞数
分类专栏： 数据结构 舞蹈链 线段树 文章标签： ACM 数据结构 舞蹈链 poj
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Dan__ge/article/details/51758690
版权
线段树 同时被 3 个专栏收录
207 篇文章 9 订阅
订阅专栏
数据结构
146 篇文章 0 订阅
订阅专栏
舞蹈链
9 篇文章 0 订阅
订阅专栏

点击打开链接

题意：完成数独

思路：直接用舞蹈链的不可重复覆盖直接就可以过了，建图时列为81+81+81+81，第一个为81个空都要有元素，然后第二个是9行，有9个数要添，第三个就是9列有9个元素要添，最后一个是9块有9个元素要添，然后建完模版OK

#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*400],R[maxn*400],U[maxn*400],D[maxn*400];//节点的上下左右四个方向的链表
int C[maxn*400],H[maxn],cnt[maxn],ans[maxn],Row[maxn*400];//C列H行cnt列链表中元素个数
int n,m,id,len;
void init(int lll){
    for(int i=0;i<=lll;i++){
        cnt[i]=0;U[i]=D[i]=i;
        L[i+1]=i;R[i]=i+1;
    }
    R[lll]=0;id=lll+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;
    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 c;
     if(!R[0]){
         len=k;
         return 1;
     }
     for(int Min=inf,i=R[0];i;i=R[i])
         if(Min>cnt[i]) Min=cnt[i],c=i;
     Remove(c);
     for(int i=D[c];i!=c;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(c);
    return 0;
}
char str[110];
int X[maxn],Y[maxn],VAL[maxn];
int astr[110];
int main(){
    while(scanf("%s",str)!=-1){
        if(strcmp(str,"end")==0) break;
        init(324);
        int kkk=1;
        for(int i=1;i<=9;i++){
            for(int j=1;j<=9;j++){
                if(str[(i-1)*9+j-1]=='.'){
                    int pos=(i-1)/3*3+j/3;
                    if(j%3!=0) pos++;
                    for(int l=1;l<=9;l++){
                        X[kkk]=i;Y[kkk]=j;VAL[kkk]=l;
                        Link(kkk,(i-1)*9+j);
                        Link(kkk,81+(i-1)*9+l);
                        Link(kkk,162+(j-1)*9+l);
                        Link(kkk++,243+(pos-1)*9+l);
                    }
                }else{
                    int pos=(i-1)/3*3+j/3;
                    if(j%3!=0) pos++;
                    int l=str[(i-1)*9+j-1]-'0';
                    X[kkk]=i;Y[kkk]=j;VAL[kkk]=l;
                    Link(kkk,(i-1)*9+j);
                    Link(kkk,81+(i-1)*9+l);
                    Link(kkk,162+(j-1)*9+l);
                    Link(kkk++,243+(pos-1)*9+l);
                }
            }
        }
        int fff=Dance(0);
        for(int i=0;i<len;i++){
            int lll=(X[ans[i]]-1)*9+Y[ans[i]]-1;
            astr[lll]=VAL[ans[i]];
        }
        for(int i=0;i<81;i++) printf("%d",astr[i]);
        printf("\n");
    }
    return 0;
}
点击打开链接

题意：与3074一样完成数独

思路：改一下数据范围什么的简单处理一下就可以了

#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=5010;
int L[maxn*1010],R[maxn*1010],U[maxn*1010],D[maxn*1010];//节点的上下左右四个方向的链表
int C[maxn*1010],H[maxn],cnt[maxn],ans[maxn],Row[maxn*1010];//C列H行cnt列链表中元素个数
int n,m,id,len;
void init(int lll){
    for(int i=0;i<=lll;i++){
        cnt[i]=0;U[i]=D[i]=i;
        L[i+1]=i;R[i]=i+1;
    }
    R[lll]=0;id=lll+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;
    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 c;
     if(!R[0]){
         len=k;
         return 1;
     }
     for(int Min=inf,i=R[0];i;i=R[i])
         if(Min>cnt[i]) Min=cnt[i],c=i;
     Remove(c);
     for(int i=D[c];i!=c;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(c);
    return 0;
}
char str[20][20];
int X[maxn],Y[maxn],VAL[maxn];
int astr[110];
int main(){
    while(scanf("%s",str[0])!=-1){
        for(int i=1;i<16;i++) scanf("%s",str[i]);
        init(1024);
        int kkk=1;
        for(int i=1;i<=16;i++){
            for(int j=1;j<=16;j++){
                if(str[i-1][j-1]=='-'){
                    int pos=(i-1)/4*4+j/4;
                    if(j%4!=0) pos++;
                    for(int l=1;l<=16;l++){
                        X[kkk]=i;Y[kkk]=j;VAL[kkk]=l;
                        Link(kkk,(i-1)*16+j);
                        Link(kkk,256+(i-1)*16+l);
                        Link(kkk,512+(j-1)*16+l);
                        Link(kkk++,768+(pos-1)*16+l);
                    }
                }else{
                    int pos=(i-1)/4*4+j/4;
                    if(j%4!=0) pos++;
                    int l=str[i-1][j-1]-'A'+1;
                    X[kkk]=i;Y[kkk]=j;VAL[kkk]=l;
                    Link(kkk,(i-1)*16+j);
                    Link(kkk,256+(i-1)*16+l);
                    Link(kkk,512+(j-1)*16+l);
                    Link(kkk++,768+(pos-1)*16+l);
                }
            }
        }
        int fff=Dance(0);
        for(int i=0;i<len;i++){
            int lll=(X[ans[i]]-1)*16+Y[ans[i]]-1;
            astr[lll]=VAL[ans[i]];
        }
        for(int i=0;i<256;){
            int lllll=16;
            while(lllll--){
                char ch=astr[i++]+'A'-1;
                printf("%c",ch);
            }
            printf("\n");
        }
        printf("\n");
    }
    return 0;
}

Dan__ge
关注
专栏目录
poj (3074：数独)Dancing Link算法神模版
Weiguang_123的专栏
09-02 4703
http://www.ocf.berkeley.edu/ jchu/publicportal/sudoku/0011047.pdf 题意：玩家需要根据9×9盘面上的已知数字，推理出所有剩余空格的数字，并满足每一行、每一列、每一个粗线宫内的数字均含1-9，不重复。 /* DLX解决9*9的数独问题，转化为729*324的精确覆盖问题 行： 一共9 * 9 * 9 == 729行。一共9 *
【POJ1101&OpenJudge2802】The Game 小游戏（BFS+记忆化）
Cassie_zkq的博客
05-24 685
传送门 题目： 解题思路： 注意题里的卡片位置(x,y)对应的是(列，行)，可以将坐标转换成我们熟悉的(行，列)形式。图按常规方法存，只需要将要查询的(x,y)换成(y,x)即可。 dis[x][y]表示从(x1,y1)到(x,y)所需要的最短线段数，在bfs的时候更新这个值。由于dis[x2][y2]可以由多个位置和多个方向更新来，队列中应该记录{x,y,dir,step}，并且允许元素多次入队，但每个(x,y)只处理一次。 网上的博客好多都是错的，假代码没过自己的样例，但AC了。。数据挺假的，题也
poj3074 Sudoku--舞蹈链数独
半壶老酒
09-25 942
原题链接：http://poj.org/problem?id=3074 题意：给定一个9*9的数独，求解。 #define _CRT_SECURE_NO_DEPRECATE #include #include #include #include #include #include #include #define INF 99999999 #define eps 0.000
POJ - 3076 Sudoku （舞蹈链模板）
LP_Cong
04-27 353
16*16的数独题。网上已经有很多解释了，关于解法推荐博客：https://blog.csdn.net/bl0ss0m/article/details/17918705这里用来保存舞蹈链模板！#include &lt;stdio.h&gt; using namespace std; typedef long long ll; const int MAXN = 16*16*16+10; const ...
POJ2676 Sudoku 舞蹈链 DLX
weixin_30595035的博客
08-07 113
欢迎访问~原文出处——博客园-zhouzhendong 去博客园看该题解 题目（传送门） 题意概括 给出一个残缺的数独，求解。SPJ 题解 DLX + 矩阵构建 （两个传送门） 代码 #include <cstring> #include <cstdio> #include <algorithm> #include &...
Leetcode 352. 将数据流变为多个不相交区间 (利用set维护区间，同时使用两边使用哨兵进行二分查找）
wwxy1995的博客
10-09 89
typedef long long LL; const LL INF = 1e18; typedef pair<LL, LL> PLL; class SummaryRanges { private: set<PLL> intervals; public: SummaryRanges() { intervals.insert({-INF, -INF}); intervals.insert({INF, INF}); }...
poj3074-舞蹈链
测试
02-24 290
题目链接：http://poj.org/problem?id=3074 DescriptionIn the game of Sudoku, you are given a large 9 × 9 grid divided into smaller 3 × 3 subgrids. . 2 7 3 8 . . 1 .....
P1886&&POJ2823&&ybtoj【动态规划】6章1题【滑动窗口/【模板】单调队列】
zhanglili1597895的博客
01-24 140
滑动窗口/【模板】单调队列 题目 滑动窗口/【模板】单调队列 解析 单调队列模板题： 单调队列：用双端队列实现，取值时从队头取，插入值时将队尾的比它小/大的值全部扔出去即可维护单调递增/递减 神奇的是，洛咕这题包括图都是直接抄POJ的，但是洛咕的乱写都能过，POJ就必须加快读，最关键的是，POJ给了12秒，洛咕只有1秒…… 反观程序效率，同样是快读快输，蒟蒻我的程序POJ用了750MS，某句爷的程序POJ居然用了11922MS(还差一点就T了) 《人生巅峰》 code: #include<cstdi
贪心算法：田忌赛马－HDU1052 & POJ2287 & ZOJ2397
CSDN
12-30 5783
  Tian Ji -- The Horse Racing Time Limit: 5 Seconds      Memory Limit: 32768 KB   Here is a famous story in Chinese history. "That was about 2300 years ago. General Tian Ji was a high official in...
POJ2406 & POJ1961 详细解题思路和AC程序 KMP算法变种
IceTeaSet的专栏
07-09 844
poj2406 题目大意: 给出一个字符串s，求其最短的子串a，使得s=a^n(即n个a相连)，输出n。s的长度小于10^6。 本题可以用kmp算法来做。 先算出next数组。有些KMP算法的代码中next[i]代表的是s[i+1]匹配不上
POJ 1861 && ZOJ 1542--Network 【最小生成树 && kruscal && 水题】
_阿阿阿阿欢
08-10 924
题目大意：给n个点，m条边。求一棵生成树，并且满足任意两点之间的距离的最大值最小。输出这个最大值，然后输出树的边的数量，最后输出树的每条边。 题目思路：对于一个生成树来说，它的最大边满足在所有的生成树的最大边里最小，这其实就是一颗最小生成树。 #include #include #include #include #define maxn 1010 #define maxm 200
stl(图论鬼畜题)
06-27 377
论stl的妙用#include<bits/stdc++.h> #define inf 0x7fffffff #define ll long long #define ull unsigned long long using namespace std; const int mod =1e4+7; const ll INF=1e18; typedef pair<int ,int > P; const
POJ 3074,3076,2676 数独 Dancing Links舞蹈链
liwei的专栏
08-29 1313
思路：数独正好是精确覆盖问题的最好例子，这个也算是舞蹈链的模板了吧。 这个和网络流一样，难的是建模。如果能想好怎么建立矩阵的话，那一切都容易了。 这题怎么建矩阵，可以看论文：http://www.cnblogs.com/steady/archive/2011/03/15/1984791.html #pragma comment(linker, "/STACK:1024000000,10240
const int INF=0x3f3f3f3f;——ACM中的无穷大常量
while_black的博客
04-07 4590
0x3f3f3f3f 0x开头的 是十六进制常数， 等于 十进制 1061109567 等于 二进制： 00111111 00111111 00111111 00111111 const int inf=0x3f3f3f3f -- 声明 inf 是 const int型 , const 表示 inf 一旦有了值以后就不允许（通过赋值 来）改变它的值. 在算法竞赛中，我们常常需要用到设置一个...
HUST 1017 舞蹈链之不可重复覆盖
Dan__ge的博客
06-17 2538
点击打开链接 题意：每一行的某些列给定为1，现在问是否能找某些行，使得每一列的1出现一次 思路：与可重复覆盖的区别就是选定的列的1与其它行的又重复的，那么重复的那行也需要删除，而可重复的不需要删除，那么直接改一改之前的那个模版就可以求这道题了，而且需要输出所选行，任意一个就可以，那么在查找的时候直接赋给一个数组即可#include #include #include #include
HDU 3529 舞蹈链之可重复覆盖
Dan__ge的博客
06-17 2326
点击打开链接 题意：在点的地方放炸弹，炸弹可以炸到的地方无限远并且是十字的形式，但是如果炸到#或者*则这个方向停下来，问最少需要多少炸弹同时引爆才可以将所有#炸掉 思路：一道模版舞蹈链了，对于建图我们将所有的#看做列，然后所有的点也就是放炸弹的地方看做行，然后找一下每个点可以炸的#就可以了，建图后跑个舞蹈链可重复的模版就可以了#include #include #include #inc
FZU 1686 舞蹈链之可重复覆盖
Dan__ge的博客
06-17 2269
点击打开链接 题意：中文题了 思路：做了几个舞蹈链的题目了，第一感觉总会想到去用网络流去写，哎YY一会就是错的，乖乖用舞蹈链来写把，对于这道题目我们可以将所有的1看成是列，然后行则是一个一个的小矩阵，每一行也就是小矩阵连接自己可以消灭的1，但是首先要给所有的1编一下号，再然后就是模版舞蹈链了#include #include #include #include #include #i
HDU 2295 舞蹈链可重复覆盖+二分
Dan__ge的博客
06-21 2161
点击打开链接 题意：给了n个点和m个圆心，在用不多于k个圆的情况下，使得所有的点被覆盖，最小的圆半径多大 思路：问最小的什么什么那么用二分准没错，那么二分条件是什么呢，就是用不多于K的圆能否将所有点覆盖，覆盖这部分就是裸的舞蹈链的可重复覆盖，行为m个圆，列为n个点，然后求出的最少行小于等于K，则成立，否则不成立#include #include #include #include #i
POJ编程题解资源集合
"这些题目是来自POJ (Programming Online Judge) 的一系列算法练习题，涵盖了基础到中等难度的编程问题。题目涉及了各种算法和数据结构，包括但不限于数学计算、字符串处理、逻辑推理、递归、动态规划、位运算等。...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值