BZOJ4078 WF2014Metal Processing Plant（二分答案+2-SAT）

　　题面甚至没给范围，由数据可得n<=200。容易想到二分答案，暴力枚举某集合的价值，2-SATcheck一下即可。这样是O(n⁴logn)的。

　　2-SAT复杂度已经是下界，考虑如何优化枚举。稍微改一下，不妨从大到小枚举较大的集合价值（即枚举边），另一个集合二分答案，同样O(n⁴logn)。

　　看起来没什么卵用。但注意到较大集合所不能包含的边不可以成奇环，否则肯定有一条环上边被选中。那么考虑当前边，如果形成奇环，最大值不可能比它更小了，做完这个就可以退出；如果加上这条边后形成偶环，可以直接跳过，因为如果要让其被选中，相当于将偶环缩成了奇环，剩余边肯定有一条被选中而又比它大。不形成环时则直接做。由于树边数量O(n)，复杂度O(n³logn)。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define ll long long
#define N 210
#define inf 1000000010
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
int n,a[N][N],b[N*N],dfn[N<<],low[N<<],stk[N<<],belong[N<<],p[N<<],fa[N<<],deep[N<<],size[N<<],len[N<<],t,top,tot,cnt,ans=inf,m;
bool flag[N<<];
struct data{int to,nxt;
}edge[N*N<<];
struct data2
{
    int x,y,z;
    bool operator <(const data2&a) const
    {
        return z<a.z;
    }
}e[N*N];
void addedge(int x,int y){t++;edge[t].to=y,edge[t].nxt=p[x],p[x]=t;}
void tarjan(int k)
{
    dfn[k]=low[k]=++cnt;
    stk[++top]=k;flag[k]=;
    for (int i=p[k];i;i=edge[i].nxt)
    if (!dfn[edge[i].to]) tarjan(edge[i].to),low[k]=min(low[k],low[edge[i].to]);
    else if (flag[edge[i].to]) low[k]=min(low[k],dfn[edge[i].to]);
    if (dfn[k]==low[k])
    {
        tot++;
        while (stk[top]!=k)
        {
            flag[stk[top]]=;
            belong[stk[top]]=tot;
            top--;
        }
        flag[k]=,belong[k]=tot,top--;
    }
}
bool check(int x,int y)
{
       memset(dfn,,sizeof(dfn));
       memset(p,,sizeof(p));t=tot=cnt=;
       for (int i=;i<n;i++)
        for (int j=i+;j<=n;j++)
        {
            if (a[i][j]>x) addedge(i*-,j*),addedge(j*-,i*);
            if (a[i][j]>y) addedge(i*,j*-),addedge(j*,i*-);
        }
    for (int i=;i<=n*;i++)
       if (!dfn[i]) tarjan(i);
    for (int i=;i<=n;i++) if (belong[i*-]==belong[i*]) return ;
    return ;
}
int find(int x)
{
    if (fa[x]==x) return x;
    int f=find(fa[x]);
    deep[x]=deep[fa[x]]^len[x];
    return f;
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj4078.in","r",stdin);
    freopen("bzoj4078.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read();
    for (int i=;i<n;i++)
        for (int j=i+;j<=n;j++)
        {
            a[i][j]=a[j][i]=read();
            m++;e[m].z=a[i][j],e[m].x=i,e[m].y=j;
        }
    sort(e+,e+m+);
    for (int i=;i<=n*;i++) fa[i]=i,size[i]=;
    for (int i=m;i>=;i--)
    {
        int p=find(e[i].x),q=find(e[i].y);
        if (size[p]<size[q]) swap(p,q);
        if (p!=q||(deep[e[i].x]^deep[e[i].y]^))
        {
            if (p!=q) fa[q]=p,size[p]+=size[q],len[q]=deep[e[i].x]^deep[e[i].y]^;
               int l=,r=i,u=inf;
            while (l<=r)
            {
                int mid=l+r>>;
                if (check(e[i].z,e[mid].z)) u=e[mid].z,r=mid-;
                else l=mid+;
            }
            ans=min(ans,e[i].z+u);
        }
        if (p==q&&(deep[e[i].x]^deep[e[i].y]^)) break;
    }
    if (ans==inf) ans=;cout<<ans;
    return ;
}