思路:首先将hate和friend建边求其次2-SAT问题,判断是否能有解,没解就输出-1,否则用二分枚举最大的长度,将两个barn的距离小于mid的看做是矛盾,然后建边,求2-SAT问题。找出最优解。

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<vector>
#define Maxn 3010
#define Maxm 1000000
using namespace std;
int dfn[Maxn],low[Maxn],vi[Maxn],head[Maxn],e,n,m,lab,top,Stack[Maxn],num,id[Maxn],A,B,ss;
struct Edge{
int u,v,next,l;
}edge[Maxm];
struct Point{
int x,y;
}p[Maxn],s1,s2;
struct Match{
int a,b;
}hate[Maxn],Friend[Maxn];
void init()
{
int i,j;
memset(dfn,,sizeof(dfn));
memset(low,,sizeof(low));
memset(head,-,sizeof(head));
memset(id,,sizeof(id));
memset(vi,,sizeof(vi));
e=lab=num=top=;
}
void add(int u,int v)
{
edge[e].u=u,edge[e].v=v,edge[e].next=head[u],head[u]=e++;
}
int Dis(Point a,Point b)
{
return abs(a.x-b.x)+abs(a.y-b.y);
}
void Tarjan(int u)
{
int i,j,v;
dfn[u]=low[u]=++lab;
Stack[top++]=u;
vi[u]=;
for(i=head[u];i!=-;i=edge[i].next)
{
v=edge[i].v;
if(!dfn[v])
{
Tarjan(v);
low[u]=min(low[u],low[v]);
}
if(vi[v])
low[u]=min(low[u],dfn[v]);
}
if(low[u]==dfn[u])
{
++num;
do{
i=Stack[--top];
vi[i]=;
id[i]=num;
}while(i!=u);
}
}
void build(int mid)
{
int i,j;
init();
for(i=;i<=n;i++)
for(j=;j<=n;j++){
if(j==i)
continue;
if(Dis(p[i],s1)+Dis(p[j],s1)>mid)
{
add(i,j+n);
add(j,i+n);
}
if(Dis(p[i],s2)+Dis(p[j],s2)>mid)
{
add(i+n,j);
add(j+n,i);
}
if(Dis(p[i],s1)+Dis(p[j],s2)+ss>mid)
{
add(i,j);
add(j+n,i+n);
}
if(Dis(p[i],s2)+Dis(p[j],s1)+ss>mid)
{
add(j,i);
add(i+n,j+n);
}
}
for(i=;i<=A;i++)
{
add(hate[i].a,hate[i].b+n);
add(hate[i].b,hate[i].a+n);
add(hate[i].a+n,hate[i].b);
add(hate[i].b+n,hate[i].a);
}
for(i=;i<=B;i++)
{
add(Friend[i].a,Friend[i].b);
add(Friend[i].a+n,Friend[i].b+n);
add(Friend[i].b,Friend[i].a);
add(Friend[i].b+n,Friend[i].a+n);
}
}
int cal()
{
int i;
for(i=;i<=n*;i++)
{
if(!dfn[i])
Tarjan(i);
}
for(i=;i<=n;i++)
{
if(id[i]==id[i+n])
return ;
}
return ;
}
int solve()
{
int i,j;
for(i=;i<=A;i++)
{
add(hate[i].a,hate[i].b+n);
add(hate[i].b,hate[i].a+n);
add(hate[i].a+n,hate[i].b);
add(hate[i].b+n,hate[i].a);
}
for(i=;i<=B;i++)
{
add(Friend[i].a,Friend[i].b);
add(Friend[i].a+n,Friend[i].b+n);
add(Friend[i].b,Friend[i].a);
add(Friend[i].b+n,Friend[i].a+n);
}
if(!cal())
return -;
int l,r,mid;
l=;r=;
while(l+<r)
{
mid=(l+r)>>;
build(mid);
if(cal())
r=mid;
else
l=mid;
}
return r;
}
int main()
{
int i,j,a,b,c;
while(scanf("%d%d%d",&n,&A,&B)!=EOF)
{
init();
scanf("%d%d%d%d",&s1.x,&s1.y,&s2.x,&s2.y);
for(i=;i<=n;i++)
scanf("%d%d",&p[i].x,&p[i].y);
for(i=;i<=A;i++)
scanf("%d%d",&hate[i].a,&hate[i].b);
for(i=;i<=B;i++)
scanf("%d%d",&Friend[i].a,&Friend[i].b);
ss=Dis(s1,s2);
printf("%d\n",solve());
}
return ;
}

poj 2749 2-SAT问题的更多相关文章

  1. HDU 1815, POJ 2749 Building roads(2-sat)

    HDU 1815, POJ 2749 Building roads pid=1815" target="_blank" style="">题目链 ...

  2. Java实现 POJ 2749 分解因数(计蒜客)

    POJ 2749 分解因数(计蒜客) Description 给出一个正整数a,要求分解成若干个正整数的乘积,即a = a1 * a2 * a3 * - * an,并且1 < a1 <= ...

  3. poj 2749

    Building roads Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 6091   Accepted: 2046 De ...

  4. poj 2749 Building roads (二分+拆点+2-sat)

    Building roads Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 6229   Accepted: 2093 De ...

  5. POJ 2749 2SAT判定+二分

    题意:图上n个点,使每个点都与俩个中转点的其中一个相连(二选一,典型2-sat),并使任意两点最大 距离最小(最大最小,2分答案),有些点相互hata,不能选同一个中转点,有些点相互LOVE,必需选相 ...

  6. poj 1687 Buggy Sat 简单计算几何

    暑期集训出的第一道一血 感觉自己萌萌哒…… 这道题本身并没有坑点 仅仅是翻译巨坑…… 解大腿在做B 安学长在做E 我闲着也没事 就一个词一个词翻译F…… 最后感觉…… 题干大多数都看不懂…… 也都没啥 ...

  7. [poj] 2749 building roads

    原题 2-SAT+二分答案! 最小的最大值,这肯定是二分答案.而我们要2-SATcheck是否在该情况下有可行解. 对于目前的答案limit,首先把爱和恨连边,然后我们n^2枚举每两个点通过判断距离来 ...

  8. POJ 2749 Building roads 2-sat+二分答案

    把爱恨和最大距离视为限制条件,可以知道,最大距离和限制条件多少具有单调性 所以可以二分最大距离,加边+check #include<cstdio> #include<algorith ...

  9. TTTTTTTTTTT POJ 2749 修牛棚 2-Sat + 路径限制 变形

    Building roads Time Limit: 2000MS   Memory Limit: 65536K Total Submissions: 7019   Accepted: 2387 De ...

随机推荐

  1. Oracle的回收站和闪回查询机制(一)

    实际工作中,我们经常会遇到一些情况,误删除某些表或某些表的某些记录,这时候就需要我们将这些记录重新插入进去.如何才能解决这个问题呢? Oracle的Flashback query(闪回查询)为我们解决 ...

  2. [iOS基础控件 - 6.10] Notification 通知机制

    A.定义      iOS程序都有一个NSNotificationCenter的单例对象,用来负责发布不同对象之间的通知      任何对象都能够在NSNotificationCenter发布通知,发 ...

  3. 判断滑动方向UITableView

    CGFloat lastContentOffset;  //ScoreView 滑动位置 -(void)scrollViewWillBeginDragging:(UIScrollView*)scrol ...

  4. Serializable 序列化为文件

    package test; import java.io.File; import java.io.FileInputStream; import java.io.FileNotFoundExcept ...

  5. 图文讲解:iOS App提交流程

    原文:http://www.toceansoft.com/ios/3287.jhtml 一.证书的导出 1.1.前期工作 首先你需要有一个苹果的开发者帐号,一个Mac系统. 如果没有帐号可以在打开ht ...

  6. C语言简单实现sizeof功能代码

    sizeof不是函数,而是运算符,C/C++语言编译器在预编译阶段的时候就已经处理完了sizeof的问题,也就是说sizeof类似于宏定义. 下面给出一个sizeof的一个宏定义实现版本 #defin ...

  7. 数据库MySQL常用命令复习

    -- 查看数据库 show databases; -- 创建数据库 create database '数据库名'; -- 删除数据库 drop database '数据库名'; -- 选库 use ' ...

  8. MySQL几个注意点

    1.在创建表.对表进行操作之前,必须首先选择数据库.通过 mysql_select_db() 函数选取数据库.当您创建 varchar 类型的数据库字段时,必须规定该字段的最大长度,例如:varcha ...

  9. vector 与map的下标操作

    1.vector的下标操作不会添加元素,只能针对已经存在的元素操作. 2.map的下标操作具有副作用,key不存在,会在map中添加一个具有该key的新元素,新元素的value使用默认构造方法. 3. ...

  10. Winform-控制边框显示属性