LIGHT - Hệ thống đèn
Tác giả: skyvn97
Ngôn ngữ: C++
#include<bits/stdc++.h>
#define FOR(i,a,b) for (int i=(a);i<=(b);i=i+1)
#define REP(i,n) for (int i=0;i<(n);i=i+1)
#define FORE(i,v) for (__typeof((v).begin()) i=(v).begin();i!=(v).end();i++)
using namespace std;
const int INF=(int)1e9+7;
class DinicFlow {
private:
vector<int> dist,head,q,work;
vector<int> capa,flow,next,point;
int n,m;
public:
DinicFlow() {
n=0;m=0;
}
DinicFlow(int n,int _m) {
this->n=n;
this->m=0;
dist=vector<int>(n+7);
head=vector<int>(n+7,-1);
q=vector<int>(n+7);
work=vector<int>(n+7);
int m=2*_m;
capa=vector<int>(m+7);
flow=vector<int>(m+7,0);
next=vector<int>(m+7);
point=vector<int>(m+7);
}
void addedge(int u,int v,int c1,int c2) {
point[m]=v;capa[m]=c1;flow[m]=0;next[m]=head[u];head[u]=m;m++;
point[m]=u;capa[m]=c2;flow[m]=0;next[m]=head[v];head[v]=m;m++;
}
bool bfs(int s,int t) {
FORE(it,dist) *it=-1;
int sz=1;
q[0]=s;dist[s]=0;
REP(x,sz) {
int u=q[x];
for (int i=head[u];i>=0;i=next[i])
if (dist[point[i]]<0 && flow[i]<capa[i]) {
dist[point[i]]=dist[u]+1;
q[sz]=point[i];
sz++;
}
}
return (dist[t]>=0);
}
int dfs(int s,int t,int f) {
if (s==t) return (f);
for (int &i=work[s];i>=0;i=next[i])
if (dist[point[i]]==dist[s]+1 && flow[i]<capa[i]) {
int d=dfs(point[i],t,min(f,capa[i]-flow[i])) ;
if (d>0) {
flow[i]+=d;
flow[i^1]-=d;
return (d);
}
}
return (0);
}
int maxflow(int s,int t) {
int ret=0;
while (bfs(s,t)) {
FOR(i,1,n) work[i]=head[i];
while (true) {
int d=dfs(s,t,INF);
if (d<=0) break;
ret+=d;
}
}
return (ret);
}
};
int m,n,k;
DinicFlow g;
void init(void) {
scanf("%d",&m);
scanf("%d",&n);
scanf("%d",&k);
g=DinicFlow(m+n+2,m+n+k);
FOR(i,1,m) {
int t;
scanf("%d",&t);
g.addedge(m+n+1,i,t,0);
}
FOR(i,1,n) {
int t;
scanf("%d",&t);
g.addedge(m+i,m+n+2,t,0);
}
REP(i,k) {
int x,y;
scanf("%d",&x);
scanf("%d",&y);
g.addedge(x,m+y,INF,0);
}
printf("%d",g.maxflow(m+n+1,m+n+2));
}
int main(void) {
init();
return 0;
}