大致题意: 有(n)种药,每种药有一个权值,且使用了若干种药材。让你选择若干种药,使得药的数量与所使用的药材并集大小相等,求最小权值总和。
网络流
(hl666):这种数据范围,一眼网络流!
然后我们两个画了半天硬是想不出一种合适的建图方式。。。
于是就默默打开了题解。。。
可惜题解里讲得都不是很详细,蒟蒻我表示难以理解,最后想了很久才搞明白为什么要这样做,且这样做为什么对。
因此,我会讲得较为详细一些,把所有让我思考了很久的地方都解释清楚。
建图
考虑从源点向每种药分别连一条容量为(INF)减去其权值的边,然后从每种药向对应的药材,从每种药材向汇点,分别连一条流量为(INF)的边。
然后,我们求出最小割(即跑一遍最大流),并用其减去从源点向外连出的边的容量总和即为答案。
看完是不是一脸懵逼?为什么边权都要设成(INF)级别?为什么这样求出最小割就能求出答案?
让我们先来借助样例,来大致感受一下按照这样建图有什么神奇效果:
此时,原图中的最小割为割去(s->1,s->2,3'->t),共为(3*INF-10)。
而从源点向外连出的边的容量总和为(3*INF-7),因此相减得到(-3),即答案。
是不是感到很神奇?
关于建图的具体解释
接下来,我们来具体解释一下这样建图的好处及作用,也就是这样建图的一些性质。
连接药与药材的边绝不会被选作最小割中的边
考虑到你一个药会向不少于(1)种药材连边,所以连接药与药材的边数显然是大于或等于药材数的。
则你设想如果一条连接药与药材的边被选作了最小割中的边,则显然,割去它没有割去它所连接的药材与汇点的连边更优。
原因是这条边只不过是连向该药材的一条边,你割去了这条边,其他边可能依然能流通。而你割去了该药材与汇点的连边,则所有连向该药材的边都断流了。
退一步说,即使没有其他连向该药材的边,考虑到这条边与该药材连向汇点的边流量相同,因此割去这两条边就是等价的,因此依然可以看作是割去了该药材与汇点的连边。
该性质得证。
则我们可由此推知,被选作最小割中的边一定是从源点向某种药的连边或从某种药材向汇点的连边。
割去一条边的实际意义
接下来,我们来考虑一下割去一条边的实际意义。
由于前面已经证过不会割药与药材的连边,因此这种边略过不提。
对于割掉一条从源点向药的连边,我们考虑到最后要将最小割减去从源点向外连出的边的容量总和,而在这一过程中,这条边的容量就抵消掉了。
因此,割掉一条从源点向药的连边,就相当于不选择这种药。
而对于割掉一条从药材连向汇点的连边则恰好相反,表示选择这种药材。
选择的药的数量与药材的数量必定相等
考虑从源点向药连出了(n)条边,从药材向汇点连入了(n)条边,则显然至少需要割(n)条边,且割(n)条边一定可行。
而割去(n)条边的代价必然可以表示成(n*INF+x),其中(x)为某个未知值。
由于割小于(n)条边不可行,而割大于(n)条边时代价可以表示为(n'*INF+y),其中(n'>n),(y)为某个未知值。
将两式相减得到:
由于(INF)是个很大的数,所以上述式子必然小于(0),因此若要求最小割,割去的边数不可能大于(n)。(这也是我们要把每条边的权值设成(INF)级别的原因)
综上所述,割去的边数必然为(n)。
又由于前面证过的连接药与药材的边绝不会被选作最小割中的边,所以这(n)条边必然由(t)条从源点向药的连边和(n-t)条从药材向汇点的连边组成。
再考虑割去一条边的实际意义,也就是说不选择的药的数量与选择的药材的数量相加为(n)。
显然不选择的药的数量与选择的药的数量相加也为(n)。
也就是说选择的药的数量与选择的药材的数量必定相等,性质得证。
以此方法能够得到合法答案
考虑前面已经证明过割去的边为(t)条从源点向药的连边和(n-t)条从药材向汇点的连边。
则代价可以表示为(n*INF)减去若干种药的权值。
而从源点向外连出的边的容量总和为(n*INF)减去所有药的权值和。
用代价减去容量总和,(n*INF)抵消,而代价中减去的那些药的权值与容量总和中这部分权值抵消(前面提到过割去一条从源点向药的连边表示不选这种药),剩下的权值负负得正成为其原本的权值。
然后求和即为答案。
以此方法选出的答案必为最优解
这是一个比较显然的性质。
我们求的是最小割,而最后答案是最小割减去容量总和。
既然是最小割,求的肯定是最小解,而答案就是最优解了。
代码
#include<bits/stdc++.h>
#define Tp template<typename Ty>
#define Ts template<typename Ty,typename... Ar>
#define Reg register
#define RI Reg int
#define Con const
#define CI Con int&
#define I inline
#define W while
#define N 300
#define INF (int)1e9
#define LL long long
using namespace std;
int n;
class FastIO
{
private:
#define FS 100000
#define tc() (A==B&&(B=(A=FI)+fread(FI,1,FS,stdin),A==B)?EOF:*A++)
#define tn (x<<3)+(x<<1)
#define D isdigit(c=tc())
int f;char c,*A,*B,FI[FS];
public:
I FastIO() {A=B=FI;}
Tp I void read(Ty& x) {x=0,f=1;W(!D) f=c^'-'?1:-1;W(x=tn+(c&15),D);x*=f;}
Ts I void read(Ty& x,Ar&... y) {read(x),read(y...);}
}F;
class MaxFlow//最大流
{
private:
#define Psz (N<<1)+2
#define Lsz (N*N<<1)+(N<<2)
#define add(x,y,c) (e[++ee].nxt=lnk[x],e[lnk[x]=ee].to=y,e[ee].Cap=c)
LL Ctot;int ee,lnk[Psz+5],cur[Psz+5],q[Psz+5],dep[Psz+5];
struct edge {int to,nxt,Cap;}e[Lsz+5];
I bool BFS()//BFS找增广路
{
RI i,k,H=1,T=1;memset(dep,0,sizeof(dep)),dep[q[1]=s]=1;W(H<=T&&!dep[t])
for(i=lnk[k=q[H++]];i;i=e[i].nxt) e[i].Cap&&!dep[e[i].to]&&(dep[q[++T]=e[i].to]=dep[k]+1);
return dep[t]?(memcpy(cur,lnk,sizeof(lnk)),true):false;
}
I int DFS(CI x,RI f)//DFS统计流量
{
if(!(x^t)||!f) return f;RI i,p,res=0;
for(i=cur[x];i;i=e[i].nxt)
{
if(cur[x]=i,(dep[x]+1)^dep[e[i].to]||!(p=DFS(e[i].to,min(f,e[i].Cap)))) continue;
if(e[i].Cap-=p,e[((i-1)^1)+1].Cap+=p,res+=p,!(f-=p)) break;
}return !res&&(dep[x]=-1),res;
}
public:
int s,t;I MaxFlow() {s=1,t=2;}
I int P1(CI x) {return x+2;}I int P2(CI x) {return x+n+2;}
I void Add(CI x,CI y,CI c) {add(x,y,c),add(y,x,0),Ctot+=c;}
I LL GetAns() {Reg LL t=0;W(BFS()) t+=DFS(s,2*INF);return t;}
}M;
int main()
{
RI i,x,y;Reg LL sum=0;for(F.read(n),i=1;i<=n;++i) for(F.read(x);x;--x) F.read(y),M.Add(M.P1(i),M.P2(y),INF);//建边
for(i=1;i<=n;++i) F.read(x),M.Add(M.s,M.P1(i),INF-x),M.Add(M.P2(i),M.t,INF),sum+=INF-x;//建边
return printf("%lld",M.GetAns()-sum),0;
}