题目链接:ヾ(≧∇≦*)ゝ

大致题意:

给出一个有向图D=(V,E).对于每个点U,定义两种操作a(u),b(u)

操作a(u):删除点U的所有出边,即属于E,操作花费为Ca(u).

操作b(u):删除点U的所有入边,即属于E,操作花费为Cb(u).

求将原图的边集的边全部删除的最小代价,总操作数和具体操作

Solution:

第一问很简单,首先,对于每一个点,把它分成出点和入点。

把每个点的出点与S相连,入点与T相连。边容量分别为删除该点所有入边和出边的花费。

然后对于每条边 a -> b,就把a的出点与b的入点连一条容量为inf的边。

根据最大流=最小割,跑一遍dinic就能得到答案了。

对于第二、三问,我们分别统计a操作和b操作。

我们先对剩余网络进行bfs(),把能够扫到的点都标记为1,不能的标记为0。

对于一个点u,如果要使用a(u),那么显然,需要至少存在一个点v,满足u -> v &&

vis[u]vis[v]0。

而对于点u,如果要使用b(u),只需要满足vis[u]==1就行了。

为了防止重复输出,在定义一个apr数组记录每个数是否加入到答案中就行了。

详见代码

Code:

#include<queue>
#include<ctype.h>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 1001
#define M 20001
#define inf 1926081700
using namespace std;
int S,T,head[N];
int n,m,cnt=1;
int ru[N],cu[N];
int ans,vis[N],apr[N];
int t1,t2,fst[N],sec[N];
struct Edge{int nxt,to,val;}edge[M];
void ins(int x,int y,int z){
edge[++cnt].nxt=head[x];
edge[cnt].to=y;edge[cnt].val=z;
head[x]=cnt;
}
namespace Network_Flow{
queue<int> q;
int dep[N];
int bfs(){
memset(dep,0,sizeof(dep));
q.push(S);dep[S]=1;
while(!q.empty()){
int x=q.front();q.pop();
for(int i=head[x];i;i=edge[i].nxt){
int y=edge[i].to,v=edge[i].val;
if(!dep[y]&&v){
q.push(y);
dep[y]=dep[x]+1;
}
}
}
return dep[T];
}
int dfs(int x,int rest){
if(x==T||rest<=0) return rest;
int flow=0;
for(int i=head[x];i;i=edge[i].nxt){
int y=edge[i].to,v=edge[i].val;
if(dep[y]==dep[x]+1&&v){
int now=dfs(y,min(rest,v));
edge[i].val-=now;
edge[i^1].val+=now;
flow+=now;rest-=now;
if(!rest) break;
}
}
return flow;
}
int dinic(){
int maxflow=0;
while(bfs()) maxflow+=dfs(S,inf);
return maxflow;
}
}
void getspj(){
queue<int> s;
s.push(S);vis[S]=1;
while(!s.empty()){
int x=s.front();s.pop();
for(int i=head[x];i;i=edge[i].nxt)
if(!vis[edge[i].to]&&edge[i].val){
s.push(edge[i].to);
vis[edge[i].to]=1;
}
}
apr[S]=apr[T]=1;
}
int read(){
int x=0,f=1;char ch=getchar();
while(!isdigit(ch)){if(ch=='-')f=-f;ch=getchar();}
while(isdigit(ch)){x=x*10+ch-48;ch=getchar();}
return x*f;
}
int main(){
n=read(),m=read();
S=n*2+1,T=S+1;
for(int i=1;i<=n;i++) ru[i]=read();
for(int i=1;i<=n;i++) cu[i]=read();
for(int i=1;i<=n;i++){
ins(S,i,cu[i]);ins(i,S,0);
ins(i+n,T,ru[i]);ins(T,i+n,0);
}
for(int x,y,i=1;i<=m;i++){
x=read(),y=read();
ins(x,n+y,inf);
ins(n+y,x,0);
}
using namespace Network_Flow;
printf("%d\n",dinic());getspj();
for(int i=1;i<=n;i++)
for(int j=head[i];j;j=edge[j].nxt){
int y=edge[j].to;
if(!vis[i]&&!vis[y]&&!apr[i]){
sec[++t2]=i;
ans++;apr[i]=1;
}
if(!apr[y]&&vis[y]){
fst[++t1]=y%n;
if(!fst[t1]) fst[t1]=n;
ans++;apr[y]=1;
}
}
printf("%d\n",ans);
for(int i=1;i<=t1;i++) printf("%d +\n",fst[i]);
for(int i=1;i<=t2;i++) printf("%d -\n",sec[i]);
return 0;
}

POJ2125 Destroying The Graph的更多相关文章

  1. POJ2125 Destroying The Graph (最小点权覆盖集)(网络流最小割)

                                                          Destroying The Graph Time Limit: 2000MS   Memo ...

  2. POJ2125 Destroying The Graph(二分图最小点权覆盖集)

    最小点权覆盖就是,对于有点权的有向图,选出权值和最少的点的集合覆盖所有的边. 解二分图最小点权覆盖集可以用最小割: vs-X-Y-vt这样连边,vs和X部点的连边容量为X部点的权值,Y部和vt连边容量 ...

  3. POJ2125 Destroying The Graph 二分图 + 最小点权覆盖 + 最小割

    思路来源:http://blog.csdn.net/lenleaves/article/details/7873441 求最小点权覆盖,同样求一个最小割,但是要求出割去了那些边, 只要用最终的剩余网络 ...

  4. POJ 2125 Destroying the Graph 二分图最小点权覆盖

    Destroying The Graph Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 8198   Accepted: 2 ...

  5. 【POJ】【2125】Destroying the Graph

    网络流/二分图最小点权覆盖 果然还是应该先看下胡伯涛的论文…… orz proverbs 题意: N个点M条边的有向图,给出如下两种操作.删除点i的所有出边,代价是Ai.删除点j的所有入边,代价是Bj ...

  6. 图论(网络流,二分图最小点权覆盖):POJ 2125 Destroying The Graph

    Destroying The Graph   Description Alice and Bob play the following game. First, Alice draws some di ...

  7. POJ 2125 Destroying The Graph [最小割 打印方案]

    Destroying The Graph Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 8311   Accepted: 2 ...

  8. poj 2125 Destroying The Graph (最小点权覆盖)

    Destroying The Graph http://poj.org/problem?id=2125 Time Limit: 2000MS   Memory Limit: 65536K       ...

  9. AC日记——Destroying The Graph poj 2125

    Destroying The Graph Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 8356   Accepted: 2 ...

随机推荐

  1. RTSP server 在mips 上莫名其妙退出(PC上则无此问题)

    http://blog.csdn.net/lubing20044793/article/details/38523701 早在这篇blog以前写过,在虚拟机下调试sn9c291时,USB 数据传输出了 ...

  2. sql语句之随机查询记录和批量插入--基于mysql

    本周遇到了好几次数据库方面的问题,一个是上一篇文章提到的因为要修改数据结构引起的在表间复制字段的需求,另一个就是这篇文章要写的:1,从某个数据表中按照某个字段不重复的随机选取几百条记录:2,然后把这些 ...

  3. 20155331《网路对抗》Exp8 WEB基础实践

    20155331<网路对抗>Exp8 WEB基础实践 基础问题回答 什么是表单 表单在网页中主要负责数据采集功能.一个表单有三个基本组成部分: 表单标签,这里面包含了处理表单数据所用CGI ...

  4. 牛客练习赛44 B题 (思维)

    链接:https://ac.nowcoder.com/acm/contest/634/B 来源:牛客网 给出n条线段,第i条线段的长度为ai, 每次可以从第i条线段的j位置跳到第i + 1条线段的j+ ...

  5. WPF编程,通过Double Animation动态缩放控件的一种方法。

    原文:WPF编程,通过Double Animation动态缩放控件的一种方法. 版权声明:我不生产代码,我只是代码的搬运工. https://blog.csdn.net/qq_43307934/art ...

  6. MFC Edit控件的使用~~

    EditBox,一般用于显示数字文本,或者与用户沟通获取数字文本. 这里介绍一种将EditBox与一个变量关联起来的方法: 快捷键:Shift + Ctrl + X,进入类导向,选择成员变量属性页: ...

  7. dxp altium pcb里面如果想让重叠的两个元件不报错怎么设置?

    dxp的设置是Design Rules里面有个Placement选项,把第一个的钩去掉即可.

  8. 使用 vi/vim 时,粘贴进新创建文件或空文件的首行内容丢失的解决方法

    只需要进入插入模式后,回车空一行或几行,再粘贴,再把上面的几个空行back回去,就不会丢失首行的内容了.

  9. Java开源博客My-Blog之docker容器组件化修改

    前言 5月13号上线了自己的个人博客,<Docker+SpringBoot+Mybatis+thymeleaf的Java博客系统开源啦>,紧接着也在github上开源了博客的代码,到现在为 ...

  10. Salesforce Apex学习 : 利用Schema命名空间中的DescribeSObjectResult类型来获取sObject对象的基本信息

    DescribeSObjectResult 对象的取得: //使用getDescribe方法和sObject token Schema.DescribeSObjectResult mySObjDesc ...