BZOJ 1221 [HNOI2001] 软件开发(费用流)

【题目链接】 http://www.lydsy.com/JudgeOnline/problem.php?id=1221

【题目大意】

  每天对毛巾都有一定的需求ni,每天可以花f价值每条购买毛巾,
  当天用完的毛巾可以花费fA价值每条通过快消毒在A天之后得到一条可用的,
  也可以通过花费fB价值每条,通过慢消毒在B天之后获得可用的
  问满足每天需求所用的最小花费。

【题解】

  首先,我们建立X集合,表示每天用完之后需要消毒的毛巾,显然第i个点值为ni,
  建立Y集合表示每天需要的毛巾数量,第i个点为ni,
  对于X中每个点,向A天后对应Y中的每个点连流量为INF,费用为fA的边
  同时向B题后对应Y中的每个点连流量为INF,费用为fB的边
  源点向X中第i个点连ni流量0费用的边,Y中第i个点向汇点连ni流量0费用的边,
  对于购买新毛巾的操作,我们从源点向Y中每个点连INF流量f费用的边
  那么求该图的最小费用最大流就是答案。
  但是我们发现按照以上方式建图边数量非常的庞大。
  考虑建图优化,我们对于X中每个点i向i+1连边,
  而对于X向Y的连边,我们只要连到有效区间的起点即可,
  这样就等价于X中每个点,向A(B)天后对应Y中的每个点连流量为INF,费用为fA(fB)的边

  对于这道题,有一种效率更高的做法,链接:三分+贪心

【代码】

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int INF=0x7fffffff,N=2010;
int S,T,cnt,ans,d[N],q[N],from[N],g[N],flow;
bool in[N];
struct edge{int from,to,nxt,c,v;}e[100010];
void add(int u,int v,int w,int c){
e[++cnt].from=u;e[cnt].to=v;
e[cnt].nxt=g[u];g[u]=cnt;
e[cnt].c=c;e[cnt].v=w;
}void add_edge(int u,int v,int w,int c){add(u,v,w,c);add(v,u,0,-c);}
bool spfa(){
for(int i=S;i<=T;i++)d[i]=INF;
int t=0,w=1;d[S]=0;in[S]=1;q[0]=S;
while(t!=w){
int now=q[t];t++;if(t==T)t=0;
for(int i=g[now];i;i=e[i].nxt)
if(e[i].v&&d[e[i].to]>d[now]+e[i].c){
d[e[i].to]=d[now]+e[i].c;from[e[i].to]=i;
if(!in[e[i].to]){in[e[i].to]=1;q[w++]=e[i].to;if(w==T)w=0;}
}in[now]=0;
}return(d[T]!=INF);
}
void mcf(){
int x=INF;
for(int i=from[T];i;i=from[e[i].from])x=min(x,e[i].v);flow+=x;
for(int i=from[T];i;i=from[e[i].from]){e[i].v-=x;e[i^1].v+=x;ans+=e[i].c*x;}
}
int n,A,B,f,fA,fB;
void solve(){
S=0; T=2*n+1;
memset(g,0,sizeof(g));
memset(e,0,sizeof(e));
ans=flow=0; cnt=1;
for(int i=1;i<=n;i++){
int x;
scanf("%d",&x);
add_edge(S,i,x,0);
add_edge(i+n,T,x,0);
}for(int i=1;i<n;i++)add_edge(i,i+1,INF,0);
for(int i=1;i<=n;i++)add_edge(S,i+n,INF,f);
for(int i=1;i+A+1<=n;i++)add_edge(i,i+A+1+n,INF,fA);
for(int i=1;i+B+1<=n;i++)add_edge(i,i+B+1+n,INF,fB);
while(spfa())mcf();
printf("%d\n",ans);
}
int main(){
while(~scanf("%d%d%d%d%d%d",&n,&A,&B,&f,&fA,&fB))solve();
return 0;
}
上一篇:npm yarn安装完成后,查不到版本号


下一篇:【BZOJ 1221】 1221: [HNOI2001] 软件开发 (最小费用流)