#include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using Int = long long; using Real = long double; using CP = complex; using P = pair; const Int MOD = 1000000007; const Int MOD2 = 998244353; const Int LINF = (1LL << 60); const int INF = (1000000007); const Real EPS = 1e-10; const long double PI = 3.141592653589793238462643383279502884L; using Matrix = vector>; Int mod; Matrix operator*(Matrix &a, Matrix &b){ int n = a.size(); Matrix c(n, vector(n, 0)); for(int i = 0;i < n;i++) for(int j = 0;j < n;j++) for(int k = 0;k < n;k++) (c[i][k] += a[i][j] * b[j][k] % mod) %= mod; return c; } Matrix matPow(Matrix &a, Int n){ if(n == 1)return a; Matrix aa = a * a; Matrix half = matPow(aa, n / 2); if(n % 2 == 0){ return half; } else{ return a * half; } } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); Int N, M, L, K, B; cin >> N >> M >> L >> K >> B; mod = B; Matrix mat(L, vector(L)); for(int i = 0;i < L;i++){ mat[i][i] = 1; if(i + 1 < L) mat[i+1][i] = 1; else mat[0][i] = M % B; } mat = matPow(mat, N); cout << mat[K][0] << endl; return 0; }