题解 P1361 【小M的作物】
ButterflyDew
2018-06-28 20:03:37
题目很棒。
从我个人的角度来说,可能有那么不太准确的一句话加深了我的理解。
最小割在数值上与最大流相等,但在本身性质上与网络流无任何关系。
它的不那么准确的定义是:**在一个图中,割去权值和最小的边集,使这个图分成两个部分,切下来的那一刀叫做最小割。最小割并不唯一,但它在数值上是等于最大流的**。
对于这个题的模型来说,它有一个名字叫做**二者取一式问题**,在具体题目中的体现即为把点集一分为二。
类似的题目有[善意的投票](https://www.luogu.org/problemnew/show/P2057)
虽然说了最小割与网络流无关系,但只是在定义或者说是意义上,真正要到建模时,还需要将两者综合考虑。
先考虑没有额外收益情况下的建图。
![](http://m.qpic.cn/psb?/V14VFGnz4fbJRr/hj1ajou0Rsm281kSiFj7z6*UXkaX0Uboq17fnL.ahB0!/b/dFYBAAAAAAAA&bo=pgJnAaYCZwEDCSw!&rf=viewer_4)
设源点S属于集合A,汇点T属于集合B,则显然点与源点所连的有向边为这个点归属于集合A所产生的收益,与汇点相连则同理。则原图边权之和-最小割就是本题的答案了。
在这个图的基础上,考虑如何把点集给加入。
首先明确一点,点集的贡献有三种情况,对集合A贡献,对集合B贡献或者不贡献。这意味着只划分出一种状态是无法描述的,至少要把A与B的情况分开描述。
讨论一个对$A$有贡献的点集$\{c,d\}$。
依据题目,我们对这个点集的要求是,只要$c,d$有一个点被割到了集合$T$,这个点集都无法产生贡献。换而言之,只要$c$或$d$在集合$B$,代表点集贡献的边必须要断开。
先尝试着连接这条贡献边,因为这条边不可能直接连到图中代表作物的点上,所以连接一个虚点上去。
![](http://m.qpic.cn/psb?/V14VFGnz4fbJRr/NNAGO3JptVlvEnKYZvkHHtQp6bXKTuZT9s.3xX2zgRM!/b/dDABAAAAAAAA&bo=fgJUAX4CVAEDGTw!&rf=viewer_4)
X为点集所产生的虚点。
如果$c$被割到了集合$B$,则所有从$S$到$c$的路径都得被断开(具体只路径的一条边断掉),如果我们想要黄边断掉,那么虚点$X$连$c$的那条路径不能被断掉。而容量为正无穷的边不可能被断掉,于是我们这样建图。
![](http://m.qpic.cn/psb?/V14VFGnz4fbJRr/ua*Uh62JYJMw246lyEHMygXlJLY2NhT3K5ZsGKq6NwQ!/b/dC0BAAAAAAAA&bo=fQJFAX0CRQEDGTw!&rf=viewer_4)
两条蓝色的边的边权为$inf$
对待贡献$B$集合的点,同理
![](http://m.qpic.cn/psb?/V14VFGnz4fbJRr/E8sYvzooansm*Ys3MhIvwFMJPcm4teMjJSPqoE8UE0U!/b/dDABAAAAAAAA&bo=cAJTAXACUwEDGTw!&rf=viewer_4)
建图跑网络流即可。
要注意的一点是,这题卡常,~~吸氧才水过去的~~
[露迭月喵~](http://www.cnblogs.com/ppprseter/)
------------
**Code**
```
#include <cstdio>
#include <cstring>
#include <queue>
using namespace std;
const int N=3010;
const int M=2000100;
const int inf=0x3f3f3f3f;
int head[N],edge[M],to[M],next[M],cnt=1;
void add(int u,int v,int w)
{
to[++cnt]=v;next[cnt]=head[u];edge[cnt]=w;head[u]=cnt;
to[++cnt]=u;next[cnt]=head[v];edge[cnt]=0;head[v]=cnt;
}
int dep[N],used[N],pre[N],tot,s[N],ans,m,n,sum;
queue <int > q;
bool bfs()
{
while(!q.empty()) q.pop();
q.push(0);
memset(dep,0,sizeof(dep));
dep[0]=1;
while(!q.empty()&&q.front()!=n+1)
{
int u=q.front();
q.pop();
for(int i=head[u];i;i=next[i])
{
int v=to[i],w=edge[i];
if(!dep[v]&&w)
{
dep[v]=dep[u]+1;
q.push(v);
}
}
}
return !q.empty();
}
int main()
{
scanf("%d",&n);
int w,v,k,c1,c2;
for(int i=1;i<=n;i++)
{
scanf("%d",&w);
sum+=w;
add(0,i,w);
}
for(int i=1;i<=n;i++)
{
scanf("%d",&w);
sum+=w;
add(i,n+1,w);
}
scanf("%d",&m);
for(int i=1;i<=m;i++)
{
scanf("%d%d%d",&k,&c1,&c2);
add(0,i+n+1,c1);sum+=c1;
add(i+n+m+1,n+1,c2);sum+=c2;
for(int j=1;j<=k;j++)
{
scanf("%d",&v);
add(i+n+1,v,inf);
add(v,i+n+m+1,inf);
}
}
while(bfs())
{
memset(used,0,sizeof(used));
s[++tot]=0;
while(tot)
{
int u=s[tot];
if(u==n+1)
{
int mi=inf,id;
for(int i=tot;i>1;i--)
if(mi>=edge[pre[s[i]]])
{
mi=edge[pre[s[i]]];
id=i;
}
ans+=mi;
for(int i=tot;i>1;i--)
{
edge[pre[s[i]]]-=mi;
edge[pre[s[i]]^1]+=mi;
}
tot=id-1;
used[n+1]=0;
}
else
{
for(int i=head[u];i;i=next[i])
{
int v=to[i],w=edge[i];
if(!used[v]&&dep[v]==dep[u]+1&&w)
{
used[v]=1;
s[++tot]=v;
pre[v]=i;
break;
}
}
if(u==s[tot]) tot--;
}
}
}
printf("%d\n",sum-ans);
return 0;
}
```