結果
問題 | 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_pairtypedef 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 templateint main() {cin.tie(0);// source code// calc strling numbervector<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;}