#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; long long MOD; template struct Matrix { int row, column; vector> val; // val[i][j]: i行j列成分 Matrix(const int &_row, const int &_column) { assert(_row > 0 && _column > 0); row = _row; column = _column; val = vector>(row, vector(column)); } Matrix(const vector> &_val) { row = _val.size(); column = _val[0].size(); assert(row > 0 && column > 0); val = _val; } Matrix operator-() const { Matrix res{row, column}; for (int i = 0; i < res.row; ++i) { for (int j = 0; j < res.column; ++j) { res.val[i][j] = -val[i][j]; } } return res; } }; template Matrix zero(const int &m, const int &n) { Matrix res{m, n}; return res; } template Matrix zero(const int &n) { Matrix res{n, n}; return res; } template Matrix unit(const int &n) { Matrix res{n, n}; for (int k = 0; k < n; ++k) { res.val[k][k] = 1; } return res; } // スカラー倍 template Matrix operator*(const T &k, const Matrix &A) { Matrix res{A.row, A.column}; for (int i = 0; i < res.row; ++i) { for (int j = 0; j < res.column; ++j) { res.val[i][j] = k * A.val[i][j]; } } 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 &left, const Matrix &right) { assert(left.row == right.row && left.column == right.column); Matrix res{left.row, left.column}; for (int i = 0; i < res.row; ++i) { for (int j = 0; j < res.column; ++j) { res.val[i][j] = left.val[i][j] + right.val[i][j]; } } return res; } template Matrix operator+=(Matrix &left, const Matrix &right) { left = left + right; } template Matrix operator-(const Matrix &left, const Matrix &right) { return left + (-right); } template Matrix operator-=(Matrix &left, const Matrix &right) { left = left - right; } // 行列積 template Matrix operator*(const Matrix &left, const Matrix &right) { assert (left.column == right.row); Matrix res{left.row, right.column}; for (int i = 0; i < res.row; ++i) { for (int j = 0; j < res.column; ++j) { for (int k = 0; k < left.column; ++k) { res.val[i][j] += left.val[i][k] * right.val[k][j]; } } } return res; } template Matrix operator*=(Matrix &left, const Matrix &right) { left = left * right; } // 行列累乗 template Matrix pow(Matrix A, long long m) { assert(A.row == A.column); assert(m >= 0); Matrix res = unit(A.row), X = A; while (m > 0) { if (m & 1) { res = res * X; } X = X * X; m >>= 1; } return res; } // 特別な積 template Matrix my_prod(Matrix left, Matrix right) { assert(left.column == right.row); Matrix res{left.row, right.column}; for (int i = 0; i < res.row; ++i) { for (int j = 0; j < res.column; ++j) { for (int k = 0; k < left.column; ++k) { // 演算例 long long temp = left.val[i][k] * right.val[k][j] % MOD; res.val[i][j] = (res.val[i][j] + temp) % MOD; } } } return res; } // 特別な積についての行列累乗 ※単位行列はE? template Matrix my_pow(Matrix A, long long m) { assert(A.row == A.column); assert(m >= 0); Matrix res = unit(A.row), X = A; while (m > 0) { if (m & 1) { res = my_prod(res, X); } X = my_prod(X, X); m >>= 1; } return res; } // debug // void show(Matrix A) { // for (int i = 0; i < A.row; ++i) { // cout << "("; // for (int j = 0; j < A.column; ++j) { // cout << A.val[i][j].val() << ", "; // if (j < A.val[0].size() - 1) { // cout << ", "; // } // } // cout << ")\n"; // } // cout << "\n"; // } template void show(Matrix A) { for (int i = 0; i < A.row; ++i) { cout << "("; for (int j = 0; j < A.column; ++j) { cout << A.val[i][j]; if (j < A.val[0].size() - 1) { cout << ", "; } } cout << ")\n"; } cout << "\n"; } int main() { int N; ll M; cin >> N >> M; MOD = M; Matrix A{{{1, 1}, {1, 0}}}; Matrix B = my_pow(A, N - 1); Matrix F = {{{1}, {0}}}; Matrix X = my_prod(B, F); cout << X.val[1][0] << endl; }