#include #include using namespace std; using ll = long long; using pii = pair; using pll = pair; using vi = vector; using vl = vector; #define rep3(i, a, b, c) for (ll i = (a); i < (b); i += (c)) #define rep2(i, a, b) rep3(i, a, b, 1) #define rep1(i, n) rep2(i, 0, n) #define rep0(n) rep1(aaaaa, n) #define ov4(a, b, c, d, name, ...) name #define rep(...) ov4(__VA_ARGS__, rep3, rep2, rep1, rep0)(__VA_ARGS__) #define per(i, a, b) for (ll i = (a) - 1; i >= (b); i--) #define fore(e, v) for (auto &&e : v) #define all(a) begin(a), end(a) #define sz(a) (int)(size(a)) #define lb(v, x) (lower_bound(all(v), x) - begin(v)) #define eb emplace_back template bool chmin(T &a, const S &b) { return a > b ? a = b, 1 : 0; } template bool chmax(T &a, const S &b) { return a < b ? a = b, 1 : 0; } const int INF = 1e9 + 100; const ll INFL = 3e18 + 100; #define i128 __int128_t struct _ { _() { cin.tie(0)->sync_with_stdio(0), cout.tie(0); } } __; using mint = atcoder::modint998244353; int main() { ll N; cin >> N; const int D = 61; mint ans = 0; rep(b, D) { vector dp0(D, mint(0)), dp1(D, mint(0)); dp0[0] = 1; per(i, D, 0) { vector ndp0(D, mint(0)), ndp1(D, mint(0)); if (i == b) { if (N >> i & 1) { rep(p, D - 1) ndp0[p + 1] += dp0[p], ndp1[p + 1] += dp1[p]; } else { rep(p, D - 1) ndp1[p + 1] += dp1[p]; } } else { if (N >> i & 1) { rep(p, D - 1) ndp0[p + 1] += dp0[p], ndp1[p + 1] += dp1[p]; rep(p, D) ndp1[p] += dp0[p] + dp1[p]; } else { rep(p, D - 1) ndp1[p + 1] += dp1[p]; rep(p, D) ndp0[p] += dp0[p], ndp1[p] += dp1[p]; } } dp0 = ndp0, dp1 = ndp1; } mint c = mint(2).pow(b); rep(i, D) { mint v = dp0[i] + dp1[i]; ans += c * v * (v + 1) / 2; } } cout << ans.val() << '\n'; }