結果

問題 No.737 PopCount
ユーザー 0w10w1
提出日時 2018-10-09 14:40:01
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 3 ms / 1,000 ms
コード長 1,210 bytes
コンパイル時間 2,216 ms
コンパイル使用メモリ 203,836 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-04-20 19:06:55
合計ジャッジ時間 2,789 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 1 ms
5,376 KB
testcase_03 AC 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 3 ms
5,376 KB
testcase_07 AC 3 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 3 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 3 ms
5,376 KB
testcase_13 AC 3 ms
5,376 KB
testcase_14 AC 3 ms
5,376 KB
testcase_15 AC 2 ms
5,376 KB
testcase_16 AC 1 ms
5,376 KB
testcase_17 AC 2 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

const int M = int(1e9 + 7);

int invmod(int v) {
  int r = 1;
  for (int n = M - 2; n; n /= 2) {
    if (n & 1) r = 1LL * r * v % M;
    v = 1LL * v * v % M;
  }
  return r;
}

signed main() {
  ios::sync_with_stdio(false);

  int64_t N;
  cin >> N;

  vector<int> s;
  for (int64_t n = N; n; n /= 2) s.emplace_back(n & 1);
  reverse(s.begin(), s.end());

  vector<vector<vector<int>>> dp(s.size() + 1, vector<vector<int>>(2, vector<int>(s.size() + 1)));
  dp[0][0][0] = 1;
  vector<vector<vector<int>>> dpdp(s.size() + 1, vector<vector<int>>(2, vector<int>(s.size() + 1)));
  for (int i = 0; i < s.size(); ++i) {
    for (int j = 0; j < 2; ++j) {
      for (int k = 0; k <= i; ++k) {
        for (int d = 0; d <= (j ? 1 : s[i]); ++d) {
          (dp[i + 1][j | d < s[i]][k + d] += dp[i][j][k]) %= M;
          (dpdp[i + 1][j | d < s[i]][k + d] += (1LL * dpdp[i][j][k] * (k + d) % M * invmod(k) % M + (1LL * d << s.size() - 1 - i) % M * (k + d) % M * dp[i][j][k] % M) % M) %= M;
        }
      }
    }
  }

  int ans = 0;
  for (int i = 1; i <= s.size(); ++i)
    for (int j = 0; j < 2; ++j)
      (ans += dpdp[s.size()][j][i]) %= M;
  
  cout << ans << endl;
}
0