結果
問題 | No.243 出席番号(2) |
ユーザー |
![]() |
提出日時 | 2018-08-26 17:34:37 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 21 ms / 2,000 ms |
コード長 | 1,444 bytes |
コンパイル時間 | 1,725 ms |
コンパイル使用メモリ | 178,008 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-01 03:57:07 |
合計ジャッジ時間 | 2,798 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include <bits/stdc++.h>using namespace std;using lint = long long int;struct fast_ios { fast_ios(){ cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(20); }; } fast_ios_;#define FOR(i, begin, end) for(int i=(begin);i<(end);i++)#define IFOR(i, begin, end) for(int i=(end)-1;i>=(begin);i--)#define REP(i, n) FOR(i,0,n)#define IREP(i, n) IFOR(i,0,n)int N;constexpr int Nmax = 5000;map<int, int> A;constexpr lint MOD = 1000000007;vector<lint> fac, facInv, inv;void facInit(int nmax){fac = facInv = inv = vector<lint>(nmax + 1, 1);for (int i = 2; i <= nmax; i++){fac[i] = fac[i-1] * i % MOD;inv[i] = MOD - inv[MOD%i] * (MOD/i) % MOD;facInv[i] = facInv[i-1] * inv[i] % MOD;}}lint nCr(int n, int r){if (n<r || r<0) return 0;return (fac[n] * facInv[r] % MOD) * facInv[n-r] % MOD;}signed main(){cin >> N;REP(i, N) {int tmp;cin >> tmp;if (tmp < N) A[tmp]++;}lint dp[Nmax + 1];memset(dp, 0, sizeof(dp));dp[0] = 1;for (auto v : A){lint n = v.second;IREP(i, Nmax) if (dp[i]) (dp[i + 1] += n * dp[i]) %= MOD;}lint ans = 0;facInit(Nmax);REP(i, N + 1) // i人が悪い番号になる{lint tmp = dp[i] * fac[N - i] % MOD;if (i & 1) (ans += (MOD - tmp)) %= MOD;else (ans += tmp) %= MOD;}cout << ans << endl;}