結果

問題 No.737 PopCount
コンテスト
ユーザー hos.lyric
提出日時 2019-01-31 13:42:08
言語 D
(dmd 2.109.1)
結果
AC  
実行時間 2 ms / 1,000 ms
コード長 2,339 bytes
コンパイル時間 738 ms
コンパイル使用メモリ 121,240 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-13 04:08:33
合計ジャッジ時間 1,586 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

import std.conv, std.functional, std.stdio, std.string;
import std.algorithm, std.array, std.bigint, std.container, std.math, std.numeric, std.range, std.regex, std.typecons;
import core.bitop;

class EOFException : Throwable { this() { super("EOF"); } }
string[] tokens;
string readToken() { for (; tokens.empty; ) { if (stdin.eof) { throw new EOFException; } tokens = readln.split; } auto token = tokens.front; tokens.popFront; return token; }
int readInt() { return readToken.to!int; }
long readLong() { return readToken.to!long; }
real readReal() { return readToken.to!real; }

bool chmin(T)(ref T t, in T f) { if (t > f) { t = f; return true; } else { return false; } }
bool chmax(T)(ref T t, in T f) { if (t < f) { t = f; return true; } else { return false; } }

int binarySearch(alias pred, T)(in T[] as) { int lo = -1, hi = cast(int)(as.length); for (; lo + 1 < hi; ) { const mid = (lo + hi) >> 1; (unaryFun!pred(as[mid]) ? hi : lo) = mid; } return hi; }
int lowerBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a >= val)); }
int upperBound(T)(in T[] as, T val) { return as.binarySearch!(a => (a > val)); }


enum MO = 10L^^9 + 7;

// sum of [a, a + d)
long calc(long a, long d) {
  long x = a + a + d - 1;
  long y = d;
  if (x % 2 == 0) {
    x /= 2;
  } else {
    y /= 2;
  }
  x %= MO;
  y %= MO;
  return (x * y) % MO;
}

void main() {
  try {
    for (; ; ) {
      const N = readLong();
      
      long ans;
      foreach (k; 0 .. 60) {
        if (((N + 1) >> k) & 1) {
          const a = (N + 1) & ~((1L << (k + 1)) - 1);
          debug {
            writefln("%s: [%s, %s)", k, a, a + 2L^^k);
          }
          ans += popcnt(a) * calc(a, 2L^^k);
          ans %= MO;
          foreach (l; 0 .. k) {
            // [2^l, 2^(l+1)) mod 2^(l+1)
            long tmp;
            tmp += (a % MO) * ((2L^^k / 2) % MO) % MO;
            tmp += calc(0, 2L^^(k - l - 1)) * (2L^^(l + 1) % MO) % MO * (2L^^l % MO) % MO;
            tmp += calc(2L^^l, 2L^^l) * (2L^^(k - l - 1) % MO) % MO;
            tmp %= MO;
            debug {
              if (N <= 2L^^10) {
                writeln("  ", l, ": ", tmp);
              }
            }
            ans += tmp;
            ans %= MO;
          }
        }
      }
      ans = (ans % MO + MO) % MO;
      writeln(ans);
    }
  } catch (EOFException e) {
  }
}
0