#include using namespace std; #include using namespace atcoder; using mint = modint; //【畳込み】O(n m) /* * a[0..n) と b[0..m) を畳み込んだ数列 c[0..n+m-1) を返す. */ vector naive_convolution(const vector& a, const vector& b) { int n = (int)a.size(), m = (int)b.size(); // c[i] = Σj∈[0..i] a[j] b[i-j] (∀i∈[0..n+m-1)) vector c(n + m - 1); for (int i = 0; i < n + m - 1; i++) { for (int j = max(i - (m - 1), 0); j <= min(i, n - 1); j++) { c[i] += a[j] * b[i - j]; } } return c; } //【テプリッツ行列の累乗】O(n^2 log d) /* * 左下から右上までの成分が順に a(-n..n) であるテプリッツ行列を d 乗したテプリッツ行列を返す. */ vector toeplitz_pow(const vector& a, long long d) { int n = ((int)a.size() + 1) / 2; // テプリッツ行列 a, b の積を返す(制約 : 結果もテプリッツ行列) auto mul = [&](const vector& a, const vector& b) { vector bl(n), br(n); for (int i = 0; i < n; i++) bl[i] = b[i]; for (int i = 0; i < n; i++) br[i] = b[n - 1 + i]; auto resl = naive_convolution(a, bl); auto resr = naive_convolution(a, br); vector res(2 * n - 1); for (int i = 0; i < n; i++) res[i] = resl[n - 1 + i]; for (int i = 1; i < n; i++) res[n - 1 + i] = resr[n - 1 + i]; return res; }; vector res(2 * n - 1), pow2 = a; res[n - 1] = 1; // ダブリングでテプリッツ行列を累乗する. while (d > 0) { if ((d & 1) != 0) { res = mul(res, pow2); } pow2 = mul(pow2, pow2); d /= 2; } return res; } int main() { // input_from_file("input.txt"); // output_to_file("output.txt"); long long n, m; int l, k, b; cin >> n >> m >> l >> k >> b; mint::set_mod(b); // l = 1 のときは二項定理で OK if (l == 1) { cout << mint(1 + m).pow(n).val() << endl; return 0; } // mat : 遷移行列(テプリッツ行列) vector mat(2 * l - 1); mat[0] = m; mat[l - 1] = mat[l] = 1; // n 回遷移する mat = toeplitz_pow(mat, n); cout << mat[l - 1 + k].val() << endl; }