#pragma GCC target("avx2") #pragma GCC optimize("unroll-loops") #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define ALL(x) (x).begin(), (x).end() #define PC(x) __builtin_popcount(x) #define PCL(x) __builtin_popcountll(x) using namespace std; typedef long long ll; typedef unsigned long long ull; typedef pair pii; typedef pair pll; template bool chmax(T& a, const T& b) { if (a < b) { a = b; return true; } return false; } template bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } return false; } const double PI = 3.14159265358979323846; const double PI2 = PI * 2.0; const double EPS = 1E-09; const ll MOD = 998244353; const ll INFL = 1E18; const int INFI = 1E09; const int MAX_N = 2E+05; struct edge { int to, cost, id; }; ll di[4] = { 0, -1, 0, 1 }, dj[4] = { 1, 0, -1, 0 }; ll N; template struct mat { // 行列m vector> m; // コンストラクタ:第1引数:行数、第2引数:列数、第3引数:初期値 mat() : m(vector>(static_cast(0))) { } mat(int h, int w) : m(vector>(h, vector(w, static_cast(0)))) { } mat(int h, int w, T d) : m(vector>(h, vector(w, d))) { } //単位行列を返すコンストラクタ:第1引数:n*n正方行列のn mat(int n) : m(vector>(n, vector(n))) { for (int i = 0; i < n; i++) { m[i][i] = static_cast(1); } } // 添字演算子 vector operator[](const int i) const { return m[i]; } //読み取り vector& operator[](const int i) { return m[i]; } //書き込み // 行数・列数 int nrow = m.size(); //行数 int ncol = m[0].size(); //列数 // 行列同士の演算 mat& operator=(const mat& a) { m.resize(a.nrow); for (int i = 0; i < a.nrow; i++) { m[i].resize(a.ncol); } for (int i = 0; i < a.nrow; i++) { for (int j = 0; j < a.ncol; j++) { m[i][j] = a[i][j]; } } nrow = a.nrow; ncol = a.ncol; return *this; } mat& operator+=(const mat& a) { for (int i = 0; i < nrow; i++) { for (int j = 0; j < ncol; j++) { m[i][j] += a[i][j]; } } return *this; } mat& operator-=(const mat& a) { for (int i = 0; i < nrow; i++) { for (int j = 0; j < ncol; j++) { m[i][j] -= a[i][j]; } } return *this; } mat& operator*=(const mat& a) { mat m2(nrow, a.ncol, static_cast(0)); for (int i = 0; i < nrow; i++) { for (int j = 0; j < a.ncol; j++) { for (int k = 0; k < nrow; k++) { m2[i][j] += m[i][k] * a[k][j] % MOD; } } } for (int i = 0; i < nrow; i++) { for (int j = 0; j < a.ncol; j++) { m[i][j] = m2[i][j]; } } ncol = a.ncol; return *this; } mat operator+(const mat& a) const { return static_cast(*this) += a; } mat operator-(const mat& a) const { return static_cast(*this) -= a; } mat operator*(const mat& a) const { return static_cast(*this) *= a; } bool operator==(const mat& a) { for (int i = 0; i < nrow; i++) { for (int j = 0; j < ncol; j++) { if (m[i][j] != a[i][j]) return false; } } return true; } // 行列とスカラの演算 mat& operator+=(const T& a) { for (int i = 0; i < nrow; i++) { for (int j = 0; j < ncol; j++) { m[i][j] += a; } } return *this; } mat& operator-=(const T& a) { for (int i = 0; i < nrow; i++) { for (int j = 0; j < ncol; j++) { m[i][j] -= a; } } return *this; } mat& operator*=(const T& a) { for (int i = 0; i < nrow; i++) { for (int j = 0; j < ncol; j++) { m[i][j] *= a; } } return *this; } mat& operator/=(const T& a) { for (int i = 0; i < nrow; i++) { for (int j = 0; j < ncol; j++) { m[i][j] /= a; } } return *this; } mat operator+(const T& a) const { return static_cast(*this) += a; } mat operator-(const T& a) const { return static_cast(*this) -= a; } mat operator*(const T& a) const { return static_cast(*this) *= a; } mat operator/(const T& a) const { return static_cast(*this) /= a; } // 標準出力 void show() { for (int i = 0; i < nrow; i++) { for (int j = 0; j < ncol; j++) { cout << m[i][j] << " "; } cout << "\n"; } return; } }; template mat powmat(mat x, ll a) { mat ans(x.nrow); while (a > 0) { if (a & 1) { ans = ans * x; } x = x * x; a >>= 1; } return ans; } int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cin >> N; mat dp(2, 1, 0); dp[0][0] = 1; mat m(2, 2, 1); m[1][1] = 0; m = powmat(m, N - 1); dp = m * dp; // for (int i = 0; i < N; i++) { // for (int j = 0; j < N; j++) { // cout << "i = " << i << ", j = " << j << ", dp = " << dp[i][j] << "\n"; // } // } cout << (dp[0][0] + dp[1][0] - 1 + MOD) % MOD << "\n"; return 0; }