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

标题:生命游戏

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

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

具体来说:

  1. 当前细胞为存活状态时,当周围低于2个(不包含2个)存活细胞时, 该细胞变成死亡状态。(模拟生命数量稀少)
  2. 当前细胞为存活状态时,当周围有2个或3个存活细胞时, 该细胞保持原样。
  3. 当前细胞为存活状态时,当周围有3个以上的存活细胞时,该细胞变成死亡状态。(模拟生命数量过多)
  4. 当前细胞为死亡状态时,当周围有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
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

大风车滴呀滴溜溜地转

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值