Coding Contest-青岛区域赛网络流

HDU5988
建立原点与汇点,将原点与(人比食物多)的点相连,流量为(人-食物)。将(食物比人多)的点与汇点相连接,流量为(食物-人)。
将每个桌子相连,假设概率为p,对p取log,可以将原来相乘的计算转化为相加,加n-1条费用为(logp)的边,加一条费用为0的边。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iomanip>
#include<set>
#include<algorithm>
#include<sstream>
#include<string>
#include<iostream>
#include<queue>
#include<map>
using namespace std;
const int Vmax=205; //需要拆点的话记得加倍
#define inf 1e9
#define eps 1e-8
namespace MCMF
{
struct Edge
{
    int from,to,cap,flow;
    double cost;
    Edge() {}
    Edge(int u,int v,int c,int f,double w):from(u),to(v),cap(c),flow(f),cost(w) {}
};
int n,m;
vector<Edge>edges;
vector<int>G[Vmax];
int inq[Vmax]; //是否在队列中
double d[Vmax]; //Bellman-Ford
int p[Vmax]; //上一条弧
int a[Vmax]; //可改进量
void init(int _Vsz)
{
    n=_Vsz;
    for(int i=0; i<=n; i++) G[i].clear();
    edges.clear();
}
void adde(int from,int to,int cap,double cost)
{
    edges.push_back(Edge(from, to, cap, 0, cost));
    edges.push_back(Edge(to, from, 0, 0, -cost));
    m=edges.size();
    G[from].push_back(m-2);
    G[to].push_back(m-1);
}
bool SPFA(int s,int t,int& flow,double& cost)
{
    for(int i=0; i<=n; i++) d[i]=inf;
    memset(inq, 0, sizeof(inq));
    d[s]=0;
    inq[s]=1;
    p[s]=0;
    a[s]=inf;
    queue<int>q;
    q.push(s);
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        inq[u]=0;
        for(int i=0; i<G[u].size(); i++)
        {
            Edge& e=edges[G[u][i]];
            if(e.cap>e.flow&&d[e.to]>d[u]+e.cost+eps)
            {
                d[e.to]=d[u]+e.cost; //松弛操作
                p[e.to]=G[u][i]; //记录上一条边信息
                a[e.to]=min(a[u], e.cap-e.flow);
                if(!inq[e.to])
                {
                    q.push(e.to);
                    inq[e.to]=1;
                }
            }
        }
    }
    if(d[t]==inf) return false; //s-t 不联通,失败退出
    flow+=a[t];
    cost+=d[t]*a[t];
    for(int u=t; u!=s; u=edges[p[u]].from)
    {
        edges[p[u]].flow+=a[t];
        edges[p[u]^1].flow-=a[t];
    }
    return true;
}
double MincostMaxflow(int s,int t,double& cost)
{
    int flow=0;
    cost=0.0;
    while(SPFA(s, t, flow, cost));
    return cost;
}
}
int main()
{
#ifdef ONLINE_JUDGE
#else
    freopen("C:\\Users\\Administrator\\Desktop\\tmp.txt","r",stdin);
#endif

    int t;
    scanf("%d",&t);
    for(int kc=0; kc<t; ++kc)
    {
        int n,m,ta,tb,tc;
        double td;
        scanf("%d%d",&n,&m);
        MCMF::init(n+2);
        for(int i=1; i<=n; ++i)
        {
            scanf("%d%d",&ta,&tb);
            if(ta>tb)
            {
                MCMF::adde(n+1,i,ta-tb,0.0);
            }
            else if(tb>ta)
            {
                MCMF::adde(i,n+2,tb-ta,0.0);
            }
        }
        for(int i=0; i<m; ++i)
        {
            scanf("%d%d%d%lf",&ta,&tb,&tc,&td);
            if(tc>0) MCMF::adde(ta,tb,1,0.0);
            if(tc>1) MCMF::adde(ta,tb,tc-1,-log2(1-td));
        }
        double ans;
        ans=MCMF::MincostMaxflow(n+1,n+2,ans);
        ans=pow(2,-ans);
        printf("%.2f\n",1-ans);
    }
    return 0;
}

最小费用最大流!最小费用最大流!最小费用最大流!
题意很明显的提示了端点之间的流量,原点与汇点如何的建立。自然而然的想到了网络流。但是以前做过得题目都是费用相加,第一次见到相乘的一脸懵逼。花了3个小时去解决,一直在想怎么样修改代码实现相乘的最小费用流。回来后看了别人的讨论说取log转化为相加,一下子恍然大悟。实力坑队友。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值