题目大意:
有m(1<=m<=1000)个猪圈,n(1<=n<=100)名客人,第i个猪圈有ci头猪.
第i人第i天去买猪,有ki把钥匙,最多买bi头猪.
第i天的时候,钥匙对应的猪圈会打开,这个时候可以对猪圈任意分配猪的数量(和不会超过打开的猪圈的猪的数量之和.)
问最多能卖多少猪.
题目思路:
(1)建立一个超级源点和一个超级汇点.
(2)源点向所有猪圈建边,容量为猪圈初始的猪的数量,这样保证了网络中的流量不会超过猪的总数量.
(3)所有客人向汇点建边,容量为其需求量,这样保证了到汇点的流量不会超过总需求量.
(3)猪圈向拥有其钥匙的人建边,容量为无穷大,能流过去即可.www.2cto.com
(4)先买的顾客向之后和他拥有相同钥匙的顾客建边,容量为无穷大,这样使得前面猪圈交换的信息可以向下传递,以达到最优值.
代码:
[cpp]
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;
#define ll long long
#define ls rt<<1
#define rs ls|1
#define lson l,mid,ls
#define rson mid+1,r,rs
#define middle (l+r)>>1
#define eps (1e-8)
#define clr_all(x,c) memset(x,c,sizeof(x))
#define clr(x,c,n) memset(x,c,sizeof(x[0])*(n+1))
#define MOD 1000000009
#define INF 0x3f3f3f3f
#define PI acos(-1.0)
#define _max(x,y) (((x)>(y)) (x):(y))
#define _min(x,y) (((x)<(y)) (x):(y))
#define _abs(x) ((x)<0 (-(x)):(x))
#define getmin(x,y) (x= ((x)<0 || (y)<(x)) (y):(x))
#define getmax(x,y) (x= ((y)>(x)) (y):(x))
template void _swap(T &x,T &y){T t=x;x=y;y=t;}
int TS,cas=1;
const int M=1100+5;
int n,m;
struct sap{
typedef int type;
struct edge{
int to,next;
type flow;
}edg[999999];
int n,m;
int g[M],cur[M],pre[M];
int dis[M],gap[M];
int q[M],head,tail;
void init(int _n){n=_n,m=0,clr_all(g,-1);}
void insert(int u,int v,type f,type c=0){
edg[m].to=v,edg[m].flow=f,edg[m].next=g[u],g[u]=m++;
edg[m].to=u,edg[m].flow=0,edg[m].next=g[v],g[v]=m++;
}
bool bfs(int s,int t){
clr(gap,0,n),clr(dis,-1,n);
gap[dis[t]=0]++,dis[s]=-1;
for(q[head=tail=0]=t;head<=tail;){
int u=q[head++];
for(int i=g[u];i!=-1;i=edg[i].next){
edge& e=edg[i],ee=edg[i^1];
if(dis[e.to]==-1 && ee.flow>0){
gap[dis[e.to]=dis[u]+1]++;
q[++tail]=e.to;
}
}
}
return dis[s]!=-1;
}
type maxFlow(int s,int t){
if(!bfs(s,t)) return 0;
type res=0,a;
int i;
for(i=0;i
pre[s]=s,cur[s]=g[s],cur[t]=g[t];
for(int u=s;dis[s]
if(u==t){
for(a=-1;u!=s;u=pre[u])
getmin(a,edg[cur[pre[u]]].flow);
for(u=t;u!=s;u=pre[u]){
edg[cur[pre[u]]].flow-=a;
edg[cur[pre[u]]^1].flow+=a;
}
res+=a;
}
bool ok=0;
for(i=cur[u];i!=-1;i=edg[i].next){
edge& e=edg[i];
if(dis[u]==dis[e.to]+1 && e.flow>0){
cur[u]=i,pre[e.to]=u,u=e.to;
ok=1;break;
}
}
if(ok) continue;
int mindis=n-1;
for(i=g[u];i!=-1;i=edg[i].next){
edge& e=edg[i];
if(mindis>dis[e.to] && e.flow>0)
mindis=dis[e.to],cur[u]=i;
}
if(--gap[dis[u]]==0) break;
gap[dis[u]=mindis+1]++,u=pre[u];
}
return res;
}
}p;
int tmp[M][111],sz[M];
bool vis[111][111];
void run(){
int i,j;
int