#include #define FOR(i, a, n) for(ll i = (ll)a; i < (ll)n; i++) #define rep(i, n) FOR(i, 0, n) #define ALL(x) (x).begin(), (x).end() using namespace std; typedef long long ll; constexpr ll mod = 1e9 + 7; template inline bool chmax(T &a, const T b) { if(a >= b) return false; a = b; return true; } template inline bool chmin(T &a, const T b) { if(a <= b) return false; a = b; return true; } /*-------------------------------------------*/ ll N, M; template struct Matrix { vector> A; Matrix() {} Matrix(size_t n, size_t m) : A(n, vector(m, 0)) {} Matrix(size_t n) : A(n, vector(n, 0)) {} size_t height() const { return A.size(); } size_t width() const { return A[0].size(); } inline const vector &operator[](int k) const { return A.at(k); } inline vector &operator[](int k) { return A.at(k); } static Matrix I(size_t n) { Matrix mat(n); for(int i = 0; i < n; i++) mat[i][i] = 1; return mat; } Matrix &operator+=(const Matrix &B) { size_t n = height(), m = width(); assert(n == B.height() && m == B.width()); for(int i = 0; i < n; i++) for(int j = 0; j < m; j++) (*this)[i][j] += B[i][j]; return (*this); } Matrix &operator-=(const Matrix &B) { size_t n = height(), m = width(); assert(n == B.height() && m == B.width()); for(int i = 0; i < n; i++) for(int j = 0; j < m; j++) (*this)[i][j] += M - B[i][j]; return (*this); } Matrix &operator*=(const Matrix &B) { size_t n = height(), m = B.width(), p = width(); assert(p == B.height()); vector> C(n, vector(m, 0)); for(int i = 0; i < n; i++) for(int j = 0; j < m; j++) for(int k = 0; k < p; k++) (C[i][j] += (*this)[i][k] * B[k][j]) %= M; A.swap(C); return (*this); } Matrix &operator^=(int64_t k) { Matrix B = Matrix::I(height()); while(k > 0) { if(k & 1) B *= (*this); (*this) *= (*this); k >>= 1LL; } A.swap(B.A); return (*this); } Matrix &operator+(const Matrix &B) const { return Matrix(*this) += B; } Matrix &operator-(const Matrix &B) const { return Matrix(*this) -= B; } Matrix &operator*(const Matrix &B) const { return Matrix(*this) *= B; } Matrix &operator^(const int64_t k) const { return Matrix(*this) ^= k; } }; int main() { cin.tie(0); ios::sync_with_stdio(0); Matrix A(2, 1), B(2); A.A = {{0}, {1}}; B.A = {{0, 1}, {1, 1}}; cin >> N >> M; B ^= (N - 1); B *= A; cout << B[0][0] << endl; return 0; }