二分解法及题目详解
解析:
这里只讨论这道题中 D i n k e l B a c h DinkelBach DinkelBach迭代法的好处,题目解析请参考上面的链接。
这道题不管是 D i n k e l B a c h DinkelBach DinkelBach还是二分都能过,但是 D i n k e l B a c h DinkelBach DinkelBach所用的时间不到二分法的 1 / 6 1/6 1/6。
这里讲一下这道题 D i n k e l B a c h DinkelBach DinkelBach怎么做。首先建图还是原来的方案建,跑费用流还是原来的方案跑。
最后怎么求出下一个迭代的答案?
考虑回到流网络里面,每条边有流量说明什么?
与源点和汇点有关系的边里面有流量什么都说明不了,因为我们求网络流的时候肯定会找到一个最大流,源汇的边肯定一直都是流满了的。
然后就是男生女生之间流的边了:
如果男生
i
i
i和女生
j
j
j之间的边有流量,说明他们构成一个匹配,那么他们之间的
a
i
,
j
a_{i,j}
ai,j和
b
i
,
j
b_{i,j}
bi,j就需要加入到最终迭代的统计答案中。
于是就可以开心地迭代下去而不用考虑二分调整上下界的问题啦 (^ _ ^)
代码:
#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const
inline int getint(){
re int num;
re char c;
while(!isdigit(c=gc()));num=c^48;
while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48);
return num;
}
cs int N=203,M=100205;
cs int S=0,T=N-1;
cs int INF=0x3f3f3f3f;
cs double eps=1e-6;
int last[N],nxt[M<<1],to[M<<1],ecnt=1;
int cap[M<<1];
double w[M<<1];
inline void addedge(int u,int v){
nxt[++ecnt]=last[u],last[u]=ecnt,to[ecnt]=v;
nxt[++ecnt]=last[v],last[v]=ecnt,to[ecnt]=u;
}
int cur[N];
double dist[N];
bool vis[N];
queue<int>q;
inline bool SPFA(){
memset(dist,127,sizeof dist);
memset(vis,0,sizeof vis);
q.push(S);cur[S]=last[S];
dist[S]=0;
while(!q.empty()){
int u=q.front();
q.pop();
vis[u]=false;
for(int re e=last[u],v=to[e];e;v=to[e=nxt[e]]){
if(cap[e]&&dist[v]>dist[u]+w[e]){
dist[v]=dist[u]+w[e];
if(!vis[v])vis[v]=true,q.push(v);
cur[v]=last[v];
}
}
}
return dist[T]<INF;
}
inline int dfs(cs int &u,cs int &flow,double &cost){
if(u==T){
cost+=flow*dist[T];
return flow;
}
int ans=0;
vis[u]=true;
for(int &e=cur[u],v=to[e];e;v=to[e=nxt[e]]){
if(cap[e]&&!vis[v]&&fabs(dist[u]+w[e]-dist[v])<eps){
int delta=dfs(v,min(flow-ans,cap[e]),cost);
if(delta){
cap[e]-=delta;
cap[e^1]+=delta;
ans+=delta;
if(ans==flow)return ans;
}
}
}
return ans;
}
inline double MCMF(){
double cost=0;
while(SPFA())dfs(S,INF,cost);
return cost;
}
int n;
int a[N][N];
int b[N][N];
int edge[N][N];
inline void build(){
for(int re i=1;i<=n;++i)
for(int re j=1;j<=n;++j)
addedge(i,j+100),edge[i][j]=ecnt;
for(int re i=1;i<=n;++i)addedge(S,i),addedge(i+100,T);
}
inline double DinkelBach(double kkk){
for(int re i=2;i<=ecnt;i+=2)cap[i]=1,cap[i^1]=0;
for(int re i=1;i<=n;++i)
for(int re j=1;j<=n;++j){
int e=edge[i][j];
w[e]=a[i][j]*1.0-b[i][j]*kkk;
w[e^1]=-w[e];
}
MCMF();
double x=0,y=0;
for(int re i=1;i<=n;++i)
for(int re j=1;j<=n;++j)
if(cap[edge[i][j]])x+=a[i][j],y+=b[i][j];
return x/y;
}
signed main(){
n=getint();
for(int re i=1;i<=n;++i)
for(int re j=1;j<=n;++j)a[i][j]=getint();
for(int re i=1;i<=n;++i)
for(int re j=1;j<=n;++j)b[i][j]=getint();
build();
double L=0,ans=10000.0;
while(fabs(L-ans)>1e-7){
ans=L;
L=DinkelBach(ans);
}
printf("%.6f",ans);
return 0;
}