QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#782007 | #9561. 树数叔术 | Milkcat2009 | 100 ✓ | 120ms | 10984kb | C++14 | 3.3kb | 2024-11-25 18:28:34 | 2024-11-25 18:28:41 |
Judging History
answer
#include <bits/stdc++.h>
#define REP(i, l, r) for (int i = (l); i <= (r); ++ i)
#define DEP(i, r, l) for (int i = (r); i >= (l); -- i)
#define fi first
#define se second
#define pb emplace_back
#define mems(x, v) memset((x), (v), sizeof(x))
#define SZ(x) (int)(x).size()
#define ALL(x) (x).begin(), (x).end()
using namespace std;
namespace math {
typedef long long LL;
template <class T> T qpow(T a, LL b) { if (!b) return {1}; T rs = a; b --; for (; b; b >>= 1, a = a * a) if (b & 1) rs = rs * a; return rs; }
LL mul(LL a, LL b, LL p) { LL rs = a * b - LL(1.L * a * b / p) * p; rs %= p; if (rs < 0) rs += p; return rs; }
template <unsigned P = 0> struct mint {
unsigned v; static unsigned mod; mint() = default;
template <class T> mint(T x) { x %= getmod(), v = x < 0 ? x + getmod() : x; }
static int getmod() { if (P > 0) return P; else return mod; }
static void setmod(unsigned m) { mod = m; }
mint operator + () const { return *this; }
mint operator - () const { return mint(0) - *this; }
mint inv() const { return assert(v), qpow(*this, getmod() - 2); }
int val() const { return v; }
mint &operator += (const mint &q) { if (v += q.v, v >= getmod()) v -= getmod(); return *this; }
mint &operator -= (const mint &q) { if (v -= q.v, v >= getmod()) v += getmod(); return *this; }
mint &operator *= (const mint &q) { v = 1ull * v * q.v % getmod(); return *this; }
mint &operator /= (const mint &q) { return *this *= q.inv(); }
friend mint operator + (mint p, const mint &q) { return p += q; }
friend mint operator - (mint p, const mint &q) { return p -= q; }
friend mint operator * (mint p, const mint &q) { return p *= q; }
friend mint operator / (mint p, const mint &q) { return p /= q; }
friend bool operator == (const mint &p, const mint &q) { return p.v == q.v; }
friend bool operator != (const mint &p, const mint &q) { return p.v != q.v; }
friend bool operator < (const mint &p, const mint &q) { return p.v < q.v; }
friend bool operator > (const mint &p, const mint &q) { return p.v > q.v; }
friend bool operator <= (const mint &p, const mint &q) { return p.v <= q.v; }
friend bool operator >= (const mint &p, const mint &q) { return p.v >= q.v; }
friend istream &operator >> (istream &is, mint &a) { is >> a.v; return is; }
friend ostream &operator << (ostream &os, const mint &a) { os << a.v; return os; }
};
template <> unsigned mint<0>::mod = 998244353;
}
namespace Milkcat {
typedef long long LL;
typedef pair<LL, LL> pii;
const int N = 155;
typedef math::mint<> MI;
int n, m, mod; MI fac[N], C[N][N], f[N][N][N];
int main() {
cin >> n >> m >> mod, MI::setmod(mod);
if (m >= n) { cout << "0\n"; return 0; }
fac[0] = C[0][0] = f[0][1][0] = 1;
REP(i, 1, n) {
C[i][0] = 1, fac[i] = fac[i - 1] * i;
REP(j, 1, i) C[i][j] = C[i - 1][j] + C[i - 1][j - 1];
}
REP(i, 1, m) REP(j, 1, n) REP(k, 0, j) {
if (f[i - 1][j][k] == 0) continue;
REP(x, 0, k - 1) f[i][j][x] += f[i - 1][j][k] * C[k][x];
REP(x, 1, n - j) f[i][j + x][k + x - 1] += f[i - 1][j][k] * j;
}
cout << f[m][n][0] * fac[n] << '\n';
return 0;
}
}
int main() {
// freopen("c.in", "r", stdin);
// freopen("c.out", "w", stdout);
ios::sync_with_stdio(0);
cin.tie(0), cout.tie(0);
int T = 1;
while (T --) Milkcat::main();
return 0;
}
詳細信息
Subtask #1:
score: 5
Accepted
Test #1:
score: 5
Accepted
time: 0ms
memory: 3664kb
input:
1 1 624295285
output:
0
result:
ok single line: '0'
Test #2:
score: 5
Accepted
time: 0ms
memory: 3668kb
input:
4 684813415 564954712
output:
0
result:
ok single line: '0'
Test #3:
score: 5
Accepted
time: 0ms
memory: 3684kb
input:
4 2 844826878
output:
24
result:
ok single line: '24'
Test #4:
score: 5
Accepted
time: 0ms
memory: 3664kb
input:
4 17724073 252218682
output:
0
result:
ok single line: '0'
Test #5:
score: 5
Accepted
time: 0ms
memory: 3600kb
input:
4 3 697681963
output:
384
result:
ok single line: '384'
Subtask #2:
score: 15
Accepted
Test #6:
score: 15
Accepted
time: 0ms
memory: 3612kb
input:
6 4 956647977
output:
238320
result:
ok single line: '238320'
Test #7:
score: 15
Accepted
time: 0ms
memory: 3580kb
input:
6 450615260 491361886
output:
0
result:
ok single line: '0'
Test #8:
score: 15
Accepted
time: 0ms
memory: 3624kb
input:
6 5 339344353
output:
933120
result:
ok single line: '933120'
Test #9:
score: 15
Accepted
time: 0ms
memory: 3684kb
input:
6 3 842228619
output:
23760
result:
ok single line: '23760'
Test #10:
score: 15
Accepted
time: 0ms
memory: 3660kb
input:
6 5 331699652
output:
933120
result:
ok single line: '933120'
Subtask #3:
score: 30
Accepted
Test #11:
score: 30
Accepted
time: 2ms
memory: 4292kb
input:
48 26 424594716
output:
362283012
result:
ok single line: '362283012'
Test #12:
score: 30
Accepted
time: 0ms
memory: 3720kb
input:
49 1000000000 738885247
output:
0
result:
ok single line: '0'
Test #13:
score: 30
Accepted
time: 0ms
memory: 4420kb
input:
48 39 688951620
output:
598399200
result:
ok single line: '598399200'
Test #14:
score: 30
Accepted
time: 0ms
memory: 3592kb
input:
50 476039414 292870080
output:
0
result:
ok single line: '0'
Test #15:
score: 30
Accepted
time: 2ms
memory: 4524kb
input:
50 48 245196368
output:
123576912
result:
ok single line: '123576912'
Subtask #4:
score: 50
Accepted
Test #16:
score: 50
Accepted
time: 0ms
memory: 3660kb
input:
150 526250070 197316869
output:
0
result:
ok single line: '0'
Test #17:
score: 50
Accepted
time: 114ms
memory: 10592kb
input:
149 116 671784452
output:
18945228
result:
ok single line: '18945228'
Test #18:
score: 50
Accepted
time: 104ms
memory: 10672kb
input:
146 144 906402626
output:
438777234
result:
ok single line: '438777234'
Test #19:
score: 50
Accepted
time: 111ms
memory: 10780kb
input:
147 143 705666477
output:
70408701
result:
ok single line: '70408701'
Test #20:
score: 50
Accepted
time: 120ms
memory: 10984kb
input:
150 147 453481175
output:
336290325
result:
ok single line: '336290325'