結果

問題 No.243 出席番号(2)
ユーザー motimoti
提出日時 2015-12-20 02:53:37
言語 C++11
(gcc 11.4.0)
結果
RE  
実行時間 -
コード長 3,734 bytes
コンパイル時間 904 ms
コンパイル使用メモリ 104,676 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-09-17 11:18:11
合計ジャッジ時間 4,817 ms
ジャッジサーバーID
（参考情報） 		judge6 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,816 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 9 ms
6,944 KB
testcase_04 AC 9 ms
6,940 KB
testcase_05 AC 9 ms
6,944 KB
testcase_06 AC 10 ms
6,940 KB
testcase_07 RE -
testcase_08 AC 32 ms
6,940 KB
testcase_09 AC 32 ms
6,940 KB
testcase_10 AC 31 ms
6,944 KB
testcase_11 AC 32 ms
6,940 KB
testcase_12 RE -
testcase_13 AC 69 ms
6,940 KB
testcase_14 AC 69 ms
6,940 KB
testcase_15 AC 69 ms
6,944 KB
testcase_16 AC 69 ms
6,940 KB
testcase_17 RE -
testcase_18 AC 120 ms
6,940 KB
testcase_19 AC 120 ms
6,940 KB
testcase_20 AC 122 ms
6,940 KB
testcase_21 AC 121 ms
6,940 KB
testcase_22 RE -
testcase_23 AC 188 ms
6,940 KB
testcase_24 AC 188 ms
6,944 KB
testcase_25 AC 186 ms
6,944 KB
testcase_26 AC 187 ms
6,940 KB
testcase_27 AC 187 ms
6,940 KB
testcase_28 AC 2 ms
6,940 KB
testcase_29 AC 1 ms
6,944 KB
testcase_30 AC 2 ms
6,940 KB
testcase_31 AC 2 ms
6,940 KB
testcase_32 AC 1 ms
6,940 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <cmath>
#include <vector>
#include <complex>
#include <queue>
#include <deque>
#include <set>
#include <map>
#include <unordered_set>
#include <unordered_map>
#include <iomanip>
#include <assert.h>
#include <array>
#include <cstdio>
#include <cstring>
#include <random>
#include <functional>
#include <numeric>
#include <bitset>

using namespace std;

#define REP(i,a,b) for(int i=a;i<(int)b;i++)
#define rep(i,n) REP(i,0,n)
#define all(c) (c).begin(), (c).end()
#define zero(a) memset(a, 0, sizeof a)
#define minus(a) memset(a, -1, sizeof a)
#define minimize(a, x) a = std::min(a, x)
#define maximize(a, x) a = std::max(a, x)

typedef long long ll;
int const inf = 1<<29;

template<class value_type, int MOD>
class ModInt {
private:
  value_type val_;
  static value_type mod_pow(value_type x, value_type n, value_type mo) { value_type ret = 1; while(n > 0) { if(n & 1) { ret = ret * x % mo; } x = x * x % mo; n >>= 1; } return ret; }

public:

  ModInt() { val_ = 0; }
  ModInt(value_type x) { val_ = (x % MOD + MOD) % MOD; }

  ModInt const operator + (ModInt const& rhs) const {
    return std::move(ModInt(val_+rhs.get()));
  }

  ModInt const operator - (ModInt const& rhs) const {
    return std::move(ModInt(val_-rhs.get()));
  }

  ModInt const operator * (ModInt const& rhs) const {
    return std::move(ModInt(val_*rhs.get()));
  }

  ModInt const operator / (ModInt const& rhs) const {
    return std::move(ModInt(val_*mod_pow(rhs.get(), MOD-2, MOD)));  // fermat theorem
  }

  friend ModInt const operator + (value_type lhs, ModInt const& rhs) {
    return std::move(ModInt(lhs+rhs.get()));
  }

  friend ModInt const operator - (value_type lhs, ModInt const& rhs) {
    return std::move(ModInt(lhs-rhs.get()));
  }

  friend ModInt const operator * (value_type lhs, ModInt const& rhs) {
    return std::move(ModInt(lhs*rhs.get()));
  }

  friend ModInt const operator / (value_type lhs, ModInt const& rhs) {
    return std::move(ModInt(lhs*mod_pow(rhs.get(), MOD-2, MOD)));  // fermat theorem
  }

  ModInt operator += (ModInt const& rhs) {
    return *this = ModInt(val_+rhs.get()).get();
  }

  ModInt operator -= (ModInt const& rhs) {
    return *this = ModInt(val_-rhs.get()).get();
  }

  ModInt operator *= (ModInt const& rhs) {
    return *this = ModInt(val_*rhs.get()).get();
  }

  ModInt operator /= (ModInt const& rhs) {
    return *this = ModInt(val_*mod_pow(rhs.get(), MOD-2, MOD)).get();
  }

  bool operator == (ModInt const& rhs) const {
    return val_ == rhs.get();
  }

  value_type const get() const { return val_; }
  value_type &     get() { return val_; }

  friend ostream& operator << (ostream& ost, ModInt const& x) { 
    return ost << x.get();
  }

  friend istream& operator >> (istream& ist, ModInt& x) { 
    string s; ist >> s;
    int size = s.size();
    x.get() = 0;
    rep(i, size) {
      x.get() *= 10;
      x.get() += s[i]-'0';
      x.get() %= MOD;
    }
    return ist;
  }
};

typedef ModInt<long long, 1000000007> mint;


int main() {

  int N; cin >> N;
  int ngnums[N]; zero(ngnums); rep(i, N) { int a; cin >> a; ngnums[a]++; }

  vector<mint> dp(N+1);
  dp[0] = 1;
  rep(i, N) for(int j=N-1; j>=0; j--)
    dp[j + 1] += dp[j] * ngnums[i]; // i1がダメな生徒 ∨ i2がダメな生徒　∨ ...

  vector<mint> P(N+1); P[0] = 1; REP(i, 1, N+1) P[i] = P[i-1] * i;

  mint ans = 0;

  // ダメな割当のない自由な割当 = 完全に自由な割当 - 1人だけダメな割当 + 2人だけダメな割当 - 3人だけダメな割当 + ... + (-1)^k k人ダメな割当 + ... (-1)^N N人ダメな割当

  rep(i, N+1) {
    ans += (i % 2 == 0 ? +1 : -1) * dp[i] * P[N - i];
  }

  cout << ans << endl;

  return 0;
}
0