結果

問題 No.737 PopCount
コンテスト
ユーザー ikd
提出日時 2018-10-02 18:41:46
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 3 ms / 1,000 ms
コード長 1,324 bytes
コンパイル時間 940 ms
コンパイル使用メモリ 97,792 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-13 01:42:04
合計ジャッジ時間 1,328 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

void main() {
  import std.stdio, std.string, std.conv, std.algorithm;

  long n;
  rd(n);
  const long mod = 10 ^^ 9 + 7;

  int get(long nn) {
    int ret = 0;
    foreach (i; 0 .. 62) {
      if (nn & (1L << i))
        ret = i;
    }
    return ret + 1;
  }

  int len = get(n);
  auto dp = new long[][][](len + 1, len + 1, 2);
  auto o = new long[][][](len + 1, len + 1, 2);
  dp[0][0][0] = 1; // ??
  foreach (i; 0 .. len) {
    for (int j = 0; j <= len; j++) {
      foreach (f; 0 .. 2) {
        foreach (b; 0 .. 2) {
          auto x = (n >> (len - i - 1)) & 1;
          if (f == 0 && x == 0 && b == 1) {
            continue;
          }
          auto nj = j + b, nf = f | (b < x);
          if (nj <= len) {
            (dp[i + 1][nj][nf] += dp[i][j][f]) %= mod;
            (o[i + 1][nj][nf] += o[i][j][f] * 2) %= mod;
            if (b) {
              (o[i + 1][nj][nf] += dp[i][j][f]) %= mod;
            }
          }
        }
      }
    }
  }
  long s = 0;
  for (int j = 1; j <= len; j++) {
    foreach (f; 0 .. 2) {
      (s += o[len][j][f] * j) %= mod;
    }
  }
  writeln(s);
}

void rd(T...)(ref T x) {
  import std.stdio : readln;
  import std.string : split;
  import std.conv : to;

  auto l = readln.split;
  assert(l.length == x.length);
  foreach (i, ref e; x)
    e = l[i].to!(typeof(e));
}
0