結果
問題 | No.329 全射 |
ユーザー |
![]() |
提出日時 | 2015-12-24 21:29:43 |
言語 | PyPy2 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 852 bytes |
コンパイル時間 | 1,397 ms |
コンパイル使用メモリ | 76,436 KB |
実行使用メモリ | 86,400 KB |
最終ジャッジ日時 | 2024-09-18 22:59:46 |
合計ジャッジ時間 | 12,054 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 11 WA * 29 |
ソースコード
import syssys.setrecursionlimit(1000000)p = 10**9+7N,M = map(int,raw_input().split())w = map(int,raw_input().split())w_max = max(w)C = [[0]*(i+1) for i in xrange(w_max+1)]C[0][0] = 1for i in xrange(w_max):for j in xrange(i+1):C[i][j] %= pC[i+1][j] += C[i][j]C[i+1][j+1] += C[i][j]for i in xrange(w_max+1): C[w_max][i] %= pE = [[] for i in xrange(N)]for t in xrange(M):i,j = map(lambda x:int(x)-1,raw_input().split())E[i].append(j)proj = set()visited = set()def dfs(s,h,w_min):if h in visited: returnvisited.add(h)if w[s] >= w[h] and w[h] <= w_min:proj.add((s,h))w_min = min(w_min,w[h])for nxt in E[h]:dfs(s,nxt,w_min)for i in xrange(N):visited = set()dfs(i,i,w[i])ans = 0for i,j in proj:wi = w[i]; wj = w[j]ans += sum(C[wj][k]*pow(wj-k,wi,p)*((-1)**(k%2)) for k in xrange(wj))ans %= pprint ans