結果

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

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
ll const m = 998244353;
ll mpow(ll a, ll n) {
	n %= (m - 1);
	ll ret = 1;
	while (n) {
		if (n & 1) ret = (ret * a) % m;
		a = (a * a) % m;
		n >>= 1;
	}
	return ret;
}
int main() {
	ll N;
	cin >> N;
	if (N < 3) {
		if (N == 1) {
			puts("1");
		}
		if (N == 2) {
			puts("3");
		}
		return 0;
	}
	N --;
	ll n = 1; int lg = -1;
	while (n <= N) {
		n <<= 1;
		lg ++;
	}
	n >>= 1;
	ll ans = mpow(2, n) - 1;
	ans += (mpow(2, n) + m - mpow(2, n - (N ^ n) - 1)) % m;
	ans %= m;
	// printf("n : %lld, ans : %lld\n", n, ans);
	ll dp[60][3];
	if (N & 1) {
		dp[0][1] = 2 + (m + 1) / 2; dp[0][0] = 2; dp[0][2] = 0;
	} else {
		dp[0][0] = 0; dp[0][1] = 2;dp[0][2] = 2 + (m + 1) / 2;
	}
	for (int i = 1; i < lg; i ++) {
		ll* pre = dp[i - 1];
		ll* now = dp[i];
		ll b = mpow(2, (1ll << i));
		ll iv = mpow(b, m - 2);
		if ((N >> i) & 1) {
			now[0] = (pre[0] + pre[1] + pre[2]) * 2; now[0] %= m;
			now[1] = 0; now[2] = 0;
			for (int j = 0; j < 3; j ++) {
				now[j] += pre[j] * (b + iv) % m;
				now[j] %= m;
			}
		} else {
			for (int j = 0; j < 3; j ++) {
				now[j] = pre[j] * 2;
				now[j] %= m;
			}
			now[2] += (pre[0] + pre[1] + pre[2]) * ((b + iv) % m);
			now[2] %= m;
		}
	}
	ll kj = dp[lg - 1][0] + dp[lg - 1][1];
	kj %= m;
	kj = (kj * mpow(2, (1ll << lg) - 1)) % m;
	cout << (ans + kj) % m << endl;
}
0