[BZOJ]4819: [Sdoi2017]新生舞会 01分数规划+费用流

Description

学校组织了一次新生舞会，Cathy作为经验丰富的老学姐，负责为同学们安排舞伴。有n个男生和n个女生参加舞会
买一个男生和一个女生一起跳舞，互为舞伴。Cathy收集了这些同学之间的关系，比如两个人之前认识没计算得出
a[i][j] ，表示第i个男生和第j个女生一起跳舞时他们的喜悦程度。Cathy还需要考虑两个人一起跳舞是否方便，
比如身高体重差别会不会太大，计算得出 b[i][j]，表示第i个男生和第j个女生一起跳舞时的不协调程度。当然，
还需要考虑很多其他问题。Cathy想先用一个程序通过a[i][j]和b[i][j]求出一种方案，再手动对方案进行微调。C
athy找到你，希望你帮她写那个程序。一个方案中有n对舞伴，假设没对舞伴的喜悦程度分别是a’1,a’2,…,a’n，
假设每对舞伴的不协调程度分别是b’1,b’2,…,b’n。令
C=(a’1+a’2+…+a’n)/(b’1+b’2+…+b’n),Cathy希望C值最大。

题解：

裸的01分数规划+费用流……不过我10s卡时过……难道是二分次数太多？好像还有一种很快的迭代方法？以后再学。

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
const int Maxn=110;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return x*f;
}
struct Edge{int x,y,d,next;double c;}e[2110000];
int last[Maxn*2],len;
void ins(int x,int y,int d,double c)
{
    int t=++len;
    e[t].x=x;e[t].y=y;e[t].d=d;e[t].c=c;
    e[t].next=last[x];last[x]=t;
}
int n,a[Maxn][Maxn],b[Maxn][Maxn];
void addedge(int x,int y,int d,double c){ins(x,y,d,c);ins(y,x,0,-c);}
int st,ed;
double f[Maxn*2];
int pre[Maxn*2];
bool in[Maxn*2];
queue<int>q;
bool spfa()
{
    for(int i=1;i<=205;i++)f[i]=100000000000000.0;f[st]=0;
    memset(in,false,sizeof(in));
    q.push(st);pre[st]=-1;
    while(!q.empty())
    {
        int x=q.front();q.pop();in[x]=false;
        for(int i=last[x];i;i=e[i].next)
        {
            int y=e[i].y;
            if(e[i].d&&f[x]+e[i].c<f[y])
            {
                pre[y]=i;
                f[y]=f[x]+e[i].c;
                if(!in[y])in[y]=true,q.push(y);
            }
        }
    }
    return(f[ed]!=100000000000000.0);
}
double ans;
void work()
{
    int mind=2147483647,now=ed;
    while(pre[now]!=-1)
    {
        mind=min(mind,e[pre[now]].d);
        now=e[pre[now]].x;
    }
    now=ed;
    while(pre[now]!=-1)
    {
        e[pre[now]].d-=mind;e[pre[now]^1].d+=mind;
        ans+=(double)(mind)*e[pre[now]].c;
        now=e[pre[now]].x;
    }
}
bool check(double c)
{
    len=1;ans=0.0;
    memset(last,0,sizeof(last));
    st=2*n+1;ed=2*n+2;
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
    addedge(i,j+n,1,-((double)(a[i][j])-(double)(b[i][j])*c));
    for(int i=1;i<=n;i++)addedge(st,i,1,0),addedge(i+n,ed,1,0);
    while(spfa())work();
    return (ans<=0.0);
}
int main()
{
    n=read();
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
    a[i][j]=read();
    for(int i=1;i<=n;i++)
    for(int j=1;j<=n;j++)
    b[i][j]=read();
    double l=0,r=100000000;
    int c=60;
    while(c--)
    {
        double mid=(l+r)/2.0;
        if(check(mid))l=mid;
        else r=mid;
    }
    printf("%.6lf",l);
}