結果

問題 No.329 全射
ユーザー motimoti
提出日時 2015-12-28 02:10:38
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 25 ms / 2,000 ms
コード長 3,217 bytes
コンパイル時間 848 ms
コンパイル使用メモリ 106,872 KB
実行使用メモリ 12,288 KB
最終ジャッジ日時 2024-09-19 07:33:49
合計ジャッジ時間 2,370 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
11,768 KB
testcase_01 AC 5 ms
11,972 KB
testcase_02 AC 6 ms
11,912 KB
testcase_03 AC 5 ms
11,672 KB
testcase_04 AC 5 ms
11,564 KB
testcase_05 AC 6 ms
11,996 KB
testcase_06 AC 5 ms
11,628 KB
testcase_07 AC 4 ms
12,136 KB
testcase_08 AC 5 ms
11,668 KB
testcase_09 AC 4 ms
12,040 KB
testcase_10 AC 5 ms
11,424 KB
testcase_11 AC 5 ms
11,744 KB
testcase_12 AC 5 ms
11,972 KB
testcase_13 AC 23 ms
11,716 KB
testcase_14 AC 13 ms
12,024 KB
testcase_15 AC 15 ms
11,776 KB
testcase_16 AC 18 ms
12,288 KB
testcase_17 AC 18 ms
11,852 KB
testcase_18 AC 23 ms
11,872 KB
testcase_19 AC 25 ms
11,932 KB
testcase_20 AC 23 ms
12,000 KB
testcase_21 AC 18 ms
11,868 KB
testcase_22 AC 22 ms
11,764 KB
testcase_23 AC 9 ms
12,020 KB
testcase_24 AC 9 ms
11,576 KB
testcase_25 AC 13 ms
11,772 KB
testcase_26 AC 9 ms
11,688 KB
testcase_27 AC 6 ms
11,888 KB
testcase_28 AC 4 ms
12,064 KB
testcase_29 AC 5 ms
11,976 KB
testcase_30 AC 6 ms
11,600 KB
testcase_31 AC 5 ms
11,584 KB
testcase_32 AC 5 ms
11,900 KB
testcase_33 AC 12 ms
11,972 KB
testcase_34 AC 9 ms
11,860 KB
testcase_35 AC 12 ms
11,712 KB
testcase_36 AC 11 ms
11,704 KB
testcase_37 AC 9 ms
11,396 KB
testcase_38 AC 10 ms
11,672 KB
testcase_39 AC 11 ms
11,764 KB
testcase_40 AC 13 ms
11,864 KB
testcase_41 AC 9 ms
11,916 KB
testcase_42 AC 12 ms
11,468 KB
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:119:42: warning: iteration 1009 invokes undefined behavior [-Waggressive-loop-optimizations]
  119 |   rep(i, MATH_COMBINATION_MAX) fact[i+1] = fact[i] * (i+1) % MOD;
      |                                ~~~~~~~~~~^~~~~~~~~~~~~~~~~~~~~~~
main.cpp:24:33: note: within this loop
   24 | #define REP(i,a,b) for(int i=a;i<(int)b;i++)
      |                                 ^
main.cpp:25:18: note: in expansion of macro ‘REP’
   25 | #define rep(i,n) REP(i,0,n)
      |                  ^~~
main.cpp:119:3: note: in expansion of macro ‘rep’
  119 |   rep(i, MATH_COMBINATION_MAX) fact[i+1] = fact[i] * (i+1) % MOD;
      |   ^~~

ソースコード

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)

constexpr int inf = 1<<29;

int const MOD = 1e9+7;
typedef long long ll;


constexpr int MATH_COMBINATION_MAX = 1010;

/***********************************************************************************/
#define DEF_MAX_BOUND(assignment) enum { assignment };
#define MEMO_BEGIN(Identifier)  \
namespace memorized { namespace Identifier {  \
  bool __computed = false; inline bool computed() {if(__computed) { return true; } else { __computed = true; return false; }}
#define MEMO_END }}

#define USE_MEMO(Identifier)        using memorized::Identifier::
#define COMPUTED_RESULT(Identifier) if(memorized::Identifier::computed()) return
#define USE_CONSTANT(Identifier)    using memorized::Identifier::
/***********************************************************************************/


/*
  スターリング数
  異なるN個をちょうどKグループに分けるときの方法の数
  高々Kグループに分ける場合は、グループi(1<=i<=K)をN回選ぶ方法の数となり、これは単なる重複順列でK^Nとなる
  Verified: http://yukicoder.me/problems/663
 */

namespace math { namespace combination {

MEMO_BEGIN(stirling_number)
  DEF_MAX_BOUND(max_value = MATH_COMBINATION_MAX)
  ll dp[max_value][max_value];
MEMO_END

int stirling_number(int N, int K) {
  USE_MEMO(stirling_number)         dp;
  COMPUTED_RESULT(stirling_number)  dp[N][K];

  USE_CONSTANT(stirling_number) max_value;

  REP(i, 1, max_value) {  // undefined when i = 0
    dp[i][1] = dp[i][i] = 1;
    REP(j, 2, i) {
      dp[i][j] = dp[i-1][j-1] + j * dp[i-1][j];
      dp[i][j] %= MOD;
    }
  }

  return dp[N][K];
}

}}

int N, M;
vector<int> g[222];
ll fact[MATH_COMBINATION_MAX];
ll dp[300][300];

int main() {

  cin >> N >> M;

  vector<ll> ws(N);
  rep(i, N) {
    cin >> ws[i];
  }

  rep(i, M) {
    int a, b; cin >> a >> b; a--, b--;
    g[a].push_back(b);
  }

  zero(dp);
  rep(i, N) {
    dp[i][i] = ws[i];
    rep(j, g[i].size()) {
      int tar = g[i][j];
      dp[i][tar] = min(dp[i][i], ws[tar]);
    }
  }

  rep(k, N) rep(i, N) rep(j, N) {
    maximize(dp[i][j], min(dp[i][k], dp[k][j]));
  }

  fact[0] = 1;
  rep(i, MATH_COMBINATION_MAX) fact[i+1] = fact[i] * (i+1) % MOD;

  // 異なるN個を異なるK個に振り分ける方法の数(全射の数) := S(N, K) * K!
  using math::combination::stirling_number;

  ll ans = 0;

  rep(i, N) rep(j, N) {
    if(dp[i][j] < ws[j]) { continue; }
    ans += stirling_number(ws[i], ws[j]) * fact[ws[j]] % MOD;
    ans %= MOD;
  }

  cout << ans << endl;

  return 0;
}
0