結果

問題 No.243 出席番号(2)
ユーザー ぷら
提出日時 2022-08-26 20:26:04
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 113 ms / 2,000 ms
コード長 970 bytes
コンパイル時間 2,207 ms
コンパイル使用メモリ 196,608 KB
最終ジャッジ日時 2025-01-31 04:08:30
ジャッジサーバーID
（参考情報） 		judge2 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 30
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

constexpr int mod = 1000000007;

int dp[2][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[1][j],dp[0][j]);
            mpl(dp[1][j+1],1ll*dp[0][j]*cnt[i]%mod);
        }
        for(int j = 0; j <= N; j++) {
            swap(dp[0][j],dp[1][j]);
            dp[1][j] = 0;
        }
    }
    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[0][i]*f[N-i]%mod);
        }
        else {
            mpl(ans,mod-1ll*dp[0][i]*f[N-i]%mod);
        }
    }
    cout << ans << endl;
}
0