BZOJ 1449/2895: [JSOI2009]球队收益|网络流

最新推荐文章于 2019-12-05 20:06:24 发布

ws_yzy

最新推荐文章于 2019-12-05 20:06:24 发布

阅读量512

点赞数

分类专栏：网络流文章标签：费用流

本文链接：https://blog.csdn.net/ws_yzy/article/details/50551041

版权

网络流专栏收录该内容

23 篇文章 0 订阅

订阅专栏

好流弊的一道网络流题啊，涨姿势
首先假设让参加比赛的两个队都输
然后调整每场比赛赢的队伍来获得最小答案
这不就是最小费用流吗！
S连每场比赛流量1费用0
每场比赛连参赛队流量1费用0
关键就是参赛队连T 需要连参赛次数条边
我们发现调整一次由win，lose变为 win+1，lose-1的费用为
(C*(win+1)^2+D*(lose-1)^2) - (C*win^2+D*lose^2)=C*(2*win+1)-D*(2*lose-1)
暴力连边就可以了
最后的答案=最想费用流+最初假设所有队伍都输的收益

#include<set>
#include<map>
#include<ctime>
#include<queue>
#include<cmath>
#include<cstdio>
#include<vector>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
#define T 6006
#define MX 1e9
using namespace std;
int sc()
{
    int i=0,f=1; char c=getchar();
    while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
    return i*f;
}
int dis[T],q[T],vis[T],from[T],pre[T];
int head[T],lst[T*10],nxt[T*10],v[T*10],c[T*10];
int win[T],lose[T],C[T],D[T],sum[T];
int n,m,S,W,tot=1,ans;
void insert(int x,int y,int _c,int _v)
{
    lst[++tot]=y; nxt[tot]=head[x]; c[tot]=_c; v[tot]=_v; head[x]=tot;
    lst[++tot]=x; nxt[tot]=head[y]; c[tot]=0; v[tot]=-_v; head[y]=tot;
}
bool spfa()
{
    for(int i=1;i<=W;i++) dis[i]=MX; dis[S]=0;
    int l=1,r=2; q[1]=S;
    while(l!=r)
    {
        int x=q[l++]; l%=T; vis[x]=0;
        for(int i=head[x];i;i=nxt[i])
            if(c[i]&&dis[lst[i]]>dis[x]+v[i])
            {
                dis[lst[i]]=dis[x]+v[i];
                from[lst[i]]=x;
                pre[lst[i]]=i;
                if(!vis[lst[i]])
                {
                    vis[lst[i]]=1;
                    q[r++]=lst[i];
                    r%=T;
                }
            }
    }
    return dis[W]!=MX;
}
void mcf()
{
    int cc=MX;
    for(int i=W;i!=S;i=from[i])cc=min(cc,c[pre[i]]);
    for(int i=W;i!=S;i=from[i])
    {
        int k=pre[i];
        ans+=v[k]*cc;
        c[k]-=cc;
        c[k^1]+=cc;
    }
}
int main()
{
    n=sc();m=sc();S=n+m+1,W=S+1;
    for(int i=1;i<=n;i++)
    {
        win[i]=sc();
        lose[i]=sc();
        C[i]=sc();
        D[i]=sc();
    }
    for(int i=1;i<=m;i++)
    {
        int x=sc(),y=sc();
        insert(S,n+i,1,0);
        insert(n+i,x,1,0);
        insert(n+i,y,1,0);
        sum[x]++,sum[y]++;
    }
    for(int i=1;i<=n;i++)lose[i]+=sum[i];
    for(int i=1;i<=n;i++)ans+=C[i]*win[i]*win[i]+D[i]*lose[i]*lose[i];
    //cout << ans<<endl;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=sum[i];j++)
        {
            insert(i,W,1,C[i]*(2*win[i]+1)-D[i]*(2*lose[i]-1));
            win[i]++; lose[i]--;
        }
    while(spfa())mcf();
    cout<<ans;
    return 0;
}