疑似乱数２

**デフォルトの名無しさん** · 2007/10/17(水) 22:34:59

擬似乱数発生器について語ろうか。その２

前スレ
擬似乱数
http://pc11.2ch.net/test/read.cgi/tech/1146071975/

関連スレ
【危険】とんでもプログラム告発スレッド【悪質】
http://pc11.2ch.net/test/read.cgi/tech/1191860116/

SIMD-oriented Fast Mersenne Twister (SFMT):
http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/SFMT/index.html

**デフォルトの名無しさん** · 2017/02/24(金) 14:00:51.79

メルセンヌツイスタ使ってなかったんだろ

**デフォルトの名無しさん** · 2017/02/24(金) 19:02:56.20

用途が用途だから暗号論的擬似乱数の適応だろ
/dev/urandomでもCryptGenRandomでも

**デフォルトの名無しさん** · 2017/02/25(土) 15:09:25.02

状態変数が同じになったから同じ数列が出たんじゃねーよ。
二セット分繰り返す数列を出力する状態変数だっただけだろ。

http://science3.2ch.net/test/read.cgi/math/1088351567/424-425
円周率の 93299341 桁目から 07214545 という数列が出るそうだ。
まあ大体その辺で出そうな桁数だ。
おそらく無限に調べれば 1/10,0000,000 より少しだけ少ない頻度で出現するだろう。
そしてそのうちの 1/10,0000,000 は二回続く、つまり 1/10+E16 程度で 0721454507214545 という数列が出現するであろう。

目立つ組み合わせだから話題だが、どの二セットを取っても、その組み合わせが発生する確率は同じだ。

**デフォルトの名無しさん** · 2017/02/25(土) 16:57:29.92

e なら割と早く出て来る

**デフォルトの名無しさん** · 2017/03/10(金) 21:31:43.44

PCGってどうなん？
ttp://www.pcg-random.org/

説明読んだりビデオ観たりした限りでは優れているみたいなんだけども
アルゴリズム的には線形合同法とXorShiftのハイブリッドの様な感じで
結構単純な構造なんだけど、高速で良質な乱数を生成するとか

因みに2年前に発表されているみたいなんだけど知らんかったわ

**デフォルトの名無しさん** · 2017/04/23(日) 03:06:26.50

周期性と分布バランスが目的用途に特化できるならそれが最強だろ

**デフォルトの名無しさん** · 2017/05/11(木) 11:22:40.04

最強(笑)

**デフォルトの名無しさん** · 2018/01/01(月) 10:34:52.74

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

**デフォルトの名無しさん** · 2018/05/23(水) 20:53:42.20

僕の知り合いの知り合いができたパソコン一台でお金持ちになれるやり方
役に立つかもしれません
グーグルで検索するといいかも『ネットで稼ぐ方法　モニアレフヌノ』

BZKM3

**デフォルトの名無しさん** · 2018/07/05(木) 01:00:09.16

5E7

**Goldwasser** · 2019/01/03(木) 13:41:14.69

ご覧あれー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;
}
}

**デフォルトの名無しさん** · 2019/01/03(木) 13:41:50.45

後半

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;
}

**Goldwasser** · 2019/01/03(木) 13:42:22.37

最後

int main() {

rp(x1);
rp(x2);
data();

return 0;
}

**デフォルトの名無しさん** · 2019/01/06(日) 08:57:50.46

誰かTestU01の使い方を教えてください

**デフォルトの名無しさん** · 2019/02/24(日) 02:35:20.31

>>543
http://www.pcg-random.org/posts/how-to-test-with-testu01.html