心路历程

根据题面的描述,我们面临的问题无非是,每次将色块更新成什么颜色。又因为是从左上角开始更新,所以我的有了第一个想法。

将左上角的色块命名为“原色块”。

对于每个色块,定义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 方法记录的更多相关文章

  1. EF里查看/修改实体的当前值、原始值和数据库值以及重写SaveChanges方法记录实体状态

    本文目录 查看实体当前.原始和数据库值:DbEntityEntry 查看实体的某个属性值:GetValue<TValue>方法 拷贝DbPropertyValues到实体:ToObject ...

  2. 64位 SQL Server2008链接访问Oracle 过程汇总解决方法记录

    64位 SQL Server2008链接访问Oracle 过程汇总解决方法记录 经过几天不停的网上找资料,实验,终于联通了. 环境:系统:win 2008 ,SqlServer2008 R2, 连接O ...

  3. js实用方法记录-js动态加载css、js脚本文件

    js实用方法记录-动态加载css/js 附送一个加载iframe,h5打开app代码 1. 动态加载js文件到head标签并执行回调 方法调用:dynamicLoadJs('http://www.yi ...

  4. js实用方法记录-简单cookie操作

    js实用方法记录-简单cookie操作 设置cookie:setCookie(名称,值,保存时间,保存域); 获取cookie:setCookie(名称); 移除cookie:setCookie(名称 ...

  5. js实用方法记录-指不定哪天就会用到的js方法

    js实用方法记录-指不定哪天就会用到的js方法 常用或者不常用都有 判断是否在微信浏览器中 测试代码:isWeiXin()==false /** * 是否在微信中 */ function isWeix ...

  6. Java给各个方法记录执行时间

    Java给各个方法记录执行时间 long startTime = System.currentTimeMillis();...//要测试时间的方法LoggerFactory.getLogger(Bas ...

  7. make menuconfig error 解决方法记录

    新建的一个虚拟机,发现make menuconfig 后会出错,查了一下是缺少一些库. 这个错误已经错了两次了,我不希望第三次出现了还想不起来,所以特此记录. # 错误信息: make[2]: *** ...

  8. 简单一键CENTOS6 安装PPTP VPN方法记录

    申明:我们使用PPTP VPN仅仅只能用在查阅资料等正规渠道,不要用在不良用途上.方法收集于网上,这里我用在搬瓦工VPS(VPS方案直达),采用的是CENTOS6 64位系统.我们需要预先将VPS服务 ...

  9. ASP.NET页面优化性能提升方法记录

    今天与大家分享:一种优化页面执行速度的方法.采用这个方法,可以使用页面的执行速度获得[8倍]的提升效果. 为了让您对优化的效果有个直观的了解,我准备了下面的测试结果截图: 测试环境:1. Window ...

随机推荐

  1. 2511-Druid监控功能的深入使用与配置-如何记录监控数据(基于logback)

    Druid的监控很强大,但可惜的是监控数据是存在内存中的,需求就是定时把监控数据记录下来,以日志文件的形式或者数据库入库. 记录两种方式: 数据库入库 logback形式记录 原理(重点) 如果仅仅想 ...

  2. php apache 和mysql查看版本常用方法收集

    php: 1.命令行查询,下图是因为添加php进系统环境变量了 2.预定义常量PHP_VERSION查询 3.phpversion()函数查询 4.phpinfo()查询 apache: mysql: ...

  3. 使用python3.7+Vue.js2.0+Django2.0.4异步前端通过api上传文件到七牛云云端存储

    原文转载自「刘悦的技术博客」https://v3u.cn/a_id_130 之前一篇文章是通过普通js+tornado来上传七牛云:使用Tornado配合七牛云存储api来异步切分上传文件,本次使用v ...

  4. 蔚来杯2022牛客暑期多校训练营6 ABGJM

    比赛链接 A 题解 知识点:数学,构造. 题目要求构造一个长为 \(m\) 的序列 \(c\) ,\(m\) 自选,使得 \(c\) 的无限循环序列 \(b\) 中任意连续 \(a_i\) 个数中都存 ...

  5. Apache Dolphinscheduler 1.3.x 系列配置文件指南

    前言 本文档为dolphinscheduler配置文件指南,针对版本为 dolphinscheduler-1.3.x 版本. 考虑公众号对markdown文件格式支持不那么友好的问题,建议大家在PC端 ...

  6. Excel 统计函数(四):AVERAGEIF 和 AVERAGEIFS

    AVERAGEIF [语法]AVERAGEIF(range, criteria, [average_range]) [参数] range:要计算平均值的一个或多个单元格: criteria:筛选条件: ...

  7. 【java】学习路径37-练习:任意文件的复制

    使用字节完成复制 import java.io.FileInputStream; import java.io.FileOutputStream; import java.io.IOException ...

  8. Windows 系统 PostgreSQL 手工安装配置方法

    自从2020年底开始接触 PostgreSQL 以来就喜欢上了这个数据库,个人感觉比 MySQL 好用,多表联合查询性能好很多,同时也不存在 SQLServer 的版权授权费用问题.搭配 .NET 开 ...

  9. MFRC522学习笔记

    MFRC522主要特性 容量为8K位(bits)=1K字节(bytes)EEPROM 分为16个扇区,每个扇区为4块,每块16个字节,以块为存取单位 每个扇区有独立的一组密码及访问控制 每张卡有唯一序 ...

  10. Altium Designer 18学习

    目录 目录 快捷键 通孔 敷铜 修改铜皮与导线之间的间隔 去除指定敷铜区域 DRC设计规则检查问题: 快捷键 EJC 快速跳转到器件 M 移动 CTRL+M 测量距离 通孔 敷铜 放置多边形平面 -- ...