結果

問題 No.329 全射
ユーザー rickythetarickytheta
提出日時 2015-12-22 22:26:25
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 22 ms / 2,000 ms
コード長 2,373 bytes
コンパイル時間 1,753 ms
コンパイル使用メモリ 171,404 KB
実行使用メモリ 11,480 KB
最終ジャッジ日時 2024-09-18 18:49:23
合計ジャッジ時間 3,365 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 6 ms
10,852 KB
testcase_01 AC 5 ms
10,736 KB
testcase_02 AC 5 ms
10,820 KB
testcase_03 AC 5 ms
10,632 KB
testcase_04 AC 5 ms
10,816 KB
testcase_05 AC 5 ms
10,848 KB
testcase_06 AC 5 ms
10,740 KB
testcase_07 AC 6 ms
10,668 KB
testcase_08 AC 5 ms
10,864 KB
testcase_09 AC 5 ms
10,712 KB
testcase_10 AC 5 ms
10,856 KB
testcase_11 AC 6 ms
10,768 KB
testcase_12 AC 5 ms
10,816 KB
testcase_13 AC 21 ms
10,948 KB
testcase_14 AC 11 ms
11,160 KB
testcase_15 AC 15 ms
10,996 KB
testcase_16 AC 20 ms
11,180 KB
testcase_17 AC 17 ms
10,804 KB
testcase_18 AC 22 ms
11,480 KB
testcase_19 AC 22 ms
10,832 KB
testcase_20 AC 21 ms
10,936 KB
testcase_21 AC 18 ms
11,012 KB
testcase_22 AC 18 ms
10,712 KB
testcase_23 AC 12 ms
11,312 KB
testcase_24 AC 10 ms
10,688 KB
testcase_25 AC 13 ms
11,064 KB
testcase_26 AC 10 ms
10,788 KB
testcase_27 AC 6 ms
10,932 KB
testcase_28 AC 5 ms
10,944 KB
testcase_29 AC 6 ms
10,836 KB
testcase_30 AC 5 ms
10,924 KB
testcase_31 AC 5 ms
10,836 KB
testcase_32 AC 6 ms
10,892 KB
testcase_33 AC 5 ms
11,000 KB
testcase_34 AC 5 ms
10,808 KB
testcase_35 AC 6 ms
10,836 KB
testcase_36 AC 6 ms
10,868 KB
testcase_37 AC 5 ms
10,992 KB
testcase_38 AC 6 ms
10,780 KB
testcase_39 AC 7 ms
10,664 KB
testcase_40 AC 5 ms
10,880 KB
testcase_41 AC 5 ms
10,776 KB
testcase_42 AC 6 ms
10,968 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef complex<double> P;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
#define REP(i,n) for(ll i=0;i<n;++i)
#define REPR(i,n) for(ll i=1;i<n;++i)
#define FOR(i,a,b) for(ll i=a;i<b;++i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define MOD (ll)(1e9+7)
#define ADD(a,b) a=((a)+(b))%MOD
#define FIX(a) ((a)%MOD+MOD)%MOD

// O(logMOD)
ll inv(ll x){
  ll t = MOD-2;
  ll r = 1;
  while(t){
    if(t&1)r=r*x%MOD;
    x=x*x%MOD;
    t>>=1;
  }
  return r;
}

// O(logx)
ll mypow(ll a,ll x){
  if(x==0)return 1;
  ll p = mypow(a,x/2);
  p = p*p%MOD;
  if(x&1) p = p*a%MOD;
  return p;
}

ll fact[1010];
ll ifact[1010];

ll comb(ll a,ll b){
  if(a-b<0)return 0;
  return fact[a]*ifact[a-b]%MOD*ifact[b]%MOD;
}

ll zensha[1010][1010];

// // O(k)
// ll zensha(ll n,ll k){
//   if(zmemo[n][k])return zmemo[n][k];
//   ll ret = 0;
//   REPR(i,k+1){
//     ll tmp = 1;
//     if((k-i)%2==1)tmp = MOD-1;
//     tmp = tmp*comb(k,i)%MOD;
//     tmp = tmp*mypow(i,n)%MOD;
//     ret = (ret+tmp)%MOD;
//   }
//   return zmemo[n][k]=ret;
// }

int main(){
  fact[0] = 1;
  REPR(i,1010) fact[i] = fact[i-1]*i%MOD;
  REP(i,1010) ifact[i] = inv(fact[i]);

  zensha[1][1] = 1;
  FOR(i,2,1010)FOR(j,1,1010){
    if(j>i)break;
    zensha[i][j]=j*(zensha[i-1][j]+zensha[i-1][j-1])%MOD;
  }

  ll n,m;
  cin>>n>>m;
  vl w(n);
  REP(i,n) cin>>w[i];
  vl G[n];
  REP(i,m){
    ll x,y;
    cin>>x>>y;
    --x;--y;
    if(x==y)continue;
    G[x].push_back(y);
  }
  ll result = 0;
  REP(src,n){
    vl flow(n,0);
    flow[src] = w[src];
    priority_queue<pll> Q;
    Q.push(make_pair(flow[src],src));
    vector<bool> used(n,false);
    while(!Q.empty()){
      pll P = Q.top(); Q.pop();
      ll f = P.first;
      ll u = P.second;
      if(used[u])continue;
      used[u] = true;
      REP(ed,G[u].size()){
        ll v = G[u][ed];
        if(flow[v]>=min(f,w[v]))continue;
        flow[v]=min(f,w[v]);
        Q.push(make_pair(flow[v],v));
      }
    }
    REP(to,n){
      if(flow[to]==w[to]){
        // result += zensha(w[src],w[to]);
        result += zensha[w[src]][w[to]];
        result %= MOD;
      }
    }
  }
  cout << result << endl;
  return 0;
}
0