結果
| 問題 |
No.329 全射
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 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;
}