#pragma GCC optimize("Ofast") #include using namespace std; #define rep(i, n) for(int i = 0; i < (n); ++i) #define all(x) (x).begin(),(x).end() #define ln '\n' const long long MOD = 1000000007LL; //const long long MOD = 998244353LL; typedef long long ll; typedef unsigned long long ull; typedef pair pii; typedef pair pll; template inline bool chmax(T &a, T b) { if (a < b) { a = b; return true;} return false; } template inline bool chmin(T &a, T b) { if (a > b) { a = b; return true;} return false; } /////////////////////////////////////////////////////////////////////////////////////////////////// //modint::set(M) struct modint { using u64 = std::uint_fast64_t; public: u64 a; modint(u64 x = 0) : a(x >= 0 ? x%Modulus() : (Modulus() - (-x)%Modulus()) % Modulus()) {} static u64 &Modulus() {static u64 Modulus = 0; return Modulus;} static void set (u64 M) {Modulus() = M;} modint operator+(const modint rhs) { return modint(*this) += rhs; } modint operator-(const modint rhs) { return modint(*this) -= rhs; } modint operator*(const modint rhs) { return modint(*this) *= rhs; } modint operator/(const modint rhs) { return modint(*this) /= rhs; } modint &operator+=(const modint rhs) { a += rhs.a; if (a >= Modulus()) { a -= Modulus(); } return *this; } modint &operator-=(const modint rhs) { if (a < rhs.a) { a += Modulus(); } a -= rhs.a; return *this; } modint &operator*=(modint rhs) { a = a * rhs.a % Modulus(); return *this; } modint &operator/=(modint rhs) noexcept { u64 exp = Modulus() - 2; while (exp) { if (exp % 2) { *this *= rhs; } rhs *= rhs; exp /= 2; } return *this; } }; using mint = modint; template struct Matrix { array, 500> A; int N; Matrix() {} Matrix(int N) : N(N) {} inline const array &operator[](int k) const {return A[k];} inline array &operator[](int k) {return A[k];} static Matrix I(int N) { Matrix mat(N); for(int i = 0; i < N; ++i) mat[i][i] = 1; return mat; } Matrix &operator+=(const Matrix &B) { for(int i = 0; i < N; i++) { for(int j = 0; j < N; j++) { (*this)[i][j] += B[i][j]; } } return (*this); } Matrix &operator-=(const Matrix &B) { for(int i = 0; i < N; i++) { for(int j = 0; j < N; j++) { (*this)[i][j] -= B[i][j]; } } return (*this); } Matrix &operator*=(const Matrix &B) { Matrix C(N); for(int i = 0; i < N; i++) { for(int j = 0; j < N; j++) { for(int k = 0; k < N; k++) { C[i][j] = (C[i][j] + (*this)[i][k] * B[k][j]); } } } A.swap(C.A); return *this; } Matrix &operator^=(long long k) { Matrix B = Matrix::I(N); 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 long long k) const { return (Matrix(*this) ^= k); } }; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); ll N,M; cin >> N >> M; modint::set(M); Matrix mat(2); mat[0][0] = mat[0][1] = mat[1][0] = 1; mat ^= N-1; cout << mat[1][0].a << ln; }