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;
}