bzoj1221: [HNOI2001] 软件开发

传送门
把每天分为二分图两个集合中的顶点Xi,Yi，建立附加源S汇T。
1、从S向每个Xi连一条容量为ri，费用为0的有向边。
2、从每个Yi向T连一条容量为ri，费用为0的有向边。
3、从S向每个Yi连一条容量为无穷大，费用为p的有向边。
4、从每个Xi向Xi+1(i+1<=N)连一条容量为无穷大，费用为0的有向边。
5、从每个Xi向Yi+m+1(i+m+1<=N)连一条容量为无穷大，费用为f的有向边。
6、从每个Xi向Yi+n+1(i+n+1<=N)连一条容量为无穷大，费用为s的有向边。
求网络最小费用最大流，费用流值就是要求的最小总花费。
感觉费用流好神奇

#include<iostream>#include<cstdio>#include<queue>#include<cstring>using namespace std;const int maxn=2005,inf=0x3f3f3f3f;int Begin[maxn],Next[maxn*10],to[maxn*10],c[maxn*10],w[maxn*10],id[maxn][maxn],e;void add(int u,int v,int p,int val){    to[++e]=v,c[e]=p,Next[e]=Begin[u],id[u][v]=Begin[u]=e,w[e]=val;    to[++e]=u,c[e]=0,Next[e]=Begin[v],id[v][u]=Begin[v]=e,w[e]=-val;}int d[maxn],pre[maxn];bool vis[maxn];queue<int>q;int spfa(int s,int t){    memset(d,inf,sizeof(d));    memset(pre,-1,sizeof(pre));    q.push(s),d[s]=0;    while(!q.empty()){        int u=q.front();        q.pop(),vis[u]=0;        for(int i=Begin[u];i;i=Next[i]){            int v=to[i];            if((d[v]>d[u]+w[i])&&(c[i])){                d[v]=d[u]+w[i];                pre[v]=u;                if(!vis[v]){                    vis[v]=1;                    q.push(v);                }            }        }    }    if(pre[t]==-1)return 0;    int u,v=t,Min=inf;    while((u=pre[v])!=-1){        if(Min>c[id[u][v]])            Min=c[id[u][v]];        v=u;    }    int ans=Min*d[v=t];    while((u=pre[v])!=-1){        c[id[u][v]]-=Min;        c[id[v][u]]+=Min;        v=u;    }    return ans;}int main(){    int n,A,B,fF,fA,fB,x;    scanf("%d%d%d%d%d%d",&n,&A,&B,&fF,&fA,&fB);    for(int i=1;i<=n;i++){        scanf("%d",&x);        add(0,i,x,0);        add(i+n,2*n+1,x,0);        add(0,i+n,inf,fF);        if(i<n)add(i,i+1,inf,0);        if(i+A+1<=n)add(i,i+A+1+n,inf,fA);        if(i+B+1<=n)add(i,i+B+1+n,inf,fB);    }    int res,ans=0;    while(res=spfa(0,2*n+1))ans+=res;    printf("%d\n",ans);    return 0;}