結果
問題 | No.329 全射 |
ユーザー | ctyl_0 |
提出日時 | 2015-12-23 02:03:41 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 25 ms / 2,000 ms |
コード長 | 1,537 bytes |
コンパイル時間 | 899 ms |
コンパイル使用メモリ | 97,532 KB |
実行使用メモリ | 11,520 KB |
最終ジャッジ日時 | 2024-09-18 19:25:15 |
合計ジャッジ時間 | 2,681 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 40 |
ソースコード
#include <iostream> #include <iomanip> #include <vector> #include <algorithm> #include <numeric> #include <functional> #include <cmath> #include <queue> #include <stack> #include <set> #include <map> #include <sstream> #include <string> #define repd(i,a,b) for (int i=(a);i<(b);i++) #define rep(i,n) repd(i,0,n) #define var auto #define mod 1000000007 #define inf 2147483647 #define nil -1 #define mp make_pair typedef long long ll; using namespace std; template <typename T> inline void output(T a, int p) { if(p){ cout << fixed << setprecision(p) << a << "\n"; } else{ cout << a << "\n"; } } // end of template int main() { cin.tie(0); // source code // calc strling number vector<vector<ll>> S(1001, vector<ll>(1001, 0)); repd(i, 1, 1001) S[i][1] = 1; repd(i, 2, 1001) repd(j, 2, i + 1){ S[i][j] = j * (S[i - 1][j] + S[i - 1][j - 1]); S[i][j] %= mod; } int N, M; cin >> N >> M; vector<vector<int>> E(N, vector<int>(N, 0)); vector<int> W(N); rep(i, N) cin >> W[i], E[i][i] = W[i]; vector<pair<int, int>> A(M); rep(i, M) cin >> A[i].first >> A[i].second, E[A[i].first - 1][A[i].second - 1] = W[A[i].second - 1]; rep(k, N) rep(i, N) rep(j, N){ E[i][j] = max(E[i][j], min(E[i][k], E[k][j])); } ll ret = 0; rep(i, N) rep(j, N){ if (E[i][j] >= W[j]) { ret += S[W[i]][W[j]]; ret %= mod; } } output(ret, 0); return 0; }