P2156 [SDOI2009]细胞探索

最新推荐文章于 2024-05-31 09:44:20 发布
最新推荐文章于 2024-05-31 09:44:20 发布
阅读量140 收藏
点赞数
文章标签： 数据结构与算法
原文链接：http://www.cnblogs.com/olinr/p/10445464.html
版权

$ \color{#0066ff}{ 题目描述 }$

生物课上，老师开始为同学们介绍细胞。为了加深同学们的印象，老师在一张N×M的矩阵中定义了一种细胞，矩阵中仅有井号“#”和点“.”：

细胞由细胞核、细胞质及细胞膜构成。细胞核是一个4连通（上下左右相连）的全为“#”的连通块，它必须实心，即不能存在一个4连通的“.”连通块被其完全包围（所谓完全包围指的是，这个“.”连通块不能位于矩阵边界相邻，且它的4相邻格子均属于包含它的“#”连通块）。细胞膜是一个8连通（上下左右，以及4个对角方向）的全为“#”的非实心连通块。细胞膜仅包围一个4连通的区域，且这个区域内有且仅有一个细胞核，这个区域剩下的位置全为“.”。

所有连通块必须极大化，即一个8连通块周围不能找到一个“#”与这个连通块的任意一个“#”8连通；同样，对于一个4连通块周围不能找到一个“#”与这个连通块的任意一个“#”4连通。

现在，老师画了一幅图画，并让小E回答图画中一共有几个细胞，并把图画中不属于任何一个细胞的“#”改成“.”。

\(\color{#0066ff}{输入格式}\)

输入文件explore.in的第一行包含两个用空格分隔的正整数N和M，表示矩阵的高和长。

接下来一个N行M列的矩阵，矩阵中仅含井号“#”和点“.”，保证没有多余字符。

\(\color{#0066ff}{输出格式}\)

输出文件explore.out第一行包含一个整数，表示输入的矩阵中的细胞数。

接下来一个N行M列的矩阵，矩阵中仅含井号“#”和点“.”，表示更改后的图画。

\(\color{#0066ff}{输入样例}\)

12 13
.###..#####..
#...#.#....#.
#.#.#.#..#.#.
#...#..#...#.
.###.#..###..
....#..##...#
..........###
##########..#
#...........#
#.###...###.#
#...........#
#############


9 14
#########.....
#.......#....#
#.#####.#...#.
#.#...#.#..#..
#.#.#.#.#.#..#
#.#...#.#..#..
#.#####.#...#.
#.......#....#
#########.....


7 15
#######.#######
#.....#.#.....#
#.###.#.#.###.#
#.#.#.#.#.#...#
#.###.#.#.###.#
#.....#.#.....#
#######.#######

\(\color{#0066ff}{输出样例}\)

1
......#####..
......#....#.
......#..#.#.
.......#...#.
........###..
.......##....
.............
.............
.............
.............
.............
.............

    
1
..............
..............
..#####.......
..#...#.......
..#.#.#.......
..#...#.......
..#####.......
..............
..............

    
1
........#######
........#.....#
........#.###.#
........#.#...#
........#.###.#
........#.....#
........#######

\(\color{#0066ff}{数据范围与提示}\)

对于20%的数据，满足1 ≤ N, M ≤ 20。

另有20%的数据，满足所有“#”都属于某一个正确的细胞。

对于100%的数据，满足1 ≤ N, M ≤ 1,000。

\(\color{#0066ff}{题解}\)

看起来很难下手，如果搜井号的话，还要考虑内外，肯定很麻烦的
那么我们不如从点入手
不难发现，一个合法的细胞，细胞质连接且仅连接两个联通块
所以可以先预处理出所有井号联通块，然后。。。过程见代码。。
#include<bits/stdc++.h>
#define LL long long
LL in() {
    char ch; LL x = 0, f = 1;
    while(!isdigit(ch = getchar()))(ch == '-') && (f = -f);
    for(x = ch ^ 48; isdigit(ch = getchar()); x = (x << 1) + (x << 3) + (ch ^ 48));
    return x * f;
}
const int maxn = 1050;
const int inf = 0x7fffffff;
int bel[maxn][maxn], n, m, mp[maxn][maxn], num, siz[maxn * maxn], cnt[maxn * maxn];
int up[maxn * maxn], bj[maxn * maxn], ls, fa[maxn * maxn], must[maxn * maxn];
std::pair<int, int> spe[maxn * maxn];
char sss[maxn];
int dx[] = {-1, -1, -1, 0, 0, 1, 1, 1};
int dy[] = {-1, 0, 1, -1, 1, -1, 0, 1};
int rx[] = {1, -1, 0, 0};
int ry[] = {0, 0, 1, -1};
bool vis[maxn][maxn], flag, have;
std::set<int> v;
void dfs(int x, int y, int id) {
    siz[id]++;
    up[id] = std::min(up[id], x);
    bel[x][y] = id;
    vis[x][y] = true;
    for(int i = 0; i < 8; i++) {
        int xx = x + dx[i];
        int yy = y + dy[i];
        if(xx >= 1 && xx <= n && yy >= 1 && yy <= m && !vis[xx][yy] && mp[xx][yy]) dfs(xx, yy, id);
    }
}
void expand(int x, int y) {
    if(x == 1 || x == n || y == 1 || y == m) flag = true;
    vis[x][y] = true;
    for(int i = 0; i < 4; i++) {
        int xx = x + rx[i];
        int yy = y + ry[i];
        if(xx >= 1 && yy >= 1 && xx <= n && yy <= m) {
            if(mp[xx][yy]) v.insert(bel[xx][yy]);
            if(!mp[xx][yy] && !vis[xx][yy]) expand(xx, yy);
        }
    }
}
void init(int x, int y) {
#ifdef olinr
    printf("init %d %d\n", x, y);
#endif
    ls++;
    vis[x][y] = true;
    for(int i = 0; i < 4; i++) {
        int xx = x + rx[i];
        int yy = y + ry[i];
        if(xx >= 1 && xx <= n && yy >= 1 && yy <= m) {
#ifdef olinr
            printf("in");
#endif
            if(!mp[xx][yy] && !vis[xx][yy]) have = true;
#ifdef olinr
            if(mp[xx][yy]) printf("!");
            if(!vis[xx][yy]) printf("$");
#endif
            if(mp[xx][yy] && !vis[xx][yy]) init(xx, yy);
        }
    }
}


void work(int x, int y) {
    v.clear();
    flag = false;
    have = false;
    ls = 0;
    expand(x, y);
    if(flag) {
#ifdef olinr
        //搜索到边界一定不是细胞质
    printf("break at 1\n");
#endif
        return;
    }
    if(v.size() != 2) {
        //一个膜里面一大堆不知道什么玩意
        int pos = 0, mx = inf;
        for(std::set<int>::iterator it = v.begin(); it != v.end(); it++)
            if(mx > up[*it]) mx = up[*it], pos = *it;
        //这个膜肯定没用了
        bj[pos] = true, cnt[pos] += v.size() - 1, must[pos] = true;
        for(std::set<int>::iterator it = v.begin(); it != v.end(); it++)
            if(*it != pos) fa[*it] = pos;
#ifdef olinr
        printf("vis: ");
        for(std::set<int>::iterator it = v.begin(); it != v.end(); it++) printf("%d ", *it);
        puts("");
        printf("break at 2\n");
#endif
        return;
    }
    std::set<int>::iterator it = v.begin();
    int aa = *it;
    it++;
    int bb = *it; 
    if(up[aa] > up[bb]) std::swap(aa, bb);
    if(bj[aa]) {
#ifdef olinr
        printf("break at 3\n");
#endif
        //已经用过了，就是两个粘在一起的细胞，然而在本题是不合法的
        must[aa] = true;
        return;
    }
    init(spe[bb].first, spe[bb].second);
    if(have) {
        //细胞核里有TM细胞质！！！我去
        must[aa] = true;
#ifdef olilnr
        printf("break at 4\n");
#endif
        return;
    }
    if(ls != siz[bb]) {
        //因为搜井号联通块是8联通的，这里要4联通搜一下，看是不是完整
#ifdef olinr
        printf("%d %d\n", ls, siz[bb]), printf("break at 5\n");
#endif
        must[aa] = true;
        return;
    }
    fa[aa] = 0, fa[bb] = aa;
    cnt[aa]++;
}
//must：是否一定不合法，bj是否曾经被作为细胞膜，fa，细胞核属于哪个细胞膜， cnt，作为细胞膜的次数
int main() {
    n = in(), m = in();
    for(int i = 1; i <= n; i++) {
        scanf("%s", sss);
        for(int j = 1; j <= m; j++)
            mp[i][j] = sss[j - 1] == '#';
    }
    //处理井号联通块
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
            if(mp[i][j] && !vis[i][j]) {
                num++;
                up[num] = inf;
                spe[num] = std::make_pair(i, j);
                dfs(i, j, num);
#ifdef olinr
                printf("id: %d, siz = %d\n", num, siz[num]);
#endif
            }
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++) 
            vis[i][j] = 0;
    //对于每个点的联通块，处理
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
            if(!mp[i][j] && !vis[i][j]) {
#ifdef olinr
                printf("now dfs %d %d\n", i, j);
#endif
                work(i, j);
            }
    int ans = 0;
    for(int i = 1; i <= num; i++) {
        bj[i] = false;
        if(!must[i] && !fa[i] && cnt[i] == 1) ans++, bj[i] = true;
    }
    printf("%d\n", ans);
    for(int i = 1; i <= n; i++) {
        for(int j = 1; j <= m; j++) {
            if(!mp[i][j]) putchar('.');
            else {
                if(bj[bel[i][j]] || bj[fa[bel[i][j]]]) putchar('#');
                else putchar('.');
            }
        }
        puts("");
    }
    return 0;
}

转载于:https://www.cnblogs.com/olinr/p/10445464.html

weixin_30376453
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
LUOGU P2156 [SDOI2009]细胞探索
REWINE
03-19 261
题目描述 生物课上，老师开始为同学们介绍细胞。为了加深同学们的印象，老师在一张N×M的矩阵中定义了一种细胞，矩阵中仅有井号“#”和点“.”： 细胞由细胞核、细胞质及细胞膜构成。细胞核是一个4连通（上下左右相连）的全为“#”的连通块，它必须实心，即不能存在一个4连通的“.”连通块被其完全包围（所谓完全包围指的是，这个“.”连通块不能位于矩阵边界相邻，且它的4相邻格子均属于包含它的“#”连通块）...
sdoi2017R1D1
04-16
sdoi2017R1D1
五分钟成为记忆王
热门推荐
burbuly的专栏
06-23 2万+
一、记忆的面纱1、记忆的含义（1）就在我嘴边上　　有多少次你这样说过，就在我嘴边上，又有过多少次在你需要什么时候，任凭你如何拼命地想，就是想不起来。　　当然，这问题不是你一个人才有，几乎所有的人都受到过记忆力差的困扰。这也是人类的一个最常见的不幸。（2）改变你的记忆力　　自身内部，就蕴藏着一种由于记忆力差而产生的烦恼的能力。如果你真想利用这一能力的话，这能力就能使你的记忆力在几天内提高
【SDOI2009】解题汇总
weixin_30813225的博客
03-20 100
又开了波专题，感觉就和炉石开冒险一样...（说的好像我有金币开冒险似的） /—————————————————————————————————————————————/ BZOJ-1226 【SDOI2009】学校食堂Dining 状态压缩DP f【i】【j】【k】表示前i-1人都吃过饭，j表示i与i之后7人的吃饭情况，k表示上一个吃饭的人与i的相对位置 转移如...
P2156 [SDOI2009] 细胞探索
m0_71227570的博客
03-11 407
细胞核是一个4连通（上下左右相连）的全为“#”的连通块，它必须实心，即不能存在一个4连通的“.”连通块被其完全包围（所谓完全包围指的是，这个“.”连通块不能位于矩阵边界相邻，且它的4相邻格子均属于包含它的“#”连通块）。所有连通块必须极大化，即一个8连通块周围不能找到一个“#”与这个连通块的任意一个“#”8连通；同样，对于一个4连通块周围不能找到一个“#”与这个连通块的任意一个“#”4连通。接下来一个N行M列的矩阵，矩阵中仅含井号“#”和点“.”，表示更改后的图画。### 样例输入 #1。
P1972 [SDOI2009]HH的项链
Ajwallet
07-09 434
莫队
P2152 [SDOI2009]SuperGCD
qq_35975367的博客
08-11 154
P2152 [SDOI2009]SuperGCD 题意： 求a和b的最大公约数 a,b<=101000010^{10000}1010000 题解： 高精度，java高精度直接有模板，c++高精度结合gcd的辗转相减法 对于：a,b的gcd(a,b)有： 若a为奇数，b为偶数，gcd(a,b)=gcd(a,b/2) 若a为偶数，b为奇数，gcd(a,b)=gcd(a/2,b) 若a为偶数，b为偶数，gcd(a,b)=2gcd(a,b) 若a为奇数，b为奇数，gcd(a,b)=gcd(a-b,b)（a&
洛谷P2152 [SDOI2009] SuperGCD
puvakyOK的博客
08-26 232
P2152 [SDOI2009] SuperGCD
P2152-[SDOI2009]SuperGCD
weixin_30836759的博客
09-03 106
1 // luogu-judger-enable-o2 2 #include<bits/stdc++.h> 3 using namespace std; 4 typedef long long ll; 5 typedef double db; 6 #define INF 0x3f3f3f3f 7 #define _for(i,a,b...
P3303 [SDOI2013]淘金
Ajwallet
11-01 231
数位dp，堆维护最优解
sdoi2017R1D2
04-16
sdoi2017R1D2
SDOI2018省选第一试DAY1试题
04-27
本文档为2018年noi山东省选拔第一轮的第一天试题，有需要欢迎下载
SDOI2018省选一轮day1试题数据
04-27
SDOI2018省选一轮day1试题，这是数据部分，题面在另一个文档文件里
2017年山东省选(SDOI)第一轮题目及解析
04-11
2017年山东省选(SDOI)第一轮题目及解析 2017年山东省选(SDOI)第一轮题目及解析
159.二叉树：二叉树的层平均值（力扣）
最新发布
weixin_63779802的博客
05-31 110
代码同样使用了广度优先搜索（BFS）的方法。主要思路是使用一个队列来进行层次遍历，在每一层遍历的时候，计算所有节点的值的总和以及节点的数量，然后用总和除以节点数量得到平均值。这个算法的时间复杂度是 O(n)，因为每个节点都会被访问一次，其中 n 是树中节点的数量。空间复杂度也是 O(n)，因为需要存储整个树的节点。代码是解决“二叉树的层平均值”问题的C++代码，这个问题要求计算二叉树中每一层的节点值的平均值，并返回一个包含这些平均值的列表。
算法训练营day41
stmfresher的博客
05-27 322
动态规划理论基础（主要就是确定动态规划的几个步骤）
一篇文章讲透排序算法之堆排序
duanku111的博客
05-27 451
在学习这篇文章之前，请大家先学习堆这一数据结构中堆的概念，向下调整算法，向下调整建堆。堆的实现堆排序就是利用堆里面学习过的知识点进行排序，如何进行排序呢？
二叉树间奏
chestnut_orenge的博客
05-30 359
已经学习过了二叉树的基础知识，那么间奏部分就来感受一下二叉树部分的OJ题吧。
【Java面试】四、MySQL篇（上）
llg___的博客
05-27 947
如果去维护一个类似二叉树的结构，再查age=45的数据，则直接从根节点开始⇒ 45 > 36，去右侧 ⇒ 45 < 48 ⇒去左侧 ⇒ 查找完毕，如此，查找效率提升，这即索引的思想。是覆盖索引，虽然select * ，但其where是根据id过滤的，即用的是主键索引、聚簇索引，索引的叶子节点存了整行数据，需要返回的字段，在索引中能够全部找到。给name字段加了非聚簇索引，因此，执行如上SQL，先根据name的非聚簇索引的B+树 ⇒ A小于L，走左边，到G和J，再走左边，找到Arm ⇒ 因为是。
【贪心算法】C++ 解决算法题：买卖股票 / K次取反 / 按身高排序 / 优势洗牌
卜及中的博客
05-28 928
贪心算法在每一步选择中都采取当前状态下最优决策的算法。贪心算法通常用来解决最优化问题，其在贪心算法中，我们并不总是考虑到未来可能发生的情况，而是只关注当前的最优选择。这种贪心选择性质使得贪心算法特别适合解决那些具有最优子结构性质的问题，即局部最优解能够推导出全局最优解的问题。需要注意的是，贪心算法并不适用于所有的问题，因为并非所有问题都具有最优子结构性质。在某些情况下，贪心算法得到的结果可能并不是全局最优解，而只是一个较好的解。因此，在应用贪心算法时，需要仔细分析问题的特性，以确定贪心算法是否适用于该问题。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值