bzoj4819: [Sdoi2017]新生舞会

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值最大。
Input

第一行一个整数n。
接下来n行，每行n个整数，第i行第j个数表示a[i][j]。
接下来n行，每行n个整数，第i行第j个数表示b[i][j]。
1<=n<=100,1<=a[i][j],b[i][j]<=10^4
Output

一行一个数，表示C的最大值。四舍五入保留6位小数，选手输出的小数需要与标准输出相等
Sample Input

3

19 17 16

25 24 23

35 36 31

9 5 6

3 4 2

7 8 9
Sample Output

5.357143

题解

今天CF打挂了，不开心
本来不想写这题的。。但是觉得自己总是觉得..就不写代码不是很好。。
比较切..还是很有必要的。。
于是这题你就分数规划一下，费用流盘可行就好了

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<queue>
using namespace std;
const int MAX=1<<28;
const int N=205;
const double eps=1e-8;
int n;
int a[N][N];
int b[N][N];
int st,ed;
struct qt
{
    int x,y,z,last;//流量    费用 
    double z1;
}s[105*105*2];int num,last[N];
void init (int x,int y,int z,double z1)
{
    num++;
    s[num].x=x;s[num].y=y;s[num].z=z;s[num].z1=z1;
    s[num].last=last[x];
    last[x]=num;
    swap(x,y);z=0;z1=-z1;
    num++;
    s[num].x=x;s[num].y=y;s[num].z=z;s[num].z1=z1;
    s[num].last=last[x];
    last[x]=num;
}
int from[N];
double f[N];
bool in[N];
bool SPFA ()
{
    for (int u=1;u<=ed;u++) 
    {
        f[u]=MAX;
        from[u]=-1;
        in[u]=false;
    }
    queue<int> q;
    q.push(st);
    f[st]=0;in[st]=true;
    while (!q.empty())
    {
        int x=q.front();q.pop();
        for (int u=last[x];u!=-1;u=s[u].last)
        {
            int y=s[u].y;
            if (s[u].z>0&&f[y]>f[x]+s[u].z1)
            {
                f[y]=f[x]+s[u].z1;
                from[y]=u;
                if (in[y]==false)
                {
                    in[y]=true;
                    q.push(y);
                }
            }
        }
        in[x]=false;
    }
    return from[ed]!=-1;
}
double ans=0;
void lalal ()
{
    int x=from[ed];
    int t=MAX;
    while (x!=-1)
    {
        t=min(t,s[x].z);
        x=from[s[x].x];
    }

    x=from[ed];
    while (x!=-1)
    {
        ans=ans+t*s[x].z1;
        s[x].z-=t;
        s[x^1].z+=t;
        x=from[s[x].x];
    }
    return ;
}
bool check (double k)
{
    num=1;memset(last,-1,sizeof(last));
    for (int u=1;u<=n;u++)
    {
        init(st,u,1,0);
        for (int i=1;i<=n;i++)
            init(u,i+n,1,k*(double)b[u][i]-(double)a[u][i]);
        init(u+n,ed,1,0);
    }
    ans=0;
    while (SPFA()) lalal();
    return ans<=0.0;
}
int main()
{
    scanf("%d",&n);st=n*2+1;ed=st+1;
    for (int u=1;u<=n;u++)
        for (int i=1;i<=n;i++)
            scanf("%d",&a[u][i]);
    for (int u=1;u<=n;u++)
        for (int i=1;i<=n;i++)
            scanf("%d",&b[u][i]);
    double l=0.0,r=2000000;
    while (l+eps<r)
    {
        double mid=(l+r)/2;
        if (check(mid))l=mid;
        else r=mid;
    }
    printf("%.6lf",l);
    return 0;
}