BZOJ1449: [JSOI2009]球队收益【胜负收益转化 费用流】

题目描述：

题目分析：

这个题困难在于胜者负者都会有收益，没办法用一个流量来表示。

那么就先假设所有人在比赛中都是负，算出一个总收益，然后每一场必须有一个胜者，也就是说对于其中的一个人要将负变成胜，这样就可以用一个流量表示了。

假设第$i$个人当前胜了$wi{n}_i$场，败了$los{e}_i$场，现在将一个lose变成一个win，收益变化为：$$C_i\ast (2\ast win+1)-D_i\ast (2\ast lose-1)$$
然后$win+=1,lose-=1$
随着胜的场次越多，收益变化也越大。对于收益不固定这种情况，通常我们采用拆点或拆边解决，此处我们采用拆边，总共建出第$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());
}