LA 3211 Now or later

这是一典型的2—SAT问题 ， 关键在于建图 还有就是“最小值尽量大”的二分法求法 。
1、建图：
     这个题目可以转化为求：是否存在一种调度方案 ， 使相邻之间的两个着陆时间差总是不小于p ；
又可以进一步转化为：任意两个着陆时间差不小于p ；
   这样就好做了 ， 把每个飞机的两个着陆时间点 ， 看成是有向图中的两个点 ， 那么假如a的早着陆时间和b的早着陆时间小于p ， 那么如果a 在早着陆时间着陆 ，  那么b 就只能在晚着陆时间着陆 ， 相反也是 ， 因此这样就可以连出两条有向边 ， 而图就是这样建立的。
2、最小值尽量大
     这个问题可以用二分法来做 ， 先求出所给出时间中的最大时间 为maxime， 然后从0 到 maxtime 之间有二分法一次求， 就能得到“最小值尽量大”的值。

代码： 
#include
#include
#include
#include
#include
#include
#include
using namespace std;

const int MAXN = 2100 ;
vectorgrap[MAXN*2] ;
bool mark[MAXN*2] ;
int s[MAXN*2] , c = 0 , n = 0;
int xy[2100][2] ;

void init()
{
    int i ;
    for(i = 0; i < n*2 ; i++)
        grap[i].clear();
    memset(mark , 0 , sizeof(mark));
    c = 0;
}

void add(int x , int a , int y , int b)
{
    x = x*2 + a;
    y = y*2 + b;
    grap[x^1].push_back(y);
    grap[y^1].push_back(x);
}

bool dfs(int x)
{
    if(mark[x^1])  return false;
    if(mark[x])  return true ;
    mark[x] =  true ;
    s[c++] = x;
    for(int i = 0 ; i < grap[x].size() ; i++)
        if(!dfs(grap[x][i]))  return false;
    return true ;
}

bool solve()
{
    for(int i = 0 ; i < n*2 ; i += 2)
        if(!mark[i] && !mark[i+1])
        {
            c = 0;
            if(!dfs(i))
            {
                while( c > 0 ) 
                    mark[s[--c]] = false ;
                if(!dfs(i+1))  return false;
            }
        }
    return true ;
}

bool test(int diff)
{
    init();
    for(int i = 0 ; i < n; i++)
        for(int a = 0 ; a < 2; a++)
            for(int j = i+1 ; j < n; j++)
                for(int b = 0 ; b < 2; b++)
                    if(abs(xy[i][a] - xy[j][b]) < diff)
                    {
                    //    cout<<i<<","<<a<<","<<j<<","<<b<<endl;
                        add(i , a^1 , j , b^1) ;
                    }
    return solve();
}

int main()
{
    while(scanf("%d" , &n) != EOF)
    {
        int l = 0 , r = 0 , i;
        for(i = 0 ; i < n; i++)
            for(int a = 0 ; a < 2; a++)
            {
                scanf("%d" , &xy[i][a]);
                if(xy[i][a] > r)  r = xy[i][a] ;
            }
        while(l < r)
        {
            int m = l + (r-l+1)/2 ;
        //    cout<<r<<endl;
            if(test(m))  l = m;
            else  r = m-1 ;
        //    cout<<m<<endl;
        }
        printf("%d\n" , l);
    }
    return 0;
}