• 费用流


    https://www.luogu.org/problemnew/show/P3980

    源点连第一天

    汇点连最后一天

    容量为INF费用为0

    这样跑网络流是沿时间流的(就是依次解决每一天的问题)

    然后每一天向后一天连一条容量为INF-a[i]

    费用为0的边

    为什么容量为INF-a[i]

    这就相当于少了a[i]

    得用带权边补全INF

    这就是志愿者连续干时沿这条边跑

    因为连续干不花钱

    所以优先选这种边

    然后将每一类志愿者s[i]与t[i]+1连一条容量为

    INF花费为c[i]的边

    当连续干的人不够时

    就得 使劲往里塞人

    补全INF

    #include<bits/stdc++.h>
    using namespace std;
    const int M=2e3+3;
    const int inf=0x3f3f3f3f;
    struct node{
        int u,v,w,cost,nextt;
    }e[30000];
    int mincost;
    int head[M],cur[M],dis[M],vis[M],a[M],tot,s,t;
    void addedge(int u,int v,int w,int cost){
        e[tot].u=u;
        e[tot].v=v;
        e[tot].w=w;
        e[tot].cost=cost;
        e[tot].nextt=head[u];
        head[u]=tot++;
        e[tot].u=v;
        e[tot].v=u;
        e[tot].w=0;
        e[tot].cost=-cost;
        e[tot].nextt=head[v];
        head[v]=tot++;
    }
    bool bfs(){
        for(int i=0;i<=t;i++)
            dis[i]=inf;
        queue<int>que;
        que.push(s);
        dis[s]=0;
        while(!que.empty()){
            int u=que.front();
            que.pop();
            vis[u]=0;
            for(int i=head[u];~i;i=e[i].nextt){
                int v=e[i].v;
                if(e[i].w&&dis[u]+e[i].cost<dis[v]){
                    dis[v]=dis[u]+e[i].cost;
                    if(!vis[v]){
                        vis[v]=1;
                        que.push(v);
                    }
                }
            }
        }
        return dis[t]!=inf;
    }
    int dfs(int u,int fl){
        if(u==t)
            return fl;
        int ans=0;
        vis[u]=1;
        for(int i=cur[u];~i;i=e[i].nextt){
            int v=e[i].v;
            if(dis[v]==dis[u]+e[i].cost&&!vis[v]&&e[i].w){
                cur[u]=i;
                int x=dfs(v,min(e[i].w,fl-ans));
                e[i].w-=x;
                e[i^1].w+=x;
                ans+=x;
                mincost+=x*e[i].cost;
                if(ans==fl)
                    break;
            }
        }
        vis[u]=0;
        return ans;
    }
    void MCMF(){
        while(bfs()){
            for(int i=0;i<=t;i++)
                cur[i]=head[i];
            dfs(s,inf);
        }
    }
    int main(){
        int n,m;
        scanf("%d%d",&n,&m);
        memset(head,-1,sizeof(head));
        s=0,t=n+2;
        for(int i=1;i<=n;i++)
            scanf("%d",&a[i]),addedge(i,i+1,inf-a[i],0);
        addedge(0,1,inf,0);
        addedge(n+1,t,inf,0);
        while(m--){
            int l,r,cost;
            scanf("%d%d%d",&l,&r,&cost);
            addedge(l,r+1,inf ,cost);
        }
        MCMF();
        printf("%d",mincost);
        return 0;
    }
    View Code
  • 相关阅读:
    codeforces707B:Bakery
    codeforces707A:Brain's Photos
    BZOJ1084 [SCOI2005]最大子矩阵
    BZOJ1264 [AHOI2006]基因匹配Match
    BZOJ2764 [JLOI2011]基因补全
    codevs1257 打砖块
    BZOJ1079 [SCOI2008]着色方案
    BZOJ1026 [SCOI2009]windy数
    菜鸟学自动化测试(一)----selenium IDE
    关于w3school的html5部分output 元素实例代码(点亲自试一试进去)的问题纠正
  • 原文地址:https://www.cnblogs.com/starve/p/11177965.html
Copyright © 2020-2023  润新知