poj2032Square Carpets(IDA* + dancing links)
题目大意:给一个H行W列的01矩阵,求最少用多少个正方形框住所有的1.
题目分析:又是一个红果果的重复覆盖模型.DLX搞之!
枚举矩阵所有的子正方形,全1的话建图.判断全1的时候,用了一个递推,dp[i][j][w][h]表示左上角(i,j)的位置开始长h宽w的矩形中1的个数,这样后面可以迅速判断某个正方形是否全1.
不过此题直接搜一直TLE,然后改成迭代加深就比较愉快啦
详情请见代码:
#include <iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std; const int N = 11;
const int M = 50005;
int dp[N][N][N][N];
int n,m,num,ans;
int mp[N][N];
bool vis[105];
int s[M],h[M],col[M],u[M],d[M],l[M],r[M]; void init()
{
memset(h,0,sizeof(h));
memset(s,0,sizeof(s));
int i,c;
c = m * n;
for(i = 0;i <= c;i ++)
{
u[i] = d[i] = i;
l[i] = (i + c)%(c + 1);
r[i] = (i + 1)%(c + 1);
}
num = c + 1;
} void ins(int i,int j)
{
if(h[i])
{
r[num] = h[i];
l[num] = l[h[i]];
r[l[num]] = l[r[num]] = num;
}
else
h[i] = l[num] = r[num] = num;
s[j] ++;
u[num] = u[j];
d[num] = j;
d[u[num]] = num;
u[j] = num;
col[num] = j;
num ++;
}
void del(int c)
{
for(int i = u[c];i != c;i = u[i])
l[r[i]] = l[i],r[l[i]] = r[i];
}
void resume(int c)
{
for(int i = d[c];i != c;i = d[i])
l[r[i]] = r[l[i]] = i;
} int A()
{
int i,j,k,ret;
ret = 0;
memset(vis,false,sizeof(vis));
for(i = r[0];i;i = r[i])
{
if(vis[i] == false)
{
ret ++;
vis[i] = true;
for(j = d[i];j != i;j = d[j])
for(k = r[j];k != j;k = r[k])
vis[col[k]] = true;
}
}
return ret;
} bool dfs(int k,int lim)
{
if(k + A() > lim)
return false;
if(!r[0])
{
// ans = min(ans,k);
return true;
}
int i,j,c,mn = M;
for(i = r[0];i;i = r[i])
{
if(s[i] < mn)
{
mn = s[i];
c = i;
}
}
for(i = d[c];i != c;i = d[i])
{
del(i);
for(j = l[i];j != i;j = l[j])
del(j);
if(dfs(k + 1,lim)) return true;
for(j = r[i];j != i;j = r[j])
resume(j);
resume(i);
}
return false;
} bool canfuck(int x,int y,int z)
{
return dp[x][y][z][z] == z * z;
int i,j;
for(i = x;i <= x + z - 1;i ++)
for(j = y;j <= y + z - 1;j ++)
if(mp[i][j] == 0)
return false;
return true;
} void fuckba(int x,int y,int z,int id)
{
int i,j;
for(i = x;i <= x + z - 1;i ++)
for(j = y;j <= y + z - 1;j ++)
ins(id,(i - 1) * m + j);
}
void build()
{
int i,j,k;
memset(dp,0,sizeof(dp));
for(i = 1;i <= n;i ++)
for(j = 1;j <= m;j ++)
scanf("%d",&mp[i][j]),dp[i][j][1][1] = mp[i][j];
for(int ii = 1;ii <= n;ii ++)
for(int jj = 1;jj <= m;jj ++)
{
for(i = 1;i + ii <= n + 1;i ++)
{
for(j = 1;j + jj <= 1 + m;j ++)
{
dp[i][j][ii][jj] = dp[i][j][ii - 1][jj - 1] + dp[i + ii - 1][j][1][jj - 1] + dp[i][j + jj - 1][ii - 1][1] + dp[i + ii - 1][j + jj - 1][1][1];
}
}
}
init();
int rownum = 1;
for(k = 1;k <= min(n,m);k ++)//枚举边长
{
for(i = 1;i <= n - k + 1;i ++)
{
for(j = 1;j <= m - k + 1;j ++)
{
if(canfuck(i,j,k))
fuckba(i,j,k,rownum);
rownum ++;
}
}
}
for(i = 1;i <= n * m;i ++)
if(s[i] == 0)
l[r[i]] = l[i],r[l[i]] = r[i];
}
void fuck()
{
// ans = M;TLE....
// dfs(0);
// printf("%d\n",ans);
ans = 0;
while(!dfs(0,ans ++));
printf("%d\n",ans - 1);
}
int main()
{
while(scanf("%d%d",&m,&n),(m + n))
{
build();
fuck();
}
return 0;
}
poj2032Square Carpets(IDA* + dancing links)的更多相关文章
- Dancing Links and Exact Cover
1. Exact Cover Problem DLX是用来解决精确覆盖问题行之有效的算法. 在讲解DLX之前,我们先了解一下什么是精确覆盖问题(Exact Cover Problem)? 1.1 Po ...
- 跳跃的舞者,舞蹈链(Dancing Links)算法——求解精确覆盖问题
精确覆盖问题的定义:给定一个由0-1组成的矩阵,是否能找到一个行的集合,使得集合中每一列都恰好包含一个1 例如:如下的矩阵 就包含了这样一个集合(第1.4.5行) 如何利用给定的矩阵求出相应的行的集合 ...
- ZOJ 3209 Treasure Map (Dancing Links)
Treasure Map Time Limit: 2 Seconds Memory Limit: 32768 KB Your boss once had got many copies of ...
- HUST 1017 - Exact cover (Dancing Links 模板题)
1017 - Exact cover 时间限制:15秒 内存限制:128兆 自定评测 5584 次提交 2975 次通过 题目描述 There is an N*M matrix with only 0 ...
- Dancing Links初学记
记得原来备战OI的时候,WCX大神就研究过Dancing Links算法并写了一篇blog.后来我还写了个搜索策略的小文章( http://www.cnblogs.com/pdev/p/3952279 ...
- Dancing Links
Dancing Links用来解决如下精确匹配的问题: 选择若干行使得每一列恰好有一个1.Dancing Links通过对非零元素建立双向十字循环链表.上面的例子建立的链表如下所示: 计算的时候使用搜 ...
- 【转】Dancing Links题集
转自:http://blog.csdn.net/shahdza/article/details/7986037 POJ3740 Easy Finding [精确覆盖基础题]HUST1017 Exact ...
- 【转】Dancing Links精确覆盖问题
原文链接:http://sqybi.com/works/dlxcn/ (只转载过来一部分,全文请看原文,感觉讲得很好~)正文 精确覆盖问题 解决精确覆盖问题 舞蹈步骤 效率分析 ...
- POJ 3074 Sudoku (Dancing Links)
传送门:http://poj.org/problem?id=3074 DLX 数独的9*9的模板题. 具体建模详见下面这篇论文.其中9*9的数独怎么转化到精确覆盖问题,以及相关矩阵行列的定义都在下文中 ...
随机推荐
- poj2000---和1969一样的模板
#include <stdio.h> #include <stdlib.h> int main() { int d; while(scanf("%d",&a ...
- Gdal 1.11.0 添加 Postgresql 9.1 sqlite3 支持
OS环境Ubuntu12.04 32bit 因为公司一个功能要用到gdal 的ogr2ogr命令转换shp数据,需要能往postgis和sqlite 中插入数据. 用gdal1.11.0的源码默认安装 ...
- 通过crash了解linux页表
目的: 通过一个c语言实例,了解linux页表的组织结果和mmu的工作原理. 通过页表找到一个物理地址, 对比物理地址与虚拟地址的内容是否一致. 运行环境: $ uname -r3.15.6-200. ...
- Android版xx助手之天天酷跑外挂具体分析
Android版xx助手之天天酷跑外挂具体分析 图/文 莫灰灰 背景 近些年来,移动互联网的大肆崛起,潜移默化中影响着人们的生活和工作习惯.当腾讯的微信平台接入手机游戏之后,移动端的游戏也開 ...
- 统计分析SQL Server Profiler 跟踪的SQL
--跟踪文件读入到表中分析 SELECT * INTO ZGSJY FROM fn_trace_gettable('E:\wxxcdbprofiler.trc', default); --某时间内,最 ...
- 【Android进阶学习】shape和selector的结合使用
shape和selector是Android UI设计中经常用到的,比如我们要自定义一个圆角Button,点击Button有些效果的变化,就要用到shape和selector.可以这样说,shape和 ...
- 自学JavaScript的几个例子
学习了广泛使用的浏览器脚本JavaScript和HTML的DOM模型(也是用JavaScript实现),下面给出两个自己学习时的例子,具体JavaScript语法请参考W3C相关网页(http://w ...
- 我的学习笔记之API-Guides翻译------AppComponent_Activites
10.26第一天开始:贵在坚持,边看遍整理 一个应用程序通常由多个Activity组成,它们之间是松耦合的关系.特别的,有一个Activity作为app的主Activity,当app首次启动时呈现给用 ...
- 给定数组A,大小为n,现给定数X,判断A中是否存在两数之和等于X
题目:给定数组A,大小为n,现给定数X,判断A中是否存在两数之和等于X 思路一: 1,先采用归并排序对这个数组排序, 2,然后寻找相邻<k,i>的两数之和sum,找到恰好sum>x的 ...
- 提供一段Excel获取Title的标题,类似于A、AA、AAA,我们操作Excel的时候通常根据次标题来获取一定的操作范围。
/******************************************** FormatExcelColumTitle Purpose Get excel title like &qu ...