HDU1815 Building roads(二分+2-SAT)
Clever John just had another good idea. He first builds two transferring point S1 and S2, and then builds a road connecting S1 and S2 and N roads connecting each barn with S1 or S2, namely every barn will connect with S1 or S2, but not both. So that every pair of barns will be connected by the roads. To make the cows don't spend too much time while dropping around, John wants to minimize the maximum of distances between every pair of barns.
That's not the whole story because there is another troublesome problem. The cows of some barns hate each other, and John can't connect their barns to the same transferring point. The cows of some barns are friends with each other, and John must connect their barns to the same transferring point. What a headache! Now John turns to you for help. Your task is to find a feasible optimal road-building scheme to make the maximum of distances between every pair of barns as short as possible, which means that you must decide which transferring point each barn should connect to.
We have known the coordinates of S1, S2 and the N barns, the pairs of barns in which the cows hate each other, and the pairs of barns in which the cows are friends with each other.
Note that John always builds roads vertically and horizontally, so the length of road between two places is their Manhattan distance. For example, saying two points with coordinates (x1, y1) and (x2, y2), the Manhattan distance between them is |x1 - x2| + |y1 - y2|.
Next line contains 4 integer sx1, sy1, sx2, sy2, which are the coordinates of two different transferring point S1 and S2 respectively.
Each of the following N line contains two integer x and y. They are coordinates of the barns from the first barn to the last one.
Each of the following A lines contains two different integers i and j(1 <= i < j <= N), which represent the i-th and j-th barns in which the cows hate each other.
The same pair of barns never appears more than once.
Each of the following B lines contains two different integers i and j(1 <= i < j <= N), which represent the i-th and j-th barns in which the cows are friends with each other. The same pair of barns never appears more than once.
You should note that all the coordinates are in the range [-1000000, 1000000].
思路是没有问题的,但是写的有点丑,很多地方还可以合并。一开始思路就没有问题的,但是这种题就是找bug很头疼。QwQ!
结果是maxn开了510,忘记了两倍,晕死了。
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
const int maxn=;
const int B=;
const int R=;
const int W=;
vector<int>G[maxn];
vector<int>G2[maxn];
int dis[][maxn],Dis;//Dis是S点,T点
int x,y,x1,y1,x2,y2,a,b,n,ans;
int col[maxn],q[maxn],num;
void init()
{
for(int i=;i<=*n;i++) G[i].clear();
for(int i=;i<=*n;i++) G2[i].clear();
ans=-;
}
void scan()
{
int i;
scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
Dis=abs(x1-x2)+abs(y1-y2);
for(i=;i<=n;i++){
scanf("%d%d",&x,&y);
dis[][i]=abs(x-x1)+abs(y-y1);
dis[][i]=abs(x-x2)+abs(y-y2);
}
for(i=;i<=a;i++){
scanf("%d%d",&x,&y);
G[x].push_back(y+n);
G[x+n].push_back(y);
G[y].push_back(x+n);
G[y+n].push_back(x);
}
for(i=;i<=b;i++){
scanf("%d%d",&x,&y);
G[x].push_back(y);
G[y].push_back(x);
G[x+n].push_back(y+n);
G[y+n].push_back(x+n);
} }
bool dfs(int u)
{
if(col[u]==R) return false;
if(col[u]==B) return true;
col[u]=B;
col[u>n?u-n:u+n]=R;
q[++num]=u;
for(int i=;i<G[u].size();i++) if(!dfs(G[u][i])) return false;
for(int i=;i<G2[u].size();i++) if(!dfs(G2[u][i])) return false;
return true;
}
bool check(int x)
{
int i,j;
for(i=;i<=n;i++) if(dis[][i]>x&&dis[][i]>x) return false;
for(i=;i<=*n;i++) G2[i].clear();
for(i=;i<=*n;i++) col[i]=;
for(i=;i<=n;i++)
for(j=i+;j<=n;j++){
int d1=dis[][i]+dis[][j];
int d2=dis[][i]+dis[][j]+Dis;
int d3=dis[][i]+dis[][j];
int d4=dis[][i]+dis[][j]+Dis;
if(d1>x&&d2>x&&d3>x&&d4>x) return false;
if(d1>x){
G2[i].push_back(j+n);
G2[j].push_back(i+n);
}
if(d2>x){
G2[i].push_back(j);
G2[j+n].push_back(i+n);
}
if(d3>x){
G2[i+n].push_back(j);
G2[j+n].push_back(i);
}
if(d4>x){
G2[i+n].push_back(j+n);
G2[j].push_back(i);
}
}
for(i=;i<=*n;i++){
if(col[i]) continue;
num=;
if(!dfs(i)){
for(j=;j<=num;j++) {
col[q[j]>n?q[j]-n:q[j]+n]=W;
col[q[j]]=W;
}
if(!dfs(i>n?i-n:i+n)) return false;
}
}
return true;
}
int main()
{
while(~scanf("%d%d%d",&n,&a,&b)){
init();
int L=,R=;
scan();
while(L<=R){
int mid=(L+R)>>;
if(check(mid)){ ans=mid;R=mid-;}
else L=mid+;
}
printf("%d\n",ans);
}
return ;
}
下面这样建图有些问题,读者可以思考一下。
if(d2<=x&&d1>x&&d3>x){
G2[i].push_back(j+n);
G2[j+n].push_back(i);
}
if(d1<=x&&d2>x&&d4>x){
G2[i].push_back(j);
G2[j].push_back(i);
}
if(d4<=x&&d3>x&&d1>x){
G2[i+n].push_back(j);
G2[j].push_back(i+n);
}
if(d3<=x&&d4>x&&d2>x){
G2[i+n].push_back(j+n);
G2[j+n].push_back(i+n);
}
HDU1815 Building roads(二分+2-SAT)的更多相关文章
- POJ Building roads [二分答案 2SAT]
睡觉啦 #include <iostream> #include <cstdio> #include <cstring> #include <algorith ...
- poj 2749 Building roads (二分+拆点+2-sat)
Building roads Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 6229 Accepted: 2093 De ...
- HDU 1815, POJ 2749 Building roads(2-sat)
HDU 1815, POJ 2749 Building roads pid=1815" target="_blank" style="">题目链 ...
- Building roads
Building roads Time Limit: 10000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Tot ...
- [POJ2749]Building roads(2-SAT)
Building roads Time Limit: 2000MS Memory Limit: 65536K Total Submissions: 8153 Accepted: 2772 De ...
- poj 3625 Building Roads
题目连接 http://poj.org/problem?id=3625 Building Roads Description Farmer John had just acquired several ...
- BZOJ 1626: [Usaco2007 Dec]Building Roads 修建道路( MST )
计算距离时平方爆了int结果就WA了一次...... ------------------------------------------------------------------------- ...
- bzoj1626 / P2872 [USACO07DEC]道路建设Building Roads
P2872 [USACO07DEC]道路建设Building Roads kruskal求最小生成树. #include<iostream> #include<cstdio> ...
- bzoj 1626: [Usaco2007 Dec]Building Roads 修建道路 -- 最小生成树
1626: [Usaco2007 Dec]Building Roads 修建道路 Time Limit: 5 Sec Memory Limit: 64 MB Description Farmer J ...
随机推荐
- nfs共享存储
1.下载软件包 yum install nfs-utils nfs-utils-lib -y 2.编辑/etc/exports文件: 1.创建目录:mkdir -p /home/glance2.编辑e ...
- loadrunder之脚本篇——加密解密
密码加密 可以给密码加密,意在把结果字符串作为脚本的参数或者参数值.例如,完整可能有一个用户密码填写的表单,你想测试网站针对不同密码的反应,但是你又想保护密码的安全.Password Encoder允 ...
- Funq之Lambda表达式2
Last month I started a series of posts covering some of the new VB and C# language features that are ...
- 【LeetCode】最大子阵列 Maximum Subarray(贪婪&分治)
描述: Given an integer array nums, find the contiguous subarray (containing at least one number) which ...
- Python编程-继承和接口
一.继承 1.什么是继承 继承是一种创建新类的方式,在python中,新建的类可以继承一个或多个父类,父类又可称为基类或超类,新建的类称为派生类或子类. 继承的好处: 可以使用现有类的所有功能,并在无 ...
- redis 第二篇 系统命令简介 上
一.系统级命令(不归属任何一个数据结构) 从keys命令说起 1. 可以模式匹配 也就是可以用正则表达式的方式来匹配 格式为KEYS pattern glob风格通配符规则 ? 匹配一个字符 ...
- 第四节课-反向传播&&神经网络1
2017-08-14 这节课的主要内容是反向传播的介绍,非常的详细,还有神经网络的部分介绍,比较简短. 首先是对求导,梯度的求解.反向传播的核心就是将函数进行分解,分段求导,前向计算损失,反向计算各个 ...
- php-fpm docker 容器 搭建
继续上一篇文章(centos 7 容器的搭建) 下面构建一个php-fpm镜像: dockerfile 文件如下: # # MAINTAINER # DOCKER-VERSION # # Docker ...
- MapReduce-排序(全部排序、辅助排序)
排序 排序是MapReduce的核心技术. 1.准备 示例:按照气温字段对天气数据集排序.由于气温字段是有符号的整数,所以不能将该字段视为Text对象并以字典顺序排序.反之,用顺序文件存储数据,其In ...
- 批处理之IPC连接工具
由于自己需要重复一些工作,闲来没事发个搞个批处理玩玩. 功能介绍:建立IPC连接,IPC磁盘映射,断开IPC连接等 @echo off color 0A echo Welcome!!! echo 本机 ...