标签:== for ros ret min cos front lin 中间
题解:直接连边跑最小费用最大流即可。
#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#define P(X,Y) ((X-1)*B+Y)
#define Q(X,Y) ((X-1)*B+Y+A*B)
using namespace std;
int A,B,n,m,sum,ans,cnt,S,T;
int cost[100000],flow[100000],next[100000],head[100000],to[100000],dis[100000],inq[100000],pv[100000],pe[100000];
queue<int> q;
inline int rd()
{
int ret=0,f=1; char gc=getchar();
while(gc<‘0‘||gc>‘9‘) {if(gc==‘-‘) f=-f; gc=getchar();}
while(gc>=‘0‘&&gc<=‘9‘) ret=ret*10+gc-‘0‘,gc=getchar();
return ret*f;
}
inline void add(int a,int b,int c,int d)
{
to[cnt]=b,cost[cnt]=c,flow[cnt]=d,next[cnt]=head[a],head[a]=cnt++;
to[cnt]=a,cost[cnt]=-c,flow[cnt]=0,next[cnt]=head[b],head[b]=cnt++;
}
inline int bfs()
{
memset(dis,0x3f,sizeof(dis));
dis[S]=0,q.push(S);
int u,i;
while(!q.empty())
{
u=q.front(),q.pop(),inq[u]=0;
for(i=head[u];i!=-1;i=next[i]) if(dis[to[i]]>dis[u]+cost[i]&&flow[i])
{
dis[to[i]]=dis[u]+cost[i],pe[to[i]]=i,pv[to[i]]=u;
if(!inq[to[i]]) inq[to[i]]=1,q.push(to[i]);
}
}
return dis[T]<0x3f3f3f3f;
}
int main()
{
A=rd(),B=rd(),S=0,T=2*A*B+1;
int i,j,a,b,c,d,mf;
memset(head,-1,sizeof(head));
for(i=1;i<=A;i++) for(j=1;j<=B;j++) add(P(i,j),Q(i,j),rd(),1);
m=rd();
for(i=1;i<=m;i++)
{
a=rd(),b=rd(),c=rd(),d=rd();
add(Q(a,b),P(c,d),0,1),add(Q(c,d),P(a,b),0,1);
}
n=rd();
for(i=1;i<=n;i++) a=rd(),b=rd(),add(S,P(a,b),0,1);
for(i=1;i<=n;i++) a=rd(),b=rd(),add(Q(a,b),T,0,1);
while(bfs())
{
mf=1<<30;
for(i=T;i!=S;i=pv[i]) mf=min(mf,flow[pe[i]]);
sum+=mf,ans+=dis[T]*mf;
for(i=T;i!=S;i=pv[i]) flow[pe[i]]-=mf,flow[pe[i]^1]+=mf;
}
if(sum<n) printf("-1");
else printf("%d",ans);
return 0;
}
标签:== for ros ret min cos front lin 中间
原文地址:http://www.cnblogs.com/CQzhangyu/p/7747194.html