【BZOJ-1449&2895】球队收益&球队预算最小费用最大流

1449: [JSOI2009]球队收益

Time Limit: 5 Sec Memory Limit: 64 MB
Submit: 648 Solved: 364
[Submit][Status][Discuss]

Description

Input

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

HINT

Source

2895: 球队预算

Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 51 Solved: 50
[Submit][Status][Discuss]

Description

在一个篮球联赛里，有n支球队，球队的支出是和他们的胜负场次有关系的，具体来说，第i支球队的赛季总支出是Ci*x^2+Di*y^2，Di<=Ci。(赢得多，给球员的奖金就多嘛)

其中x,y分别表示这只球队本赛季的胜负场次。现在赛季进行到了一半，每只球队分别取得了a[i]场胜利和b[i]场失利。而接下来还有m场比赛要进行。问联盟球队的最小总支出是多少。

Input

第一行n，m

接下来n行每行4个整数a[i],b[i],Ci,Di

再接下来m行每行两个整数s，t表示第s支队伍和第t支队伍之间将有一场比赛，注意两只队间可能有多场比赛。

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
Data Limit
对于20%的数据2<=n<=10,0<=m<=20
对于100%的数据2<=n<=5000,0<=m<=1000,0<=di<=ci<=10,0<=a[i],b[i]<=50.

HINT

Source

Solution

idea很好的建图，不能考虑直接建，而应该利用增量的思想

个人的想法是这样的：

源到各比赛连边，容量为1，费用为0

建每个球队分成两个，由个比赛连向两只球队a,b，再由a连b'，b连a'，分别表示a赢了这场比赛或b赢了这场比赛

然后由各个点i'连向汇容量为inf，费用为0

然后预处理出来已比完的比赛胜负的价值，加上最小费用即可；

思路是考虑一场比赛的胜负对于结果的贡献，从这场比赛结束后对结果的增量来考虑，但显然这种方法并不是很优，但理论上可A

那么同样的思想，换一种方法建图：

不妨先假设后面的M场比赛中双方都是输家，这样我们只要在模型中表示一方成为赢家即可。

此应该有一个初步的模型了：对于N支球队和M场比赛各建一个点，从源向每场比赛连流量1费用0的边，从比赛向参与这场比赛的两支队伍各连一条流量1费用0的边。剩下的就是队伍收益的费用表示了。

多赢一场比赛产生的收益。即(C*(w+1)^2+D*(l-1)^2)-(C*w^2+D*l^2)=2w*C-2l*D+C+D。对于第i支队伍，假设后M场中i参加的有x场，那么最初w=win,l=lose+x，之后每赢一场w++,l--。我们从第i支队伍的点向汇连x条边，分别代表第i支队伍赢了j场比赛时相对赢j-1场时收益的增量。由于增量一定越来越大（平方嘛），所以流量最先流过的一定是费用较小的边，即j最小的边。

　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　-----by huzecong

Code

#include<iostream>

#include<cstdio>

#include<algorithm>

#include<cmath>

#include<cstring>

#include<queue>

using namespace std;

int read()

{

    int x=,f=; char ch=getchar();

    while (ch<'' || ch>'') {if (ch=='-')f=-; ch=getchar();}

    while (ch>='' && ch<='') {x=x*+ch-''; ch=getchar();}

    return x*f;

}

#define maxn 5100+1100

#define maxm 100100

int n,m,ans,S,T,num[];

struct Teamnode{int win,lose,C,D;}team[];

struct Edgenode{int to,next,cap,cost,from;}edge[maxm];

int head[maxn],cnt=;

void add(int u,int v,int w,int c)

{cnt++;edge[cnt].from=u;edge[cnt].to=v;

edge[cnt].next=head[u];head[u]=cnt;

edge[cnt].cost=c;edge[cnt].cap=w;}

void insert(int u,int v,int w,int c){add(u,v,w,c);add(v,u,,-c);}

#define inf 0x7fffffff

int q[maxn],h,t; int dis[maxn];bool mark[maxn];

bool spfa()

{

    memset(mark,,sizeof(mark));

    for (int i=S; i<=T; i++) dis[i]=inf;

    h=,t=;

    q[]=T; mark[T]=; dis[T]=;

    while (h<t)

        {

            int now=q[h++]; mark[now]=;

            for (int i=head[now]; i; i=edge[i].next)

                if (edge[i^].cap && dis[now]+edge[i^].cost<dis[edge[i].to])

                    {

                        dis[edge[i].to]=dis[now]+edge[i^].cost;

                        if (!mark[edge[i].to])

                            mark[edge[i].to]=,q[t++]=edge[i].to;

                    }

        }

    return dis[S]!=inf;

}

int dfs(int loc,int low)

{

    mark[loc]=;

    if (loc==T) return low;

    int w,used=;

    for (int i=head[loc]; i; i=edge[i].next)

        if (!mark[edge[i].to] && edge[i].cap && dis[edge[i].to]==dis[loc]-edge[i].cost)

            {

                w=dfs(edge[i].to,min(low-used,edge[i].cap));

                edge[i].cap-=w; edge[i^].cap+=w;

                ans+=w*edge[i].cost; used+=w;

                if (used==low) return low;

            }

    return used;

}

int zkw()

{

    int tmp=;

    while (spfa())

        {

            mark[T]=;

            while (mark[T])

                {

                    memset(mark,,sizeof(mark));

                    tmp+=dfs(S,inf);

                }

        }

    return tmp;

}

int power(int x){return x*x;}

void build()

{

    S=,T=n+m+;

    for(int a,b,i=; i<=m; i++)

        insert(S,i,,),

        a=read(),b=read(),

        insert(i,a+m,,),insert(i,b+m,,),

        num[a]++,num[b]++;

    for(int i=; i<=n; i++)

        team[i].lose+=num[i];

    for(int i=; i<=n; i++)

        ans+=power(team[i].win)*team[i].C+power(team[i].lose)*team[i].D;

    for(int i=; i<=n; i++)

         for(int j=;j<=num[i];j++)

            insert(i+m,T,,*team[i].C*team[i].win+team[i].C+team[i].D-*team[i].D*team[i].lose),

            team[i].lose--,team[i].win++;

}

int main()

{

    n=read(),m=read();

    for (int i=; i<=n; i++)

        team[i].win=read(),team[i].lose=read(),

        team[i].C=read(),team[i].D=read();

    build();

    int Maxflow=zkw();

    printf("%d\n",ans);

    return ;

}

吐槽一句..增广的过程中，一开始使用STL中的queue然后莫名的WA了？手写了个q...怒A...跑得飞快