/* -*- coding: utf-8 -*- * * 3364.cc: No.3364 Push_back Operation - yukicoder */ #include #include using namespace std; /* constant */ const int MOD = 998244353; /* typedef */ using ll = long long; template struct MI { int v; MI(): v() {} MI(int _v): v(_v % MOD) { if (v < 0) v += MOD; } MI(long long _v): v(_v % MOD) { if (v < 0) v += MOD; } explicit operator int() const { return v; } MI operator+(const MI m) const { return MI(v + m.v); } MI operator-(const MI m) const { return MI(v + MOD - m.v); } MI operator-() const { return MI(MOD - v); } MI operator*(const MI m) const { return MI((long long)v * m.v); } MI &operator+=(const MI m) { return (*this = *this + m); } MI &operator-=(const MI m) { return (*this = *this - m); } MI &operator*=(const MI m) { return (*this = *this * m); } bool operator==(const MI m) const { return v == m.v; } bool operator!=(const MI m) const { return v != m.v; } MI pow(long long n) const { // a^n % MOD MI pm = 1, a = *this; while (n > 0) { if (n & 1) pm *= a; a *= a; n >>= 1; } return pm; } MI inv() const { return pow(MOD - 2); } MI operator/(const MI m) const { return *this * m.inv(); } MI &operator/=(const MI m) { return (*this = *this / m); } }; using mi = MI; /* global variables */ /* subroutines */ /* main */ int main() { ll n; scanf("%lld", &n); // f(l)=[n/l]^l mi sum = 0; for (ll i = 1;; i++) { ll j = n / i; if (j * j < n) break; sum += mi(j).pow(i); //printf("%lld\n", i); } for (ll j = 1; j * j < n; j++) { ll i0 = n / j, i1 = n / (j + 1); //printf(" %lld-%lld(j=%lld)\n", i1 + 1, i0, j); // j^(i1+1)+j^(i1+2)+..+j^i0 // = sum_{k<=i0} j^k - sum_{k<=i1} j^k // = (j^(i0+1)-1)/(j-1) - (j^(i1+1)-1)/(j-1) if (j == 1) sum += i0 - i1; else { mi e0 = (mi(j).pow(i0 + 1) - 1); mi e1 = (mi(j).pow(i1 + 1) - 1); sum += (e0 - e1) / (j - 1); } } printf("%d\n", (int)sum); return 0; }