BZOJ3232: 圈地游戏【最小割解决环限制】

题目描述：

DZY家的后院有一块地，由N行M列的方格组成，格子内种的菜有一定的价值，并且每一条单位长度的格线有一定的费用。
DZY喜欢在地里散步。他总是从任意一个格点出发，沿着格线行走直到回到出发点，且在行走途中不允许与已走过的路线有任何相交或触碰（出发点除外）。记这条封闭路线内部的格子总价值为V，路线上的费用总和为C，DZY想知道V/C的最大值是多少。

题目分析：

首先分数规划二分一个$\text{mid}$，转化为求$V-\text{mid}\ast C$的最大值。

圈住一个格子有收益，走一条线有代价。数据范围提示我们可以用最小割解决。可以把格子看作一个点，与源点相连表示选，流量为格子价值；与汇点相连表示不选，容量为0。问题是怎么解决走出的线能够形成一个环。

分析用最小割解决问题的一般性质，其中一点便是把点分成两类，可以发现如果一条线的两边的格子一个选一个不选，那么这条线就是要选的，也就是说要付出这个格子的代价，所以可以从每个格子向相邻的格子连流量为对应的线*mid的双向边。格子外圈的一圈格子必定不选，所以向汇点连流量为inf的边。

这样建边可以保证选出来的线一定是环，但是可能会割中可能有多个环。注意到答案是分数形式，如果是两个独立的环，那么一定选其中一个比两个都选更优，所以最小割可以保证正确性。

这道题似乎有些卡精度（个人感觉出题人的标程精度可能有点问题，我WA了好几遍，把网上的代码拿来对拍改了改eps才过）。

Code：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 3005
#define maxm 18005
using namespace std;
const double inf = 1e9, eps = 1e-5;
int n,m,S,T,sum;
int fir[maxn],cur[maxn],nxt[maxm],to[maxm],tot=1;
double c[maxm];
inline void line(int x,int y,double z,double rz=0){
    nxt[++tot]=fir[x],fir[x]=tot,to[tot]=y,c[tot]=z;
    nxt[++tot]=fir[y],fir[y]=tot,to[tot]=x,c[tot]=rz;
}
namespace Maxflow{
    int d[maxn],vd[maxn],sz;
    double aug(int u,double augco){
        if(u==T) return augco;
        double need=augco,delta;
        for(int &i=cur[u];i;i=nxt[i]) if(c[i]>eps&&d[u]==d[to[i]]+1){
            delta=aug(to[i],min(need,c[i]));
            c[i]-=delta,c[i^1]+=delta;
            if((need-=delta)<eps||d[S]==sz) return augco-need; 
        }
        cur[u]=fir[u];
        if(!--vd[d[u]]) d[S]=sz;
        else vd[++d[u]]++;
        return augco-need;
    }
    double SAP(){
        memset(d,0,(T+1)<<2),memset(vd,0,(T+1)<<2),memset(cur,0,(T+1)<<2);
        double flow=0; sz=T+1;
        while(d[S]<sz) flow+=aug(S,inf);
        return flow;
    }
}
#define P(i,j) ((i)*(m+2)+(j))
int R[55][55],C[55][55],a[55][55];
double check(double mid){
    memset(fir,0,(T+1)<<2),tot=1;
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) line(S,P(i,j),a[i][j]);
    for(int i=1;i<=n+1;i++) for(int j=1;j<=m;j++) line(P(i-1,j),P(i,j),mid*R[i][j],mid*R[i][j]);
    for(int i=1;i<=n;i++) for(int j=1;j<=m+1;j++) line(P(i,j-1),P(i,j),mid*C[i][j],mid*C[i][j]);
    for(int i=1;i<=n;i++) line(P(i,0),T,inf),line(P(i,m+1),T,inf);
    for(int j=1;j<=m;j++) line(P(0,j),T,inf),line(P(n+1,j),T,inf);
    return sum-Maxflow::SAP();
}
int main()
{
    scanf("%d%d",&n,&m),S=0,T=P(n+1,m+1);
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) scanf("%d",&a[i][j]),sum+=a[i][j];
    for(int i=1;i<=n+1;i++) for(int j=1;j<=m;j++) scanf("%d",&R[i][j]);
    for(int i=1;i<=n;i++) for(int j=1;j<=m+1;j++) scanf("%d",&C[i][j]);
    double l=0,r=sum,mid;
    while(r-l>1e-6){
        mid=(l+r)*0.5;
        if(check(mid)>eps) l=mid;
        else r=mid;
    }
    printf("%.3f\n",l);
}