http://blog.tenhou.net/article/30503297.html
より転載

2009年07月13日
牌山生成の公開方法
牌山の生成方法を完全公開にする予定です。

公開する情報は「配山の生成手法」と「乱数シード」になります。

これによって
１．作為的に配牌を選択していない
２．特定の局面で山を操作していない
ことが牌譜と対戦ログから検証可能になります。
※現状でも２．は検証可能です。

以下試験的に作成したサンプルコードです。
(MTやSHAについては省略)
MT : http://www.google.co.jp/search?q=mt19937ar.c
SHA : http://www.google.co.jp/search?q=sha512.c

サンプルコードについてのご意見・アドバイスなどなどは support@c-egg.comまでお願いします。

公開開始スケジュールは準備が整い次第告知していきます。


■サンプルコード
void SampleYamaShuffle(){
    static const char *haiDisp[34]={
        "一","二","三","四","五","六","七","八","九",
        "①","②","③","④","⑤","⑥","⑦","⑧","⑨",
        "１","２","３","４","５","６","７","８","９",
        "東","南","西","北","白","發","中"
    };

    WORD wShuffleVersionMajor=1;
    int i;
    _MTRAND mtRoot;
    TCHAR szSeedSeqName[32];
    {
        DWORD seed[MTRAND_N]; // これは公開されない
        if (1){
            // 実際のサーバ動作(予定)
            HCRYPTPROV hCP; // for Win32Only
            if (!CryptAcquireContext(&hCP,NULL,NULL,PROV_RSA_FULL,0)) throw 0;
            if (!CryptGenRandom(hCP,sizeof(seed),(BYTE*)seed)) throw 0;
            if (!CryptReleaseContext(hCP,0)) throw 0;

            SYSTEMTIME st; // 日付以外を使用する可能性アリ
            GetLocalTime(&st);
            wsprintf(szSeedSeqName,_T("%d.%04d.%02d%02d.%02d%02d"),
                wShuffleVersionMajor,
                st.wYear,st.wMonth,st.wDay,st.wHour,st.wMinute);
            printf("seedDailyPublic=%s",szSeedSeqName);
            if (1) for(i=0;i<sizeof(seed)/sizeof(*seed);++i) printf(",%08X",seed[i]);
            printf("\r\n");
        }else{
            // 検証用の入力 // szSeedSeqNameをキーに検索
            for(i=0;i<sizeof(seed)/sizeof(*seed);++i) seed[0]=/* REPLACE HERE */0;
        }
        mtRoot.init_by_array(seed,sizeof(seed)/sizeof(*seed));
    }

    int nGame=0;
    for(;nGame<10;++nGame){ // １０回目の対戦まで
        mtRoot.genrand_int32(); // 席順決定などで３つ消費
        mtRoot.genrand_int32();
        mtRoot.genrand_int32(); // サンマは未使用
        mtRoot.genrand_int32(); // 未使用

        _MTRAND mtLocal;
        {
            DWORD seed[MTRAND_N]; // 136!より十分に大きく
            for(i=0;i<sizeof(seed)/sizeof(*seed);++i) seed[i]=mtRoot.genrand_int32();
            mtLocal.init_by_array(seed,sizeof(seed)/sizeof(*seed));
            printf("mt.seed=%s,%d",szSeedSeqName,nGame);
            if (1) for(i=0;i<sizeof(seed)/sizeof(*seed);++i) printf(",%08X",seed[i]);
            printf("\r\n");
        }
        // ここまででmtRootは合計(MTRAND_N+4)tick

        // ここで牌譜にszSeedSeqName,nGame,seedなどを出力

        int nKyoku=0;
        for(;nKyoku<10;++nKyoku){ // 配牌が１０回行われた場合
            DWORD rnd[SHA512_DIGEST_SIZE/sizeof(DWORD)*9]; // 135+2以上を確保
            {
                DWORD src[sizeof(rnd)/sizeof(*rnd)*2]; // 1024bit単位で512bitへhash
                for(i=0;i<sizeof(src)/sizeof(*src);++i) src[i]=mtLocal.genrand_int32();
                for(i=0;i<sizeof(rnd)/SHA512_DIGEST_SIZE;++i){
                    SHA2::sha512_ctx ctx;
                    SHA2::sha512_init(&ctx);
                    SHA2::sha512_update(&ctx,(BYTE*)src+i*SHA512_DIGEST_SIZE*2,SHA512_DIGEST_SIZE*2); // in=1024bit
                    SHA2::sha512_final (&ctx,(BYTE*)rnd+i*SHA512_DIGEST_SIZE); // out=512bit
                }
            }

            BYTE yama[136]; // サンマは108
            for(i=0;i<136;++i) yama[i]=i;
            for(i=0;i<136-1;++i) swap(yama[i],yama[i + (rnd[i]%(136-i))]); // 1/2^32以下の誤差は許容

            printf("nGame=%d nKyoku=%d yama=",nGame,nKyoku);
            for(i=0;i<136;++i) printf("%s",haiDisp[yama[i]/4]);
            printf("\r\n");

            int dice0=rnd[135]%6;
            int dice1=rnd[136]%6;
            // rnd[137]～rnd[143]は未使用
        }
    }
}