bzoj3663/4660CrazyRabbit && bzoj4206最大团
题意
给出平面上N个点的坐标,和一个半径为R的圆心在原点的圆.对于两个点,它们之间有连边,当且仅当它们的连线与圆不相交.求此图的最大团.
点数<=2000,坐标的绝对值和半径<=5000.
注意这里的连线指的是直线而不是线段.保证所有点都在圆外(不会落在圆的边界上).保证任意两个点的连线不会与圆相切.
分析
半年前写过题解.bzoj3663&&bzoj4660 Crazy Rabbit
因为最近联考考了这道题,然后发现我当时的题解写得一塌糊涂...于是重写一遍,顺便给当时看过我的题解的朋友道歉.
(捂脸熊.jpg)
首先对于圆外的点A,我们可以过它作两条圆的切线,找出两个切点之间的圆弧,这代表着从A放一个点光源时光线能照到的圆上的区间(也就是从点A能够"看到"的一段圆弧).我们可以用原点到切点的向量的极角来表示区间的端点.比如从点(0,sqrt(2))向半径为1的单位圆引切线,对应的区间就是[pi/4,3*pi/4].
关于"从某个点能够看到的一段圆弧",我们有一个至关重要的性质:
两个点所连的直线和圆没有交点,当且仅当这两个点分别"看到"的两端圆弧有公共部分但并不相互包含.
为了理解这个性质,可以自己在纸上画个圆,观察以下3条性质.
- 如果两个点"看到"的圆弧没有公共部分,那么它们之间连一条线段会和圆相交.
- 如果两个点"看到"的圆弧完全相互包含,那么它们之间连一条直线会和圆相交.
- 如果两个点"看到"的圆弧不满足(1)(2),那么它们之间连一条直线不和圆相交.
这里本来想画个圆,不过比较懒就没有画
求出"看见"的区间只需要用基本的计算几何和一些三角函数.这一部分细节见代码.
那么问题至此转化为:在圆上选出一些区间使得它们两两有公共部分且不能相互包含.
我们发现,两个区间是否"有公共部分且不相互包含"只与这区间的端点有关,与区间代表的是端点划分出的优弧还是劣弧无关.
上一句话可能比较难理解.我画画图.
为了看起来方便,我们把圆周画得比较有宽度,把一些区间画在圆周内侧和外侧,实际上这些区间都是放在没有宽度的圆周上.
上面的红色区间和绿色区间满足"相交且不包含",但是红色区间和紫色区间也满足,绿色区间和蓝色区间也满足,蓝色区间和紫色区间也满足.
也就是说:我们随便把一些区间变成每个区间在圆上对应的互补区间,"相交且不包含"的关系是不变的.
如果这个结论比较难理解,可以考虑我们是怎样判断两个区间是否"相交且不包含".因为任意两点连线和圆不相切所以不会出现区间端点相同的情况.那么第一个区间的两个端点中必然一个被第二个区间包含,一个不被第二个区间包含.显然,这里我们把哪一个区间变成其在圆上的互补区间都不会对答案有任何影响.
这个性质的用处在于把环上的问题变成链上的问题.我们随便找一个点,从那个地方把圆断开为链(方便起见不如取x轴负半轴和圆的交点).如果有的区间跨过来这个断点,只需把那个区间变成它在圆上的互补区间.
现在问题转化为直线上有一些区间,选出尽量多的区间使得它们两两之间"相交但不包含"
考虑选出来的任意两个区间[l1,r1],[l2,r2],假如l1<l2,那么由于"相交且不包含"的限制,必然有r1<r2.也就是说,把所有的区间按照左端点排序之后,选出来的区间一定是右端点递增的.因此我们naive地猜想,答案是不是对所有区间按照左端点排个序,再按照排序之后的顺序把右端点拿出来跑LIS呢?通过看数据范围我们发现不是这样.因为这样选出来的一些区间虽然满足不会相互包含,却不一定满足两两有公共部分.考虑最终方案中选出来的所有区间中左端点最小的一个,因为区间两两之间有公共部分所以这个区间内必须包含其他所有区间的左端点.那么我们枚举这个左端点最小的区间,只把其他区间中左端点在这个区间内的区间拿出来,按左端点排序后对右端点跑LIS,就可以保证选出来的区间有公共部分而且不包含.复杂度为O(n^2logn).
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
const int maxn=1505;
const double pi=acos(-1);
const double eps=1e-10;
int cmp(double x){
return x<-eps?-1:x>eps;
}
struct point{
double x,y;
point(){}
point(double a,double b){
x=a;y=b;
}
void read(){
scanf("%lf%lf",&x,&y);
}
}P[maxn];
int n,r;
double L[maxn],R[maxn],ori[maxn];
double dot(const point &A,const point &B){
return A.x*B.x+A.y*B.y;
}
point operator -(const point &A,const point &B){
return point(A.x-B.x,A.y-B.y);
}
double length(const point &A){
return sqrt(dot(A,A));
}
double arg(double x,double y){
return atan2(y,x);
}//the angle of vector(x,y),from(-pi,pi] ,arg(-1,0)=pi,arg(1,0)=0 arg(0,1)=pi/2
int tot,seq[maxn];
bool inside(double l,double r,double x){
if(l<r){
return cmp(l-x)*cmp(r-x)==-1;
}else{
return cmp(l-x)*cmp(r-x)==1;
}
}
bool inter(int i,int j){
return inside(L[i],R[i],L[j])^inside(L[i],R[i],R[j]);
}
bool cmp1(const int &a,const int &b){
return cmp(L[a]-L[b])==-1;
}
double m[maxn];
double v[maxn];
int f[maxn];
int binary(double x){
int l=0,r=tot;
while(l<=r){
int mid=(l+r)>>1;
if(m[mid]<x)l=mid+1;
else r=mid-1;
}
return l-1;
}
int lis(){
for(int i=0;i<=tot;++i)m[i]=1e10;
m[0]=-1e10;
for(int i=1;i<=tot;++i)v[i]=R[seq[i]];
int ans=0;
for(int i=1;i<=tot;++i){
f[i]=binary(v[i])+1;
if(f[i]>ans)ans=f[i];
if(v[i]<m[f[i]])m[f[i]]=v[i];
}
return ans;
}
int main(){
int tests;scanf("%d",&tests);
while(tests--){
scanf("%d%d",&n,&r);
for(int i=1;i<=n;++i)P[i].read();
for(int i=1;i<=n;++i)ori[i]=arg(P[i].x,P[i].y);
for(int i=1;i<=n;++i){
double dlt=acos(r/length(P[i]));
L[i]=ori[i]-dlt;R[i]=ori[i]+dlt;
if(L[i]<-pi)L[i]+=2*pi;
if(R[i]>pi)R[i]-=2*pi;
if(L[i]>R[i])swap(L[i],R[i]);
}
// for(int i=1;i<=n;++i)printf("%.3f %.3f\n",L[i],R[i]);
int ans=0;
for(int i=1;i<=n;++i){
//enum start point
tot=0;
for(int j=1;j<=n;++j){
if(j==i)continue;
if(cmp(L[i]-L[j])==-1&&inter(i,j)){
seq[++tot]=j;
}
}
sort(seq+1,seq+tot+1,cmp1);
int tmp=lis()+1;
if(tmp>ans)ans=tmp;
}
printf("%d\n",ans);
}
return 0;
}
bzoj3663/4660CrazyRabbit && bzoj4206最大团的更多相关文章
- 三倍经验——bzoj3663、4660、4206 Crazy Rabbit/最大团
题目描述: 3663 4660 4206 题解: 第一眼:不成立的互相连边,然后用网络流求解无向图最小点覆盖! 好吧我不会. 正解: 每个点对应圆上的一段圆弧,长这样: 设对应圆弧$(l,r)$. 若 ...
- 代码的坏味道(5)——数据泥团(Data Clumps)
坏味道--数据泥团(Data Clumps) 特征 有时,代码的不同部分包含相同的变量组(例如用于连接到数据库的参数).这些绑在一起出现的数据应该拥有自己的对象. 问题原因 通常,数据泥团的出现时因为 ...
- 软件海贼团 OnePiece (版权所有)
最近迷上了“海贼王”这部动画片,不仅仅是因为其中的人物个个性格鲜明,剧情跌宕起伏扣人心弦,各种耍宝搞笑,还感觉到这个团队很像理想中的敏捷软件团队. 作为一直带团队的我,感觉“海贼王”这个动画片给了我很 ...
- Golang友团无闻Go语言Web基础视频教程
教程内容:GO语言资料Golang友团无闻Go语言编程基础Golang友团无闻Go语言Web基础教程 Go语言Web基础教程列表:[Go Web基础]12Go Web 扩展学习.mp4[Go Web基 ...
- 【深度分享】千团大战:看今天商业WiFi乱局及其破解之道
不知道还有没有人记得起始于2010年的千团大战.从2010年初开始,第一家团购网站上线以来,到2011年底,团购网站的数量超过了5000家.当时就有很多媒体预言,2013年,团购的泡沫就将褪去,将有9 ...
- 猿团YTFCloud生态系统,全面服务创业者
9月15日,YTFCloud已正式开启了内测. 创业者翘首以待的YTFCloud,虽然让部分创业者感受到了它的神奇,但对于更多暂时无法尝试的创业者来说,它依然有一层神秘的面纱. 今天小编就来带你近距离 ...
- 猿团YTFCloud--5分钟自制APP,开发从未如此简单
9月15日,YTFCloud将正式开启内测, 这意味着猿团YTF框架产品线全面升级.同时,公测过后,YTFCloud的APP线上DIY服务将面向所有用户,让人人都能成为APP“开发商”. 什么是YTF ...
- zoj1492 最大团
Maximum Clique Time Limit: 10 Seconds Memory Limit: 32768 KB Given a graph G(V, E), a clique is ...
- ecshop 团购点击价格变动
前提:价格阶梯只能设置一级 需要用到: jquery,transport.js(transport_jquery.js),Ajax.call html页面 js代码,还需要插入jquery,trans ...
随机推荐
- 20155320 2016-2017-2 《Java程序设计》第二周学习总结
20155320 2016-2017-2 <Java程序设计>第二周学习总结 教材学习内容总结 本周主要学习了java语言的基础语法,我发现与以前学过的c语言有许多共同点,无论是从一些基本 ...
- [2016北京集训测试赛7]isn-[树状数组+dp+容斥]
Description Solution 定义dp[i][j]为在1到i个数中选了j个数,并且保证选了i的选法总数. dp[i][j]为所有满足A[k]>A[i]的k(k<i)的dp[k] ...
- 北京Uber优步司机奖励政策(4月17日)
滴快车单单2.5倍,注册地址:http://www.udache.com/ 如何注册Uber司机(全国版最新最详细注册流程)/月入2万/不用抢单:http://www.cnblogs.com/mfry ...
- Python socket网络模块
一.基于TCP协议的socket通信 以打电话为理解方式进行TCP的通信. Server端代码: import socket phone = socket.socket(socket.AF_INET, ...
- apache-kylin-2.5.2-bin-cdh57与cdh-5.13.0集群整合运用
1.下载kylin最新版apache-kylin-2.5.2-bin-cdh57: 2.解压配置环境变量: export BASE_PATH="/opt/cloudera/parcels/C ...
- python根据正则表达式的简单爬虫
今天根据正则表达式简单的爬了一下大众点评,把北京的美食爬了爬,(店铺名,人均消费,地址) import re import urllib.request from urllib.request imp ...
- Windows下MongoDB优化及问题处理
1.MongoDB 服务器CPU占用100% 给Mongodb对应数据库中的表建立索引,这里我采用使用工具:NoSQL Manager for MongoDB 直接在表的属性栏,选择Indexes,右 ...
- Kibana TypeError : Object #<GlobalState> has no method 'setDefaults'
在windows server中装完elasticsearch和kibana后,elasticsearch能正常访问(http://localhost:9200): 而访问kibana的地址(http ...
- Vue 项目在其他电脑 npm run dev 运行报错的解决方法
一个 Vue 项目从一台电脑上传到 github 上之后,再另外一台电脑上 git clone .并使用 npm run dev 或 npm run start 发生以下报错的解决方法. 报错原因 ...
- 《机器学习实战》6.2小节,KKT条件代码理解
<机器学习实战>6.2小节 #这句是检测 当前样本点i 是否满足KKT条件的 if (alphas[i, :] < C and E_i * labelMat[i, :] < - ...