2015北京现场赛的金牌题,同步赛的时候最后一小时没开出来,受到了DAG的题目背景的误导,完全没有往网络流的方向去考虑。。T_T
题意:
给定一个描述游戏技能约束关系的DAG,约束关系u->v表示获得技能V的必要条件是获得技能u。获得第i个技能需要费用P[i]。每一个约束关系可以通过花费费用Ci而将其取消。另外,可以无视一切技能约束关系,而直接以费用Di强行获得技能i。问题是求解获得技能s所需的最小费用。
很好的一个题目~第一感觉是直接在DAG上进行DP,然而很快就会发现状态无法转移= =。正解是构建成最小割的模型。那么我们应该考虑,如何把获得技能u的任一方式,表为从源点到点u的一个割。从这个方向思考的话,便可以得到以下的建图策略:
(1)将每个点i,拆成i和i',另外添加一个源点source,以及汇点sink。
(2)对于原来DAG中的边,u->v,连边u'->v,容量为DAG中对应边的费用。
(3)对于每个点u,连边source->u,容量为P[u];连边u->u',容量为D[u]。
(4)对于目标技能s,连边s'->sink,流量INF。
然后答案就是source->sink的最小割。以上的建图方式的正确性仔细想想便能明白。
借助这个题也重新复习了一遍最大流的代码= = 希望能够在4天后的上海站有裸敲最大流的能力~ (T_T感觉自己好弱~一个最大流还要对着模板写)
#define N (1<<10)
#define INFP (1<<30)
#include <bits/stdc++.h>
using namespace std;
int T,n,m,s,x,y,z,cnt;
int head[N],pre[N],gap[N],level[N],cur[N];
struct Node
{
int v,cap,next;
Node(){}
Node(int v,int cap,int next):v(v),cap(cap),next(next){}
}Edge[N<<6];
void AddEdge(int u,int v,int cap)
{
Edge[cnt++]=Node(v,cap,head[u]);
Edge[cnt++]=Node(u,0,head[v]);
head[u]=cnt-2,head[v]=cnt-1;
}
int ISAP(int vs,int vt,int n)
{
int res=0;
memset(gap,0,sizeof(gap));
memset(pre,-1,sizeof(pre));
memset(level,0,sizeof(level));
for(int i=0;i<n;i++)
cur[i]=head[i];
pre[vs]=vs,gap[0]=n;
for(int u=vs,aug=INFP;level[vs]<n;)
{
bool flag=0;
for(int i=cur[u];i!=-1;i=Edge[i].next)
{
int v=Edge[i].v;
if(Edge[i].cap && level[u]==level[v]+1)
{
flag=1,pre[v]=u,cur[u]=i,u=v;
aug=min(aug,Edge[i].cap);
break;
}
}
if(!flag)
{
int mh=n;
for(int i=head[u];i!=-1;i=Edge[i].next)
{
int v=Edge[i].v;
if(Edge[i].cap && level[v]<mh) mh=level[v],cur[u]=i;
}
if(!--gap[level[u]]) break;
gap[level[u]=mh+1]++,u=pre[u];
}
else if(u==vt)
{
for(;u!=vs;u=pre[u])
{
int v=pre[u];
Edge[cur[v]].cap-=aug;
Edge[cur[v]^1].cap+=aug;
}
res+=aug,aug=INFP;
}
}
return res;
}
int main()
{
for(cin>>T;T--;)
{
cin>>n>>m>>s;
memset(head,-1,sizeof(head)),cnt=0;
for(int i=0;i<m;i++)
scanf("%d%d%d",&x,&y,&z),AddEdge(n+x,y,z);
for(int i=1;i<=n;i++)
scanf("%d",&x),AddEdge(0,i,x);
for(int i=1;i<=n;i++)
scanf("%d",&x),AddEdge(i,n+i,x);
AddEdge(n+s,n+n+1,INFP);
cout<<ISAP(0,n+n+1,n+n+2)<<endl;
}
return 0;
}