zoukankan      html  css  js  c++  java
  • 洛谷P4174 [NOI2006]最大获利(最大流)

    题目描述

    新的技术正冲击着手机通讯市场,对于各大运营商来说,这既是机遇,更是挑战。THU 集团旗下的 CS&T 通讯公司在新一代通讯技术血战的前夜,需要做太多的准备工作,仅就站址选择一项,就需要完成前期市场研究、站址勘测、最优化等项目。

    在前期市场调查和站址勘测之后,公司得到了一共 N 个可以作为通讯信号中转站的地址,而由于这些地址的地理位置差异,在不同的地方建造通讯中转站需要投入的成本也是不一样的,所幸在前期调查之后这些都是已知数据:建立第 i个通讯中转站需要的成本为 PiP_iPi (1≤i≤N)。

    另外公司调查得出了所有期望中的用户群,一共 M 个。关于第 i 个用户群的信息概括为 AiA_iAi , BiB_iBiCiC_iCi :这些用户会使用中转站 A i 和中转站 B i 进行通讯,公司可以获益 CiC_iCi 。(1≤i≤M, 1≤AiA_iAi , BiB_iBi ≤N)

    THU 集团的 CS&T 公司可以有选择的建立一些中转站(投入成本),为一些用户提供服务并获得收益(获益之和)。那么如何选择最终建立的中转站才能让公司的净获利最大呢?(净获利 = 获益之和 – 投入成本之和)

    输入输出格式

    输入格式:

    输入文件中第一行有两个正整数 N 和 M 。

    第二行中有 N 个整数描述每一个通讯中转站的建立成本,依次为 P1,P2,…,PNP_1 , P_2 , …,P_NP1,P2,,PN

    以下 M 行,第(i + 2)行的三个数 Ai,BiA_i , B_iAi,BiCiC_iCi 描述第 i 个用户群的信息。

    所有变量的含义可以参见题目描述。

    输出格式:

    你的程序只要向输出文件输出一个整数,表示公司可以得到的最大净获利。

    输入输出样例

    输入样例#1: 复制
    5 5
    1 2 3 4 5
    1 2 3
    2 3 4
    1 3 3
    1 4 2
    4 5 3
    输出样例#1: 复制
    4

    说明

    样例:选择建立 1、2、3 号中转站,则需要投入成本 6,获利为 10,因此得到最大收益 4。

    题解:这据说是一道最大权闭合子图...

    但是我根本不知道这是什么东西qwq

    大概的思路是先将所有客户的贡献加起来,然后此时我们要减去的有两个东西:届不到的用户和届到的中转站

    根据这类题的思路

    将源点向中转站建cost的边,中转站向用户建inf的边,用户向汇点建cost的边

    接着跑一遍dinic,得到的就是要减掉的数

    话说当前弧优化跑的真是飞快

    代码如下:

    #pragma GCC optimize(3,"inline")
    #include<queue>
    #include<cstdio>
    #include<cstring>
    #include<iostream>
    #include<algorithm>
    #define N 500050
    #define inf 0x3f3f3f3f
    using namespace std;
    
    int v[N],w[N],pos[N],nxt[N],cnt=-1,deep[55555],cur[N];
    int n,m,tmp;
    
    int addedge(int from,int to,int cost)
    {
        v[++cnt]=to;
        w[cnt]=cost;
        nxt[cnt]=pos[from];
        pos[from]=cnt;
    }
    
    int add(int from,int to,int cost)
    {
        addedge(from,to,cost);
        addedge(to,from,0);
    }
    
    int bfs(int s,int t)
    {
        queue<int> q;
        memset(deep,0,sizeof(deep));
        deep[s]=1;
        q.push(s);
        while(!q.empty())
        {
            int u=q.front();
            q.pop();
            for(int i=pos[u];~i;i=nxt[i])
            {
                if(w[i]>0&&(!deep[v[i]]))
                {
                    deep[v[i]]=deep[u]+1;
                    q.push(v[i]);
                }
            }
        }
        return deep[t]!=0;
    }
    
    int dfs(int s,int t,int dist)
    {
        if(s==t) return dist;
        for(int &i=cur[s];~i;i=nxt[i])
        {
            if(w[i]!=0&&deep[v[i]]==deep[s]+1)
            {
                int di=dfs(v[i],t,min(dist,w[i]));
                if(di>0)
                {
                    w[i]-=di;
                    w[i^1]+=di;
                    return di;
                }
            }
        }
        return 0;
    }
    
    int dinic(int s,int t)
    {
        int ans=0;
        while(bfs(s,t))
        {
            for(int i=0;i<=n+m+1;i++)
            {
                cur[i]=pos[i];
            }
            while(int di=dfs(s,t,inf))
            {
                ans+=di;
            }
        }
        return ans;
    }
    
    long long anss;
    
    int main()
    {
        memset(nxt,-1,sizeof(nxt));
        memset(pos,-1,sizeof(pos));
        scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++)
        {
            scanf("%d",&tmp);
            add(0,i,tmp);
        }
        int to1,to2,cost;
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d%d",&to1,&to2,&cost);
            add(to1,i+n,inf);
            add(to2,i+n,inf);
            add(i+n,n+m+1,cost);
            anss+=cost;
        }
        anss-=dinic(0,n+m+1);
        printf("%lld
    ",anss);
    }

    100%的数据中:N≤5 000,M≤50 000,0≤CiC_iCi ≤100,0≤PiP_iPi ≤100。

  • 相关阅读:
    Django测试开发-4-django3.0.3报错:/mysql ImproperlyConfigured: mysqlclient 1.3.13 or newer is required; you have 0.9.3
    Django测试开发-3-新建一个Django工程
    Django测试开发-2-创建虚拟环境
    Django测试开发-1-MVC/MVT的概念
    sqlserver中复合索引和include索引到底有多大区别?
    PLSQL开发笔记和小结(转载)
    PL/sql语法单元
    jquery 获取下拉框值与select text
    jquery 事件
    jquery dom操作
  • 原文地址:https://www.cnblogs.com/stxy-ferryman/p/9600596.html
Copyright © 2011-2022 走看看