#define _USE_MATH_DEFINES #include using namespace std; using ll = long long; using ull = unsigned long long; using pii = pair; using vi = vector; using vvi = vector; using vvvi = vector; using pll = pair; using vl = vector; using vvl = vector; using vvvl = vector; using pdd = pair; using vd = vector; using vvd = vector; using vvvd = vector; using pbb = pair; using vb = vector; using vvb = vector; using vvvb = vector; using pss = pair; using vs = vector; using vvs = vector; using vvvs = vector; #define repeat(counter, start, goal, step) for (int counter = (int)(start); counter < (int)(goal); counter += step) #define rep(counter, times) repeat(counter, 0, times, 1) #define rep2(counter, start, goal) repeat(counter, start, goal, 1) #define rep0(times) repeat(counter, 0, times, 1) #define rrepeat(counter, start, goal, step) for (int counter = (int)(start); counter > (int)(goal); counter -= step) #define rep3(counter, start, goal) rrepeat(counter, start, goal, 1) #define rrep(counter, start, goal) rrepeat(counter, start - 1, goal - 1, 1) #define all(container) begin(container), end(container) template bool chmax(T &x, const T &y) { if (x < y) { x = y; return true; } return false; } template bool chmin(T &x, const T &y) { if (x > y) { x = y; return true; } return false; } void yn(bool flag) { cout << (flag ? "yes" : "no") << "\n"; } void Yn(bool flag) { cout << (flag ? "Yes" : "No") << "\n"; } void YN(bool flag) { cout << (flag ? "YES" : "NO") << "\n"; } void set_prec(const int &digits) { cout << fixed << setprecision(digits); cerr << fixed << setprecision(digits); } template void pr(const T &obj) { cerr << obj; } template void pr(const T &first, const Ts &...rest) { pr(first); pr(", "); pr(rest...); } template void pr(const pair &pair) { pr("("); pr(pair.first); pr(", "); pr(pair.second); pr(")"); } template void pr(const vector &vec) { pr("{"); for (T obj : vec) { pr(obj); pr(", "); } pr("}"); } template void pr(const vector> &vv) { pr("\n"); rep(index, vv.size()) { pr("["); pr(index); pr("]: "); pr(vv[index]); pr("\n"); } } template void pr(const set &vec) { pr("{"); for (T obj : vec) { pr(obj); pr(", "); } pr("}"); } template void pr(const map &map) { pr("{"); for (pair pair : map) { pr("("); pr(pair.first); pr(": "); pr(pair.second); pr("), "); } pr("}"); } #define db(obj) cerr << #obj << ": "; pr(obj); cerr << " " #define dl(obj) db(obj); cerr << "\n"; #define dm(...) cerr << "(" << #__VA_ARGS__ << "): ("; pr(__VA_ARGS__); cerr << ") " #define dml(...) dm(__VA_ARGS__); cerr << "\n" // #include // using namespace atcoder; // const ll MOD = 1000000007; // const ll MOD = 998244353; ll M; template T zero() { return 0; } // 加法の単位元 template T one() { return 1; } // 乗法の単位元(零元) template T sum(T x, T y) { return (x + y) % M; } // 加法 template T neg(T x) { return (-x + M) % M; } // 加法の逆元 template T prod(T x, T y) { return x * y % M; } // 乗法 template // 乗法の逆元 T inv(T x) { assert(x != zero()); return 1 / x; } // 行列 template struct Matrix { int row, col; vector> a; Matrix(const int &_row, const int &_col) { assert(_row > 0); assert(_col > 0); row = _row; col = _col; a = vector>(row, vector(col, zero())); } Matrix(const int &_n) : Matrix(_n, _n) {} Matrix(const vector> &_a) : Matrix(_a.size(), _a[0].size()) { a = _a; } // 加法の逆行列 Matrix operator-() const { Matrix res = this; for (int i = 0; i < res.row; ++i) { for (int j = 0; j < res.col; ++j) { res.a[i][j] = neg(a[i][j]); } } return res; } // debug void show() { for (int i = 0; i < row; ++i) { cout << "("; for (int j = 0; j < col; ++j) { cout << a[i][j]; if (j < col - 1) cout << ", "; } cout << ")\n"; } cout << "\n"; } }; // 零行列(加法の単位行列) template Matrix zero(const int &m, const int &n) { return new Matrix(m, n); } template Matrix zero(const int &n) { return zero(n, n); } // 単位行列 template Matrix one(const int &m, const int &n) { Matrix res = Matrix(m, n); for (int i = 0; i < min(m, n); ++i) { res.a[i][i] = one(); } return res; } template Matrix one(const int &n) { return one(n, n); } // 行列のスカラー倍 template Matrix operator*(const T &k, const Matrix &A) { Matrix res = A; for (int i = 0; i < res.row; ++i) { for (int j = 0; j < res.col; ++j) { res.a[i][j] = prod(res.a[i][j], k); } } return res; } template Matrix operator*(const Matrix &A, const T &k) { return k * A; } template void operator*=(Matrix &A, const T &k) { A = A * k; } // 行列の加法・減法 template Matrix operator+(const Matrix &A, const Matrix &B) { assert(A.row == B.row); assert(A.col == B.col); Matrix res = Matrix(A.row, A.col); for (int i = 0; i < res.row; ++i) { for (int j = 0; j < res.col; ++j) { res.a[i][j] = sum(A.a[i][j], B.a[i][j]); } } return res; } template void operator+=(Matrix &A, const Matrix &B) { A = A + B; } template Matrix operator-(const Matrix &A, const Matrix &B) { return A + (-B); } template void operator-=(Matrix &A, const Matrix &B) { A = A - B; } // 行列の乗法 template Matrix operator*(const Matrix &A, const Matrix &B) { assert(A.col == B.row); Matrix res = Matrix(A.row, B.col); for (int i = 0; i < res.row; ++i) { for (int j = 0; j < res.col; ++j) { for (int k = 0; k < A.col; ++k) { res.a[i][j] = sum(res.a[i][j], prod(A.a[i][k], B.a[k][j])); } } } return res; } template void operator*=(Matrix &A, const Matrix &B) { A = A * B; } // 行列累乗 template Matrix pow(const Matrix &A, long long d) { assert(A.row == A.col); assert(d >= 0); Matrix res = one(A.row), X = A; while (d > 0) { if (d & 1) res *= X; X *= X; d >>= 1; } return res; } // mintのMatrixのdebug // void show(const Matrix &A) { // for (int i = 0; i < A.row; ++i) { // cout << "("; // for (int j = 0; j < A.col; ++j) { // cout << A.a[i][j].val(); // if (j < A.col - 1) cout << ", "; // } // cout << ")\n"; // } // cout << "\n"; // } // debug template void show(const Matrix &A) { for (int i = 0; i < A.row; ++i) { cout << "("; for (int j = 0; j < A.col; ++j) { cout << A.a[i][j]; if (j < A.col - 1) cout << ", "; } cout << ")\n"; } cout << "\n"; } int main() { int N; cin >> N >> M; Matrix A{{{1, 1}, {1, 0}}}, x{{{1}, {0}}}; Matrix ans = pow(A, N - 1) * x; cout << ans.a[1][0] << endl; }