疑似乱数２ - 暇つぶし2ch

疑似乱数２ at TECH

534:デフォルトの名無しさん
17/02/23 13:59:17.36 7wYQlXII.net
一つずつ乱数生成してないのかもね

535:デフォルトの名無しさん
17/02/23 16:45:07.24 eG9NqkD/.net
>>523
真じゃない乱数は
ループの長さが長いか短いかだけで
いずれはループするんだよな

536:デフォルトの名無しさん
17/02/23 16:46:44.91 WFbVQHcp.net
>>525
なんらかのバグで内部状態が変化しなくなって
２回連続で同じ値を出してしまったら
再起動がかかるシステムだったりするかもしれない

537:デフォルトの名無しさん
17/02/23 21:50:11.33 5geihC2S.net
>>522
「前回(の乱数取得時)と同じ」という意味ではなく
「前回試合予想をした時と同じ」という意味ね
どちらも、同じ値から同じ試合の予想をして、
結果、前週と全く同じ予想が14試合x5口分で起きた

538:デフォルトの名無しさん
17/02/24 14:00:51.79 xRGcfmim.net
メルセンヌツイスタ使ってなかったんだろ

539:デフォルトの名無しさん
17/02/24 19:02:56.20 3z/u0Cfz.net
用途が用途だから暗号論的擬似乱数の適応だろ
/dev/urandomでもCryptGenRandomでも

540:デフォルトの名無しさん
17/02/25 15:09:25.02 u8Ry8YI7.net
状態変数が同じになったから同じ数列が出たんじゃねーよ。
二セット分繰り返す数列を出力する状態変数だっただけだろ。
ｽﾚﾘﾝｸ(math板:424-425番)
円周率の 93299341 桁目から 07214545 という数列が出るそうだ。
まあ大体その辺で出そうな桁数だ。
おそらく無限に調べれば 1/10,0000,000 より少しだけ少ない頻度で出現するだろう。
そしてそのうちの 1/10,0000,000 は二回続く、つまり 1/10+E16 程度で 0721454507214545 という数列が出現するであろう。
目立つ組み合わせだから話題だが、どの二セットを取っても、その組み合わせが発生する確率は同じだ。

541:デフォルトの名無しさん
17/02/25 16:57:29.92 usTDxsWv.net
e なら割と早く出て来る

542:デフォルトの名無しさん
17/03/10 21:31:43.44 ARsP6JfD.net
PCGってどうなん？
URLﾘﾝｸ(www.pcg-random.org)
説明読んだりビデオ観たりした限りでは優れているみたいなんだけども
アルゴリズム的には線形合同法とXorShiftのハイブリッドの様な感じで
結構単純な構造なんだけど、高速で良質な乱数を生成するとか
因みに2年前に発表されているみたいなんだけど知らんかったわ

543:デフォルトの名無しさん
17/04/23 03:06:26.50 H8Cvp+NU.net
周期性と分布バランスが目的用途に特化できるならそれが最強だろ

544:デフォルトの名無しさん
17/05/11 11:22:40.04 yPF7Zec4.net
最強(笑)

545:デフォルトの名無しさん
18/01/01 10:34:52.74 Adreh4LC.net
Xn+

546:1=Xn*201+8+1/1 Xn+1=Xn*(2018>>1)+1 あけおめ

547:デフォルトの名無しさん
18/05/23 20:53:42.20 Au5e7VGg.net
僕の知り合いの知り合いができたパソコン一台でお金持ちになれるやり方
役に立つかもしれません
グーグルで検索するといいかも『ネットで稼ぐ方法　モニアレフヌノ』
BZKM3

548:デフォルトの名無しさん
18/07/05 01:00:09.16 RfoszcD2.net
5E7

549:Goldwasser
19/01/03 13:41:14.69 r8gE2use.net
ご覧あれーw
#include "pch.h"
#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#define N 32
unsigned char x1[N], x2[N], x3[N];
void rp(unsigned char* a) {
int i, j, x;
for (i = 0; i < N; i++) {
a[i] = i;
}
for (i = 0; i < N - 2; i++) {
// rand from i+1 to N-1
j = (rand() % (N - 1 - i)) + i + 1;
// swap a[i] and a[j]
x = a[j];
a[j] = a[i];
a[i] = x;
}
if (a[N - 1] == N - 1) {
a[N - 1] = a[N - 2];
a[N - 2] = N - 1;
}
}

550:デフォルトの名無しさん
19/01/03 13:41:50.45 r8gE2use.net
後半
int data() {
int i, j = 0, k = 0;
unsigned int a[N];
unsigned int z[N];
unsigned char w[N];
for (i = 0; i < N; i++)
a[i] = rand()%256;
for (i = 0; i < N; i++)
z[i] = 0;
k = 0;
while (k< 4000) {
for (i = 0; i < N; i++)
z[i] ^= a[x2[i]];
for (i = 0; i < N; i++)
a[i] ^= z[i];
for (i = 0; i < N; i++)
w[i] = x1[x2[x3[i]]];
*x2 = *w;
k++;
for(i=0;i<N;i++)
printf("%u,",a[i]);
}
printf("\n");
return 0;
}

551:Goldwasser
19/01/03 13:42:22.37 r8gE2use.net
最後
int main() {
rp(x1);
rp(x2);
data();
return 0;
}

552:デフォルトの名無しさん
19/01/06 08:57:50.46 kWOVO8kw.net
誰かTestU01の使い方を教えてください

553:デフォルトの名無しさん
19/02/24 02:35:20.31 LCfXPkf7.net
>>542
URLﾘﾝｸ(www.pcg-random.org)

554:過去ログ ★
[過去ログ]
■ このスレッドは過去ログ倉庫に格納されています