結果

問題 No.654 Air E869120
コンテスト
ユーザー chocorusk
提出日時 2019-05-20 23:55:32
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 56 ms / 2,000 ms
コード長 2,579 bytes
コンパイル時間 954 ms
コンパイル使用メモリ 115,520 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-17 06:44:37
合計ジャッジ時間 2,910 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 5
other AC * 35
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#define popcount __builtin_popcount
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
const ll INF=1e16;
struct edge {int to; long long int cap; int rev;} ;
 
vector<edge> G[2010];
int level[2010];
int iter[2010];
 
void add_edge(int from, int to, long long int cap){
	edge e;
	e.to=to, e.cap=cap, e.rev=G[to].size();
	G[from].push_back(e);
	e.to=from, e.cap=0, e.rev=G[from].size()-1;
	G[to].push_back(e);
}
 
void bfs(int s){
	memset(level, -1, sizeof(level));
	queue<int> que;
	level[s]=0;
	que.push(s);
	while(!que.empty()){
		int v=que.front();
		que.pop();
		for(int 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);
			}
		}
	}
}
 
long long int dfs(int v, int t, long long int f){
	if(v==t) return f;
	for(int &i=iter[v]; i<G[v].size(); i++){
		edge &e=G[v][i];
		if(e.cap>0 && level[v]<level[e.to]){
			long long int d=dfs(e.to, t, min(f, e.cap));
			if(d>0){
				e.cap-=d;
				G[e.to][e.rev].cap+=d;
				return d;
			}
		}
	}
	return 0;
}
 
long long int max_flow(int s, int t){
	long long int flow=0;
	while(1){
		bfs(s);
		if(level[t]<0) return flow;
		memset(iter, 0, sizeof(iter));
		long long int f;
		while((f=dfs(s, t, INF))>0){
			flow+=f;
		}
	}
}

int main()
{
    int n, m, d; cin>>n>>m>>d;
    map<P, int> mp;
    int c=0;
    for(int i=0; i<m; i++){
        int u, v, p, q; ll w;
        cin>>u>>v>>p>>q>>w;
        if(mp.find({2*u, p})==mp.end()){
            mp[{2*u, p}]=c; c++;
        }
        if(mp.find({2*v+1, q})==mp.end()){
            mp[{2*v+1, q}]=c; c++;
        }
        add_edge(mp[{2*u, p}], mp[{2*v+1, q}], w);
    }
    for(auto p:mp){
        for(auto q:mp){
            int u=p.first.first, v=q.first.first;
            if(u%2==0 || v%2==1) continue;
            if(u/2!=v/2) continue;
            int x=p.first.second, y=q.first.second;
            if(x+d<=y){
                add_edge(p.second, q.second, INF);
            }
        }
    }
    int s=c, t=c+1;
    for(auto p:mp){
        if(p.first.first==2){
            add_edge(s, p.second, INF);
        }else if(p.first.first==2*n+1){
            add_edge(p.second, t, INF);
        }
    }
    cout<<max_flow(s, t)<<endl;
    return 0;
}
0