LibreOJ #2003. 「SDOI2017」新生舞会

　　　　　　　　　　　　　　　　　　　　　　　　　　内存限制：256 MiB 时间限制：1500 ms 标准输入输出

　　　　　　　　　　　　　　　　　　　　　　　　　　　　题目类型：传统 评测方式：文本比较

　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　上传者： 匿名

01分数规划(并不知道这是啥。。)

km或费用流(并不会)验证

屠龙宝刀点击就送

#include <cstring>
#include <cstdio>
#include <queue>
#include <cmath>
#define N 100100
#define inf 0x3f3f3f3f
#define eps 1e-8
using namespace std;
bool vis[N];
int nextt[N],to[N],fa[N],flow[N],head[N],cnt=,n,ans1,ans2,a[][],b[][];
double cost[N],dis[N];
inline void ins(int u,int v,int f,double c)
{
    nextt[++cnt]=head[u];
    to[cnt]=v;
    flow[cnt]=f;
    cost[cnt]=c;
    head[u]=cnt;
}
bool spfa(int s,int t)
{
    for(int i=s;i<=t;++i) vis[i]=,dis[i]=-1e9;
    dis[s]=fa[s]=;
    queue<int>q;
    q.push(s);
    for(int now;!q.empty();)
    {
        now=q.front();q.pop();
        vis[now]=;
        for(int i=head[now];i;i=nextt[i])
        {
            int v=to[i];
            if(dis[v]<dis[now]+cost[i]&&flow[i]>)
            {
                dis[v]=dis[now]+cost[i];
                fa[v]=i;
                if(!vis[v])
                {
                    vis[v]=;
                    q.push(v);
                }
            }
        }
    }
    return dis[t]!=-1e9;
}
bool dinic(int s,int t,double k)
{
    cnt=;
    memset(head,,sizeof(head));
    for(int i=;i<=n;++i) ins(s,i,,),ins(i,s,,);
    for(int i=n+;i<=n*;++i) ins(i,t,,),ins(t,i,,);
    for(int i=;i<=n;++i)
     for(int j=;j<=n;++j)
      ins(i,n+j,,(double)a[i][j]-k*b[i][j]),ins(n+j,i,,-((double)a[i][j]-k*b[i][j]));
    double tmp=;
    bool flag=false;
    for(;spfa(s,t);)
    {
        if(tmp+dis[t]>=)
        {
            flag=true;
            for(int i=t;i!=s&&i;i=to[fa[i]^])
            {
                flow[fa[i]]-=;
                flow[fa[i]^]+=;
            }
            tmp+=dis[t];
        }
        else return false;
    }
    return flag;
}
int Main()
{
    scanf("%d",&n);
    for(int i=;i<=n;++i)
     for(int j=;j<=n;++j)
      scanf("%d",&a[i][j]);
    for(int i=;i<=n;++i)
     for(int j=;j<=n;++j)
      scanf("%d",&b[i][j]);
    int s=,t=n*+;
    double l=,r=,mid,ans;
    for(;abs(r-l)>eps;)
    {
        mid=(l+r)/;
        if(dinic(s,t,mid)) ans=mid,l=mid+eps;
        else r=mid-eps;
    }
    printf("%.6lf",ans);
    return ;
}
int sb=Main();
int main(int argc,char *argv[]) {;}