hdu 3879 hdu 3917 构造最大权闭合图 俩经典题

hdu3879  base station : 各一个无向图,点的权是负的,边的权是正的。自己建一个子图,使得获利最大。

一看,就感觉按最大密度子图的构想:选了边那么连接的俩端点必需选,于是就以边做点,轻轻松松构造了最大权闭合图。简单题。分分钟搞定。

hdu3917 :road  constructions :这题题目看了半天没理解。。。感觉描述的不好。。。一个有向图,每条路有响应公司承保,若选了该公司,那么该公司的路必需全部选,还有,该公司的承保的路的下面的一条路对应公司也要选,求最大获利。构图:开始在原图上瞎折腾,其实不然,用公司做结点,在原图上跑出关系,(若a公司选,则对应哪些公司 选),构造新图:最大权闭合之即可。。。。中午没睡醒。。迷迷糊糊的,因为1A的,数组开小了。。最讨厌用俩套链前来构图了。晕。。

再重说一次最大权闭合图解法:源点向正权点连边对应点权值,原图边容量inf,负权点向汇点连边权值绝对值,ans=正权和-最大流。

#include<iostream> //hdu3879
#include<queue>
#include<cstdio>
#include<cstring>
#include<set>
#include<vector>
using namespace std;
const int inf=0x3f3f3f3f;
const int maxv=61000,maxe=300000;
int nume=0;int head[maxv];int e[maxe][3];
void inline adde(int i,int j,int c)
{
    e[nume][0]=j;e[nume][1]=head[i];head[i]=nume;
    e[nume++][2]=c;
    e[nume][0]=i;e[nume][1]=head[j];head[j]=nume;
    e[nume++][2]=0;
}
int ss,tt,n,m;
int vis[maxv];int lev[maxv];
bool bfs()
{
    for(int i=0;i<maxv;i++)
      vis[i]=lev[i]=0;
    queue<int>q;
    q.push(ss);
    vis[ss]=1;
    while(!q.empty())
    {
        int cur=q.front();
        q.pop();
        for(int i=head[cur];i!=-1;i=e[i][1])
        {
            int v=e[i][0];
            if(!vis[v]&&e[i][2]>0)
            {
                lev[v]=lev[cur]+1;
                vis[v]=1;
                q.push(v);
            }
        }
    }
    return vis[tt];
}
int dfs(int u,int minf)
{
    if(u==tt||minf==0)return minf;
    int sumf=0,f;
    for(int i=head[u];i!=-1&&minf;i=e[i][1])
    {
        int v=e[i][0];
        if(lev[v]==lev[u]+1&&e[i][2]>0)
        {
            f=dfs(v,minf<e[i][2]?minf:e[i][2]);
            e[i][2]-=f;e[i^1][2]+=f;
            sumf+=f;minf-=f;
        }
    }
    if(!sumf) lev[u]=-1;
    return sumf;
}
int dinic()
{
    int sum=0;
    while(bfs())sum+=dfs(ss,inf);
    return sum;
};

int cpy[maxv];int sumz=0;
void read_build()
{
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&cpy[i]);
        adde(i,tt,cpy[i]);
    }
    int aa,bb,cc;
    int numv=n+1;
    for(int i=0;i<m;i++)
     {
        scanf("%d%d%d",&aa,&bb,&cc);
        sumz+=cc;
        adde(numv,aa,inf);
        adde(numv,bb,inf);
        adde(ss,numv++,cc);
     }
  /*  for(int i=1;i<=m+2;i++)
      for(int j=head[i];j!=-1;j=e[j][1])
      {
           if(j%2==0)
          printf("%d->%d:%d\n",i,e[j][0],e[j][2]);
      }*/
}
void init()
{
    nume=0;sumz=0;
    ss=n+m+1;tt=ss+1;
    for(int i=0;i<maxv;i++)
      {
       head[i]=-1;
      }
}
int main()
{
      while(~scanf("%d%d",&n,&m))
     {
         init();
        read_build();
       int ans;
       ans=sumz-dinic();
      printf("%d\n",ans);
    }
    return 0;
}

#include<iostream>//3917
#include<queue>
#include<cstdio>
#include<cstring>
#include<set>
#include<vector>
using namespace std;
const int inf=0x3f3f3f3f;
const int maxv=5100,maxe=100000;
int nume=0;int head[maxv];int e[maxe][3];
void inline adde(int i,int j,int c)
{
    e[nume][0]=j;e[nume][1]=head[i];head[i]=nume;
    e[nume++][2]=c;
    e[nume][0]=i;e[nume][1]=head[j];head[j]=nume;
    e[nume++][2]=0;
}
int ss,tt,n,m,k;
int vis[maxv];int lev[maxv];
bool bfs()
{
    for(int i=0;i<maxv;i++)
      vis[i]=lev[i]=0;
    queue<int>q;
    q.push(ss);
    vis[ss]=1;
    while(!q.empty())
    {
        int cur=q.front();
        q.pop();
        for(int i=head[cur];i!=-1;i=e[i][1])
        {
            int v=e[i][0];
            if(!vis[v]&&e[i][2]>0)
            {
                lev[v]=lev[cur]+1;
                vis[v]=1;
                q.push(v);
            }
        }
    }
    return vis[tt];
}
int dfs(int u,int minf)
{
    if(u==tt||minf==0)return minf;
    int sumf=0,f;
    for(int i=head[u];i!=-1&&minf;i=e[i][1])
    {
        int v=e[i][0];
        if(lev[v]==lev[u]+1&&e[i][2]>0)
        {
            f=dfs(v,minf<e[i][2]?minf:e[i][2]);
            e[i][2]-=f;e[i^1][2]+=f;
            sumf+=f;minf-=f;
        }
    }
    if(!sumf) lev[u]=-1;
    return sumf;
}
int dinic()
{
    int sum=0;
    while(bfs())sum+=dfs(ss,inf);
    return sum;
};
int nume2=0;int head2[maxv];int e2[maxe][3];
void inline adde2(int i,int j,int c)
{
    e2[nume2][0]=j; e2[nume2][1]=head2[i]; head2[i]=nume2;
    e2[nume2++][2]=c;
}
int cpy[maxv];int sumz=0;
void read_build()
{
    for(int i=1;i<=m;i++)
    {
        scanf("%d",&cpy[i]);
    }
    scanf("%d",&k);
    int aa,bb,cc,dd;
    for(int i=0;i<k;i++)
     {
        scanf("%d%d%d%d",&aa,&bb,&cc,&dd);
        cpy[cc]-=dd;
        adde2(aa,bb,cc);
     }
    for(int i=1;i<=n;i++)
      for(int j=head2[i];j!=-1;j=e2[j][1])
      {
          for(int vj=head2[e2[j][0]];vj!=-1;vj=e2[vj][1])
          {
              if(e2[vj][2]!=e2[j][2])
                 adde(e2[j][2],e2[vj][2],inf);
          }
      }
    for(int i=1;i<=m;i++)
    {
        if(cpy[i]>0)
          {
              adde(ss,i,cpy[i]);
              sumz+=cpy[i];
          }
        else
          adde(i,tt,-cpy[i]);
    }
  /*  for(int i=1;i<=m+2;i++)
      for(int j=head[i];j!=-1;j=e[j][1])
      {
           if(j%2==0)
          printf("%d->%d:%d\n",i,e[j][0],e[j][2]);
      }*/
}
void init()
{
    nume=0;sumz=0;nume2=0;
    ss=m+1;tt=ss+1;
    for(int i=0;i<maxv;i++)
      {
         head2[i]=head[i]=-1;
      }
}
int main()
{
      while(~scanf("%d%d",&n,&m)&&(n||m))
     {
         init();
        read_build();
       int ans;
       ans=sumz-dinic();
       if(ans<0)ans=0;
      printf("%d\n",ans);
    }
    return 0;
}

hdu 3879 hdu 3917 构造最大权闭合图 俩经典题,布布扣,bubuko.com

时间: 2024-10-07 06:38:48

hdu 3879 hdu 3917 构造最大权闭合图 俩经典题的相关文章

hdu 2987最大权闭合图模板类型题

/* 最大权闭合图模板类型的题,考验对知识概念的理解. 题意:现在要辞退一部分员工,辞退每一个员工可以的到一部分利益(可以是负的),并且辞退员工,必须辞退他的下属,求最大利益和辞退的最小人数. 最大权闭合图模板类型. 求出最大权后沿着源点s,dfs到的点就为最小的人数. 证明/* 转载:利用一个经典的trick:多关键字 > 建图前,对所有b[i],执行变换b[i]=b[i]*10000-1,然后,会惊异地发现, > 此时最大流所对应的方案就是满足辞退最少人数的了. > 为什么?显然,变

hdu 3879 最大权闭合图(裸题)

/* 裸的最大权闭合图 解:参见胡波涛的<最小割模型在信息学竞赛中的应用> */ #include<stdio.h> #include<string.h> #include<stdlib.h> #include<queue> using namespace std; #define N 55100//刚开始开的是5100一直越界应该是n+m #define NN 510000 #define inf 0x3fffffff struct node

HDU 3061:Battle(最大权闭合图)

http://acm.hdu.edu.cn/showproblem.php?pid=3061 题意:中文题意. 思路:和上一题神似啊,比上一题还简单,重新看了遍论文让我对这个理解更加深了. 闭合图:如果某个点在图中的话,那么这个点的后继点全部都要在图中. 对应至题目,这里的必须攻占b以后才能攻占a,那么是a依赖于b.如果a在图中的话,那么b必定在图中(因为a是依赖于b的),所以是a连向b(而不是b连向a). 这里总结一下做最大权闭合图的套路:把权值为正的点与超级源点S相连,容量为该权值,把权值为

HDU 3061 Battle(最小割----最大权闭合图)

题目地址:HDU 3061 多校中遇到的最小割的最大权闭合模型花了一上午时间终于看懂啦. 最大权闭合图就是将一些互相有依赖关系的点转换成图,闭合图指的是在图中的每一个点的后继点都是在图内的. 还要明白简单割的概念,就是指所有的割边都与源点或汇点相连.然后让源点与正权点相连,汇点与负权点相连,权值均为其绝对值,有依赖关系的点连一条有向边,如果a必须在b的基础上,那么就连一条a->b的有向边,权值为INF.最后用所有正权值得和减去最小割的值就是答案. 具体证明可看胡伯涛大牛的国家队集训论文<最小割

hdu 4971 多校10最大权闭合图

/* 很明显的最大权闭合图题 */ #include<stdio.h> #include<string.h> #include<queue> using namespace std; #define N 2100 #define inf 0x3fffffff struct node { int u,v,w,next; }bian[N*N*20]; int head[N],yong,dis[N],work[N]; void init(){ yong=0; memset(h

hdu 4971/ 2014多校/最大权闭合图

题意:n个项目(每一个相应获得一定价值).m个技术问题(每一个须要支出一定价值),每一个项目必须要攻克若干个技术问题.技术难题之间有拓扑关系. 关键是建图.一看,第一感觉就是最大权闭合图,马上建好了图.不难:以项目为正权点,问题为负权点.有依赖关系的点边就可以. ps:这题题目有句话有问题,按例子的来!害我贡献一次WA..... #include<iostream> #include<queue> #include<cstdio> #include<cstring

HDU 3879 Base Station(最大权闭合子图)

经典例题,好像说可以转化成maxflow(n,n+m),暂时只可以勉强理解maxflow(n+m,n+m)的做法. 题意:输入n个点,m条边的无向图.点权为负,边权为正,输出最大权闭合子图的权值. (最大权闭合子图:图中各点的后继必然也在图中) 构图攻略:将边看做点(有的题边是没有权的,建模稍微有点不同), 对原本的边e[i](u,v,w)连3条边(S,n+i,w),(n+i,u,inf),(n+i,v,inf). 对原本的点v,连1条边(v,T,p[v]). 即正权点与源点连,负权点与汇点连.

poj Firing(最大权闭合图)

Firing 题目: 要解雇一些人,而解雇的这些人如果人跟他有上下级的关系,则跟他有关系的人也要一起解雇.每个人都会创造一定的价值,要求你求出在最大的获利下,解雇的人最小. 算法分析: 在这之前要知道一个定理: 最小割 = 最大流 一道最大权闭合图的裸题,而这又可以转换成最小割来求解.证明可以看2007年胡伯涛的论文则可以直接套出模板,没看过的最好去看一下,那里解释的清楚.这里我给出他的原文的一些构造方法. 增加源s汇t 源s连接原图的正权点,容量为相应点权 原图的负权点连接汇t,容量为相应点权

太空飞行计划问题-最大权闭合图

题目描述 W 教授正在为国家航天中心计划一系列的太空飞行.每次太空飞行可进行一系列商业性实验而获取利润.现已确定了一个可供选择的实验集合E={E1,E2,…,Em},和进行这些实验需要使用的全部仪器的集合I={I1,I2,…In}.实验Ej需要用到的仪器是I的子集RjÍI.配置仪器Ik的费用为ck美元.实验Ej的赞助商已同意为该实验结果支付pj美元.W教授的任务是找出一个有效算法,确定在一次太空飞行中要进行哪些实验并因此而配置哪些仪器才能使太空飞行的净收益最大.这里净收益是指进行实验所获得的全部