BZOJ1391: [Ceoi2008]order

[Ceoi2008]order

Description

有N个工作，M种机器，每种机器你可以租或者买过来. 每个工作包括若干道工序，每道工序需要某种机器来完成,你可以通过购买或租用机器来完成。 现在给出这些参数，求最大利润。

构图方法

若不考虑租借，则是裸最大权闭合子图。那么这题就需要考虑如何修改建图方式，使得“可以选x的时候可以不选y，但要付出一定代价”。

只需将中间相邻的边权为INF的边，容量改为租借费用即可。

①：S向每个工作连边，容量为获益值；

②：每个机器向T连边，容量为购买费用；

③：每个工作向所需机器连边，容量为租借费用。

PS：要加当前弧优化！

【代码】

#include <cstdio>
#include <iostream>
#include <queue>
#include <vector>
#include <algorithm>
#include <cstring>

4000
#include <cmath>
#include <stack>
#define N 2405
#define M 2884805
#define INF 1000000001
using namespace std;
typedef long long ll;
typedef pair<ll,ll> pa;

int read()
{
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-') f=-1;ch=getchar();}
while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}
return x*f;
}

int n,m,cnt=1,S,T,ans;
int b[M],p
,nextedge[M],w[M],cur
;
int Level
;

void Add(int x,int y,int z)
{
cnt++;
b[cnt]=y;
nextedge[cnt]=p[x];
p[x]=cnt;
w[cnt]=z;
}

void Anode(int x,int y,int z){
Add(x,y,z);Add(y,x,0);
}

void Input_Init()
{
T=n+m+1;
for(int i=1;i<=n;i++)
{
static int x,y,z,num;
x=read(),num=read();
Anode(0,i,x);ans+=x;
while(num--)
{
y=read(),z=read();
Anode(i,n+y,z);
}
}
for(int i=1;i<=m;i++)
{
static int x;
x=read();
Anode(n+i,T,x);
}
}

bool Bfs()
{
queue<int>q;
q.push(S);
for(int i=0;i<=T;i++) Level[i]=0;
Level[S]=1;
while(!q.empty())
{
int k=q.front();q.pop();
for(int i=p[k];i;i=nextedge[i])
{
int v=b[i],f=w[i];
if(!Level[v]&&f)
{
Level[v]=Level[k]+1;
q.push(v);
}
}
}
return Level[T];
}

int Dfs(int x,int maxf)
{
if(x==T||!maxf) return maxf;
int rtn=0;
for(int i=cur[x];i&&maxf>rtn;i=nextedge[i])
{
int v=b[i],f=w[i];
if(Level[v]==Level[x]+1&&f)
{
f=Dfs(v,min(maxf-rtn,f));
w[i]-=f;w[i^1]+=f;
if(w[i]>0) cur[x]=i;
rtn+=f;
}
}
if(!rtn) Level[x]=0;
return rtn;
}

void Dinic()
{
while(Bfs())
{
for(int i=0;i<=T;i++) cur[i]=p[i];
ans-=Dfs(S,INF);
}
printf("%d\n",ans);
}

int main()
{
n=read();m=read();
Input_Init();
Dinic();
return 0;
}