1 名前:デフォルトの名無しさん [2007/10/17(水) 22:34:59 .net] 擬似乱数発生器について語ろうか。その2 前スレ 擬似乱数 pc11.2ch.net/test/read.cgi/tech/1146071975/ 関連スレ 【危険】とんでもプログラム告発スレッド【悪質】 pc11.2ch.net/test/read.cgi/tech/1191860116/ SIMD-oriented Fast Mersenne Twister (SFMT): www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/SFMT/index.html
549 名前:Goldwasser [2019/01/03(木) 13:41:14.69 ID: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 名前:デフォルトの名無しさん [2019/01/03(木) 13:41:50.45 ID: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 [2019/01/03(木) 13:42:22.37 ID:r8gE2use.net] 最後 int main() { rp(x1); rp(x2); data(); return 0; }
552 名前:デフォルトの名無しさん [2019/01/06(日) 08:57:50.46 ID:kWOVO8kw.net] 誰かTestU01の使い方を教えてください
553 名前:デフォルトの名無しさん mailto:sage [2019/02/24(日) 02:35:20.31 ID:LCfXPkf7.net] >>542 www.pcg-random.org/posts/how-to-test-with-testu01.html
554 名前:過去ログ ★ [[過去ログ]] ■ このスレッドは過去ログ倉庫に格納されています