結果

問題 No.702 中央値を求めよ LIMITED
コンテスト
ユーザー te-sh
提出日時 2018-07-04 10:17:01
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 134 ms / 5,000 ms
コード長 988 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 676 ms
コンパイル使用メモリ 102,604 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-13 01:22:50
合計ジャッジ時間 5,134 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 25
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import std.algorithm, std.container, std.conv, std.math, std.range, std.typecons, std.stdio, std.string;

auto rdsp(){return readln.splitter;}
void pick(R,T)(ref R r,ref T t){t=r.front.to!T;r.popFront;}
void readV(T...)(ref T t){auto r=rdsp;foreach(ref v;t)pick(r,v);}

const n = 10000001;

void main()
{
  uint seed; readV(seed);

  init(seed);
  auto a = new int[](1<<16);
  foreach (_; 0..n) {
    auto r = generate();
    ++a[r>>16];
  }

  auto i = 0u, s = 0;
  for (; i < (1<<16); ++i) {
    if (s+a[i] >= (n+1)/2) break;
    s += a[i];
  }

  init(seed);
  a[] = 0;
  foreach (_; 0..n) {
    auto r = generate();
    if ((r>>16) == i)
      ++a[r&((1<<16)-1)];
  }

  foreach (j; 0..1u<<16) {
    s += a[j];
    if (s >= (n+1)/2) {
      writeln((i<<16)|j);
      return;
    }
  }
}

uint x, y, z, w;

void init(uint seed)
{
  x = seed;
  y = 1;
  z = 2;
  w = 3;
}
uint generate() {
  uint t = (x^(x<<11));
  x = y;
  y = z;
  z = w;
  w = (w^(w>>19))^(t^(t>>8)); 
  return w;
}
0