[HNOI2001] 软件开发

题目描述：

雾

题目分析：

最小费用最大流
原题解
拆点，每个点拆成入点和出点。

显然，要满足每天的餐巾供应，于是有附加源汇ST，S向入点引一条容量为ni，费用为0的边，出点向T引一条容量为ni，费用为0的边。

若满足题目要求，这必有所有出点到T的弧满载，为了保证这一点，由S向出点引容量为无穷，费用为f的边。

又因为每天没用完的餐巾可以留到下一天，所以每天的入点向下一天的入点引流量为无穷，费用为0的边。

又由题目的条件，对于餐巾消毒，第i天的入点向第i+a(b)天的出点引一条容量为无穷，费用为fa(fb)的边。

求ST得最小费用最大流即得答案。

题目链接：

Luogu 2223
BZOJ 1221

Ac 代码：

#include <cstdio>
#include <iostream>
#include <cstring>
#include <queue>
#define il inline
using namespace std;
const int inf=0x7fffffff;
const int maxm=110000;
int head[maxm],to[maxm*2],cap[maxm*2],net[maxm*2],cost[maxm*2];
int cnt=1;
il void add(int x,int y,int c,int z){cnt++,to[cnt]=y,cost[cnt]=z,cap[cnt]=c,net[cnt]=head[x],head[x]=cnt;}
int flow[maxm],pre[maxm],id[maxm],dis[maxm],max_flow,min_cost;
bool vis[maxm];
queue <int> dl;
struct node{
    int x,y,cap,w;
}a[maxm];
il bool BFS(int s,int t)
{
    while(!dl.empty()) dl.pop();
    memset(pre,-1,sizeof(pre));
    memset(dis,127/3,sizeof(dis));
    dis[s]=0,flow[s]=inf,pre[s]=0,vis[s]=1;
    dl.push(s);
    while(!dl.empty())
    {
        int x=dl.front();
        dl.pop();
        vis[x]=0;
        for(int i=head[x];i;i=net[i])
        {
           int tmp=to[i];
           if(cap[i]>0&&dis[tmp]>dis[x]+cost[i])
           {
              dis[tmp]=dis[x]+cost[i];
              pre[tmp]=x;
              id[tmp]=i;
              flow[tmp]=min(flow[x],cap[i]);
              if(!vis[tmp]) vis[tmp]=1,dl.push(tmp);
           }
        }
    }
    return pre[t]==-1?0:1;
}
il void change_cap(int s,int t,int x)
{
    int now=t;
    while(now!=s)
    {
        cap[id[now]]-=x,cap[id[now]^1]+=x;
        now=pre[now];
    }
}
void il get_ans(int s,int t)
{
    max_flow=0,min_cost=0;
    while(BFS(s,t))
    {
        //printf("%d\n",flow[t]);
        max_flow+=flow[t],min_cost+=dis[t]*flow[t];
        change_cap(s,t,flow[t]);
    }
}
il void adx(int x,int y,int cax,int c)
{
    add(x,y,cax,c),add(y,x,0,-c);
}
int main()
{
    int n,ya,yb,f,fa,fb;  
    scanf("%d%d%d%d%d%d",&n,&ya,&yb,&f,&fa,&fb);  
    int S=2*n+1,T=S+1;  
    for(int i=1;i<=n;i++)
    {  
        int w;
        scanf("%d",&w);  
        adx(S,i,w,0);  
        adx(i+n,T,w,0);  
        if(i+ya+1<=n) adx(i,i+ya+n+1,inf,fa);  
        if(i+yb+1<=n) adx(i,i+yb+n+1,inf,fb);  
        if(i<n) adx(i,i+1,inf,0);  
        adx(S,i+n,inf,f);  
    }  
    get_ans(S,T);
    printf("%d",min_cost);  
    return 0;
}