結果

問題 No.654 Air E869120
ユーザー どららどらら
提出日時 2018-02-24 00:13:37
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 70 ms / 2,000 ms
コード長 3,454 bytes
コンパイル時間 2,945 ms
コンパイル使用メモリ 214,520 KB
実行使用メモリ 23,552 KB
最終ジャッジ日時 2024-10-10 05:30:58
合計ジャッジ時間 4,626 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 3 ms
5,248 KB
testcase_02 AC 3 ms
5,248 KB
testcase_03 AC 3 ms
5,248 KB
testcase_04 AC 2 ms
5,248 KB
testcase_05 AC 2 ms
5,248 KB
testcase_06 AC 2 ms
5,248 KB
testcase_07 AC 2 ms
5,248 KB
testcase_08 AC 2 ms
5,248 KB
testcase_09 AC 2 ms
5,248 KB
testcase_10 AC 11 ms
6,144 KB
testcase_11 AC 10 ms
5,760 KB
testcase_12 AC 10 ms
6,016 KB
testcase_13 AC 10 ms
5,888 KB
testcase_14 AC 10 ms
5,760 KB
testcase_15 AC 11 ms
5,888 KB
testcase_16 AC 69 ms
23,552 KB
testcase_17 AC 68 ms
23,424 KB
testcase_18 AC 69 ms
23,168 KB
testcase_19 AC 70 ms
23,296 KB
testcase_20 AC 36 ms
13,440 KB
testcase_21 AC 38 ms
13,568 KB
testcase_22 AC 34 ms
13,568 KB
testcase_23 AC 36 ms
13,568 KB
testcase_24 AC 36 ms
13,184 KB
testcase_25 AC 30 ms
11,776 KB
testcase_26 AC 33 ms
12,544 KB
testcase_27 AC 25 ms
9,984 KB
testcase_28 AC 23 ms
9,600 KB
testcase_29 AC 22 ms
9,344 KB
testcase_30 AC 7 ms
5,248 KB
testcase_31 AC 7 ms
5,248 KB
testcase_32 AC 7 ms
5,248 KB
testcase_33 AC 7 ms
5,248 KB
testcase_34 AC 7 ms
5,248 KB
testcase_35 AC 2 ms
5,248 KB
testcase_36 AC 2 ms
5,248 KB
testcase_37 AC 2 ms
5,248 KB
testcase_38 AC 2 ms
5,248 KB
testcase_39 AC 2 ms
5,248 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define REP(i,a,n) for(int i=(a); i<(int)(n); i++)
#define rep(i,n) REP(i,0,n)
#define FOR(it,c) for(__typeof((c).begin()) it=(c).begin(); it!=(c).end(); ++it)
#define ALLOF(c) (c).begin(), (c).end()
typedef long long ll;
typedef unsigned long long ull;

const ll INF = (ll)(1e19);

class Dinic {
  ll MAX_V;
  struct edge{ ll to, cap, rev, icap, flow; };

  vector< vector<edge> > G;
  vector<ll> level; //sからの距離
  vector<ll> iter; //どこまで調べたか

  void max_flow_bfs(ll s){
    fill(level.begin(), level.end(), -1);
    queue<ll> que;
    level[s] = 0;
    que.push(s);
    while(!que.empty()){
      ll v = que.front(); que.pop();
      for(ll i=0; i<G[v].size(); i++){
        edge &e = G[v][i];
        if(e.cap>0 && level[e.to]<0){
          level[e.to] = level[v] + 1;
          que.push(e.to);
        }
      }
    }
  }
  ll max_flow_dfs(ll v, ll t, ll f){
    if(v==t) return f;
    for(ll &i=iter[v]; i<G[v].size(); i++){
      edge &e = G[v][i];
      if(e.cap>0 && level[v]<level[e.to]){
        ll d = max_flow_dfs(e.to, t, min(f, e.cap));
        if(d>0){
          e.cap -= d;
          G[e.to][e.rev].cap += d;
          e.flow += d;
          return d;
        }
      }
    }
    return 0;
  }

public:
  Dinic(ll N):MAX_V(N),G(N),level(N),iter(N){
  }

  void add_edge(ll from, ll to, ll cap){
    G[from].push_back((edge){to, cap, (ll)G[to].size(), cap, 0});
    G[to].push_back((edge){from, 0, (ll)G[from].size()-1, 0, 0});
  }

  ll get_flow(ll from, ll to){ //untried
    rep(i,G[from].size()){
      if(G[from][i].to == to){
        return G[from][i].flow;
      }
    }
    return -1;
  }

  ll max_flow(ll s, ll t){
    ll flow = 0;
    while(true){
      max_flow_bfs(s);
      if(level[t]<0) return flow;
      fill(iter.begin(), iter.end(), 0);
      ll f;
      while((f = max_flow_dfs(s, t, INF))>0){
        flow += f;
      }
    }
  }
};

int N, M;
ll d;
vector<vector<ll>> air;


int main(){
  cin >> N >> M >> d;
  rep(i,M){
    ll u, v, p, q, w;
    cin >> u >> v >> p >> q >> w;
    air.push_back({u,v,p,q,w});
  }

  {
    map<vector<ll>,ll> check;
    rep(i,air.size()){
      vector<ll> tmp = {air[i][0], air[i][1], air[i][2], air[i][3]};
      check[tmp] += air[i][4];
    }
    air.clear();
    FOR(it,check){
      vector<ll> tmp = it->first;
      tmp.push_back(it->second);
      air.push_back(tmp);
    }
  }

  vector<set<int>> G(N+1);
  rep(i,air.size()){
    G[air[i][0]].insert(air[i][2]);
    G[air[i][1]].insert(air[i][3]);
  }
  
  map<vector<ll>,int> memo;
  int id = 0;
  rep(i,G.size()){
    FOR(it,G[i]){
      memo[{(ll)i,(*it), 0}] = id++;
      memo[{(ll)i,(*it), 1}] = id++;
    }
  }


  Dinic dinic(id+2);
  rep(i,G.size()){
    FOR(it1,G[i]){
      FOR(it2,G[i]){
        int prev = *it1;
        int next = *it2;
        if(next-prev>=d){
          dinic.add_edge(memo[{(ll)i,prev,1}], memo[{(ll)i,next,0}], INF);
        }
      }
    }
  }
  
  rep(i,air.size()){
    dinic.add_edge(memo[{air[i][0],air[i][2],0}], memo[{air[i][1],air[i][3],1}], air[i][4]);
  }

  rep(i,G.size()){
    FOR(it,G[i]){
      if(i==1) dinic.add_edge(id, memo[{(ll)i,*it,0}], INF);
      if(i==1) dinic.add_edge(id, memo[{(ll)i,*it,1}], INF);
      if(i==N) dinic.add_edge(memo[{(ll)i,*it,0}], id+1, INF);
      if(i==N) dinic.add_edge(memo[{(ll)i,*it,1}], id+1, INF);
    }
  }
  cout << dinic.max_flow(id, id+1) << endl;

  return 0;
}
0