1449: [JSOI2009]球队收益
Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 941 Solved: 544
[ Submit][ Status][ Discuss]
Description
![](http://www.lydsy.com/JudgeOnline/images/1449_1.jpg)
Input
![](http://www.lydsy.com/JudgeOnline/images/1449_2.jpg)
Output
一个整数表示联盟里所有球队收益之和的最小值。
Sample Input
3 3
1 0 2 1
1 1 10 1
0 1 3 3
1 2
2 3
3 1
1 0 2 1
1 1 10 1
0 1 3 3
1 2
2 3
3 1
Sample Output
43
题解:
费用流好题。
我们先假设所有队伍都输,然后赢得队伍加上相差的值。
这样就可以跑费用流了。
(C*(win+1)^2+D*(lose-1)^2) - (C*win^2+D*lose^2)=C*(2*win+1)-D*(2*lose-1)
这是赢一次要加上的值。
以此类推。
对于每场比赛我们,将x->比赛,y->比赛,流量为1,费用为0、
比赛->ed,流量为1,费用为0.
接着,从st->每个队伍连di条边(di等于这个队伍还未比的赛),流量为1,费用为赢(1~di)场比赛要加的值。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<queue>
#define INF 0x7fffffff
using namespace std;
const int N=6060;
struct node{
int x,y,z,next,other,cost;
}sa[1001001];int len=0,first[N];
int n,m,st,ed;
int win[N],lose[N],c[N],D[N],d[N],ans=0;
void ins(int x,int y,int z,int w)
{
len++;
sa[len].x=x;
sa[len].y=y;
sa[len].z=z;
sa[len].cost=w;
sa[len].next=first[x];
sa[len].other=len+1;
first[x]=len;
len++;
sa[len].x=y;
sa[len].y=x;
sa[len].z=0;
sa[len].cost=-w;
sa[len].next=first[y];
sa[len].other=len-1;
first[y]=len;
}
int dis[N],pre[N],vis[N],frpe[N];
bool spfa()
{
memset(pre,-1,sizeof(pre));
memset(vis,0,sizeof(vis));
for(int i=0;i<=ed;i++) dis[i]=INF;
dis[st]=0;
pre[st]=0;
vis[st]=1;
queue<int>q;
q.push(st);
while(!q.empty())
{
int x=q.front();
q.pop();
vis[x]=0;
for(int i=first[x];i!=-1;i=sa[i].next)
{
int y=sa[i].y;
if(sa[i].z&&dis[y]>dis[x]+sa[i].cost)
{
dis[y]=dis[x]+sa[i].cost;
if(!vis[y])
{
vis[y]=1;
q.push(y);
}
pre[y]=x;
frpe[y]=i;
}
}
}
if(dis[ed]>=INF) return false;
return true;
}
void MCMF()
{
while(spfa())
{
int minl=INF;
for(int i=ed;i!=st;i=pre[i])
{
if(minl>sa[frpe[i]].z)
minl=sa[frpe[i]].z;
}
ans+=minl*dis[ed];
for(int i=ed;i!=st;i=pre[i])
{
sa[frpe[i]].z-=minl;
sa[sa[frpe[i]].other].z+=minl;
}
}
printf("%d\n",ans);
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++)
scanf("%d%d%d%d",&win[i],&lose[i],&c[i],&D[i]);
memset(d,0,sizeof(d));
memset(first,-1,sizeof(first));
st=0;ed=6060-2;
for(int i=1;i<=m;i++)
{
int x,y;
scanf("%d%d",&x,&y);
ins(x,n+i,1,0);
ins(y,n+i,1,0);
ins(n+i,ed,1,0);
d[x]++;d[y]++;
}
for(int i=1;i<=n;i++)
{
ans+=c[i]*win[i]*win[i]+(d[i]+lose[i])*(lose[i]+d[i])*D[i];
for(int j=1;j<=d[i];j++)
ins(st,i,1,c[i]*(2*win[i]+2*j-1)-D[i]*(2*(d[i]+lose[i])-2*j+1));
}
MCMF();
}