题目描述:
题目分析:
这个题困难在于胜者负者都会有收益,没办法用一个流量来表示。
那么就先假设所有人在比赛中都是负,算出一个总收益,然后每一场必须有一个胜者,也就是说对于其中的一个人要将负变成胜,这样就可以用一个流量表示了。
假设第
i
i
i个人当前胜了
w
i
n
i
win_i
wini场,败了
l
o
s
e
i
lose_i
losei场,现在将一个lose变成一个win,收益变化为:
C
i
∗
(
2
∗
w
i
n
+
1
)
−
D
i
∗
(
2
∗
l
o
s
e
−
1
)
C_i*(2*win+1)-D_i*(2*lose-1)
Ci∗(2∗win+1)−Di∗(2∗lose−1)
然后
w
i
n
+
=
1
,
l
o
s
e
−
=
1
win+=1,lose-=1
win+=1,lose−=1
随着胜的场次越多,收益变化也越大。对于收益不固定这种情况,通常我们采用拆点或拆边解决,此处我们采用拆边,总共建出第
i
i
i个人的比赛场次条边,费用依次递增。
然后就是把每场比赛看成一个点,源点向其连容量为1的边,比赛向对应的两个人连容量为1,费用为0的边。然后每个人向汇点连容量为1,改变胜负状况所获收益为费用的比赛场次条边。
因为胜的场次越多,收益也越大,所以最小费用最大流一定会让胜负转化是连续的,即一定会先走赢一场的边,再走赢两场的边…
Code:
#include<cstdio>
#include<deque>
#include<cstring>
#include<algorithm>
#define maxn 6005
#define maxm 10005
using namespace std;
const int inf = 0x3f3f3f3f;
int n,m,S,T;
int fir[maxn],cur[maxn],nxt[maxm],to[maxm],c[maxm],w[maxm],tot=1;
inline void line(int x,int y,int z,int v){
nxt[++tot]=fir[x],fir[x]=tot,to[tot]=y,c[tot]=z,w[tot]=v;
nxt[++tot]=fir[y],fir[y]=tot,to[tot]=x,c[tot]=0,w[tot]=-v;
}
namespace Costflow{
int dis[maxn],ans;
bool vis[maxn],inq[maxn];
deque<int>q;
bool SPFA(){
memset(dis,0x3f,(T+1)<<2);
dis[T]=0,q.push_back(T);
while(!q.empty()){
int u=q.front();q.pop_front(),inq[u]=0;
for(int i=fir[u],v;i;i=nxt[i]) if(c[i^1]&&dis[v=to[i]]>dis[u]+w[i^1]){
dis[v]=dis[u]+w[i^1];
if(!inq[v]) {inq[v]=1;if(!q.empty()&&dis[v]<dis[q.front()]) q.push_front(v);else q.push_back(v);}
}
}
return dis[S]!=inf;
}
int aug(int u,int augco){
if(u==T) {ans+=dis[S]*augco;return augco;}
vis[u]=1;
int need=augco,delta;
for(int &i=cur[u];i;i=nxt[i]) if(c[i]&&!vis[to[i]]&&dis[u]==dis[to[i]]+w[i]){
delta=aug(to[i],min(need,c[i]));
c[i]-=delta,c[i^1]+=delta;
if(!(need-=delta)) break;
}
vis[u]=0;
return augco-need;
}
int solve(){
ans=0;
while(SPFA()) memcpy(cur,fir,(T+1)<<2),aug(S,inf);
return ans;
}
}
int W[maxn],L[maxn],C[maxn],D[maxn],X[maxn],Y[maxn],deg[maxn],ans;
int main()
{
scanf("%d%d",&n,&m),S=0,T=n+m+1;
for(int i=1;i<=n;i++) scanf("%d%d%d%d",&W[i],&L[i],&C[i],&D[i]);
for(int i=1;i<=m;i++) scanf("%d%d",&X[i],&Y[i]),deg[X[i]]++,deg[Y[i]]++,line(n+i,X[i],1,0),line(n+i,Y[i],1,0),line(S,n+i,1,0);
for(int i=1;i<=n;i++){
ans+=C[i]*W[i]*W[i]+D[i]*(L[i]+deg[i])*(L[i]+deg[i]);
for(int j=1;j<=deg[i];j++) line(i,T,1,C[i]*(2*(W[i]+j)-1)-D[i]*(2*(L[i]+deg[i]-j)+1));
}
printf("%d\n",ans+Costflow::solve());
}