結果
| 問題 | No.737 PopCount | 
| コンテスト | |
| ユーザー |  | 
| 提出日時 | 2018-09-28 23:15:17 | 
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) | 
| 結果 | 
                                AC
                                 
                             | 
| 実行時間 | 19 ms / 1,000 ms | 
| コード長 | 1,051 bytes | 
| コンパイル時間 | 1,677 ms | 
| コンパイル使用メモリ | 193,696 KB | 
| 最終ジャッジ日時 | 2025-01-06 14:00:08 | 
| ジャッジサーバーID (参考情報) | judge4 / judge2 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 3 | 
| other | AC * 15 | 
ソースコード
//=================================
// Created on: 2018/09/28 22:51:11
//=================================
#include <bits/stdc++.h>
using ll = long long;
constexpr ll MOD = 1000000007LL;
constexpr ll HALF = (MOD + 1) / 2;
int main()
{
    ll N;
    std::cin >> N;
    ll d = 1;
    ll ans = 0;
    for (; d <= N; d *= 2) {
        const ll step = d * 2;
        const ll num = (N + 1) / step;
        ll sum = 0;
        if (num >= 1) {
            const ll STEP = (d % MOD) * (step % MOD) % MOD;
            const ll S = (((3 * d - 1) % MOD) * (d % MOD) % MOD) * HALF % MOD;
            const ll T = (S + ((num - 1) % MOD) * STEP % MOD) % MOD;
            sum = (((S + T) % MOD) * (num % MOD) % MOD) * HALF % MOD;
            const ll A = step * num + d;
            if (A <= N) { (sum += (((A + N) % MOD) * ((N - A + 1) % MOD) % MOD) * HALF % MOD) %= MOD; }
        } else {
            sum = (((N + d) % MOD) * ((N - d + 1) % MOD) % MOD) * HALF % MOD;
        }
        (ans += sum) %= MOD;
    }
    std::cout << ans << std::endl;
    return 0;
}
            
            
            
        