bzoj 4078 [Wf2014]Metal Processing Plant 2SAT 二分图

忍住想骂人和打人的冲动。。。。。。。。。

bz没写数据范围，n<=200
从大到小枚举一块的权值，另一块的权值可以二分。用2SAT验证合法性。

草细节实在不想写了，自己看代码吧。。。

#include <bits/stdc++.h>
using namespace std;
#define N 410
#define M 410000
int n,num,cnt,ans,scc,top,tot;
int v[N][N],p[N][2],fa[N],d[N];
int head[N],nex[M],to[M];
int low[N],deep[N],bel[N],st[N],ins[N];
struct edge
{
    int x,y,v;
    edge(){}
    edge(int x,int y,int v):x(x),y(y),v(v){}
    friend bool operator < (const edge &r1,const edge &r2)
    {return r1.v>r2.v;};
}e[M];
int find(int x)
{
    if(fa[x]==x)return x;
    int ret=find(fa[x]);
    d[x]^=d[fa[x]];
    return fa[x]=ret;
}
void add(int x,int y)
{
    tot++;
    nex[tot]=head[x];head[x]=tot;
    to[tot]=y;
}
void init()
{
    cnt=0;tot=0;scc=0;
    memset(head,0,sizeof(head));
    memset(deep,0,sizeof(deep));
    memset(low,0,sizeof(low));
}
void tarjan(int x)
{
    low[x]=deep[x]=++cnt;
    st[++top]=x;ins[x]=1;
    for(int i=head[x];i;i=nex[i])
    {
        if(!deep[to[i]])
            tarjan(to[i]),low[x]=min(low[x],low[to[i]]);
        else if(ins[to[i]])
            low[x]=min(low[x],deep[to[i]]);
    }
    if(deep[x]==low[x])
    {
        int tmp;
        bel[tmp=st[top--]]=++scc;
        ins[tmp]=0;
        while(tmp!=x)
        {
            bel[tmp=st[top--]]=scc;
            ins[tmp]=0;
        }
    }
}
int check(int x,int y)
{   
    if(x>y)swap(x,y);
    init();
    for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++)
        {
            if(v[i][j]>y)
            {
                add(p[i][0],p[j][1]);
                add(p[j][1],p[i][0]);
                add(p[i][1],p[j][0]);
                add(p[j][0],p[i][1]);
            }
            else if(v[i][j]>x)
            {
                add(p[i][0],p[j][1]);
                add(p[j][0],p[i][1]);
            }
        }
    for(int i=1;i<=p[n][1];i++)
        if(!deep[i])tarjan(i);
    for(int i=1;i<=n;i++)
        if(bel[p[i][0]]==bel[p[i][1]])
            return 0;
    return 1;
}
int search(int x)
{
    int l=0,r=x;
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(check(mid,x))r=mid-1;
        else l=mid+1;
    }
    return l;
}
int main()
{
    scanf("%d",&n);
    if(n<=2)return puts("0"),0;
    for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++)
            scanf("%d",&v[i][j]),e[++num]=edge(i,j,v[i][j]);
    sort(e+1,e+1+num);
    for(int i=1;i<=n;i++)fa[i]=i;
    for(int i=1;i<=n;i++)p[i][0]=++cnt,p[i][1]=++cnt;
    ans=2e9+10;
    for(int i=1;i<=num;i++)
    {
        int x=e[i].x,y=e[i].y,v=e[i].v;
        int x1=find(x),y1=find(y);
        if(x1!=y1)
        {
            ans=min(ans,v+search(v));
            d[x1]=d[x]^d[y]^1;
            fa[x1]=y1;
        }
        else if(d[x]==d[y])
        {
            ans=min(ans,v+search(v));
            break;
        }
    }
    printf("%d\n",ans);
    return 0;
}