結果

問題 No.329 全射
コンテスト
ユーザー Pulmn
提出日時 2018-07-17 16:43:56
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 62 ms / 2,000 ms
コード長 1,431 bytes
コンパイル時間 1,332 ms
コンパイル使用メモリ 174,160 KB
実行使用メモリ 11,468 KB
最終ジャッジ日時 2024-11-25 08:39:57
合計ジャッジ時間 3,450 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 40
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <bits/stdc++.h>
#define syosu(x) fixed<<setprecision(x)
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
typedef pair<int,int> P;
typedef pair<double,double> pdd;
typedef pair<ll,ll> pll;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<double> vd;
typedef vector<vd> vvd;
typedef vector<ll> vl;
typedef vector<vl> vvl;
typedef vector<string> vs;
typedef vector<P> vp;
typedef vector<vp> vvp;
typedef vector<pll> vpll;
typedef pair<int,P> pip;
typedef vector<pip> vip;
const int inf=1<<30;
const ll INF=1ll<<60;
const double pi=acos(-1);
const double eps=1e-8;
const ull mod=1e9+7;
const int dx[4]={-1,0,1,0},dy[4]={0,-1,0,1};

const int M=1005;
ll n,m,dp[M][M];
vi a,b;
vip e;
vvi g;

void dfs(int v,int t){
	b[v]=t;
	for(auto u:g[v]) if(b[u]==inf) dfs(u,t);
}

int main(){
	cin>>n>>m;
	a=vi(n);
	e=vip(m);
	for(int i=0;i<n;i++) cin>>a[i];
	for(int i=0;i<m;i++){
		int u,v;
		cin>>u>>v;
		u--;v--;
		e[i]={a[v],{u,v}};
	}
	sort(e.rbegin(),e.rend());
	dp[1][1]++;
	for(int i=2;i<M;i++) for(int j=1;j<=i;j++) dp[i][j]=(dp[i-1][j]+(j?dp[i-1][j-1]:0))*j%mod;
	ll res=0;
	for(int i=0;i<n;i++){
		int A=a[i];
		b=vi(n,inf);
		b[i]=A;
		g=vvi(n);
		for(int j=0;j<m;j++){
			P p=e[j].second;
			int u=p.first,v=p.second;
			g[u].push_back(v);
			if(b[u]<inf&&b[v]==inf) dfs(v,min(A,a[v]));
		}
		for(int j=0;j<n;j++) if(a[j]<=b[j]&&b[j]<inf) (res+=dp[A][a[j]])%=mod;
	}
	cout<<res<<endl;
}
0