题目分析:
由于属于一个点的正方形,点必须在该正方形上边或者下边的中间,因为只有两个朝向,所以可以按正方形的朝向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 ;
}