結果

問題 No.329 全射
コンテスト
ユーザー piyoko_212
提出日時 2015-12-23 12:18:56
言語 C++11
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 149 ms / 2,000 ms
コード長 1,088 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 511 ms
コンパイル使用メモリ 55,460 KB
実行使用メモリ 22,512 KB
最終ジャッジ日時 2026-04-06 23:06:30
合計ジャッジ時間 5,188 ms
ジャッジサーバーID
（参考情報） 		judge1_1 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 40
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<stdio.h>
#include<algorithm>
using namespace std;
long long mod=1000000007;
long long rj[1100][1100];
long long C[1100][1100];
int c[110000];
int g[210][210];
int main(){
	for(int i=0;i<1100;i++){
		rj[i][0]=1;
		for(int j=1;j<1100;j++)rj[i][j]=rj[i][j-1]*i%mod;
	}
	C[0][0]=1;
	for(int i=0;i<1050;i++){
		for(int j=0;j<=i;j++){
			C[i+1][j]=(C[i+1][j]+C[i][j])%mod;
			C[i+1][j+1]=(C[i+1][j+1]+C[i][j])%mod;
		}
	}
	int a,b;scanf("%d%d",&a,&b);
	for(int i=0;i<a;i++)scanf("%d",c+i);
	for(int i=0;i<a;i++)for(int j=0;j<a;j++)g[i][j]=-999999999;
	for(int i=0;i<a;i++)g[i][i]=999999999;
	
	for(int i=0;i<b;i++){
		int p,q;scanf("%d%d",&p,&q);p--;q--;
		g[p][q]=min(c[q],c[p]);
	}
	for(int k=0;k<a;k++)for(int i=0;i<a;i++)for(int j=0;j<a;j++)g[i][j]=max(g[i][j],min(g[i][k],g[k][j]));
	long long ret=0;
	for(int i=0;i<a;i++)for(int j=0;j<a;j++){
		if(g[i][j]>=c[j]){
			for(int k=0;k<c[j];k++){
				long long tmp=rj[c[j]-k][c[i]]*C[c[j]][k]%mod;
				if(k%2)ret=(ret+mod-tmp)%mod;
				else ret=(ret+tmp)%mod;
			}
		//	printf("%d %d: %lld\n",i,j,ret);
		}
	}
	printf("%lld\n",ret);
}
0