結果

問題 No.737 PopCount
コンテスト
ユーザー tsutaj
提出日時 2018-09-29 18:46:00
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3 ms / 1,000 ms
コード長 997 bytes
コンパイル時間 668 ms
コンパイル使用メモリ 69,248 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-10-12 08:19:21
合計ジャッジ時間 1,369 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 15
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <algorithm>
#include <string>
using namespace std;

// 場合の数、総和
long long int dp[65][65][2][2];
const long long int MOD = 1000000007LL;

int main() {
    long long int N, ans = 0; cin >> N;
    string s = "";
    for(int i=0; i<60; i++, N>>=1) {
        s += (char)('0' + (N & 1));
    }
    reverse(s.begin(), s.end());

    fill(dp[0][0][0], dp[60][0][0], 0);
    dp[0][0][0][0] = 1;

    for(int i=0; i<60; i++) {
        for(int j=0; j<=i; j++) {
            for(int f=0; f<2; f++) {
                int lim = f ? 1 : (s[i] - '0');
                for(int x=0; x<=lim; x++) {
                    (dp[i+1][j+x][f || x < lim][0] += dp[i][j][f][0]) %= MOD;
                    (dp[i+1][j+x][f || x < lim][1] += dp[i][j][f][1] * 2 + dp[i][j][f][0] * x) %= MOD;
                }
            }
        }
    }
    
    for(int i=0; i<60; i++) {
        (ans += i * (dp[60][i][0][1] + dp[60][i][1][1])) %= MOD;
    }
    cout << ans << endl;
    return 0;
}
0