元胞自动机 C实现 + OpenCV界面

前面

元胞自动机
建模的时候一直听到但是不知道干啥的
然后找了点代码还是不知道他干嘛的
反正蛮好看就是了orz

https://zh.wikipedia.org/zh-hans/%E7%B4%B0%E8%83%9E%E8%87%AA%E5%8B%95%E6%A9%9F
https://zh.wikipedia.org/wiki/%E5%BA%B7%E5%A8%81%E7%94%9F%E5%91%BD%E6%B8%B8%E6%88%8F

类似这种的自然演变的模型都一副（我就笑笑不说话）的样子
或许有着某些规律把

http://www.cnblogs.com/Firefly727/articles/1856328.html
这是一个很详细（但是排版很糟糕）的介绍

http://www.cnblogs.com/Firefly727/articles/1856328.html
这是一个比较成形的proj。

代码

C++（？）的代码
OpenCV(3.3.1)的界面，顺便取随机数用的也是OpenCV的RNG模块。

main.cpp

//Birdy 2017/12/7#include "opencv2/core.hpp"#include "opencv2/imgproc.hpp"#include "opencv2/highgui.hpp"#include "app.h"#include <stdio.h>using namespace cv;int main(){    char wndname[] = "Demo - Birdy";    int i, width = 1000, height = 800;    RNG rng(0xFFFFFFFF);    Mat image = Mat::zeros(height, width, CV_8UC3);    class ca rand(rng);    rand.nextstep();    for (int t = 0; t < 200; t++)    {        image = Mat(height, width, CV_8UC3, Scalar(0, 0, 0));        for (int i = 0; i < 60; i++)        {            for (int j = 0; j < 60; j++)            {                if (rand.a[i][j] == 1)                {                    circle(image, Point(200 + 10 * i, 100 + 10 * j), 5, Scalar(200 + rng.uniform(0, 50), 100, 0),                        1, LINE_AA);                }            }        }        imshow(wndname, image);        waitKey(50);        rand.nextstep();    }    waitKey();    return 0;}

ca.h

#pragma once#include "opencv2/core.hpp"#include "opencv2/imgproc.hpp"#include <stdio.h>using namespace cv;class ca{public:    bool a[60][60];    ca(RNG &rng);    void nextstep();};

ca.cpp

#include "app.h"#include <math.h>ca::ca(RNG & rng){    for (int i = 1; i < 59; i++)    {        for (int j = 1; j < 59; j++)        {            if (rng.uniform(0, 20) < 2)            {                a[i][j] = 1;            }            else            {                a[i][j] = 0;            }        }    }}void ca::nextstep(){        bool c[60][60];        for (int i = 1; i < 59; i++)        {            for (int j = 1; j < 59; j++)            {                int count = a[i - 1][j - 1] + a[i - 1][j] + a[i - 1][j + 1] + a[i][j - 1] + a[i][j + 1] + a[i + 1][j - 1] + a[i + 1][j] + a[i + 1][j + 1];                if (count == 2)                {                    c[i][j] = a[i][j];                }                else if (count == 3)                {                    c[i][j] = 1;                }                else                {                    c[i][j] = 0;                }            }        }        memcpy(*a, *c, sizeof(bool) * 3600);}

后面

效果

加一幅图（我找不到出处了反正不是我画的orz）、或许是找资料过程中印象最深的东西了
肯定有一天能画小绿小蓝的