【POJ】2296 Map Labeler 2-sat

传送门：【POJ】2296 Map Labeler

题目分析：

由于属于一个点的正方形，点必须在该正方形上边或者下边的中间，因为只有两个朝向，所以可以按正方形的朝向2-sat。

我们只需要找冲突然后建边。

二分正方形的边长mid，

对于两个点i、j，如果abs(xi-xj)>=mid，则一定怎么放正方形都不会冲突。

如果abs(yi-yj)>=2*mid，同样不会冲突。

所以我们只需要考虑abs(xi-xj)<mid且abs(yi-yj)<2*mid时的情况。

当abs(yi-yj)>=mid时：

yi<yj:<i,j>,<~i,~j>，同向且朝上

yj<yi:<j,i>,<~j,~i>，同向且朝下

当abs(yi-yj)<mid时：

yi<yj:<i,~i>,<~j,j>，i必须朝下，j必须朝上

yj<yi:<~i,i>,<j,~j>，i必须朝上，j必须朝下

yi=yj:<i,~j>,<j,~i>,<~i,j>,<~j,i>，i，j朝向必须不同

代码如下：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;

#define REP( i , a , b ) for ( int i = ( a ) ; i < ( b ) ; ++ i )
#define FOR( i , a , b ) for ( int i = ( a ) ; i <= ( b ) ; ++ i )
#define REV( i , a , b ) for ( int i = ( a ) ; i >= ( b ) ; -- i )
#define CLR( a , x ) memset ( a , x , sizeof a )

const int MAXN = 205 ;
const int MAXE = 40000 ;

struct Edge {
	int v ;
	Edge* next ;
} ;

struct Line {
	int x , y ;
} ;

Edge E[MAXE] , *H[MAXN] , *cur ;
int dfn[MAXN] , low[MAXN] , scc[MAXN] , scc_cnt ;
int S[MAXN] , top , dfs_clock ;
int n ;

Line l[MAXN] ;

void init () {
	cur = E ;
	top = 0 ;
	scc_cnt = 0 ;
	dfs_clock = 0 ;
	CLR ( H , 0 ) ;
	CLR ( dfn , 0 ) ;
	CLR ( scc , 0 ) ;
}

void addedge ( int u , int v ) {
	cur -> v = v ;
	cur -> next = H[u] ;
	H[u] = cur ++ ;
}

void tarjan ( int u ) {
	dfn[u] = low[u] = ++ dfs_clock ;
	S[top ++] = u ;
	for ( Edge* e = H[u] ; e ; e = e -> next ) {
		int v = e -> v ;
		if ( !dfn[v] ) {
			tarjan ( v ) ;
			low[u] = min ( low[u] , low[v] ) ;
		} else if ( !scc[v] ) low[u] = min ( low[u] , dfn[v] ) ;
	}
	if ( low[u] == dfn[u] ) {
		++ scc_cnt ;
		do {
			scc[S[-- top]] = scc_cnt ;
		} while ( u != S[top] ) ;
	}
}

/*
void scanf ( int& x , char c = 0 ) {
	while ( ( c = getchar () ) < '0' || c > '9' ) ;
	x = c - '0' ;
	while ( ( c = getchar () ) >= '0' && c <= '9' ) x = x * 10 + c - '0' ;
}
*/

bool ok () {
	REP ( i , 0 , n << 1 ) if ( !dfn[i] ) tarjan ( i ) ;
	REP ( i , 0 , n ) if ( scc[i << 1] == scc[i << 1 | 1] ) return 0 ;
	return 1 ;
}

void solve () {
	scanf ( "%d" , &n ) ;
	REP ( i , 0 , n ) scanf ( "%d%d" , &l[i].x , &l[i].y ) ;
	int low = 0 , high = 20000 ;
	while ( low < high ) {
		int m = ( low + high ) >> 1 ;
		init () ;
		REP ( i , 0 , n ) REP ( j , i + 1 , n ) {
			if ( abs ( l[i].x - l[j].x ) >= m ) continue ;
			if ( abs ( l[i].y - l[j].y ) >= 2 * m ) continue ;
			if ( abs ( l[i].y - l[j].y ) >= m ) {
				if ( l[i].y < l[j].y ) {
					addedge ( i << 1 , j << 1 ) ;
					addedge ( j << 1 | 1 , i << 1 | 1 ) ;
				} else {
					addedge ( i << 1 | 1 , j << 1 | 1 ) ;
					addedge ( j << 1 , i << 1 ) ;
				}
			} else {
				if ( l[i].y < l[j].y ) {
					addedge ( i << 1 , i << 1 | 1 ) ;
					addedge ( j << 1 | 1 , j << 1 ) ;
				} else if ( l[i].y > l[j].y ) {
					addedge ( i << 1 | 1 , i << 1 ) ;
					addedge ( j << 1 , j << 1 | 1 ) ;
				} else {
					addedge ( i << 1 , j << 1 | 1 ) ;
					addedge ( j << 1 , i << 1 | 1 ) ;
					addedge ( i << 1 | 1 , j << 1 ) ;
					addedge ( j << 1 | 1 , i << 1 ) ;
				}
			}
		}
		if ( ok () ) low = m + 1 ;
		else high = m ;
	}
	printf ( "%d\n" , low - 1 ) ;
}

int main () {
	int T ;
	scanf ( "%d" , &T ) ;
	while ( T -- ) solve () ;
	return 0 ;
}