题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6598
题目大意:
你有很多个士兵,每个士兵有两种选择:成为法师或者成为战士。
现在有m种关系,表示两个士兵u和v,他们有三种选择:
- 两者都是战士,贡献为a
- 两者都是法师,贡献为c
- 两者一个是战士,一个是法师,贡献为a/4+c/3
问你该如何选择你的士兵的职业,让你的团队贡献最大?
题解:
官方题解很详细
(据说这是常用建图方式,反正我是没想到)
自己写的时候还有一个最小费用流的思路
一个点,分为Mi和Wi两个点,如果它与点j有关联,则
Mi连接Mj,费用为-a;
Mi连接Wj,费用为-b;
Wi链接Mj,费用为-b,
Wi连接Wj,费用为-c
这样下来,在起点跑流量为1的最小费用流即可。
不过我写的T了,如果有时间再来好好写验证一下吧
#include <bits/stdc++.h>
#define LL long long
#define INF 0x3f3f3f3f3f3f3f3fll
#define sc(x) scanf("%d",&x)
#define scc(x,y) scanf("%d%d",&x,&y)
#define sccc(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define file(x) freopen(#x".in","w",stdout);
using namespace std;
using std::default_random_engine;
const int mod = 998244353;
const int N = 3010;
const int M = 2e5 + 10;
namespace ISAP
{
LL H[M],d[M],cur[M],pre[M],gap[M],Q[M];
struct Edge{int u,v,c,n;} E[M];
LL nv,flow,head,tail,cntE,f;
void init(){cntE=0; memset(H,-1,sizeof(H));}
void addedge(int u,int v,int c)
{
E[cntE]=Edge{u,v,c,H[u]}; H[u]=cntE++;
E[cntE]=Edge{v,u,0,H[v]}; H[v]=cntE++;
}
void revbfs(int s,int t)
{
head=tail=0 ;
memset(d,-1,sizeof(d));
memset(gap,0,sizeof(gap));
Q[tail++]=t;d[t]=0;gap[d[t]]=1;
while (head!=tail)
{
int u=Q[head++];
for (int i=H[u];~i;i=E[i].n)
{
int v=E[i].v; if (~d[v]) continue;
d[v]=d[u]+1; gap[d[v]]++; Q[tail++]=v;
}
}
}
LL isap(int s,int t)
{
memcpy(cur,H,sizeof(cur)); nv=t;
flow=0; revbfs(s,t); int u=pre[s]=s,i;
while (d[s]<nv)
{
if (u==t)
{
f=INF;
for (i=s;i!=t;i=E[cur[i]].v)
if (f>E[cur[i]].c) f=E[cur[i]].c,u=i;
flow += f;
for (i=s;i!=t;i=E[cur[i]].v)
E[cur[i]].c-=f,E[cur[i]^1].c+=f;
}
for (i=cur[u];~i;i=E[i].n)
if (E[i].c&&d[u]==d[E[i].v]+1) break ;
if (~i) cur[u]=i,pre[E[i].v]=u,u=E[i].v;
else
{
if (0==--gap[d[u]]) break ;
LL minv=nv,v;
for (int i=H[u];~i;i=E[i].n)
{
v=E[i].v;
if (E[i].c&&minv>d[v]) minv=d[v],cur[u]=i;
}
d[u]=minv+1; gap[d[u]]++; u=pre[u];
}
}
return flow ;
}
}
int main()
{
int n,m;
while(~scc(n,m))
{
ISAP::init();
LL ans=0;
for(int i=1;i<=m;i++)
{
int u,v,a,b,c;
scc(u,v); sccc(a,b,c);
ans+=(a+b+c)*2;
ISAP::addedge(n+1,u,5*a/4+c/3);
ISAP::addedge(n+1,v,5*a/4+c/3);
ISAP::addedge(u,n+2,a/4+4*c/3);
ISAP::addedge(v,n+2,a/4+4*c/3);
ISAP::addedge(u,v,a/2+c/3);
ISAP::addedge(v,u,a/2+c/3);
}
LL tmp = ans-ISAP::isap(n+1,n+2);
assert(tmp%2==0);
printf("%lld\n",tmp/2);
}
return 0;
}