結果

問題 No.2483 Yet Another Increasing XOR Problem
ユーザー SSRSSSRS
提出日時 2023-09-22 23:43:46
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 3 ms / 2,000 ms
コード長 1,878 bytes
コンパイル時間 2,228 ms
コンパイル使用メモリ 206,456 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-16 00:19:32
合計ジャッジ時間 3,312 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

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

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
const long long MOD = 998244353;
const long long HALF = 499122177;
int main(){
  long long M;
  cin >> M;
  if (M == 1){
    cout << 1 << endl;
  } else {
    int p = 64 - __builtin_clzll(M - 1);
    vector<long long> POW(p + 1);
    vector<long long> POW_HALF(p + 1);
    POW[0] = 2;
    POW_HALF[0] = HALF;
    for (int i = 0; i < p; i++){
      POW[i + 1] = POW[i] * POW[i] % MOD;
      POW_HALF[i + 1] = POW_HALF[i] * POW_HALF[i] % MOD;
    }
    long long ans = POW[p];
    long long tmp = 1;
    for (int i = 0; i < p; i++){
      if (((((long long) 1 << p) - M) >> i & 1) == 1){
        tmp *= POW[i];
        tmp %= MOD;
      }
    }
    ans += MOD - tmp;
    vector<vector<long long>> dp(p + 1, vector<long long>(4, 0));
    dp[p][0] = POW[p] * HALF % MOD;
    for (int i = p - 1; i >= 0; i--){
      for (int j = 0; j < 4; j++){
        for (int c = 0; c < 2; c++){
          for (int d = 0; d < 2; d++){
            bool ok = true;
            int j2 = j;
            if ((j & 1) == 0 && c == 1 && d == 0){
              ok = false;
            }
            if (c == 0 && d == 1){
              j2 |= 1;
            }
            if ((j & 2) == 0 && ((M - 1) >> i & 1) == 1 && c == d){
              ok = false;
            }
            if (((M - 1) >> i & 1) == 0 && c != d){
              j2 |= 2;
            }
            if (ok){
              if (c == d){
                dp[i][j2] += dp[i + 1][j];
              }
              if (c == 1 && d == 0){
                dp[i][j2] += dp[i + 1][j] * POW[i];
              }
              if (c == 0 && d == 1){
                dp[i][j2] += dp[i + 1][j] * POW_HALF[i];
              }
              dp[i][j2] %= MOD;
            }
          }
        }
      }
    }
    ans += MOD - dp[0][2];
    ans += MOD - dp[0][3];
    ans %= MOD;
    cout << ans << endl;
  }
}
0