結果
問題 | No.329 全射 |
ユーザー |
![]() |
提出日時 | 2016-05-24 22:08:44 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 24 ms / 2,000 ms |
コード長 | 1,387 bytes |
コンパイル時間 | 731 ms |
コンパイル使用メモリ | 83,300 KB |
実行使用メモリ | 7,604 KB |
最終ジャッジ日時 | 2024-10-07 06:49:50 |
合計ジャッジ時間 | 2,328 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 40 |
ソースコード
#define _CRT_SECURE_NO_WARNINGS#define _USE_MATH_DEFINES#include<cstdio>#include<algorithm>#include<vector>#include<iostream>#include<queue>#include<stack>#include<cstring>#include<string>#include<map>#include<set>#include<sstream>#include<bitset>#include<numeric>#include<climits>#include<complex>using namespace std;typedef long long ll;template<typename A,typename B>inline void chmin(A &a, B b) { if (a > b)a = b; }template<typename A,typename B>inline void chmax(A &a, B b) { if (a < b)a = b; }const int mod = 1000000007;int N, M;int W[200];int f[200][200];int dp[1010][1010];int F[1010];signed main() {cin >> N >> M;for (int i = 0; i < N; i++)cin >> W[i];for (int i = 0; i < N; i++)f[i][i] = W[i];for (int i = 0; i < M; i++) {int a, b;cin >> a >> b;a--; b--;f[a][b] = min(W[a], W[b]);}for (int k = 0; k < N; k++) {for (int i = 0; i < N; i++) {for (int j = 0; j < N; j++) {chmax(f[i][j], min(f[i][k], f[k][j]));}}}dp[0][0] = 1;F[0] = 1;for (int i = 1; i <= 1000; i++) {F[i] = (ll)F[i - 1] * i%mod;for (int j = 1; j <= i; j++) {dp[i][j] = (dp[i - 1][j - 1] + (ll)j*dp[i - 1][j]) % mod;}}int ans = 0;for (int i = 0; i < N; i++) {for (int j = 0; j < N; j++) {if (f[i][j] >= W[j])ans = (ans + (ll)dp[W[i]][W[j]] * F[W[j]]) % mod;}}cout << ans << endl;return 0;}