2017/National _C_C++_A/2/生命游戏

最新推荐文章于 2024-01-20 15:21:19 发布

大风车滴呀滴溜溜地转

最新推荐文章于 2024-01-20 15:21:19 发布

阅读量488

点赞数 1

分类专栏： # 蓝桥杯

本文链接：https://blog.csdn.net/weixin_43336281/article/details/89889377

版权

蓝桥杯专栏收录该内容

202 篇文章 25 订阅

订阅专栏

标题：生命游戏

康威生命游戏是英国数学家约翰·何顿·康威在1970年发明的细胞自动机。
这个游戏在一个无限大的2D网格上进行。

初始时，每个小方格中居住着一个活着或死了的细胞。
下一时刻每个细胞的状态都由它周围八个格子的细胞状态决定。

具体来说：

当前细胞为存活状态时，当周围低于2个（不包含2个）存活细胞时，该细胞变成死亡状态。（模拟生命数量稀少）
当前细胞为存活状态时，当周围有2个或3个存活细胞时，该细胞保持原样。
当前细胞为存活状态时，当周围有3个以上的存活细胞时，该细胞变成死亡状态。（模拟生命数量过多）
当前细胞为死亡状态时，当周围有3个存活细胞时，该细胞变成存活状态。（模拟繁殖）

当前代所有细胞同时被以上规则处理后, 可以得到下一代细胞图。按规则继续处理这一代的细胞图，可以得到再下一代的细胞图，周而复始。

例如假设初始是:(X代表活细胞，.代表死细胞)

.....
.....
.XXX.
.....

下一代会变为:

.....
..X..
..X..
..X..
.....

康威生命游戏中会出现一些有趣的模式。例如稳定不变的模式：

....
.XX.
.XX.
....

还有会循环的模式：

......      ......       ......
.XX...      .XX...       .XX...
.XX...      .X....       .XX...
...XX.   -> ....X.  ->   ...XX.
...XX.      ...XX.       ...XX.
......      ......       ......

本题中我们要讨论的是一个非常特殊的模式，被称作"Gosper glider gun"：

......................................
.........................X............
.......................X.X............
.............XX......XX............XX.
............X...X....XX............XX.
.XX........X.....X...XX...............
.XX........X...X.XX....X.X............
...........X.....X.......X............
............X...X.....................
.............XX.......................
......................................

假设以上初始状态是第0代，请问第1000000000(十亿)代一共有多少活着的细胞？

注意：我们假定细胞机在无限的2D网格上推演，并非只有题目中画出的那点空间。
当然，对于遥远的位置，其初始状态一概为死细胞。

注意：需要提交的是一个整数，不要填写多余内容。

Ideas

这一看1000000000(十亿)就知道，暴力迭代肯定run不出来，必有规律，模拟几次找到规律就可以了。
在这里插入图片描述
模拟了前101代，可以看出来大致有规律，接下来就是具体分析规律是啥。

Code

C++

#include <map>
#include <iostream>

#define MAP map<Node,bool>

using namespace std;

int Left=0,Top=0,Right=40,Bottom=15;
int index=1,count=36;

typedef struct Node{
	int x,y;
	Node(int x,int y):x(x),y(y){
	}
	bool operator < (const Node& o) const{
		if(x!=o.x) return x<o.x;
		return y<o.y;
	}
}Node;

MAP pre,now;

int refresh_bound(int x,int y){
	Top=min(Top,x);
	Bottom=max(Bottom,x);
	Left=min(Left,y);
	Right=max(Right,y);
}

int get_neighbor(MAP &mp,int x,int y){
	int ans=0;
	for(int i=-1;i<2;i++)
		for(int j=-1;j<2;j++)
			if((i||j)&&(mp[Node(x+i,y+j)]))
				ans++;
	return ans;
}

void epoch(){
	pre=now;
	now=MAP();
	
	int tmp_count=0;
	for(int i=Top-1;i<Bottom+1;i++){
		for(int j=Left-1;j<Right+1;j++){
			int nums=get_neighbor(pre,i,j);
			if(pre[Node(i,j)]){		//活细胞 
				if(nums>1&&nums<4){		//周围细胞适合 
					now[Node(i,j)]=1;
					tmp_count++;
				}
			}else{		//死细胞 
				if(nums==3){	//繁殖 
					now[Node(i,j)]=1;
					refresh_bound(i,j);
					tmp_count++;
				}
			}
		}
	}
	cout<<index<<' '<<tmp_count<<' '<<tmp_count-count<<endl;
	count=tmp_count;
	index++;
}

void display(MAP &mp){
	for(int i=Top;i<Bottom+1;i++){
		for(int j=Left;j<Right+1;j++){
			if(mp[Node(i,j)]) putchar('X');
			else putchar(' ');
		}
		puts("");
	} 
	puts("");
}

int main(){
	freopen("input.txt","r",stdin);
	freopen("output.txt","w",stdout);
	
	char buf[100];
	for(int i=0;i<Bottom+1;i++){
		scanf("%s",buf);
		for(int j=0;j<Right+1;j++){
			if(buf[j]=='X'){
				now[Node(i,j)]=1;
			}
		}
	}
	for(int i=0;i<300;i++)
		epoch();
	return 0;
}

Python

from copy import deepcopy


def countAround(m, x, y):
    res = 0
    for i in [-1, 0, 1]:
        for j in [-1, 0, 1]:
            try:
                res += 1 if not (i == 0 and j == 0) and m[x + i][y + j] == 'X' else 0
            except IndexError:
                continue
    return res


def countTotalX(m):
    res = 0
    for row in m:
        for col in row:
            if col == 'X':
                res += 1
    return res


def simulate(maps):
    count = []
    newMaps = [['.' for _ in range(1000)] for _ in range(1000)]
    for i in range(len(maps)):
        for j in range(len(maps[i])):
            if maps[i][j] == 'X':
                newMaps[500 + i][500 + j] = maps[i][j]
    # print(countNumber(newMaps))

    for _ in range(100):
        maps = deepcopy(newMaps)
        count.append(countTotalX(newMaps))
        print(count)
        for i in range(1000):
            for j in range(1000):
                cnt = countAround(maps, i, j)
                if maps[i][j] == 'X':
                    if cnt < 2 or cnt > 3:
                        newMaps[i][j] = '.'
                    else:
                        newMaps[i][j] = 'X'
                else:
                    if cnt == 3:
                        newMaps[i][j] = 'X'


if __name__ == '__main__':
    simulate([
        ['.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.',
         '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', ],
        ['.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.',
         '.', '.', '.', 'X', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', ],
        ['.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.',
         '.', 'X', '.', 'X', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', ],
        ['.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', 'X', 'X', '.', '.', '.', '.', '.', '.', 'X',
         'X', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', 'X', 'X', '.', ],
        ['.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', 'X', '.', '.', '.', 'X', '.', '.', '.', '.', 'X',
         'X', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', 'X', 'X', '.', ],
        ['.', 'X', 'X', '.', '.', '.', '.', '.', '.', '.', '.', 'X', '.', '.', '.', '.', '.', 'X', '.', '.', '.', 'X',
         'X', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', ],
        ['.', 'X', 'X', '.', '.', '.', '.', '.', '.', '.', '.', 'X', '.', '.', '.', 'X', '.', 'X', 'X', '.', '.', '.',
         '.', 'X', '.', 'X', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', ],
        ['.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', 'X', '.', '.', '.', '.', '.', 'X', '.', '.', '.', '.',
         '.', '.', '.', 'X', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', ],
        ['.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', 'X', '.', '.', '.', 'X', '.', '.', '.', '.', '.',
         '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', ],
        ['.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', 'X', 'X', '.', '.', '.', '.', '.', '.', '.',
         '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', ],
        ['.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.',
         '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', '.', ],
    ])
    """
    0  -> 36
                +12
    10 -> 48
                +6
    20 -> 54
                -13
    30 -> 41
                +12
    40 -> 53
                +6
    50 -> 59
                -13
    60 -> 46
                +12
    70 -> 58
                +6
    80 -> 64
                -13
    90 -> 51
    """
    ans, pool = 36, [12, 6, -13]
    for i in range(1000000000 // 10):
        ans += pool[i % 3]
    print(ans)

大风车滴呀滴溜溜地转

关注

1
点赞
踩
1

收藏

觉得还不错? 一键收藏
打赏
0
评论
2017/National _C_C++_A/2/生命游戏

标题：生命游戏康威生命游戏是英国数学家约翰·何顿·康威在1970年发明的细胞自动机。这个游戏在一个无限大的2D网格上进行。初始时，每个小方格中居住着一个活着或死了的细胞。下一时刻每个细胞的状态都由它周围八个格子的细胞状态决定。具体来说：当前细胞为存活状态时，当周围低于2个（不包含2个）存活细胞时，该细胞变成死亡状态。（模拟生命数量稀少）当前细胞为存活状态时，当周围有2个或3个存活...
复制链接

扫一扫