【BZOJ1449】 球队收益

BZOJ1449 球队收益

Output

一个整数表示联盟里所有球队收益之和的最小值。

Sample Input

3 3

1 0 2 1

1 1 10 1

0 1 3 3

1 2

2 3

3 1

Sample Output

43

我们先假设所有的球队都赢，算出答案。然后每场比赛都要提供一个输的场次。

考虑费用流。源点向每场比赛连边，每场比赛向两只队伍连边，队伍再向汇点连边。

注意到一只队伍的得分是关于输的场次的一个二次函数，所以每增加一个输场，增加或减少的收益不一样。所以我们拆边。设\(f_{i,k}\)表示第\(i\)只队伍输\(k\)场的收益，则第\(a\)条边的权值为\(f_{i,a}-f_{i,a-1}\)。然后再跑最小费用流。

因为这是个开口向上的二次函数，它的二阶导是大于\(0\)的。 所以\(f_{i,a}-f_{i,a-1}\)单调递增，这么拆边是对的。

代码：

#include<bits/stdc++.h>
#define ll long long
#define N 15005

using namespace std;
inline int Get() {int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=getchar();}while('0'<=ch&&ch<='9') {x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f;}

int n,m;
int S,T;
int win[N],lose[N];
int c[N],d[N];
int a[N],b[N];
struct road {
	int to,next;
	int f,c;
}s[N<<4];

int h[N<<2],cnt=1;
void add(int i,int j,int f,int c) {
	s[++cnt]=(road) {j,h[i],f,c};h[i]=cnt;
	s[++cnt]=(road) {i,h[j],0,-c};h[j]=cnt;
}

int tim[N];
ll ans;
ll dis[N];
queue<int>q;
int fr[N],e[N];
bool in[N];
bool spfa(int S) {
	memset(dis,0x3f,sizeof(dis));
	dis[S]=0;
	q.push(S);
	while(!q.empty()) {
		int v=q.front();q.pop();
		in[v]=0;
		for(int i=h[v];i;i=s[i].next) {
			int to=s[i].to;
			if(s[i].f&&dis[to]>dis[v]+s[i].c) {
				dis[to]=dis[v]+s[i].c;
				fr[to]=v;
				e[to]=i;
				if(!in[to]) {
					in[to]=1;
					q.push(to);
				}
			}
		}
	}
	if(dis[T]>1e9) return 0;
	ans+=dis[T];
	for(int i=T;i;i=fr[i]) {
		s[e[i]].f--;
		s[e[i]^1].f++;
	}
	return 1;
}

int main() {
	n=Get(),m=Get();
	for(int i=1;i<=n;i++) {
		win[i]=Get(),lose[i]=Get();
		c[i]=Get(),d[i]=Get();
	}
	T=m+n+1;
	for(int i=1;i<=m;i++) {
		a[i]=Get(),b[i]=Get();
		add(S,i,1,0);
		tim[a[i]]++,tim[b[i]]++;
		add(i,a[i]+m,1,0),add(i,b[i]+m,1,0);
	}
	for(int i=1;i<=n;i++) {
		ll last=1ll*c[i]*(win[i]+tim[i])*(win[i]+tim[i])+1ll*d[i]*lose[i]*lose[i],now;
		ans+=last;
		for(int j=1;j<=tim[i];j++) {
			now=1ll*c[i]*(win[i]+tim[i]-j)*(win[i]+tim[i]-j)+1ll*d[i]*(lose[i]+j)*(lose[i]+j);
			add(i+m,T,1,now-last);
			last=now;
		}
	}
	while(spfa(S));
	cout<<ans;
	return 0;
}