4078: [Wf2014]Metal Processing Plant 2-sat+二分+二分图优化

题解：

网上大多数题解都是这样的：从大到小枚举对A集合中最大的边，然后二分B集中最大的边，用2-sat判断是否合法，这样做是比较慢的，需要枚举$n^2$条边，栋老师给出了一种只需要枚举$n$条边的方法：每次枚举一条边时，在图中加入这条边，若这样产生了一个奇环，那么计算完这条边的贡献后就不需要继续枚举，因为我们现在需要图中的任意两点是不能在同一个集合中的，若出现奇环则无法满足；若这样产生了偶环，那么则这条边一定不会产生贡献，因为当前这条边的两个端点在同一个集合中，且为图中权值最小的边，那么偶环上一定会有除这条边以外，两个端点在同一个集合中，权值大于等于它的边。利用这两点，可以将枚举的边数降为$O(n)$。
最后再补充一下如何判断奇偶环：将一个点拆成两个点，设点$v$拆出的两个点分别为$v$和$v'$，对于一条边$(x,y)$，若$x$与$y'$相连，那么连这条边后会有偶环；否则把$x,y$交叉相连，即$x$连$y'$，$x'$连$y$，然后此时若$x$与$x'$相连，则有奇环。
有几个需要注意的地方：
1、二分的下界要设为0，因为B集合有可能没有点。
2、我代码中连了重边，其实是没有必要的。

代码：

#include<bits/stdc++.h>
using namespace std;
#define LL long long
#define pa pair<int,int>
const int Maxn=410;
const int inf=2147483647;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+(ch^48),ch=getchar();
    return x*f;
}
int n;
struct Edge
{
    int x,y,next,v;
    Edge(int _x=0,int _y=0,int _v=0){x=_x,y=_y,v=_v;}
}E[Maxn*Maxn],e[Maxn*Maxn<<3];int le;
bool cmp(Edge a,Edge b){return a.v<b.v;}
int last[Maxn],len;
void ins(int x,int y)
{
    int t=++len;
    e[t].y=y;e[t].next=last[x];last[x]=t;
}
int low[Maxn],dfn[Maxn],id,sta[Maxn],top,bel[Maxn],Cnt;
bool in[Maxn];
void Tarjan(int x)
{
    sta[++top]=x;in[x]=true;
    low[x]=dfn[x]=++id;
    for(int i=last[x];i;i=e[i].next)
    {
        int y=e[i].y;
        if(!dfn[y])Tarjan(y),low[x]=min(low[x],low[y]);
        else if(in[y])low[x]=min(low[x],dfn[y]);
    }
    if(low[x]==dfn[x])
    {
        int i;Cnt++;
        do
        {
            i=sta[top--];
            in[i]=false;
            bel[i]=Cnt;
        }while(i!=x);
    }
}
int fa[Maxn];
int findfa(int x)
{
    if(fa[x]==x)return x;
    return fa[x]=findfa(fa[x]);
}
int ans;
bool check(int p,int v)
{
    id=top=Cnt=len=0;
    memset(last,0,sizeof(last));
    memset(low,0,sizeof(low));
    memset(dfn,0,sizeof(dfn));
    memset(in,false,sizeof(in));
    for(int i=p+1;i<=le;i++)
    {
        int x=E[i].x,y=E[i].y;
        if(E[i].v>v)ins(y<<1,(x<<1)-1),ins(x<<1,(y<<1)-1),ins((y<<1)-1,x<<1),ins((x<<1)-1,y<<1);
        if(E[i].v>E[p].v)ins(y<<1,(x<<1)-1),ins(x<<1,(y<<1)-1);
    }
    for(int i=1;i<=(n<<1);i++)if(!dfn[i])Tarjan(i);
    for(int i=1;i<=n;i++)
    if(bel[i<<1]==bel[(i<<1)-1])return false;
    return true;
}
void solve()
{
    for(int i=1;i<=(n<<1);i++)fa[i]=i;
    sort(E+1,E+1+le,cmp);E[le+1].v=inf;
    for(int i=le;i;i--)
    {
        int x=E[i].x,y=E[i].y;
        int fx1=findfa((x<<1)-1),fx2=findfa(x<<1);
        int fy1=findfa((y<<1)-1),fy2=findfa(y<<1);
        if(fx2==fy1)continue;//偶环
        int l=0,r=i-1;
        while(l<=r)
        {
            int mid=l+r>>1;
            if(check(mid,E[i].v))r=mid-1;
            else l=mid+1;
        }
        ans=min(ans,E[i].v+E[r+1].v);
        if(fx1!=fy1&&fx2!=fy2)fa[fx1]=fy2,fa[fx2]=fy1;
        else return;
    }
}
int main()
{
    while(scanf("%d",&n)!=EOF)
    {
        le=0;
        for(int i=1;i<n;i++)
        for(int j=i+1;j<=n;j++)
        E[++le]=Edge(i,j,read());
        if(n<=2){puts("0");continue;}
        ans=inf;
        sort(E+1,E+1+le,cmp);
        solve();
        printf("%d\n",ans);
    }
}