提出 #67551603

---

ソースコード 拡げる

#include <bits/stdc++.h>

using namespace std;

// b != c
// a == c when b > a
// b > a is bad
// a divides into b is bad

const long long MOD = 998244353;
const long long inv2 = MOD / 2 + 1;
const long long lim = 2000000;

long long mult(long long a, long long b) {
    return (a % MOD) * (b % MOD) % MOD;
}
long long sub(long long a, long long b) {
    return ((a - (b % MOD)) % MOD + MOD) % MOD;
}

int main() {
    ios::sync_with_stdio(0);
    cin.tie(0);
    long long N;
    cin >> N;
    long long ans = mult(N, N + 1);
    ans = mult(ans, inv2);
    for (int i= 1; i <= lim; i++) {
        ans = sub(ans, N / i);
    }
    for (int k = 1; k <= lim; k++) {
        long long lb = N / (k + 1);
        lb++;
        lb = max(lim + 1, lb);
        long long amt = N / k - lb + 1;
        amt = max(0LL, amt);
        ans = sub(ans, mult(k, amt));
    }
    cout << ans << "\n";
}

提出情報

ジャッジ結果