結果
問題 | No.243 出席番号(2) |
ユーザー |
![]() |
提出日時 | 2020-06-18 00:57:10 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 48 ms / 2,000 ms |
コード長 | 1,300 bytes |
コンパイル時間 | 711 ms |
コンパイル使用メモリ | 83,788 KB |
実行使用メモリ | 8,420 KB |
最終ジャッジ日時 | 2024-07-03 12:34:02 |
合計ジャッジ時間 | 2,112 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 30 |
ソースコード
#include <iostream>#include <algorithm>#include <iomanip>#include <vector>#include <queue>#include <set>#include <map>#define N_MAX 200002using namespace std;typedef long long ll;const ll MOD = 1000000007;ll inv[N_MAX],fac[N_MAX],finv[N_MAX];void init(){fac[0]=fac[1]=1;finv[0]=finv[1]=1;inv[1]=1;for(int i=2;i<N_MAX;i++){inv[i]=MOD-inv[MOD%i]*(MOD/i)%MOD;fac[i]=fac[i-1]*(ll) i%MOD;finv[i]=finv[i-1]*inv[i]%MOD;}}int N;int A[5000];int cnt[5000];ll dp[2][5005];int main(){ios::sync_with_stdio(false);cin.tie(0);cout << setprecision(10) << fixed;cin >> N;init();for(int i = 0; i < N; i++) {cin >> A[i];cnt[A[i]]++;}dp[0][0] = 1;for(int i = 0; i < N; i++){for(int j = 0; j <= i; j++) dp[(i+1)%2][j] = 0;for(int j = 0; j <= i; j++){dp[(i+1)%2][j] += dp[i%2][j];dp[(i+1)%2][j] %= MOD;dp[(i+1)%2][j+1] += dp[i%2][j]*cnt[i];dp[(i+1)%2][j+1] %= MOD;}}ll ans = 0;for(int i = 0; i <= N; i++){ll c = dp[N%2][i];ll tmp = (fac[N-i]*c)%MOD;if(i%2 == 0) ans += tmp;else ans -= tmp;ans %= MOD;}cout << (ans+MOD)%MOD << endl;}