結果

問題 No.243 出席番号(2)
ユーザー machy
提出日時 2015-07-10 23:48:36
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 1,675 bytes
コンパイル時間 572 ms
コンパイル使用メモリ 68,828 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-08 02:22:03
合計ジャッジ時間 3,341 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26 WA * 4
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <string>
#include <algorithm>

using namespace std;

typedef long long LL;
const LL MOD = 1000000007;

int main(){
    int N;
    cin >> N;
    vector<int> deny(N);
    for(int i = 0; i < N; i++){
        cin >> deny[i];
    }
    sort(deny.begin(), deny.end());
    vector<LL> cur((N+1)*2);
    vector<LL> next((N+1)*2);
    int prev_deny = -1;
    cur[0] = 1;
    for(int i = 1; i <= N; i++){
        for(int u = 0; u < 2; u++){
            for(int j = 0; j <= N; j++){
                int u2 = u;
                if(prev_deny != deny[i-1]) u2 = 0;
                next[j*2 + u2] += cur[j*2 + u];
                next[j*2 + u2] %= MOD;
            }
            if(u == 0 || prev_deny != deny[i-1]){
                for(int j = 0; j+1 <= N; j++){
                    next[(j+1)*2 + 1] += cur[j*2 + u];
                    next[(j+1)*2 + 1] %= MOD;
                }
            }
        }
        prev_deny = deny[i-1];
        next.swap(cur);
        for(int j = 0; j < (N+1)*2; j++){
            next[j] = 0;
        }
    }
    vector<LL> pat(N+1);
    for(int i = 0; i <= N; i++){
        pat[i] += cur[i*2 + 0] + cur[i*2 + 1];
        pat[i] %= MOD;
    }

    LL ans = 1;
    for(int i = 0; i < N; i++){
        ans *= (N-i);
        ans %= MOD;
    }

    for(int cnt = 1; cnt <= N; cnt++){
        LL tmp = pat[cnt];
        for(int i = cnt; i < N; i++){
            tmp *= (N-i);
            tmp %= MOD;
        }
        if(cnt % 2 == 1){
            ans += MOD - tmp;
            ans %= MOD;
        }else{
            ans += tmp;
            ans %= MOD;
        }
    }

    cout << ans << endl;

    return 0;
}
0