[随机增量法 最小圆覆盖] BZOJ 1366 [Balkan2002]Alien最小圆覆盖 & 1337 最小圆覆盖 & 2823 [AHOI2012]信号塔

最新推荐文章于 2022-02-06 21:23:22 发布
最新推荐文章于 2022-02-06 21:23:22 发布
阅读量1.2k 收藏
点赞数
分类专栏： 其他算法 随机化
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u014609452/article/details/50636384
版权


论文：顾研《浅谈随机化思想在几何问题中的应用》

以下摘录复杂度分析




期望O(n2)的张角法

实现的常数爆炸 不开O2就作死


#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;

const int N=2005;
const double eps=1e-7;
const double PI=acos(-1.0);

inline int dcmp(double a,double b){
  if (fabs(a-b)<eps) return 0;
  if (a>b) return 1; return -1;
}

inline double sqr(double x){ return x*x; }

struct Point{
  double x,y;
  Point(double x=0,double y=0):x(x),y(y) { }
  void read(){ scanf("%lf%lf",&x,&y); }
  double len(){ return sqrt(x*x+y*y); }
  friend Point operator - (Point A,Point B){ return Point(A.x-B.x,A.y-B.y); }
  friend double operator * (Point A,Point B){ return A.x*B.x+A.y*B.y; }
  friend double Cos(Point A,Point B){
    if (A.len()<eps || B.len()<eps) return -1.0; 
    return A*B/A.len()/B.len();
  }
  friend double Dis(Point A,Point B){ return (A-B).len(); }
  friend double Cross(Point A,Point B){ return A.x*B.y-A.y*B.x; }
}P[N];

inline Point Get(Point a,Point b,Point c){
  double a1,a2,b1,b2,c1,c2;
  Point ans;
  a1=2*(b.x-a.x),b1=2*(b.y-a.y),c1=sqr(b.x)-sqr(a.x)+sqr(b.y)-sqr(a.y);
  a2=2*(c.x-a.x),b2=2*(c.y-a.y),c2=sqr(c.x)-sqr(a.x)+sqr(c.y)-sqr(a.y);
  if(!dcmp(a1,0)){
    ans.y=c1/b1; 
    ans.x=(c2-ans.y*b2)/a2; 
  }else if(!dcmp(b1,0)){ 
    ans.x=c1/a1; 
    ans.y=(c2-ans.x*a2)/b2; 
  }else{ 
    ans.x=(c2*b1-c1*b2)/(a2*b1-a1*b2); 
    ans.y=(c2*a1-c1*a2)/(b2*a1-b1*a2); 
  }
  return ans; 
}

inline Point Get(Point A,Point B){
  return Point((A.x+B.x)/2,(A.y+B.y)/2);
}

int n;

int main(){
  freopen("t.in","r",stdin);
  freopen("t2.out","w",stdout);
  scanf("%d",&n);
  for (int i=1;i<=n;i++) P[i].read();
  random_shuffle(P+1,P+n+1);
  if (n==1){
    printf("%.5lf\n%.5lf %.5lf\n",0,P[1].x,P[1].y);
    return 0;
  }
  Point C,c; double R,r;
  C=Get(P[1],P[2]); R=Dis(P[1],C);
  for (int i=3;i<=n;i++)
    if (dcmp(Dis(P[i],C),R)>0){
      R=1e130;
      for (int j=1;j<i;j++){
	int k1=i,k2=i;
	for (int k=1;k<i;k++){
	  if (k==j) continue;
	  if (Cross(P[j]-P[i],P[k]-P[i])>0){
	    if (!k1 || Cos(P[j]-P[k],P[i]-P[k])>Cos(P[j]-P[k1],P[i]-P[k1]))
	      k1=k;
	  }else{
	    if (!k2 || Cos(P[j]-P[k],P[i]-P[k])>Cos(P[j]-P[k2],P[i]-P[k2]))
	      k2=k;
	  }
	}
	if (dcmp(acos(Cos(P[j]-P[k1],P[i]-P[k1]))+acos(Cos(P[j]-P[k2],P[i]-P[k2])),PI)>=0){
	  c=Get(P[i],P[j],P[k1]),r=Dis(c,P[i]);
	  if (r<R) R=r,C=c;
	  c=Get(P[i],P[j],P[k2]),r=Dis(c,P[i]);
	  if (r<R) R=r,C=c;
	  c=Get(P[i],P[j]),r=Dis(c,P[i]);
	  if (dcmp(Dis(P[k1],c),r)<=0 && dcmp(Dis(P[k2],c),r)<=0)
	    if (r<R)
	      R=r,C=c;
	}
      }
    }
  printf("%.5lf\n%.5lf %.5lf\n",R,C.x+eps,C.y+eps);
  return 0;
}


期望O(n) 棒极了


1336&1337


#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#define eps 1e-8
#define X first
#define Y second
using namespace std;
typedef pair<double,double> Point;

int n;
Point P[100005],C;
double R;

inline double sqr(double x) { return x*x; }

inline int dcmp(double a,double b)
{
	if (fabs(a-b)<eps) return 0;
	if (a<b) return -1; return 1;
}

inline double Dist(Point A,Point B)
{
	return sqrt((A.X-B.X)*(A.X-B.X)+(A.Y-B.Y)*(A.Y-B.Y));
}

inline Point Get(Point a,Point b,Point c)
{
      double a1,a2,b1,b2,c1,c2;
      Point ans;
      a1=2*(b.X-a.X),b1=2*(b.Y-a.Y),c1=sqr(b.X)-sqr(a.X)+sqr(b.Y)-sqr(a.Y);
      a2=2*(c.X-a.X),b2=2*(c.Y-a.Y),c2=sqr(c.X)-sqr(a.X)+sqr(c.Y)-sqr(a.Y);
      if(!dcmp(a1,0)) 
      { 
           ans.Y=c1/b1; 
           ans.X=(c2-ans.Y*b2)/a2; 
      } 
      else if(!dcmp(b1,0)) 
      { 
           ans.X=c1/a1; 
           ans.Y=(c2-ans.X*a2)/b2; 
      } 
      else
      { 
          ans.X=(c2*b1-c1*b2)/(a2*b1-a1*b2); 
          ans.Y=(c2*a1-c1*a2)/(b2*a1-b1*a2); 
      } 
      return ans; 
}

inline Point Get(Point A,Point B)
{
	return Point((A.X+B.X)/2,(A.Y+B.Y)/2);
}

inline void MinCover()
{
	C=P[1]; R=0;
    for(int i=2;i<=n;i++)
    	if (dcmp(Dist(P[i],C),R)>0)
		{
			C=P[i];
			R=0;
        	for (int j=1;j<=i;j++)
				if (dcmp(Dist(P[j],C),R)>0)
				{
					C=Get(P[i],P[j]);
					R=Dist(C,P[i]);
            		for (int k=1;k<=j;k++)
						if (dcmp(Dist(P[k],C),R)>0)
						{
							C=Get(P[i],P[j],P[k]);
							R=Dist(C,P[i]);
						}
				}
		}
}

int main()
{
	freopen("t.in","r",stdin);
	freopen("t.out","w",stdout);
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
		scanf("%lf%lf",&P[i].X,&P[i].Y);
	MinCover();
	printf("%.10lf\n%.10lf %.10lf\n",R,C.X,C.Y);
	return 0;
}

2823

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
#define eps 1e-8
#define X first
#define Y second
using namespace std;
typedef pair<double,double> Point;

int n;
Point P[1000005],C;
double R;

inline double sqr(double x) { return x*x; }

inline int dcmp(double a,double b)
{
	if (fabs(a-b)<eps) return 0;
	if (a<b) return -1; return 1;
}

inline double Dist(Point A,Point B)
{
	return sqrt((A.X-B.X)*(A.X-B.X)+(A.Y-B.Y)*(A.Y-B.Y));
}

inline Point Get(Point a,Point b,Point c)
{
      double a1,a2,b1,b2,c1,c2;
      Point ans;
      a1=2*(b.X-a.X),b1=2*(b.Y-a.Y),c1=sqr(b.X)-sqr(a.X)+sqr(b.Y)-sqr(a.Y);
      a2=2*(c.X-a.X),b2=2*(c.Y-a.Y),c2=sqr(c.X)-sqr(a.X)+sqr(c.Y)-sqr(a.Y);
      if(!dcmp(a1,0)) 
      { 
           ans.Y=c1/b1; 
           ans.X=(c2-ans.Y*b2)/a2; 
      } 
      else if(!dcmp(b1,0)) 
      { 
           ans.X=c1/a1; 
           ans.Y=(c2-ans.X*a2)/b2; 
      } 
      else
      { 
          ans.X=(c2*b1-c1*b2)/(a2*b1-a1*b2); 
          ans.Y=(c2*a1-c1*a2)/(b2*a1-b1*a2); 
      } 
      return ans; 
}

inline Point Get(Point A,Point B)
{
	return Point((A.X+B.X)/2,(A.Y+B.Y)/2);
}

inline void MinCover()
{
	C=P[1]; R=0;
    for(int i=2;i<=n;i++)
    	if (dcmp(Dist(P[i],C),R)>0)
		{
			C=P[i];
			R=0;
        	for (int j=1;j<=i;j++)
				if (dcmp(Dist(P[j],C),R)>0)
				{
					C=Get(P[i],P[j]);
					R=Dist(C,P[i]);
            		for (int k=1;k<=j;k++)
						if (dcmp(Dist(P[k],C),R)>0)
						{
							C=Get(P[i],P[j],P[k]);
							R=Dist(C,P[i]);
						}
				}
		}
}

int main()
{
	freopen("t.in","r",stdin);
	freopen("t.out","w",stdout);
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
		scanf("%lf%lf",&P[i].X,&P[i].Y);
	MinCover();
	printf("%.2lf %.2lf %.2lf\n",C.X,C.Y,R);
	return 0;
}

里阿奴摩西
关注
专栏目录
bzoj1336&1337 最小圆覆盖
AaronPolaris
04-27 4491
随机增量法求最小圆覆盖
【BZOJ2823】【AHOI2012】信号塔 最小圆覆盖 计算几何
辗转山河弋流歌
06-23 2256
题解之前： 首先最小圆覆盖虽然有三层 forfor 循环，但是它是期望 O(n)O(n) 的。什么？你问我为啥？那我只能呵呵了，50W的 O(n3)O(n^3) 高速跑过。 后交的是不求凸包直接跑的，先交的是求了凸包再跑的。。并没有什么差距。 题解： 这道题我们可以先写一份求凸包来缩减点的规模，如果点是随机生成的，那么期望有不到100个点在凸包上，然后就可以乱搞了（其实毛用没有23
最小圆覆盖（随机增量法&模拟退火法）
weixin_30865427的博客
09-18 239
http://acm.hdu.edu.cn/showproblem.php?pid=3007 相关题型连接： http://acm.hdu.edu.cn/showproblem.php?pid=3932 http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=450 Buried memory Time...
O(n)最小圆覆盖算法 （随机增量法）
一位蒟蒻的小博客
01-31 2051
为什么叫随机增量我也不知道 首先考虑一个点集A，设他的最小覆盖圆为g(A)。 g(A)是存在且唯一的。 （假设存在两个，则必定存在更小覆盖圆） 而且g(A)必定满足以下两者之一： 1. 圆上有三个（或以上）A中的点 2. 圆上有两个A中的点并且其连线为直径。 若上述条件不成立则显然有更小覆盖圆。（感性调整） 就像这样： 记K(A)为A中在g(A)上的点集合，这样的点下文叫做关键
【luogu P1742】最小圆覆盖（随机增量法）（计算几何）
欢迎您的光顾awa
02-06 418
给你二维平面上的一些点，要你找到一个半径最小的圆覆盖所有的点。
hdu 3932,3007 最小圆覆盖 随机增量法
_rabbit的专栏
04-22 1032
题目链接： http://acm.hdu.edu.cn/showproblem.php?pid=3932 http://acm.hdu.edu.cn/showproblem.php?pid=3007 都是最小圆覆盖的题目，学习一下随机增量法。 盗用别人的讲解，慢慢理解。 最小覆盖圆， 增量法 假设圆O是前i-1个点得最小覆盖圆，加入第i个点，如果在圆内或边上则什么也不做。否,新得到
【BZOJ 1336&1337】 [Balkan2002]Alien最小圆覆盖
Regina8023
03-23 1407
随机增量法O(n)求最小圆覆盖~
BZOJ1336&1337：[Balkan2002]Alien最小圆覆盖
commonc的博客
08-23 445
题目大意：给定n个点，求最小圆覆盖 最小圆覆盖裸题....算法讲解戳→这里 #include #include #include #include #include #define N 100010 using namespace std; struct node{double x,y;}b[N]; node O; double R; double sqr(double x){
BZOJ 1336 Balkan2002 Alien最小圆覆盖
(╯°口°)╯(┴—┴
02-26 961
题目大意：最小圆覆盖。 思路：再拍一份模板。做法见：http://blog.csdn.net/jiangyuze831/article/details/43950601 CODE： #define _CRT_SECURE_NO_WARNINGS #include #include #include #include #include #include #de
python opencv最小外接矩形_Opencv绘制最小外接矩形、最小外接圆
weixin_33641654的博客
01-12 2511
Opencv中求点集的最小外结矩使用方法minAreaRect，求点集的最小外接圆使用方法minEnclosingCircle。minAreaRect方法原型：RotatedRect minAreaRect( InputArray points );输入参数points是所要求最小外结矩的点集数组或向量；minEnclosingCircle方法原型：void minEnclosingCircle(...
BZOJ1797: [Ahoi2009]Mincut 最小割（洛谷P4126）
forezxl的博客
09-26 203
最小割 Tarjan BZOJ题目传送门 洛谷题目传送门 结论题。。。 跑一遍最小割，对残量网络进行缩点，使得剩下的边都是满流边。 一条边属于最小割的并集当这条边两个端点所在的联通块不同。 一条边属于最小割的交集当这条边两个端点所在的联通块分别与源点和汇点相同。 下面给出口胡： 如果这条边的两个端点在同一个联通块内，那么它就不是满流边，一定没有割它流入的满流边划算，也就不是最小割。 一条边如果一个...
BZOJ 2150. 部落战争（最小路径覆盖问题）【BZOJ千题计划】
繁凡さん的博客
09-16 479
【BZOJ修复计划 #15】BZOJ 2150 部落战争 【国家集训队2011】
BZOJ 1337 最小圆覆盖 随机增量法
-Claude
09-11 2180
BZOJ 1337 最小圆覆盖 随机增量法
[最大团 随机化 || 二分图最大独立集] BZOJ 4080 [Wf2014]Sensor Network
雯舞
06-17 3077
最大团这种数据范围，标算只有被随机艹的节奏了吧正解 抄自这里 我们穷举两个点，这两点距离要小于限制 然后我们分别以这两个点为圆心，两点距离为半径画圆 圆圆相交的部分被两点练成线段划分成两部分，不难发现 每个部分内点点之间的距离是小于限制的，很明显想到二分图 对于上半部分与下半部分的两点，如果距离大于限制则连边 然后我们求最大点独立集即可 #include<cstdi
[三维线性规划 空间解析几何 随机增量法 || 单纯形] HDU 2979 Expensive Drink
雯舞
12-02 1865
论文：顾研《浅谈随机化思想在几何问题中的应用》 以下摘录复杂度分析 #include #include #include #include using namespace std; typedef pair abcd; typedef double Db; const Db oo=1e130; const int N=405; const Db eps=1e-8;
[几何 模拟退火 || 随机增量法] Codeforces 442E #253 (Div. 1) E. Gena and Second Distance
雯舞
03-14 1355
怎么都在理性愉悦啊 总算做到一道有点思路的题了官方题解 题解里说的是什么呢 我们先二分答案AnsAns 然后我们考虑AnsAns合法的条件是存在一个圆OO 他只包含其中一个特殊点 我们肯定可以挪动这个圆 使得有至少一个点在边界上 枚举点pp在边界上 那么圆心OO就在以pp为圆心的圆CC上 圆心OO合法的条件是不存在其他点到他的距离 ≤Ans\leq Ans 实际上就是C的圆周上的不被任
[CDQ分治 并查集 || 随机化权值Xor] BZOJ 3237 [Ahoi2013]连通图
雯舞
05-08 911
考虑CDQ分治 把这半边对后半边没有影响的操作做了 然后分治 用并查集维护 开个栈暴力还原 #include #include using namespace std; inline char nc() { static char buf[100000],*p1=buf,*p2=buf; if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,s
最小费用流算法详解与实现
"该资源提供了一个关于最小费用流问题的参考实现，主要涉及网络流理论、最大流算法和最短增广路算法。题目来源于BZOJ1292，通过残留图和代价矩阵来存储网络信息。" 最小费用流问题是在网络流的基础上引入费用概念，...
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值