結果

問題 No.2483 Yet Another Increasing XOR Problem
コンテスト
ユーザー SSRS
提出日時 2023-09-22 23:43:46
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,878 bytes
コンパイル時間 1,815 ms
コンパイル使用メモリ 200,960 KB
実行使用メモリ 7,844 KB
最終ジャッジ日時 2025-06-20 01:52:30
合計ジャッジ時間 2,379 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 24
権限があれば一括ダウンロードができます

ソースコード

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