結果

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

ソースコード

diff #
プレゼンテーションモードにする

//Let's join Kaede Takagaki Fan Club !!
#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
using namespace std;
typedef long long ll;
typedef pair<int,int> P;
typedef pair<int,P> P1;
typedef pair<P,P> P2;
#define pu push
#define pb push_back
#define mp make_pair
#define eps 1e-7
#define INF 1000000000
#define fi first
#define sc second
#define rep(i,x) for(int i=0;i<x;i++)
#define repn(i,x) for(int i=1;i<=x;i++)
#define SORT(x) sort(x.begin(),x.end())
#define ERASE(x) x.erase(unique(x.begin(),x.end()),x.end())
#define POSL(x,v) (lower_bound(x.begin(),x.end(),v)-x.begin())
#define POSU(x,v) (upper_bound(x.begin(),x.end(),v)-x.begin())
#define all(x) x.begin(),x.end()
template<class T>
void dmp(T a){
    rep(i,a.size()) cout << a[i] << " ";
    cout << endl;
}
template<class T>
bool chmax(T&a, T b){
    if(a < b){
        a = b;
        return 1;
    }
    return 0;
}
template<class T>
bool chmin(T&a, T b){
    if(a > b){
        a = b;
        return 1;
    }
    return 0;
}
template<class T>
void g(T &a){
    cin >> a;
}
template<class T>
void o(const T &a,bool space=false){
    cout << a << (space?' ':'\n');
}
//ios::sync_with_stdio(false);
const ll mod = 1000000007;//998244353
template<class T>
void add(T&a,T b){
    a+=b;
    if(a >= mod) a-=mod;
}
int n, a[5005];
ll dp[2][5005], w[5005], F[5005];
int main(){
    F[0] = 1;
    repn(i, 5004) F[i] = F[i-1] * i % mod;
    cin >> n;
    repn(i, n) {
        cin >> a[i];
        w[a[i]]++;
    }
    dp[0][0] = 1;
    int cur = 0, nxt = 1;
    rep(j, n){
        memset(dp[nxt], 0, sizeof(dp[nxt]));
        rep(i, j+1){
            {
                add(dp[nxt][i], dp[cur][i]);
                add(dp[nxt][i+1], dp[cur][i] * w[j] % mod);
            }
        }
        swap(cur, nxt);
    }
    ll ans = 0;
    rep(j, n+1) {
        if(j%2 == 0) ans += dp[cur][j] * F[n-j] % mod;
        else ans -= dp[cur][j] * F[n-j] % mod;
    }
    printf("%lld\n", (ans%mod+mod)%mod);
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0