acwing-2325. 有向图破坏(最小割之最小权覆盖集)「建议收藏」

acwing-2325. 有向图破坏(最小割之最小权覆盖集)「建议收藏」爱丽丝和鲍勃正在玩以下游戏。首先,爱丽丝绘制一个 N 个点 M 条边的有向图。然后,鲍勃试图毁掉它。在每一步操作中,鲍勃都可以选取一个点,并将所有射入该点的边移除或者将所有从该点射出的边移除。已知,对于第 i 个点,将所有射入该点的边移除所需的花费为 W+i,将所有从该点射出的边移除所需的花费为 W−i。鲍勃需要将图中的所有边移除,并且还要使花费尽可能少。请帮助鲍勃计算最少花费。输入格式第一行包含 N 和 M。第二行包含 N 个正整数,第 i 个为 W+i。第三行包含 N 个正整数,第.

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

在这里插入图片描述
爱丽丝和鲍勃正在玩以下游戏。

首先,爱丽丝绘制一个 N 个点 M 条边的有向图。

然后,鲍勃试图毁掉它。

在每一步操作中,鲍勃都可以选取一个点,并将所有射入该点的边移除或者将所有从该点射出的边移除。

已知,对于第 i 个点,将所有射入该点的边移除所需的花费为 W+i,将所有从该点射出的边移除所需的花费为 W−i。

鲍勃需要将图中的所有边移除,并且还要使花费尽可能少。

请帮助鲍勃计算最少花费。

输入格式
第一行包含 N 和 M。

第二行包含 N 个正整数,第 i 个为 W+i。

第三行包含 N 个正整数,第 i 个为 W−i。

接下来 M 行,每行包含两个整数 a,b,表示从点 a 到点 b 存在一条有向边。

所有点编号从 1 到 N。

图中可能由重边或自环。

输出格式
第一行输出整数 W,表示鲍勃所需的最少花费。

第二行输出整数 K,表示鲍勃需要进行的操作步数。

接下来 K 行,每行输出一个鲍勃的具体操作。

如果操作为将所有射入点 i 的边移除,则输出格式为 i +。

如果操作为将所有从点 i 射出的边移除,则输出格式为 i -。

如果答案不唯一,则输出任意一种即可。

数据范围
1≤N≤100,
1≤M≤5000,
1≤W+i,W−i≤106

输入样例:
3 6
1 2 3
4 2 1
1 2
1 1
3 2
1 2
3 1
2 3
输出样例:
5
3
1 +
2 -
2 +
#include<bits/stdc++.h>
using namespace std;
const int N = 210;
const int M = 2 * (5000 + 200) + 10;
const int INF = 0x3f3f3f3f;
struct Edge{ 

int v,next,w;
}edge[M];
int head[N],cnt;
int n,m,s,e;
int st[N];
void add(int u,int v,int w){ 

edge[cnt].v = v;
edge[cnt].w = w;
edge[cnt].next = head[u];
head[u] = cnt ++;
}
int d[N],q[N],hh = 0,tt = 0,cur[N];
bool bfs(){ 

hh = tt = 0;
memset(d,-1,sizeof d);
d[s] = 0,q[tt ++] = s,cur[s] = head[s];
while(hh < tt){ 

int t = q[hh ++];
for(int i = head[t];~i;i = edge[i].next){ 

int v = edge[i].v,w = edge[i].w;
if(d[v] == -1 && w){ 

d[v] = d[t] + 1;
cur[v] = head[v];
q[tt ++] = v;
if(v == e)return true;
}
}
}
return false;
}
int dfs(int u,int limit){ 

if(u == e)return limit;
int flow = 0;
for(int i = cur[u];~i && flow < limit;i = edge[i].next){ 

int v = edge[i].v,w = edge[i].w;
cur[u] = i;
if(d[v] == d[u] + 1 && w){ 

int t = dfs(v,min(w,limit - flow));
if(!t)d[v] = -1;
flow += t,edge[i].w -= t,edge[i ^ 1].w += t; 
}
}
return flow;
}
int dinic(){ 

int maxflow = 0;
int flow = 0;
while(bfs())while(flow = dfs(s,INF))maxflow += flow;
return maxflow;
}
void dfs(int u){ 

st[u] = true;
for(int i = head[u];~i;i = edge[i].next){ 

int v = edge[i].v,w = edge[i].w;
if(w > 0  && !st[v])dfs(v);
}
}
int main(){ 

cin>>n>>m;
memset(head,-1,sizeof head);
s = 0,e = 2 * n + 1;
int x;
for(int i = 1;i <= n;i ++){ 

cin>>x;
add(s,i,x),add(i,s,0);
}
for(int i = 1;i <= n;i ++){ 

cin>>x;
add(n + i,e,x),add(e,n + i,0);
}
int a,b;
for(int i = 0;i < m;i ++){ 

cin>>a>>b;
add(b,a + n,INF),add(a + n,b,0);
}
cout<<dinic()<<endl;
dfs(s);
int tot = 0;
for(int i = 0;i < cnt;i += 2){ 

if(edge[i].w == 0 && !st[edge[i].v] && st[edge[i ^ 1].v])tot ++;  //割边一定是满流
}
cout<<tot<<endl;
for(int i = 0;i < cnt;i += 2){ 

if(!st[edge[i].v] && st[edge[i ^ 1].v]){ 

if(edge[i ^ 1].v == s){ 

cout<<edge[i].v<<" +"<<endl;
}
if(edge[i].v == e){ 

cout<<edge[i ^ 1].v - n<<" -"<<endl;
}
}
}
return 0;
}
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/169077.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)
blank

相关推荐

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号