結果
| 問題 |
No.243 出席番号(2)
|
| ユーザー |
ぷら
|
| 提出日時 | 2022-08-26 20:25:16 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
MLE
|
| 実行時間 | - |
| コード長 | 866 bytes |
| コンパイル時間 | 1,977 ms |
| コンパイル使用メモリ | 196,336 KB |
| 最終ジャッジ日時 | 2025-01-31 04:08:21 |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 20 MLE * 10 |
ソースコード
#include <bits/stdc++.h>
using namespace std;
constexpr int mod = 1000000007;
int dp[5050][5050];
void mpl(int &x,int y) {
x += y;
if(x >= mod) x -= mod;
}
int main() {
int N;
cin >> N;
vector<int>cnt(5000);
for(int i = 0; i < N; i++) {
int a;
cin >> a;
cnt[a]++;
}
dp[0][0] = 1;
for(int i = 0; i < N; i++) {
for(int j = 0; j < N; j++) {
mpl(dp[i+1][j],dp[i][j]);
mpl(dp[i+1][j+1],1ll*dp[i][j]*cnt[i]%mod);
}
}
vector<int>f(N+1);
f[0] = 1;
for(int i = 1; i <= N; i++) {
f[i] = 1ll*f[i-1]*i%mod;
}
int ans = 0;
for(int i = 0; i <= N; i++) {
if(i%2 == 0) {
mpl(ans,1ll*dp[N][i]*f[N-i]%mod);
}
else {
mpl(ans,mod-1ll*dp[N][i]*f[N-i]%mod);
}
}
cout << ans << endl;
}
ぷら