【CF331E】Biologist(网络流,最小割)
题面
翻译:
有一个长度为\(n\)的\(01\)串,将第\(i\)个位置变为另外一个数字的代价是\(v_i\)。
有\(m\)个要求
每个要求的形式是
首先确定若干位置都要是\(0\)或者\(1\)
然后给定这\(K\)个位置,如果些位置上都满足要求
那么就可以得到\(W_k\)元
某些要求如果失败了还要倒着给\(g\)元
问最终能够得到的最大利润
输入格式:
第一行是\(n,m,g\)
第二行是\(V_i\)
接下来\(m\)行
第一个数字表示这个集合都要是\(0\)还是\(1\)
第二个数字\(W_i\)表示利润,接下来\(k_i\)表示这个集合中有\(k\)个位置
接下来是这\(k\)个位置,
最后还有一个\(0/1\),如果是\(1\),表示如果失败了还要倒着给\(g\)元。
题解
不是很难的最小割
首先很明显的,源点表示\(0\),汇点表示\(1\)
与自己\(01\)相同的源或者汇连容量为\(0\)的边
然后往另外一侧连容量为\(v_i\)的边
先假设所有的要求都能够拿到利润
老套路的变成了计算最少的损失
如果要倒着给\(G\)元的不过是把这个任务失败的损失变成\(W_i+G\)
然后考虑怎么强行选任务
因为一个任务和其他的点之间不能断开
因此任务和其他的点用\(INF\)连接
然后这个任务要求是\(0/1\)就向对应的源/汇连容量为损失的边就行了
总的来说不是很难???
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 22222
#define INF 1000000000
inline int read()
{
RG int x=0,t=1;RG char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
struct Line{int v,next,w;}e[3333333];
int h[MAX],cnt=2;
inline void Add(int u,int v,int w)
{
e[cnt]=(Line){v,h[u],w};h[u]=cnt++;
e[cnt]=(Line){u,h[v],0};h[v]=cnt++;
}
int level[MAX],S,T,cur[MAX];
queue<int> Q;
bool bfs()
{
memset(level,0,sizeof(level));level[S]=1;
while(!Q.empty())Q.pop();Q.push(S);
while(!Q.empty())
{
int u=Q.front();Q.pop();
for(RG int i=h[u];i;i=e[i].next)
if(e[i].w&&!level[e[i].v])
{
level[e[i].v]=level[u]+1,Q.push(e[i].v);
if(e[i].v==T)return true;
}
}
return level[T];
}
int dfs(int u,int flow)
{
if(u==T||!flow)return flow;
int ret=0,used=0;
for(RG int &i=cur[u];i;i=e[i].next)
if(e[i].w&&level[e[i].v]==level[u]+1)
{
int d=dfs(e[i].v,min(flow-used,e[i].w));
used+=d;ret+=d;e[i].w-=d;e[i^1].w+=d;
if(used==flow)return ret;
}
if(!ret)level[u]=0;
return ret;
}
int Dinic()
{
RG int ret=0;
while(bfs())
{
for(RG int i=S;i<=T;++i)cur[i]=h[i];
while(int res=dfs(S,INF))ret+=res;
}
return ret;
}
int a[MAX],n,m,G;
int St[MAX],top=0;
int Ans;
int main()
{
n=read();m=read();G=read();
for(int i=1;i<=n;++i)a[i]=read();
S=0;T=n+m+1;
for(int i=1;i<=n;++i)
{
int x=read();
if(a[i])Add(i,T,0),Add(S,i,x);
else Add(i,T,x),Add(S,i,0);
}
for(int i=1;i<=m;++i)
{
int v=read(),W=read();
Ans+=W;top=read();
for(int j=1;j<=top;++j)
{
int x=read();
v?Add(i+n,x,INF):Add(x,i+n,INF);
}
W+=read()*G;
v?Add(S,i+n,W):Add(i+n,T,W);
}
printf("%d\n",Ans-Dinic());
return 0;
}