bzoj4670: 佛罗里达

这题直接随机化+贪心就可以爆踩过去，我加了个退火增加容错率而已。。。。其实你随机的次数够多根本不需要。。。

然后来自肉丝哥哥的正经做法：

先钦定D(A)>D(B)，那么可以枚举D(A)，然后再去二分D(B)，强行2-sat判就可以做到O(n^4logn)了

再考虑优化，把边按大到小排序，对于当前枚举D(A)枚举到的边，假如它在原图的一个偶环里面，那么在2-sat中一定也会和偶环里的一半点放在放在同一侧。而这些点的边之前已经尝试计算答案了，当前不可能成为最优解。假如处于一个奇环里面，进行完以后2-sat里面一定会有奇环导致不合法。这样就只会有O(n)条边贡献

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
inline double dbrand(){return (double)(rand()%)/10000.0;}

int n,Au,Av,Bu,Bv,ans;
int mp[][],bel[];
int atp,a[],btp,b[],au,av,bu,bv;
int calc()
{
    atp=,btp=;
    for(int i=;i<=n;i++)
        if(bel[i]==)a[++atp]=i;
        else b[++btp]=i;

    int A=-,B=-;
    for(int i=;i<=atp;i++)
        for(int j=i+;j<=atp;j++)
            if(mp[a[i]][a[j]]>A)A=mp[a[i]][a[j]],au=a[i],av=a[j];
    for(int i=;i<=btp;i++)
        for(int j=i+;j<=btp;j++)
            if(mp[b[i]][b[j]]>B)B=mp[b[i]][b[j]],bu=b[i],bv=b[j];
    if(A==-)A=;
    if(B==-)B=;
    if(A+B<ans)ans=A+B;
    return A+B;
}

void annealing()
{
    double T=1e5;
    int num=calc();Au=au,Av=av,Bu=bu,Bv=bv;
    while(T>1e-)
    {
        for(int i=;i<=;i++)
        {
            if(i%==)
            {
                bool bk=false; int cc=;
                for(int p=;p<=;p++)
                {
                    int x;
                         if(p==)x=Au;
                    else if(p==)x=Av;
                    else if(p==)x=Bu;
                    else if(p==)x=Bv;
                    bel[x]^=;
                    int ddd=calc();
                    if(ddd<num)
                        num=ddd,Au=au,Av=av,Bu=bu,Bv=bv;
                    else bel[x]^=,cc++;
                    T*=0.97;
                }
                if(cc==)break;
            }

            int p=rand()%;
            int x;
                 if(p==)x=Au;
            else if(p==)x=Av;
            else if(p==)x=Bu;
            else if(p==)x=Bv;
            bel[x]^=;
            int ddd=calc();
            if(ddd<num)
                num=ddd,Au=au,Av=av,Bu=bu,Bv=bv;
            else bel[x]^=;
            T*=0.97;
        }

        for(int i=;i<=;i++)
        {
            int x=rand()%n+;
            bel[x]^=;
            int ddd=calc();
            int delta=ddd-num;
            if(ddd<num||exp(delta/T)>=dbrand())
                num=ddd,Au=au,Av=av,Bu=bu,Bv=bv;
            else bel[x]^=;
            T*=0.97;
        }
    }
    for(int i=;i<=;i++)
    {
        int x=rand()%n+;
        bel[x]^=;
        int ddd=calc();
        bel[x]^=;
    }
}

int main()
{
    srand();
    while(scanf("%d",&n)!=EOF)
    {
        int mx=,sx;
        for(int i=;i<=n;i++)
            for(int j=i+;j<=n;j++)
            {
                scanf("%d",&mp[i][j]),mp[j][i]=mp[i][j];
                if(mx<mp[i][j])
                {
                    sx=mx;
                    mx=mp[i][j];
                }
                else sx=max(sx,mp[i][j]);
            }

        ans=sx;
        for(int i=;i<=;i++)
        {
            for(int j=;j<=n;j++)bel[j]=rand()%;
            annealing();
        }
        printf("%d\n",ans);
    }

    return ;
}

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
const int maxn=+;
const int maxp=*maxn;
int n;
struct edge{int x,y,d;}e[maxn*maxn];int elen;
bool cmp(edge e1,edge e2){return e1.d>e2.d;}

namespace DISSCUSS
{
    int fa[maxn];
    int findfa(int x)
    {
        if(x==fa[x])return x;
        fa[x]=findfa(fa[x]);return fa[x];
    }
    vector<int>vec[maxn]; int ed,dis;
    void dfs(int x,int fr,int d)
    {
        if(dis!=-)return ;
        if(x==ed){dis=d;return ;}
        for(int i=;i<vec[x].size();i++)
        {
            int y=vec[x][i];
            if(y!=fr)dfs(y,x,d^);
        }
    }
    int getflag(int x,int y)
    {
        int fx=findfa(x),fy=findfa(y);
        if(fx!=fy)
        {
            fa[fx]=fy;
            vec[x].push_back(y);
            vec[y].push_back(x);
            return ;
        }
        else
        {
            dis=-,ed=y;dfs(x,,);
            if(dis==)return ;
            else return ;
        }
    }
    void pre()
    {
        for(int i=;i<=n;i++)
            fa[i]=i,vec[i].clear();
    }
}

namespace TWOSAT
{
    struct node
    {
        int x,y,next;
    }a[**maxn*maxn];int len,last[maxp];
    void ins(int x,int y)
    {
        len++;
        a[len].x=x;a[len].y=y;
        a[len].next=last[x];last[x]=len;
    }
    int z,dfn[maxp],low[maxp];
    int top,sta[maxp];bool v[maxp];
    int cnt,bel[maxp];
    void SCC(int x)
    {
        dfn[x]=low[x]=++z;
        sta[++top]=x;v[x]=true;
        for(int k=last[x];k;k=a[k].next)
        {
            int y=a[k].y;
            if(dfn[y]==)
            {
                SCC(y);
                low[x]=min(low[x],low[y]);
            }
            else if(v[y]==true)
                low[x]=min(low[x],dfn[y]);
        }
        if(dfn[x]==low[x])
        {
            int i; cnt++;
            do
            {
                i=sta[top];top--;
                v[i]=false;
                bel[i]=cnt;
            }while(i!=x);
        }
    }
    void composition(int da,int db)
    {
        len=;memset(last,,sizeof(last));
        for(int i=;i<=elen;i++)
        {
            if(e[i].d>da)ins(e[i].x,e[i].y+n),ins(e[i].y,e[i].x+n);
            if(e[i].d>db)ins(e[i].x+n,e[i].y),ins(e[i].y+n,e[i].x);
        }
    }
    bool main(int da,int db)
    {
        composition(da,db);
        z=top=cnt=;
        memset(dfn,,sizeof(dfn));
        memset(v,false,sizeof(v));
        for(int i=;i<=*n;i++)
            if(dfn[i]==)SCC(i);

        for(int i=;i<=n;i++)
            if(bel[i]==bel[i+n])return false;
        return true;
    }
}

bool check(int da,int db){return TWOSAT::main(da,db);}
int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        elen=;
        for(int i=;i<=n;i++)
            for(int j=i+;j<=n;j++)
                scanf("%d",&e[++elen].d),e[elen].x=i,e[elen].y=j;
        sort(e+,e+elen+,cmp);
        e[elen+].d=;

        DISSCUSS::pre();
        int ans=;
        for(int i=;i<=elen;i++)
        {
            int flag=DISSCUSS::getflag(e[i].x,e[i].y);
            if(flag>=)
            {
                int el=i,er=elen+,j;
                while(el<=er)
                {
                    j=(el+er)/;
                    if(check(e[i].d,e[j].d))
                    {
                        ans=min(ans,e[i].d+e[j].d);
                        el=j+;
                    }
                    else er=j-;
                }
                if(flag==)break;
            }
        }
        printf("%d\n",ans);
    }

    return ;
}