bzoj1449: [JSOI2009]球队收益

传送门
感觉就是费用流呀。
可以发现这道题与之前的费用流题有所不同，因为一个球队不论输还是赢都会获得收益。
这里就要用到一个技巧，我们可以假装比赛双方都输，然后修改赢的就ok辣。
然后就是每个人向终点连边，这里的费用随着流量的变化而变化，所以我们要用到拆边法。考虑一个人从赢win[i]次、输lose[i]次变为赢(win[i]+1)、输(lose[i]-1)次，总费用的变化量为c[i](win[i]+1)^2+d[i](lose[i]-1)^2-c[i]win[i]^2-d[i]*lose[i]^2=c[i](2*win[i]+1)-d[i]*(2*lose[i]-1)。就这样一条一条地加边，同时修改win[i]和lose[i]。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<queue>
#define inf 0x3f3f3f3f
#define N 6005
using namespace std;
queue<int> q;
struct edge{int from,to,next,f,v;}e[N*3];
int head[N],dis[N],v[N],from[N];
int win[N],lose[N],a[N],b[N],c[N],d[N];
int tot=1,n,m,S,T,ans;
void add(int x,int y,int f,int v){
    e[++tot]=(edge){x,y,head[x],f,v};
    head[x]=tot;
    e[++tot]=(edge){y,x,head[y],0,-v};
    head[y]=tot;
}
bool bfs(){
    for (int i=1;i<=T;i++) dis[i]=inf;
    int h=dis[S]=0,t=v[S]=1,x;
    q.push(S);
    while (!q.empty()){
        int x=q.front(); q.pop(); v[x]=0;
        for (int i=head[x];i;i=e[i].next)
            if (dis[e[i].to]>dis[x]+e[i].v&&e[i].f){
                dis[e[i].to]=dis[x]+e[i].v;
                from[e[i].to]=i;
                if (!v[e[i].to]){
                    v[e[i].to]=1;
                    q.push(e[i].to);
                }
            }
    }
    return dis[T]!=inf;
}
void del(){
    int x=inf;
    for (int i=from[T];i;i=from[e[i].from]) x=min(x,e[i].f);
    for (int i=from[T];i;i=from[e[i].from]){
        e[i].f-=x; e[i^1].f+=x; ans+=x*e[i].v;
    }
}
int main(){
    scanf("%d%d",&n,&m);
    S=n+m+1; T=S+1;
    for (int i=1;i<=n;i++)
        scanf("%d%d%d%d",&win[i],&lose[i],&c[i],&d[i]);
    for (int i=1;i<=m;i++){
        scanf("%d%d",&a[i],&b[i]);
        lose[a[i]]++; lose[b[i]]++;
    }
    for (int i=1;i<=n;i++)
        ans+=c[i]*win[i]*win[i]+d[i]*lose[i]*lose[i];
    for (int i=1;i<=m;i++){
        add(S,i,1,0),add(i,a[i]+m,1,0),add(i,b[i]+m,1,0);
        int tmp=c[a[i]]*(2*win[a[i]]+1)-d[a[i]]*(2*lose[a[i]]-1);
        add(a[i]+m,T,1,tmp);
        win[a[i]]++; lose[a[i]]--;
        tmp=c[b[i]]*(2*win[b[i]]+1)-d[b[i]]*(2*lose[b[i]]-1);
        add(b[i]+m,T,1,tmp);
        win[b[i]]++; lose[b[i]]--;
    }
    while (bfs()) del();
    printf("%d",ans);
}