HDOJ 3622 - Bomb Game 2-sat+二分....细心...
题意:
有N个炸弹..每个炸弹有两个位置可以选择..把炸弹放到其中一个地方去...炸弹的爆炸范围是其为圆心的圆...两个炸弹不能有攻击范围上的重合..问要满足条件..炸弹爆炸范围的半径最长能是多少...
题解:
每个炸弹看成一类..其在两个中比选一个..符合2-sat的构图条件....那么就二分枚举炸弹的爆炸范围..枚举相互是否干扰来做边构造2-sat模型...tarjan来判断是否合法..
题目不难..但是我2B了...在初始化vector时..for (i=0;i<(n<<1);i++) T[i].clear() 写成了for (i=0;i<(1<<n);i++) T[i].clear()...搞得我查了好久的错误..还写了几个版本的代码..这道题不好自己出数据..搞得很郁闷..最后怀疑是vector使用的问题..就索性写了个链表..结果就A了~~再反过头来看..发现了错误的根源...弱啊......
Program-----链表存图..整个程序int型计算..结果再开方除2输出
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<cmath>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<algorithm>
#define ll long long
#define oo 1000000007
#define pi acos(-1.0)
#define MAXN 205
#define MAXM 100005
using namespace std;
struct node
{
int x0,y0,x1,y1;
}b[MAXN];
struct LINE
{
int x,y,next;
}line[MAXM];
int Lnum,_next[MAXN],dfn[MAXN],low[MAXN],tp[MAXN],tpnum,DfsIndex,arc[MAXN][MAXN];
bool instack[MAXN];
stack<int> mystack;
int dis(int x0,int y0,int x1,int y1)
{
return (x0-x1)*(x0-x1)+(y0-y1)*(y0-y1);
}
void addline(int x,int y)
{
line[++Lnum].next=_next[x],_next[x]=Lnum;
line[Lnum].x=x,line[Lnum].y=y;
}
void tarjan(int x)
{
int y,k;
dfn[x]=low[x]=++DfsIndex;
instack[x]=true;
mystack.push(x);
for (k=_next[x];k;k=line[k].next)
{
y=line[k].y;
if (!dfn[y])
{
tarjan(y);
low[x]=min(low[x],low[y]);
}else
if (instack[y]) low[x]=min(low[x],dfn[y]);
}
if (dfn[x]==low[x])
{
tpnum++;
do
{
x=mystack.top();
mystack.pop();
instack[x]=false;
tp[x]=tpnum;
}while (dfn[x]!=low[x]);
}
return;
}
bool _2sat(int n,int d)
{
int i,j;
Lnum=0;
memset(_next,0,sizeof(_next));
for (i=0;i<n;i++)
for (j=i+1;j<n;j++)
{
if (arc[i<<1][j<<1]<d)
addline(i<<1,j<<1|1),addline(j<<1,i<<1|1);
if (arc[i<<1][j<<1|1]<d)
addline(i<<1,j<<1),addline(j<<1|1,i<<1|1);
if (arc[i<<1|1][j<<1]<d)
addline(i<<1|1,j<<1|1),addline(j<<1,i<<1);
if (arc[i<<1|1][j<<1|1]<d)
addline(i<<1|1,j<<1),addline(j<<1|1,i<<1);
}
memset(dfn,0,sizeof(dfn));
memset(instack,false,sizeof(instack));
while (!mystack.empty()) mystack.pop();
tpnum=DfsIndex=0;
for (i=0;i<(n<<1);i++)
if (!dfn[i]) tarjan(i);
for (i=0;i<n;i++)
if (tp[i<<1]==tp[i<<1|1]) return false;
return true;
}
int main()
{
int n,i;
while (~scanf("%d",&n))
{
for (i=0;i<n;i++) scanf("%d%d%d%d",&b[i].x0,&b[i].y0,&b[i].x1,&b[i].y1);
for (i=0;i<n;i++)
for (int j=i+1;j<n;j++)
{
arc[i<<1][j<<1]=arc[j<<1][i<<1]=dis(b[i].x0,b[i].y0,b[j].x0,b[j].y0);
arc[i<<1][j<<1|1]=arc[j<<1|1][i<<1]=dis(b[i].x0,b[i].y0,b[j].x1,b[j].y1);
arc[i<<1|1][j<<1]=arc[j<<1][i<<1|1]=dis(b[i].x1,b[i].y1,b[j].x0,b[j].y0);
arc[i<<1|1][j<<1|1]=arc[j<<1|1][i<<1|1]=dis(b[i].x1,b[i].y1,b[j].x1,b[j].y1);
}
int l,r,mid;
l=0,r=2*oo;
while (r-l>1)
{
mid=(r+l)/2;
if (_2sat(n,mid)) l=mid;
else r=mid;
}
printf("%.2lf\n",sqrt(l)*0.5);
}
return 0;
}
Program-----vector存图..其他同上
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<cmath>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<algorithm>
#define ll long long
#define oo 1000000007
#define pi acos(-1.0)
#define MAXN 205
using namespace std;
struct node
{
int x0,y0,x1,y1;
}b[MAXN];
vector<int> T[MAXN];
int dfn[MAXN],low[MAXN],tp[MAXN],tpnum,DfsIndex,arc[MAXN][MAXN];
bool instack[MAXN];
stack<int> mystack;
int dis(int x0,int y0,int x1,int y1)
{
return (x0-x1)*(x0-x1)+(y0-y1)*(y0-y1);
}
void tarjan(int x)
{
int i,y,m=T[x].size();
dfn[x]=low[x]=++DfsIndex;
instack[x]=true;
mystack.push(x);
for (i=0;i<m;i++)
{
y=T[x][i];
if (!dfn[y])
{
tarjan(y);
low[x]=min(low[x],low[y]);
}else
if (instack[y]) low[x]=min(low[x],dfn[y]);
}
if (dfn[x]==low[x])
{
tpnum++;
do
{
x=mystack.top();
mystack.pop();
instack[x]=false;
tp[x]=tpnum;
}while (dfn[x]!=low[x]);
}
return;
}
bool _2sat(int n,int d)
{
int i,j;
for (i=0;i<(n<<1);i++) T[i].clear(); //小错误..搞了N久才发现..
for (i=0;i<n;i++)
for (j=i+1;j<n;j++)
{
if (arc[i<<1][j<<1]<d)
T[i<<1].push_back(j<<1|1),T[j<<1].push_back(i<<1|1);
if (arc[i<<1][j<<1|1]<d)
T[i<<1].push_back(j<<1),T[(j<<1)|1].push_back(i<<1|1);
if (arc[i<<1|1][j<<1]<d)
T[(i<<1)|1].push_back((j<<1)|1),T[j<<1].push_back(i<<1);
if (arc[i<<1|1][j<<1|1]<d)
T[(i<<1)|1].push_back(j<<1),T[(j<<1)|1].push_back(i<<1);
}
memset(dfn,0,sizeof(dfn));
memset(instack,false,sizeof(instack));
while (!mystack.empty()) mystack.pop();
tpnum=DfsIndex=0;
for (i=0;i<(n<<1);i++)
if (!dfn[i]) tarjan(i);
for (i=0;i<n;i++)
if (tp[i<<1]==tp[i<<1|1]) return false;
return true;
}
int main()
{
int n,i;
while (~scanf("%d",&n))
{
for (i=0;i<n;i++) scanf("%d%d%d%d",&b[i].x0,&b[i].y0,&b[i].x1,&b[i].y1);
for (i=0;i<n;i++)
for (int j=i+1;j<n;j++)
{
arc[i<<1][j<<1]=arc[j<<1][i<<1]=dis(b[i].x0,b[i].y0,b[j].x0,b[j].y0);
arc[i<<1][j<<1|1]=arc[j<<1|1][i<<1]=dis(b[i].x0,b[i].y0,b[j].x1,b[j].y1);
arc[i<<1|1][j<<1]=arc[j<<1][i<<1|1]=dis(b[i].x1,b[i].y1,b[j].x0,b[j].y0);
arc[i<<1|1][j<<1|1]=arc[j<<1|1][i<<1|1]=dis(b[i].x1,b[i].y1,b[j].x1,b[j].y1);
}
int l,r,mid;
l=0,r=2*oo;
while (r-l>1)
{
mid=(r+l)>>1;
if (_2sat(n,mid)) l=mid;
else r=mid;
}
printf("%.2f\n",sqrt(l)*0.5);
}
return 0;
}
Program-----vector存图..直接浮点进行运算...
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<cmath>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<algorithm>
#define ll long long
#define eps 1e-5
#define oo 1000000007
#define pi acos(-1.0)
#define MAXN 205
using namespace std;
struct node
{
int x0,y0,x1,y1;
}b[MAXN];
vector<int> T[MAXN];
int dfn[MAXN],low[MAXN],tp[MAXN],tpnum,DfsIndex;
double arc[MAXN][MAXN];
bool instack[MAXN];
stack<int> mystack;
double dis(int x0,int y0,int x1,int y1)
{
return sqrt((x0-x1)*(x0-x1)+(y0-y1)*(y0-y1));
}
void tarjan(int x)
{
int i,y,m=T[x].size();
dfn[x]=low[x]=++DfsIndex;
instack[x]=true;
mystack.push(x);
for (i=0;i<m;i++)
{
y=T[x][i];
if (!dfn[y])
{
tarjan(y);
low[x]=min(low[x],low[y]);
}else
if (instack[y]) low[x]=min(low[x],dfn[y]);
}
if (dfn[x]==low[x])
{
tpnum++;
do
{
x=mystack.top();
mystack.pop();
instack[x]=false;
tp[x]=tpnum;
}while (dfn[x]!=low[x]);
}
return;
}
bool _2sat(int n,double d)
{
int i,j;
for (i=0;i<(n<<1);i++) T[i].clear();
for (i=0;i<n;i++)
for (j=i+1;j<n;j++)
{
if (arc[i<<1][j<<1]+eps<d)
T[i<<1].push_back(j<<1|1),T[j<<1].push_back(i<<1|1);
if (arc[i<<1][j<<1|1]+eps<d)
T[i<<1].push_back(j<<1),T[(j<<1)|1].push_back(i<<1|1);
if (arc[i<<1|1][j<<1]+eps<d)
T[(i<<1)|1].push_back((j<<1)|1),T[j<<1].push_back(i<<1);
if (arc[i<<1|1][j<<1|1]+eps<d)
T[(i<<1)|1].push_back(j<<1),T[(j<<1)|1].push_back(i<<1);
}
memset(dfn,0,sizeof(dfn));
memset(instack,false,sizeof(instack));
while (!mystack.empty()) mystack.pop();
tpnum=DfsIndex=0;
for (i=0;i<(n<<1);i++)
if (!dfn[i]) tarjan(i);
for (i=0;i<n;i++)
if (tp[i<<1]==tp[i<<1|1]) return false;
return true;
}
int main()
{
int n,i;
while (~scanf("%d",&n))
{
for (i=0;i<n;i++) scanf("%d%d%d%d",&b[i].x0,&b[i].y0,&b[i].x1,&b[i].y1);
for (i=0;i<n;i++)
for (int j=i+1;j<n;j++)
{
arc[i<<1][j<<1]=arc[j<<1][i<<1]=dis(b[i].x0,b[i].y0,b[j].x0,b[j].y0);
arc[i<<1][j<<1|1]=arc[j<<1|1][i<<1]=dis(b[i].x0,b[i].y0,b[j].x1,b[j].y1);
arc[i<<1|1][j<<1]=arc[j<<1][i<<1|1]=dis(b[i].x1,b[i].y1,b[j].x0,b[j].y0);
arc[i<<1|1][j<<1|1]=arc[j<<1|1][i<<1|1]=dis(b[i].x1,b[i].y1,b[j].x1,b[j].y1);
}
double l,r,mid;
l=0,r=2*oo;
while (r-l>eps)
{
mid=(r+l)/2;
if (_2sat(n,mid)) l=mid;
else r=mid;
}
printf("%.2lf\n",l*0.5);
}
return 0;
}
HDOJ 3622 - Bomb Game 2-sat+二分....细心...的更多相关文章
- HDU 3622 Bomb Game(2-sat)
HDU 3622 Bomb Game 题目链接 题意:求一个最大半径,使得每一个二元组的点任选一个,能够得到全部圆两两不相交 思路:显然的二分半径,然后2-sat去判定就可以 代码: #include ...
- HDU 3622 Bomb Game（二分+2-SAT）
Bomb Game Time Limit: 10000/3000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total ...
- hdu 3622 Bomb Game(二分+2-SAT)
Bomb Game Time Limit: 10000/3000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total ...
- HDU 3622 Bomb Game(二分+2SAT)
题意:有一个游戏,有n个回合,每回合可以在指定的2个区域之一放炸弹,炸弹范围是一个圈,要求每回合的炸弹范围没有重合.得分是炸弹半径最小的值.求可以得到的最大分数. 思路:二分+2SAT. 二分炸弹范围 ...
- hdu 3622 Bomb Game【二分+2-SAT+tarjan】
用read()会挂 二分半径,显然最优的是所有原都用这个最小半径,然后2-SAT把相交的圆建图,跑tarjan判一下可行性即可 #include<iostream> #include< ...
- HDU - 3622 Bomb Game(二分＋2-SAT)
题目大意:玩一个放炸弹游戏,有N次放炸弹的机会,每次放炸弹时,你都有两个位置能够选择.问怎样放炸弹,能使爆炸的炸弹的半径的最小值最大(炸弹爆炸半径能够控制,可是爆炸形成的圈不能有重叠部分) 解题思路: ...
- 【HDU】3622 Bomb Game（2-SAT）
http://acm.hdu.edu.cn/showproblem.php?pid=3622 又是各种逗.. 2-SAT是一种二元约束,每个点可以置于两种状态,但只能处于一种状态,然后图是否有解就是2 ...
- HDU 3622 Bomb Game（2-sat）
Bomb Game Time Limit: 10000/3000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total ...
- HDU 3622 Bomb Game
Description \(n\) 个炸弹,每个炸弹有两个放置点,可以任选一个,问你最大的半径是多少. Sol 二分+2-SAT+Tarjan. 首先二分一下答案.然后就成了一个2-SAT问题. 建模 ...
随机推荐
- 更改Zend Studio/Eclipse代码风格主题
最近决定把几个IDE的代码样式统一一下,Visual Studio的还算好改,PHP目前用得不多,不过也打算给Zend Studio换身新装. 网上搜索的一些更改Zend Studio主题的多是修改或 ...
- 对不起，说句粗话——这个太屌了，windows1.0安装程序(附下载)
今天逛一个软件论坛发现的,仅仅有几百K.遥想当今我刚接触windows的版本号是3.1,当时记得非常清楚哦,进入windows要从dos命令行进入.如今一转眼,变成进入伪dos是执行栏里敲cmd了.唉 ...
- Linux控制台下的快捷键
Linux控制台(文本模式)下提高工作效率的快捷键 在Linux环境里,有一些按键有特殊的含意.# Ctrl-U: 擦除一行光标前面的部分.# Ctrl-H: 擦除光标前面的一个字符.# Ctrl-D ...
- Html5选择图片并及时预览图片
以往想要实现图片预览基本都是先传至服务器后等返回链接地址才能进行预览,使用Html5选择图片并及时预览图片的代码如下,使用起来更爽了. <!DOCTYPE html> <html l ...
- AS 断点调试 debug
debug面板 点击下图工具栏开启调试会话 此种调试方式是通过冻结应用运行的状态,仿佛时间停止了一般,然后我们逐一观察此时程序的各个参数是否符合我们的预期. 这种调试方法适用于对时间不敏感的程序.也就 ...
- cogs 1008 贪婪大陆
/* 不要思维定视 盯着线段树维护l r 的ans不放 显然没法区间合并 换一种思路 如果打暴力的话 O(nm) 每次询问 扫一遍之前所有的修改 有交点则说明种数++ 接下来考虑如何优化 我们把每个区 ...
- 纯 CSS 创建各种不同的图形形状
使用代码 矩形 .rectangle { width: 250px; height: 150px; background-color: #6DC75F; } <div></div&g ...
- GridView、Repeater获取当前行号
GridView: <%# Container.DataItemIndex+1 %> Repeater:<%# Container.ItemIndex+1%>
- Lesson 7: Responsive Typography In Action
Lesson 7: Responsive Typography In Action 排版一直都是设计和传达的基础.虽然现在的设计和印刷品设计差别很大,但核心原则还是不变的. Article 1: Bo ...
- WARNING OGG-01223 TCP/IP error 111 (Connection refused)
一:问题描述 GGSCI (source_pc) 64> info all Program Status Group Lag at Chkpt Time Sinc ...