BZOJ.1312.[Neerc2006]Hard Life(分数规划 最大权闭合子图)

BZOJ

最大密度子图。

二分答案\(x\)，转为求是否存在方案满足：\(边数-x*点数\geq 0\)。

选一条边就必须选两个点，所以可以转成最大权闭合子图。边有\(1\)的正权，点有\(x\)的负权。判断\(边数-最小割\)是否非负即可。

有一个结论是，任意两个密度子图，它们的密度差不超过\(\frac{1}{n^2}\)。
所以拿eps=1e-7或者更小做二分边界不对。。。
必须是\(while(l+1.0/n/n<=r)\)。

还要注意精度的问题。。
m=0要输出1。

//1300kb    236ms
#include <cstdio>
#include <cctype>
#include <cstring>
#include <algorithm>
#define gc() getchar()
#define eps 1e-8
const int N=2005,M=6005+205;
const double INF=1ll<<55;

int n,m,src,des,Ans,A[N],B[N],Enum,H[N],nxt[M],fr[M],to[M],lev[N],pre[N];
double cap[M];

inline int read()
{
    int now=0;register char c=gc();
    for(;!isdigit(c);c=gc());
    for(;isdigit(c);now=now*10+c-'0',c=gc());
    return now;
}
inline void AE(int u,int v,double w)
{
    to[++Enum]=v, fr[Enum]=u, nxt[Enum]=H[u], H[u]=Enum, cap[Enum]=w;
    to[++Enum]=u, fr[Enum]=v, nxt[Enum]=H[v], H[v]=Enum, cap[Enum]=0;
}
bool BFS()
{
    static int q[N];
    for(int i=0; i<des; ++i) lev[i]=des+1;
    int h=0,t=1; q[0]=des, lev[des]=0;
    while(h<t)
    {
        int x=q[h++];
        for(int i=H[x]; i; i=nxt[i])
            if(cap[i^1]>=eps && lev[to[i]]==des+1)
                lev[to[i]]=lev[x]+1, q[t++]=to[i];
    }
    return lev[src]<=des;
}
inline double Augment()
{
    double mn=INF;
    for(int i=des; i; i=fr[pre[i]])
        mn=std::min(mn,cap[pre[i]]);
    for(int i=des; i; i=fr[pre[i]])
        cap[pre[i]]-=mn, cap[pre[i]^1]+=mn;
    return mn;
}
double ISAP()
{
    static int cur[N],num[N];
    if(!BFS()) return 0;
    for(int i=0; i<=des; ++i) cur[i]=H[i], ++num[lev[i]];
    int x=0; double res=0;
    while(lev[0]<=des)
    {
        if(x==des) x=0, res+=Augment();
        bool can=0;
        for(int i=cur[x]; i; i=nxt[i])
            if(lev[to[i]]==lev[x]-1 && cap[i]>=eps)
            {
                can=1, cur[x]=i, pre[x=to[i]]=i;
                break;
            }
        if(!can)
        {
            int mn=des;
            for(int i=H[x]; i; i=nxt[i])
                if(cap[i]>=eps) mn=std::min(mn,lev[to[i]]);
            if(!--num[lev[x]]) break;
            ++num[lev[x]=mn+1], cur[x]=H[x];
            if(x) x=fr[pre[x]];
        }
    }
    return res;
}
bool Check(double x)
{
    Enum=1, memset(H,0,des+1<<2);
    for(int i=1; i<=m; ++i) AE(0,i+n,1), AE(i+n,A[i],INF), AE(i+n,B[i],INF);
    for(int i=1; i<=n; ++i) AE(i,des,x);
    return m-ISAP()>=eps;
}
void DFS(int x)
{
    static bool vis[N];
    vis[x]=1, Ans+=(x<=n);
    for(int i=H[x]; i; i=nxt[i])
        if(cap[i]>=eps && !vis[to[i]]) DFS(to[i]);
}

int main()
{
    n=read(),m=read(),src=0,des=n+m+1;
    if(!m) return puts("1"),0;
    for(int i=1; i<=m; ++i) A[i]=read(),B[i]=read();
    double l=0.49,r=m/2.0,mid,EPS=1.0/n/n;//l不能设0.5。虽然最优比率最小是0.5，但是因为神奇的浮点误差0.5做最优比率并不对（0.49999999403953才对）
    while(l+EPS<r)
        if(Check(mid=(l+r)*0.5)) l=mid;
        else r=mid;
    Check(l), DFS(src);
    printf("%d\n",Ans-1);

    return 0;
}

转载于:https://www.cnblogs.com/SovietPower/p/10113344.html