hdu3947 给一些已知(需费用)路径去覆盖一些边 //预先加灌法费用流

River Problem

题意:一个有向树(河流),只有一个汇点1,每条边只有一个出度。有些河道有污染指数xi,必需要治理,有m段路径,可以去覆盖这些,每被覆盖一次,xi降低响应值。

:即 给出一些边必需要覆盖的次数,用m段路径去覆盖,每次覆盖有相应费用,求最小费用。

思路:这题被誉为难题,给一个网络流,给出一些边的流量下界,以及给用某些路段流量去流满足要求。这里与正常网络流相悖,是wi>=xi,所以用:

预先加灌法:(预先灌入大流量(相应干道大小灌溉相应流量),要覆盖的减掉,使达不到预期流量,而用添加费用来补充,恰好跑最小费用最大流)

相应河道流量为U-wi(wi为要覆盖次数),再给的路径直接连边。如果没有wi,那么必然最大流为 U*主干道流量。应该有这些污染指数,所以在某些边流量变窄了!,再没有花费路径补充流量的情况下,比如达不到最大流,就无解,这样起到补充作用。

取之官方思路:

这里u取wi最大值略大(50)即可.

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
const int inf=0x3f3f3f3f;
const int maxv=301,maxe=10001;
int nume=0;int head[maxv];int e[maxe][4];
void inline adde(int i,int j,int c,int w)         //网络流图
{
    e[nume][0]=j;e[nume][1]=head[i];head[i]=nume;
    e[nume][2]=c;e[nume++][3]=w;
    e[nume][0]=i;e[nume][1]=head[j];head[j]=nume;
    e[nume][2]=0;e[nume++][3]=-w;
}
int nume2=0;int head2[maxv];int e2[maxe][4];
void inline adde2(int i,int j,int c,int w)              //原图
{
    e2[nume2][0]=j;e2[nume2][1]=head2[i];head2[i]=nume2;
    e2[nume2][2]=c;e2[nume2++][3]=w;
}
int numfa[maxv];                              //即dp[i]:i的边干道宽度/包含叶子数
int outd[maxv];
int n,m;int ss,tt;
int inq[maxv];int d[maxv]; int pre[maxv];int prv[maxv];
bool spfa(int &sum,int &sumflow)
{
    for(int i=0;i<=n;i++)
    {
        inq[i]=0;
        d[i]=inf;
    }
    queue<int>q;
    q.push(ss);
    inq[ss]=1;
    d[ss]=0;
    while(!q.empty())
    {
        int cur=q.front();
        q.pop();
        inq[cur]=0;
        for(int i=head[cur];i!=-1;i=e[i][1])
        {
            int v=e[i][0];
            if(e[i][2]>0&&d[cur]+e[i][3]<d[v])
            {
                d[v]=d[cur]+e[i][3];
                pre[v]=i;
                prv[v]=cur;
                if(!inq[v])
                {
                    q.push(v);
                    inq[v]=1;
                }
            }
        }
    }
    if(d[tt]==inf)return 0;
    int cur=tt;int minf=inf;
   while(cur!=ss)
    {
        int fe=pre[cur];
        minf=e[fe][2]<minf?e[fe][2]:minf;
        cur=prv[cur];
    }
     cur=tt;
    while(cur!=ss)
    {
        e[pre[cur]][2]-=minf;
        e[pre[cur]^1][2]+=minf;
        cur=prv[cur];
    }
    sum+=d[tt]*minf;
    sumflow+=minf;
    return 1;
}
int sumcost(int &sumflow)
{
    int sum=0;
    while(spfa(sum,sumflow));
    return sum;
}
int dfs_getnumfa(int u)                 // 记录出i所包含的叶子数(源点数)
{
    if(u==0)return 1;
    for(int i=head2[u];i!=-1;i=e2[i][1])
    {
           int v=e2[i][0];
             numfa[u]+=dfs_getnumfa(v);
    }
    return numfa[u];
}
void read_build()
{
       scanf("%d",&n);
       int aa,bb,cc,ww;
        for(int i=0;i<n-1;i++)
        {
            scanf("%d%d%d",&aa,&bb,&cc);
            adde2(bb,aa,0,cc);         //反向建树,原来的是倒的树
            outd[bb]++;
        }
         for(int i=1;i<=n;i++)
        {
            if(outd[i]==0)            //源点
            {
                adde2(i,0,0,0);
            }
        }
        numfa[0]=1;
        dfs_getnumfa(1);
      for(int i=0;i<=n;i++)
      {
          for(int j=head2[i];j!=-1;j=e2[j][1])
          {
              int v=e2[j][0];
              adde(v,i,numfa[v]*50-e2[j][3],0);
          }
      }
      scanf("%d",&m);
      for(int i=0;i<m;i++)
      {
          scanf("%d%d%d%d",&aa,&bb,&cc,&ww);
          adde(aa,bb,cc,ww);
      }
}
void init()
{
    nume=0;
    memset(head,-1,sizeof(head));
    nume2=0;
    memset(head2,-1,sizeof(head2));
    memset(numfa,0,sizeof(numfa));
    memset(outd,0,sizeof(outd));
    ss=0;tt=1;
}
int main()
{
    int T;
    scanf("%d",&T);
    for(int ii=1;ii<=T;ii++)
    {
        init();
        read_build();
        int sumflow=0;
        int ans=sumcost(sumflow);
        if(sumflow!=numfa[1]*50)
          printf("Case #%d: -1\n",ii);
        else
          printf("Case #%d: %d\n",ii,ans);
    }
    return 0;
}

hdu3947 给一些已知(需费用)路径去覆盖一些边 //预先加灌法费用流

时间: 2024-11-05 18:49:06

hdu3947 给一些已知(需费用)路径去覆盖一些边 //预先加灌法费用流的相关文章

Android百度地图 - 在地图上标注已知GPS纬度经度值的一个或一组覆盖物 - OPEN 开发经验库 - 360安全浏览器 8.1

首页   代码   文档   问答   资讯   经验   GitHub日报 登录   注册 www.open-open.com/libOPEN经验 投稿 全部经验分类  Android IOS JavaScript HTML5 CSS jQuery Python PHP NodeJS Java Spring MySQL MongoDB Redis NOSQL Vim C++ C# JSON Ruby Linux Nginx Docker 所有分类  >  开发语言与工具  >  移动开发  

[推荐]ORACLE PL/SQL编程之五:异常错误处理(知已知彼、百战不殆)

原文:[推荐]ORACLE PL/SQL编程之五:异常错误处理(知已知彼.百战不殆) [推荐]ORACLE PL/SQL编程之五: 异常错误处理(知已知彼.百战不殆) 继上三篇:ORACLE PL/SQL编程之八:把触发器说透 ORACLE PL/SQL编程之六:把过程与函数说透(穷追猛打,把根儿都拔起!) [推荐]ORACLE PL/SQL编程之四:把游标说透(不怕做不到,只怕想不到) 得到了大家的强力支持与建议,万分感谢.接下来介绍下一篇:oracle pl/sql异常处理部分,还望大家一定

【转】[推荐]ORACLE PL/SQL编程之五:异常错误处理(知已知彼、百战不殆)

[推荐]ORACLE PL/SQL编程之五: 异常错误处理(知已知彼.百战不殆) 继上三篇:ORACLE PL/SQL编程之八:把触发器说透 ORACLE PL/SQL编程之六:把过程与函数说透(穷追猛打,把根儿都拔起!) [推荐]ORACLE PL/SQL编程之四:把游标说透(不怕做不到,只怕想不到) 得到了大家的强力支持与建议,万分感谢.接下来介绍下一篇:oracle pl/sql异常处理部分,还望大家一定要支持与推荐呀~! 本篇主要内容如下: 5.1 异常处理概念 5.1.1 预定义的异常

Microsoft 修正了安全补丁 MS14-045 存在的已知问题

在2014年8月27号,微软修改了安全公告MS14-045 -内核模式驱动程序中的漏洞可能允许特权提升,以修正安装安全更新2982791所带来的一些已知问题. 概况 为了解决安全更新2982791存在的一些已知问题,微软对于所有受支持的Windows系统重新发布了安全公告MS14-045,使用更新2993651来取代存在问题的更新2982791.微软已经在2014年8月15日移除了更新2982791. 所有的客户都需安装更新2993651来取代已经过期的更新2982791. 如果想要获取更多关于

ARCgis已知线裁剪已知面

经常遇到需要在ArcGIS中,根据已知线图层(要素)切分已知面图层(要素).经过研究,利用topology拓扑菜单中的construct features可以实现.具体如下 现有用线图层A.面图层B,需使用A图层分割B图层,操作步骤如下: 1.打开Editor编辑工具条,单击Start Editing开始编辑. 2.调出topology拓扑菜单. 3.全选线层:右键单价A依次执行“选择/选择全部”. 4.点击拓扑工具条的工具构建要素(construct features). 5.在弹出的对话框中

实战分享:如何成功防护1.2T国内已知最大流量DDoS攻击

作者:腾讯云宙斯盾安全团队&腾讯安全平台部 引言: DDoS攻击势头愈演愈烈,除了攻击手法的多样化发展之外,最直接的还是攻击流量的成倍增长.3月份国内的最大规模DDoS攻击纪录还停留在数百G规模,4月,这个数据已经突破T级,未来不可期,我们唯有保持警惕之心,技术上稳打稳扎,以应对DDoS攻击卷起的血雨腥风.4月8日,腾讯云宙斯盾成功防御了1.2Tbps的超大流量攻击,也是目前国内已知的最大攻击流量,这篇文章就此次攻防事件简单地为大家做一个梳理和分析. 国内已知最大攻击流量来袭 4月8日,清明节后

已知起始点,获取每段等距离途经点的经纬度(用百度js api作)

已知两个中文地址,自动规划路径,获取路径上每个3公里的点的经纬度 <html> <head> <meta http-equiv="Content-Type" content="text/html; charset=utf-8" /> <meta name="viewport" content="initial-scale=1.0, user-scalable=no" /> <

如何过滤出已知当前目录下oldboy中的所有一级目录

第1章 linux启动过程 1.1 linux启动过程(1) 开机自检(BIOS)内存cup是否有问题(硬件检查)(2) MBR引导(3) GRUB菜单(选择不同的内核)(4) 加载内核(kernel)(5) 运行init进程(linux里面的第一个进程,初始化过程)(6) 读取/etc/inittab配置文件(运行级别配置文件)(7) 执行/etc/rc.d/rc.sysinit脚本(系统初始化脚本,设置主机名,设置IP地址)(8) 执行/etc/rc.d/rc脚本(根据系统的云心级别,再开机

Luogu-P1027 Car的旅行路线 已知三点确定矩形 + 最短路

传送门:https://www.luogu.org/problemnew/show/P1027 题意: 图中有n个城市,每个城市有4个机场在矩形的四个顶点上.一个城市间的机场可以通过高铁通达,不同城市间要通过飞机.现在问从s到t城市最少需要多少的费用. 思路: 已知矩形的三个顶点,可以用勾股定理确定斜边后,利用平行四边形原理——两对对角顶点的x之和是相同的,y之和也是相同的得到第四个顶点.然后用求最短路的dji即可. #include <algorithm> #include <iter