ГПСЧ на РНР (сугубо для самообразования)

Ммм... а почему для генерации 16-значного ты перемножаешь два 8-значных?
Это же зверски влияет на равномерность распределения результата, в частности, нечетные числа будут выдаваться с вероятностью 1/4, а не 1/2, как должно быть...

Я тут недавно почитал статью про алгоритм для генерации СЧ под названием "Mersenne Twister" (это как раз тот, который используется в mt_rand). Достаточно несложная идея, просто суперские результаты (какой-то просто невообразимый период), мне очень понравилось.

 

А разьве ПХП не может читать из "/dev/random" ?

А ещё может оказаться так, что система поддерживает аппаратный RNG.

 

А какой период данной ПСЧ (не диапазон генерируемых чисел, а именно период повторения последовательности)? Зависит ли он от начальных seed1,seed2,seed3?

 

А почему бы не использовать более простой алгоритм типа: x(n+1)=A*x(n)+B mod M ?
При грамотном выборе A и B период данной последовательности будет равен M, да и распределение можно подобрать приемлемым.

 

Пытаюсь разобраться в математике алгоритма...

Я не совсем понял, где в коде перемножение двух 8-значных?

Это не ошибка? Именно 65535, а не 65536?

 

как вам такая психоделика? :-D

PHP:

1. <?php
2.  
3. ############################
4. # Лицензионнное соглашение #
5. #--------------------------#
6. # Вам   запрещено   читать #
7. # этот  текст.  Немедленно #
8. # закройте глаза  во избе- #
9. # жание   поражения  моска #
10. ############################
11.  
12. # Desu. Desu? Desu!
13. class Desu {
14.     var $data;
15.     function __construct( $props= array() ){
16.         $this->data= (object)$props;
17.     }
18.     function __get( $name ){
19.         $method= 'get__' . $name;
20.         if( method_exists( $this, $method ) )
21.             return $this->{$method}();
22.         return $this->data->{$name};
23.     }
24.     function __set( $name, $value ){
25.         $method= 'set__' . $name;
26.         if( method_exists( $this, $method ) )
27.             return $this->{$method}( $value );
28.         return $this->data->{$name}= $value;
29.     }
30.     function __clone( ){
31.         $this->data= clone $this->data;
32.     }
33. }
34.  
35. # Бродяга, блуждающий по индексному массиву
36. class Stranger extends Desu {
37.     function __construct( $array ){
38.         $this->array= $array;
39.         $this->index= -1;
40.     }
41.     function get__haveNext( ){ # есть куда идти?
42.         return isset( $this->array[ ++$this->index ] );
43.     }
44.     function get__current( ){ # а тут чо?
45.         return $this->array[ $this->index ];
46.     }
47. }
48.  
49. # Историк, он записывает происходящие события
50. class Historian extends Desu {
51.     function __construct( $array= array() ){
52.         $this->memory= $array;
53.     }
54.     function get__tentacle( ){ # Хотим пройтись по памяти? Юзаем щупальцу!
55.         return new Stranger( $this->memory );
56.     }
57.     function push( $values= array() ){ # О! Вспомнил!
58.         $this->memory= array_merge( array_reverse( $values ), $this->memory );
59.         return $values;
60.     }
61.     function pop( ){ # Гм... забыл...
62.         return array_shift( $this->memory );
63.     }
64. }
65.  
66. # Разумное, но тупое как пробка существо
67. class Mind extends Desu {
68.     function __construct( $value= null ){
69.         $this->employers= (object)array();
70.         $this->value= $value;
71.         $this->weight= 1;
72.     }
73.     function touch( ){ # Щупаем
74.         ++$this->weight;
75.     }
76.     function grow( $years ){ # Прошло десять лет
77.         $this->weight/= $years;
78.     }
79.     function extradite( $value ){ # Выдай-ка мне своего подчинённого, занимающегося вышиванием крестиком
80.         return @$this->employers->{$value};
81.     }
82.     function employ( $memoir, $employer ){ # Вот те новый подчинённый для вышивания крестиком
83.         $this->employers->{$memoir}= $employer;
84.     }
85. }
86.  
87. # Моск Нострадамиуса
88. class Brain {
89.     function __construct( $memoirs= array() ){
90.         $this->memory= new Historian( array( 0 ) );
91.         $this->reason= new Mind( );
92.         $this->population= 0;
93.         $this->train( $memoirs );
94.     }
95.     function train( $values= array() ){ # Накачиваем моск историей
96.         foreach( $values as $value ):
97.             $tentacle= $this->memory->tentacle;
98.             $mind= $this->reason;
99.             while( $tentacle->haveNext && ( $employer= $mind->extradite( $tentacle->current ) ) ){
100.                 $mind->touch();
101.                 $mind= $employer;
102.             }
103.             $memoir= @$tentacle->current;
104.             if( !isset( $memoir ) ) $memoir= 0;
105.             if( $value !== $mind->value ) $mind->employ( $memoir, new Mind( $value ) );
106.             $this->memory->push(array( $value ));
107.         endforeach;
108.         $this->population+= count( $values );
109.     }
110.     function predict( ){ # Получаем астрологический прогноз на будущее
111.         $tentacle= $this->memory->tentacle;
112.         $mind= $this->reason;
113.         while( $tentacle->haveNext && ( $employer= $mind->extradite( $tentacle->current ) ) )
114.             $mind= $employer;
115.         return $mind->value;
116.     }
117. }
118.  
119. # Бурный Белый Шумный Поток
120. class WhiteStream extends Desu {
121.     function __construct( $initial= array() ){
122.         $this->generator= new Brain( $initial );
123.     }
124.     function get__value( ){ # а сейчас какое значение?
125.         $gen= $this->generator;
126.         $chit= 1 - $gen->predict();
127.         $gen->train(array( $chit ));
128.         return $chit;
129.     }
130. }
131.  
132. $init= @$_GET[ 'input' ];
133. $init= $init ? explode( ',', $init ) : array();
134. $random= new WhiteStream( $init );
135.  
136. $stat= array();
137. for( $i= 1000; $i; --$i ):
138.     $rnd= $random->value;
139.     @++$stat[ $rnd ];
140.     echo $rnd;
141. endfor;
142.  
143. echo '<pre>';
144. ksort( $stat );
145. print_r( $stat );
146.  
147. ?>

два вечера угробил осталось приделать старение, вымирание при перенаселённости и поддержку любых значений, а не только 0 и 1...

 

guest2013
65535+1(нулевое) значений для 16 бит, 65536 - это уже следующий разряд, 17й.

 

да ладно, не так уж с трудом и писалось

 

sword dancer, юмор местами специфичный, но веселый. (ня!)

 

я знаю! сценарии для хентая пишут на пхп!

PHP:

1. shounen::saw( $tentacle->virgin ) && shounen::comeIn();
2. $tentacle->shounen->woman;
3. Daemon::comeIn();
4. foreach( Daemon::tentacles as $tentacle ) $tentacle->{ rand(2) ? 'man' : 'woman' }
5. Daemon::force( array( 'Big Ben' => woman ) );
6. man::explode( 'tentacle', $daemon );
7. exit( 'happy endo' );

 

sword dancer, *пацталом!* А ты не виртуал dark-demon ?

 

совершенно случайно - нет.

 

Хм... сколько ж их, программирующих анимешников...

 

а что, есть программисты, которые не любят анимэ? палево! это кодеры обычные!
хороший программист всегда, прежде чем писать код, понаблюдает за классами, свыкнится с объектами, проникнится их чувствами, усвоит взаимоотношения. и только после этого начинается действие!

 

Vladson. ИМХО. Судя по asm

Код (Text):

1.  
2.    add ax,dx
3.    rcl ax,1
4.    mov seed3,ax

Последние две строки должны выглядеть так:

PHP:

1. <?
2. //    $a = ($a+$b)&65535; !!! Теряется Carry Flag
3. //    return $this->seed3 = $a = ($a*2 > 65535) ? $a*2-65535 : $a*2; // Не то условие проверки
4.       $a = $a+$b;
5.       return $this->seed3 = ( $a > 65535 ? $a*2+1 : $a*2) & 65535;
6. ?>