这是昨天周赛的题,我竟然不会怎么判断多解,后来一google,卧槽,我想复杂了。。。。。。直接看能搜出来几次就行了。

这题就是个变形,先floodfill一下,然后就是模板了

然后发现比大华的快了好几倍,然后加了输入挂后,瞬间成Best solution中的rank1了

/*
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];
int ditu[N + 5][N + 5], pos[N * 5][N * 5];
int num, ct;
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;num = 0;ct = 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; //代表数独中同一列所属的列ck = N * N * 3 + (pos[i][j] - 1) * N + k;//代表数独中同一宫所属的列
}
void search(int x, int y)
{pos[x][y] = ct;int ta = ditu[x][y];if(ta & 16) ta ^= 16;if(ta & 32) ta ^= 32;if(ta & 64) ta ^= 64;if(ta & 128) ta ^= 128;mp[x][y] = ta + '0';ta = ditu[x][y];if((ta & 16) == 0 && x >= 1 && pos[x - 1][y] == 0) search(x - 1, y);if((ta & 32) == 0 && y < N - 1 && pos[x][y + 1] == 0) search(x, y + 1);if((ta & 64) == 0 && x < N - 1 && pos[x + 1][y] == 0) search(x + 1, y);if((ta & 128) == 0 && y >= 1 && pos[x][y - 1] == 0) search(x, y - 1);
}
void readdata()
{int r, cx, cy, ck, cg;memset(pos, 0, sizeof(pos));for(int i = 0; i < N; i++)for(int j = 0; j < N; j++)scanf("%d", &ditu[i][j]);for(int i = 0; i < N; i++)for(int j = 0; j < N; j++)if(pos[i][j] == 0){++ct;search(i, j);}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){num++;  //直接用num记录有多少个解for(int i = 0; i < k; i++)ans[X[O[i]] / N] = X[O[i]] % 9 + '1';if(num >= 2) return true;else return false;}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, cas = 0;scanf("%d", &T);while(T--){init();readdata();dfs(0);printf("Case %d:\n", ++cas);if(num == 0) printf("No solution\n");else if(num == 1){for(int i = 0; i < N * N; i++){printf("%c", ans[i]);if((i + 1) % N == 0) printf("\n");}}else printf("Multiple Solutions\n");}return 0;
}

HDU 4069 Squiggly Sudoku DLX相关推荐

  1. HDU 4069 Squiggly Sudoku DLX 精确覆盖

    题意: 数独问题,给你9个连通块,每个连通块有9个位置. 现在已经有一些数字在上面,让你在空的位置上放数字. 问你是否存在方案,使得每个连通块包含1~9,并且每行每列都有1~9的数字. 输出结果参照样 ...

  2. HDU 4069 Squiggly Sudoku(DLX)(The 36th ACM/ICPC Asia Regional Fuzhou Site —— Online Contest)...

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4069 Problem Description Today we play a squiggly sud ...

  3. HDU 4069 Squiggly Sudoku 【DLX+BFS】

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4069 ★这题一开始题都看不懂,后来发现还是挺有意思的 题意: 给你一个9x9的矩阵, 矩阵里面有一些墙 ...

  4. HDU 4069 Squiggly Sudoku Dancing-Links(DLX)+Floodfill

    题目大意:..还是数独,不同的是原先的九宫格约束条件变为了给定的任意形状... 我们跑一遍floodfill 得出每一个格子属于哪一个形状 然后就是裸的数独了 这题T<=2500 绝对不能开动态 ...

  5. [DLX+bfs] hdu 4069 Squiggly Sudoku

    题意: 给你9*9的矩阵.对于每个数字,能减16代表上面有墙,能减32代表下面有墙... 最后剩下的数字是0代表这个位置数要求,不是0代表这个数已知了. 然后通过墙会被数字分成9块. 然后做数独,这里 ...

  6. HDU 4069 Squiggly Sudoku【Dancing Links精确覆盖】

    跟普通的数独有一点点不同,先预处理一下再用Dancing Links进行精确覆盖即可. #include <iostream> #include <cstdio> #inclu ...

  7. HDOJ 4069 Squiggly Sudoku 精确覆盖+搜索

    //HDOJ 4069 Squiggly Sudoku 精确覆盖+搜索 /* 题意:数独变形 9块弯弯曲曲的区域 每块有9个小格往这81个格子里面填写1-9的数字,使得每行,每列,每个区域都含有1.2 ...

  8. [kuangbin]专题三 Dancing Links Squiggly Sudoku HDU - 4069【DFS】【精确覆盖】

    [题目描述] Today we play a squiggly sudoku, The objective is to fill a 9*9 grid with digits so that each ...

  9. HDU-4069(Squiggly Sudoku)(Dancing Links + dfs)

    题目:HDU - 4069 Today we play a squiggly sudoku, The objective is to fill a 9*9 grid with digits so th ...

最新文章

  1. 利用动态加载模板,配合ajax实现无刷新操作
  2. react-router browserHistory刷新页面404问题解决
  3. Android开发--AsyncTask异步任务(一)
  4. html弹窗确认取消公告代码,js 弹出确认与取消对话框的四种方法
  5. Android下常见终端模拟器和SSH客户端感受及几个Tips
  6. 不等号属于不等式吗_初中数学——中考知识3【方程与不等式(基础)】
  7. 跟着动画学习 TCP 三次握手和四次挥手
  8. 欢迎来怼--第二十九次Scrum会议
  9. android软件游戏显示fps测试工具,fpsviewer—实时显示fps,监控Android卡顿的可视化工具...
  10. 热释电红外传感器电路
  11. 关于GPS信息在谷歌地图上实时更新位置信息的问题
  12. 数据结构设计题大题总结(非代码)
  13. PS快速修改颜色的小方法
  14. 解决:idea中文变繁体怎么办
  15. 程序员最新面试谈薪指南
  16. LaTeX期刊论文写作
  17. 台湾大学郭彦甫matlab百度云,台湾国立大学郭彦甫Matlab教程笔记(23) linear systems...
  18. JS 生成二维码实现(qrcode.js)
  19. 来 CSDN 实习技术编辑,直面改变世界的技术与人物
  20. 程设大作业之魔兽世界

热门文章

  1. 阿里云服务器购买合同怎么申请
  2. 【重磅】MIT发布2018年“全球十大突破性技术”
  3. 【blender建模功能】03 倒角工具
  4. 嵌入式linux智能小车ppt,基于FPGA嵌入式系统的智能小车全面解析
  5. 使用flownet2进行flow warping error的计算
  6. 底层嵌入式之NOR FLASH编程
  7. 康托展开详解 -csdn博客
  8. Switch中的参数
  9. Resid------set
  10. Java算法之冒泡算法