细胞分裂游戏

runintostar

回复 10# bikkuri

你需要更清楚地定义下这句话
"死亡后的.在下一年中仍保留但不做任何变化。"
是指死亡的下一年显示, 然后就彻底消失, 还是说死亡的一直都显示出来

bikkuri

一直显示啊。

runintostar 发表于 2014-04-01 11:34
回复 10# bikkuri

你需要更清楚地定义下这句话

yestreenstars

回复 10# bikkuri

script:

1. #!/usr/bin/perl
   
2. use 5.010;
   
3. my @cells;
   
4. my $cell = "o";
   
5. my $count = 0;
   
6. while ( $count < 100 ) {
   
7.         my @each = $cell =~ /./g;
   
8.         for my $i ( 0 .. $#each ) {
   
9.                 given ( rand() ) {
   
10.                         when ( $_ < 1/10 ) { $each[$i] = "." }
    
11.                         when ( $_ < 2/10 ) { $each[$i] = "o." if $each[$i] eq "8" }
    
12.                         default {
    
13.                                 given ( $each[$i] ) {
    
14.                                         when ( "o" ) { $each[$i] = "0" }
    
15.                                         when ( "0" ) { $each[$i] = "8" }
    
16.                                         when ( "8" ) { $each[$i] = "oo" }
    
17.                                 }
    
18.                         }
    
19.                 }
    
20.         }
    
21.         $cell = join "", @each;
    
22.         push @cells, $cell;
    
23.         $count++;
    
24.         say "[$count]\t$cell";
    
25.         $cell =~ /^\.+$/ && exit;
    
26.         sleep 1;
    
27. }
    
28. say "-" x 50;
    
29. say join "\t", "", "[.]", "[o]", "[0]", "[8]", "[.o08]";
    
30. for my $i ( 0 .. $#cells ) {
    
31.         my $f5 = $cells[$i] =~ s/\.o08//g;
    
32.         my $f4 = $cells[$i] =~ s/8//g;
    
33.         my $f3 = $cells[$i] =~ s/0//g;
    
34.         my $f2 = $cells[$i] =~ s/o//g;
    
35.         my $f1 = $cells[$i] =~ s/\.//g;
    
36.         say join "\t", "[".($i+1)."]:", $f1||0, $f2||0, $f3||0, $f4||0, $f5||0;
    
37. }

复制代码

result:
result.rar (155.87 KB, 下载次数: 4)

2014-04-01 12:20 上传

点击文件名下载附件

   

bikkuri

哇！这个程序很完美啊！而且运算速度很快啊，每代的计算都可以在1秒钟之内完成啊！
为什么用awk效率那么低呢？
但是后面世代细胞数的计算有一个地方错了：
[.o08]是指所有的细胞数，而不是指.o08死幼成分四个连续状态的细胞组合数。
[.o08]应该是[.]+[o]+[0]+[8]，把四种状态的细胞数量加起来就好。

1.         [.]        [o]        [0]        [8]        [.o08]
   
2. [90]:   2513    800     670     524     4507
   
3. [91]:   2769    974     701     532     4976
   
4. [92]:   3040    1013    846     565     5464
   
5. [93]:   3344    1068    898     671     5981
   
6. [94]:   3687    1243    948     704     6582
   
7. [95]:   4016    1347    1112    747     7222
   
8. [96]:   4422    1415    1179    883     7899
   
9. [97]:   4890    1610    1239    941     8680
   
10. [98]:   5381    1752    1399    984     9516
    
11. [99]:   5906    1823    1543    1128    10400
    
12. [100]:  6458    2117    1620    1222    11417

复制代码

看来即使是10%的死亡率，死去的细胞数都超过总细胞数的一半以上啊。。。
还有就是perl脚本可以在shell下直接执行，但是不能用tee或者重定向把结果重定向到文件啊？
这次执行产生的细胞数差不多是上一次的三倍：

1.         [.]        [o]        [0]        [8]        [.o08]
   
2. [90]:   6993    2179    1851    1350    12373
   
3. [91]:   7708    2479    1917    1471    13575
   
4. [92]:   8417    2783    2179    1545    14924
   
5. [93]:   9264    2839    2456    1735    16294
   
6. [94]:   10138   3268    2515    1933    17854
   
7. [95]:   11081   3602    2883    2036    19602
   
8. [96]:   12111   3815    3225    2278    21429
   
9. [97]:   13285   4278    3345    2576    23484
   
10. [98]:   14574   4765    3768    2675    25782
    
11. [99]:   15966   4979    4204    3009    28158
    
12. [100]:  17503   5642    4374    3363    30882

复制代码

yestreenstars 发表于 2014-04-01 12:21
回复 10# bikkuri

script:result:

yinyuemi

回复 14# bikkuri


    你把awk的sleep去掉试试

bikkuri

去掉sleep以后，速度没什么变化。
到80代以后速度很慢。
到95代以后计算一代需要十分钟以上。
计算第100代需要半个小时以上。

yinyuemi 发表于 2014-04-01 13:22
回复 14# bikkuri

yestreenstars

回复 14# bikkuri

你之前又没说清楚，咋知道你什么意思？
简单修改一下就好了。

1. #!/usr/bin/perl
   
2. use 5.010;
   
3. my @cells;
   
4. my $cell = "o";
   
5. my $count = 0;
   
6. while ( $count < 100 ) {
   
7.         my @each = $cell =~ /./g;
   
8.         for my $i ( 0 .. $#each ) {
   
9.                 given ( rand() ) {
   
10.                         when ( $_ < 1/10 ) { $each[$i] = "." }
    
11.                         when ( $_ < 2/10 ) { $each[$i] = "o." if $each[$i] eq "8" }
    
12.                         default {
    
13.                                 given ( $each[$i] ) {
    
14.                                         when ( "o" ) { $each[$i] = "0" }
    
15.                                         when ( "0" ) { $each[$i] = "8" }
    
16.                                         when ( "8" ) { $each[$i] = "oo" }
    
17.                                 }
    
18.                         }
    
19.                 }
    
20.         }
    
21.         $cell = join "", @each;
    
22.         push @cells, $cell;
    
23.         $count++;
    
24.         say "[$count]\t$cell";
    
25.         $cell =~ /^\.+$/ && exit;
    
26.         sleep 1;
    
27. }
    
28. say "-" x 50;
    
29. say join "\t", "", "[.]", "[o]", "[0]", "[8]", "[.o08]";
    
30. for my $i ( 0 .. $#cells ) {
    
31.         my $f1 = $cells[$i] =~ s/\.//g;
    
32.         my $f2 = $cells[$i] =~ s/o//g;
    
33.         my $f3 = $cells[$i] =~ s/0//g;
    
34.         my $f4 = $cells[$i] =~ s/8//g;
    
35.         my $f5 = $f1 + $f2 + $f3 + $f4;
    
36.         say join "\t", "[".($i+1)."]:", $f1||0, $f2||0, $f3||0, $f4||0, $f5||0;
    
37. }

复制代码

我的result.txt就是用重定向得到的。
   

bikkuri

你的重定向命令怎么写的？

yestreenstars

回复 14# bikkuri

根据我的perl脚本改写的awk脚本：

1. #!/bin/awk -f
   
2. BEGIN {
   
3.         srand()
   
4.         s="o"
   
5.         for(i=0;i++<100;){
   
6.                 l=split(s,a,"")
   
7.                 for(j=0;j++<l;){
   
8.                         r=rand()
   
9.                         if(r<1/10){
   
10.                                 a[j]="."
    
11.                         }else if(r<2/10){
    
12.                                 if(a[j]=="8")a[j]="o."
    
13.                         }else{
    
14.                                 if(a[j]=="8")a[j]="oo"
    
15.                                 if(a[j]=="0")a[j]="8"
    
16.                                 if(a[j]=="o")a[j]="0"
    
17.                         }
    
18.                 }
    
19.                 s=""
    
20.                 for(j=0;j++<=l;)s=s!=""?s""a[j]:a[j]
    
21.                 b[i]=s
    
22.                 printf "[%d]\t%s\n",i,s
    
23.                 if(s~/^\.+$/)exit
    
24.                 system("sleep 1")
    
25.         }
    
26.         for(i=0;i++<10;)printf "-----"
    
27.         print "\n\t[.]\t[o]\t[0]\t[8]\t[.o08]"
    
28.         for(i=0;i++<length(b);){
    
29.                 f5=length(b[i])
    
30.                 f4=gsub(/8/,"",b[i])
    
31.                 f3=gsub(/0/,"",b[i])
    
32.                 f2=gsub(/o/,"",b[i])
    
33.                 f1=f5-f4-f3-f2
    
34.                 printf "[%d]:\t%d\t%d\t%d\t%d\t%d\n",i,f1,f2,f3,f4,f5
    
35.         }
    
36. }

复制代码

效率测试：

1. [root@localhost ~]# time ./a > a.r
   
2. 
   
3. real    1m43.244s
   
4. user    0m1.535s
   
5. sys     0m1.153s
   
6. [root@localhost ~]# wc a.r
   
7.    202    806 297136 a.r
   
8. [root@localhost ~]# time ./p > p.r
   
9. 
   
10. real    1m40.743s
    
11. user    0m0.251s
    
12. sys     0m0.036s
    
13. [root@localhost ~]# wc p.r
    
14.   202   806 84208 p.r
    
15. [root@localhost ~]#

复制代码

从上面的数据可知，awk的效率不比perl低。
   

yestreenstars

回复 18# bikkuri

从我19楼的效率测试你应该知道我是怎么重定向的了吧？