flood_it 方法记录
心路历程
根据题面的描述,我们面临的问题无非是,每次将色块更新成什么颜色。又因为是从左上角开始更新,所以我的有了第一个想法。
将左上角的色块命名为“原色块”。
对于每个色块,定义4中状态:
0-不属于原色块势力,和原色块势力不邻接,未没有进行任何操作;
1-不属于原色块势力,和原色块势力邻接,尚不知道会不会被同化(因为还不知道下一步染什么色);
2-属于原色块势力,和非原色块势力邻接,是上一批被同化的色块;
3-属于原色块势力,和非原色块势力不邻接**;
我们可以形象地依次把这四种状态称作“未加入”,“待审核”,“新成员”,“老成员”。
以下图为例(颜色是色块,数字是状态):

第一步比较特殊,因为需要求出最初的新队员。(之后就只需要处理新成员和待审核就可以了,因为新成员处于邻接位置,需要承担审核工作。而老成员和未加入不需要做任何工作)
首先原色块自身是新成员,然后跑一遍bfs,求出所有与原色块连通且与原色块同色的色块,将它们都标记为新成员。
之后的每一步,首先统计与新成员邻接且未加入的色块,将它们标记为待审核,并统计每种颜色出现的次数。
出现次数最多的颜色就是我们即将更新的颜色。得到这个之后新成员的工作就做完了,新成员就可以标记为老成员了。
再次遍历所有待审核的色块,若颜色和即将更新的颜色相同,则标记为新成员;否则重新标记为未加入。
将所有新成员、老成员的颜色更新。
每次都需要判断是否更新完整个图。
这个算法逻辑上是行得通的,但是太慢了,因为要反复地遍历各种状态,导致大量无用的重复运算。
正解
标算:IDA*
先不谈IDA*。我们回到这道题,从模块化编程的角度出发,想想这个程序要实现什么功能。
首先我们有更新色块的操作。对于更新色块的函数,我们需要3个形参:当前色块的横坐标、纵坐标、我们即将更新的颜色。我们仍然需要一个打标记的数组:若当前位置的颜色和即将更新的颜色相同则标记为1,反之则标记为2.所以这个函数的主要功能是给需要更新颜色的色块打标记,还没有进行上色操作。
void update_flags(int i,int j,int color)//color是即将更新的颜色
{
if(board[i][j]!=color)
{
flags[i][j]=2;//和即将更新的颜色不同,打上标记2
return ;
}
flags[i][j]=1;//和即将更新的颜色相同,打上标记1
for(int k=0;k<4;k++)
{
//dir[][]用于枚举4个方向
int x=i+dir[k][0];
int y=j+dir[k][1];
if(x<0||y<0||x>=siz||y>=siz||flags[x][y]) continue;
update_flags(x,y,color);//递归更新尚未打标记的色块
}
}
然后考虑具体的上色操作。最外层循环遍历6种颜色,然后遍历全图。若当前遍历到的位置是不同颜色且正好等于要搜索的颜色,则记录存在这种颜色,并进行\(update_flags\).若全图遍历完后发现不存在这种颜色,则继续遍历下一个颜色。如果更新了棋盘,就递归搜索,成功了就返回。注意:\(flags\)数组在这个过程中可能会被改变,为了后续的使用,我们需要在操作前拷贝\(flags\)到\(tmp\)数组,操作结束后再拷贝回去。
接下来再考虑使用IDA.IDA可以理解为以迭代加深DFS的搜索框架为基础,把原来简单的深度限制加强为:若当前深度+未来估计步数>深度限制,则立即从当前分支回溯
具体到本题,我们的估价函数就应该是:如果颜色数量大于剩余步数,则行不通。所以我们还需要一个计算剩余颜色数量的函数。
int color_cnt()
{
int ret=0;
int s=0;
for(int i=0;i<siz;i++)
for(int j=0;j<siz;j++)
if(flags[i][j]!=1)
s|=(1<<board[i][j]);
while(s)//检查还有几种颜色
{
ret+=s&1;
s>>=1;
}
return ret;
}
并且,我们要在主函数中枚举每一种最大限定深度。
最终代码呈现如下:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=15;
int siz;
int board[N][N];
int dir[4][2]={{0,1},{1,0},{0,-1},{-1,0}};
int flags[N][N];
int max_depth;
void update_flags(int i,int j,int color)//color是即将更新的颜色
{
if(board[i][j]!=color)
{
flags[i][j]=2;//和即将更新的颜色不同,打上标记2
return ;
}
flags[i][j]=1;//和即将更新的颜色相同,打上标记1
for(int k=0;k<4;k++)
{
//dir[][]用于枚举4个方向
int x=i+dir[k][0];
int y=j+dir[k][1];
if(x<0||y<0||x>=siz||y>=siz||flags[x][y]) continue;
update_flags(x,y,color);//递归更新尚未打标记的色块
}
}
int color_cnt()
{
int ret=0;
int s=0;
for(int i=0;i<siz;i++)
for(int j=0;j<siz;j++)
if(flags[i][j]!=1)
s|=(1<<board[i][j]);
while(s)//检查还有几种颜色
{
ret+=s&1;
s>>=1;
}
return ret;
}
int dfs(int depth)
{
//如果颜色数量大于剩余步数,则行不通(估价函数)
if(color_cnt()>max_depth-depth)
return 0;
//如果当前深度刚好等于最大深度,说明已经找完了
if(depth==max_depth)
return 1;
int tmp[N][N];
int color_exist;
memcpy(tmp,flags,sizeof(flags));
for(int color=0;color<6;color++)
{
color_exist=0;
for(int i=0;i<siz;i++)
for(int j=0;j<siz;j++)//当前位置是不同颜色且正好等于要搜索的颜色
if(flags[i][j]==2&&board[i][j]==color)
{
color_exist=1;
update_flags(i,j,color);
}
//颜色不存在,继续搜索
if(!color_exist) continue;
//如果更新了棋盘,就递归搜索,成功则返回
if(dfs(depth+1)) return 1;
memcpy(flags,tmp,sizeof(flags));
}
return 0;
}
int main()
{
while(scanf("%d",&siz)&&siz)
{
for(int i=0;i<siz;i++)
for(int j=0;j<siz;j++)
scanf("%d",&board[i][j]);
//限定最大深度
for(max_depth=0;;max_depth++)
{
memset(flags,0,sizeof(flags));
update_flags(0,0,board[0][0]);
if(dfs(0)) break;
}
printf("%d\n",max_depth);
}
return 0;
}
flood_it 方法记录的更多相关文章
- EF里查看/修改实体的当前值、原始值和数据库值以及重写SaveChanges方法记录实体状态
本文目录 查看实体当前.原始和数据库值:DbEntityEntry 查看实体的某个属性值:GetValue<TValue>方法 拷贝DbPropertyValues到实体:ToObject ...
- 64位 SQL Server2008链接访问Oracle 过程汇总解决方法记录
64位 SQL Server2008链接访问Oracle 过程汇总解决方法记录 经过几天不停的网上找资料,实验,终于联通了. 环境:系统:win 2008 ,SqlServer2008 R2, 连接O ...
- js实用方法记录-js动态加载css、js脚本文件
js实用方法记录-动态加载css/js 附送一个加载iframe,h5打开app代码 1. 动态加载js文件到head标签并执行回调 方法调用:dynamicLoadJs('http://www.yi ...
- js实用方法记录-简单cookie操作
js实用方法记录-简单cookie操作 设置cookie:setCookie(名称,值,保存时间,保存域); 获取cookie:setCookie(名称); 移除cookie:setCookie(名称 ...
- js实用方法记录-指不定哪天就会用到的js方法
js实用方法记录-指不定哪天就会用到的js方法 常用或者不常用都有 判断是否在微信浏览器中 测试代码:isWeiXin()==false /** * 是否在微信中 */ function isWeix ...
- Java给各个方法记录执行时间
Java给各个方法记录执行时间 long startTime = System.currentTimeMillis();...//要测试时间的方法LoggerFactory.getLogger(Bas ...
- make menuconfig error 解决方法记录
新建的一个虚拟机,发现make menuconfig 后会出错,查了一下是缺少一些库. 这个错误已经错了两次了,我不希望第三次出现了还想不起来,所以特此记录. # 错误信息: make[2]: *** ...
- 简单一键CENTOS6 安装PPTP VPN方法记录
申明:我们使用PPTP VPN仅仅只能用在查阅资料等正规渠道,不要用在不良用途上.方法收集于网上,这里我用在搬瓦工VPS(VPS方案直达),采用的是CENTOS6 64位系统.我们需要预先将VPS服务 ...
- ASP.NET页面优化性能提升方法记录
今天与大家分享:一种优化页面执行速度的方法.采用这个方法,可以使用页面的执行速度获得[8倍]的提升效果. 为了让您对优化的效果有个直观的了解,我准备了下面的测试结果截图: 测试环境:1. Window ...
随机推荐
- appache ab测试高并发
ab使用范例: 要执行 1000 次的 connection, 20 次的 concurrent (并行, 同时): 语法: ab -n 1000 -c 20 http://localhost/ind ...
- python闭包函数及装饰器简介
目录: 闭包函数简介 闭包函数的实际应用 装饰器简介 装饰器初期-完整版 装饰器语法糖 闭包函数简介 1.定义在函数内部的函数(函数的嵌套) 2.内部函数运用外部函数局部名称空间中的变量名 注:函数名 ...
- Vue3.0工程创建 && setup、ref、reactive函数 && Vue3.0响应式实现原理
1 # 一.创建Vue3.0工程 2 # 1.使用vue-cli创建 3 # 官方文档: https://cli.vuejs.org/zh/guide/creating-a-project.html# ...
- Excel 数学函数(二):ISODD 和 ISEVEN
ISODD 函数用于判断一个数是否为奇数:ISEVEN 函数用于判断一个数是否为偶数. 结合文本函数 MID 和逻辑函数 IF,来判断一个身份证的性别. 身份证倒数第二位是奇数就代表男性,否则代表女性 ...
- 如何使能512个virtio_blk设备
一例virtio_blk设备中断占用分析 背景:这个是在客户的centos8.4的环境上复现的,dpu是目前很多 云服务器上的网卡标配了,在云豹的dpu产品测试中,dpu实现的virtio_blk 设 ...
- HDU4372 Count the Buildings (+题解:斯特林数)
题面 (笔者翻译) There are N buildings standing in a straight line in the City, numbered from 1 to N. The h ...
- 持久化-Powershell配置文件持久性
持久化-Powershell配置文件持久性 概述 可以使用Powershell配置文件进行持久性和/或特权升级. 执行 获取可以滥用的powershell配置文件,这取决于你拥有的权限. $PROFI ...
- 【java】非常多!学习路径24-总结目前所有知识(上)
感谢sikiedu.com的siki老师.几年前就开始看siki的课程,最近突然想写这个笔记系列,顺便回顾一下这些基础的知识,同时也希望能帮助到一些人,有问题一起交流哈. 全文共十章,大约1.5万字, ...
- 【unity游戏入门】2 使用代码编写Hello Unity游戏
作者 罗芭Remoo 2021年9月24日 第一章.许可证的安装 下载好Unity之后,我们还需要一个前置操作才可以进入Unity引擎----许可证. 当然不用担心,Unity是一个开放的引擎,一切以 ...
- 第十四章 二进制部署k8s集群的平滑升级
1.软件包下载 去github上下载较新的Kubernetes软件包https://github.com/ 2.升级说明 升级包括master节点升级和node节点的升级,本章升级至v1.15.12: ...