bzoj 1221: [HNOI2001] 软件开发 费用流

题意：

某软件公司正在规划一项n天的软件开发计划，根据开发计划第i天需要ni个软件开发人员，为了提高软件开发人员的效率，公司给软件人员提供了很多的服务，其中一项服务就是要为每个开发人员每天提供一块消毒毛巾，这种消毒毛巾使用一天后必须再做消毒处理后才能使用。消毒方式有两种，A种方式的消毒需要a天时间，B种方式的消毒需要b天（b>a），A种消毒方式的费用为每块毛巾fA, B种消毒方式的费用为每块毛巾fB，而买一块新毛巾的费用为f（新毛巾是已消毒的，当天可以使用）；而且f>fA>fB。公司经理正在规划在这n天中，每天买多少块新毛巾、每天送多少块毛巾进行A种消毒和每天送多少块毛巾进行B种消毒。当然，公司经理希望费用最低。你的任务就是：为该软件公司计划每天买多少块毛巾、每天多少块毛巾进行A种消毒和多少毛巾进行B种消毒，使公司在这项n天的软件开发中，提供毛巾服务的总费用最低。

题解：

费用流。
应该很容易看出来，是费用流。
每个点拆成两个点，消毒后得来的毛巾可以看作限量的低价毛巾。
于是大概就这样：

    ins(st,i,inf,f);
    ins(i,ed,A[i],0);
    ins(st,i+n,A[i],0);
    if(i+a+1<=n) ins(i+n,i+a+1,inf,fA);
    if(i+b+1<=n) ins(i+n,i+b+1,inf,fB);

然后很愉快的wa了。
orz，发现有个坑点，消毒完的毛巾可以放几天再用坑坑坑。
就随便向下连条边就行了。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#define LL long long
using namespace std;
struct node{
    LL x,y,c,d,next,other;
}a[310000];LL len=0,last[3100];
const LL inf=(1<<28);
LL st,ed;
LL s[3100],p[3100],q[3100];
bool u[3100];
void ins(LL x,LL y,LL c,LL d)
{
    LL k1=++len;
    a[len].x=x;a[len].y=y;a[len].c=c;a[len].d=d;
    a[len].next=last[x];last[x]=len;
    LL k2=++len;
    a[len].x=y;a[len].y=x;a[len].c=0;a[len].d=-d;
    a[len].next=last[y];last[y]=len;
    a[k1].other=k2;
    a[k2].other=k1;
}
bool spfa()
{
    memset(s,63,sizeof(s));
    memset(u,false,sizeof(u));
    LL l=1,r=2;q[l]=st;s[st]=0;u[st]=true;
    while(l!=r)
    {
        LL x=q[l];
        for(LL i=last[x];i;i=a[i].next)
        {
            LL y=a[i].y;
            if(s[y]>s[x]+a[i].d&&a[i].c>0)
            {
                s[y]=s[x]+a[i].d;
                p[y]=i;
                if(!u[y])
                {
                    u[y]=true;
                    q[r]=y;
                    r++;if(r>ed+1) r=1;
                }
            }
        }
        u[x]=false;
        l++;if(l>ed+1) l=1; 
    }
    return s[ed]<999999999;
}
LL flow() 
{
    LL x=ed;
    LL ans=0,min=-1;
    while(x!=st) 
    {
        LL i=p[x];
        if (a[i].c<min||min==-1)min=a[i].c;
        x=a[i].x;
    }
    x=ed;
    while(x!=st) 
    {
        LL i=p[x];
        a[i].c-=min;
        a[a[i].other].c+=min;
        x=a[i].x;
        ans=ans+a[i].d*min;
    }
    return ans;
}
LL n,aa,b,f,fA,fB,A[1100];
int main()
{
    scanf("%lld %lld %lld %lld %lld %lld",&n,&aa,&b,&f,&fA,&fB);
    for(LL i=1;i<=n;i++) scanf("%lld",&A[i]);
    st=0;ed=2*n+1;
    for(LL i=1;i<=n;i++)
    {
        ins(st,i,inf,f);
        ins(i,ed,A[i],0);
        ins(st,i+n,A[i],0);
        if(i+aa+1<=n) ins(i+n,i+aa+1,inf,fA);
        if(i+b+1<=n) ins(i+n,i+b+1,inf,fB);
        if(i!=n) ins(i,i+1,inf,0);
    }
    LL ans=0;
    while(spfa()) ans+=flow();
    printf("%lld",ans);
}