結果

問題 No.702 中央値を求めよ LIMITED
ユーザー yosupotyosupot
提出日時 2018-06-15 23:00:01
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 56 ms / 5,000 ms
コード長 1,329 bytes
コンパイル時間 5,357 ms
コンパイル使用メモリ 192,792 KB
最終ジャッジ日時 2025-01-05 14:03:13
ジャッジサーバーID
（参考情報） 		judge5 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
using uint = unsigned int;
using ll = long long;
using ull = unsigned long long;
constexpr ll TEN(int n) { return (n==0) ? 1 : 10*TEN(n-1); }
template<class T> using V = vector<T>;
template<class T> using VV = V<V<T>>;

uint32_t x = 0, y = 1, z = 2, w = 3;
uint32_t generate() { 
    uint32_t t = (x^(x<<11));
    x = y;
    y = z;
    z = w;
    w = (w ^ (w >> 19)) ^ (t ^ (t >> 8)); 
    return w;
}

int main(void) {
    ll seed; cin >> seed; 
    uint k = 5 * TEN(6);

    V<uint> bk(1<<16);
    fill(begin(bk), end(bk), 0);

    x = seed;
    for (int i = 0; i < 10000001; i++) {
        uint32_t a = generate();
        bk[a >> 16]++;
    }

    uint up = 0;
    for (int i = 0; i < (1<<16); i++) {
        if (k < bk[i]) {
            up = i;
            break;
        }
        k -= bk[i];
    }

    fill(begin(bk), end(bk), 0);
    x = seed; y = 1; z = 2; w = 3;
    for (int i = 0; i < 10000001; i++) {
        uint32_t a = generate();
        if ((a >> 16) != up) continue;
        bk[a % (1<<16)]++;
    }

    uint lw = 0;
    for (int i = 0; i < (1<<16); i++) {
        if (k < bk[i]) {
            lw = i;
            break;
        }
        k -= bk[i];
    }
    
//    cout << up << " " << lw << endl;
    cout << ((up << 16) + lw) << endl;
    
    return 0;
}
0