HDU5988
建立原点与汇点,将原点与(人比食物多)的点相连,流量为(人-食物)。将(食物比人多)的点与汇点相连接,流量为(食物-人)。
将每个桌子相连,假设概率为p,对p取log,可以将原来相乘的计算转化为相加,加n-1条费用为(logp)的边,加一条费用为0的边。
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
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 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; iif(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; kcint 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; iscanf("%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转化为相加,一下子恍然大悟。实力坑队友。