poj3675 求圆与多边形相交的面积_matlab 直线和多边形的交点 intersectwithline-程序员宅基地

技术标签: 计算几何  二维点线面  

每条线段与原点连起来求在圆内的有向面积,然后输出之前绝对值一下。

求三角形在圆内的面积,都在圆内,直接三角形面积,都在圆外且无交点,直接扇形面积,如果一点在圆内一点在圆外的话,求一个扇形一个三角形,两点在圆外但有2个交点,则求2个扇形中间一个三角形的面积。

使用点积解方程求线段与圆的交点过不了这题,不知道为何,但用上交小红书的解方程模板就可以

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#define maxl 510
#define eps 1e-10
const double pi=acos(-1.0);
using namespace std;

inline int sgn(double x)
{
	if(x>-eps && x<eps) return 0;
	if(x>0) return 1;
	else 	return -1;
}
inline double sqr(double x)
{
	return x*x;
}

struct point
{
	double x,y;
	point(double a=0,double b=0)
	{
		x=a;y=b;
	}
	point operator + (const point &b)const
	{
		return point(x+b.x,y+b.y);
	}
	point operator - (const point &b)const
	{
		return point(x-b.x,y-b.y);
	}
	friend point operator * (const double t,const point &p)
	{
		return point(t*p.x,t*p.y);
	}
	inline double norm()
	{
		return sqrt(x*x+y*y);
	}
	inline void transxy(double b)//逆时针旋转b弧度
	{
		double tx=x,ty=y;
		x=tx*cos(b)-ty*sin(b);
		y=tx*sin(b)+ty*cos(b);
	}
	inline void transxy(double cosb,double sinb)
	{
		double tx=x,ty=y;
		x=tx*cosb-ty*sinb;
		y=tx*sinb+ty*cosb;
	}
};
inline double dot(const point &a,const point &b)
{
	return a.x*b.x+a.y*b.y;
}
inline double det(const point &a,const point &b)
{
	return a.x*b.y-a.y*b.x;
}
inline double dist(const point &a,const point &b)
{
	return (b-a).norm();
}
struct line
{
	point s,e;
	double k;
	line(point a=point(),point b=point())
	{
		s=a;e=b;
		k=atan2(e.y-s.y,e.x-s.x);
	}
};
inline bool point_on_seg(point &p,line &a)
{
	return sgn(det(p-a.s,p-a.e))==0 && sgn(dot(p-a.s,p-a.e))<=0;
}
inline double mysqrt(double x)
{
	return sqrt(max(x,0.0));
}
/*inline void circle_cross_line(line L,point o,double r,point ret[],int &num)
{
	point v=L.e-L.s;
	double t=dot(o-L.s,v)/dot(v,v);
	point rt=L.s+t*v;
	double dis=dist(rt,o);
	num=0; 
	if(sgn(dis-r)>0) 
		return;
	if(sgn(dis-r)==0)
	{
		if(point_on_seg(rt,L)) ret[num++]=rt;
		return;
	}
	double cosval=dis/r,sinval=mysqrt(1-cosval*cosval),sita;
	if(cosval<-1.0)
		cosval+=eps;
	if(cosval>1.0)
		cosval-=eps; 
	sita=acos(cosval);
	point p1=rt,p2=rt;
	//p1.transxy(cosval,sinval);p2.transxy(cosval,-sinval);
	p1.transxy(sita);p2.transxy(-sita);
	p1=(r/dis)*p1;p2=(r/dis)*p2;
	if(point_on_seg(p1,L)) 
		ret[num++]=p1;
	if(point_on_seg(p2,L))
		ret[num++]=p2;
}*/

inline void circle_cross_line(line L,point o,double r,point ret[],int &num)
{
	point a=L.s;point b=L.e;
	double x0=o.x,y0=o.y;
	double x1=a.x,y1=a.y;
	double x2=b.x,y2=b.y;
	double dx=x2-x1,dy=y2-y1;
	double A=dx*dx+dy*dy;
	double B=2*dx*(x1-x0)+2*dy*(y1-y0);
	double C=sqr(x1-x0)+sqr(y1-y0)-sqr(r);
	double delta=B*B-4*A*C;
	num=0;
 	if(sgn(delta)>=0)
    {
		double t1=(-B-mysqrt(delta))/(2*A);
		double t2=(-B+mysqrt(delta))/(2*A);
     	if(sgn(t1-1)<=0&&sgn(t1)>=0)
        	ret[num++]=point(x1+t1*dx,y1+t1*dy);
    	if(sgn(t2-1)<=0&&sgn(t2)>=0)
        	ret[num++]=point(x1+t2*dx,y1+t2*dy);
    }
}

double R,ans;
int n;
point a[maxl];

inline void prework()
{
	scanf("%d",&n);
	for(int i=0;i<n;i++)
		scanf("%lf%lf",&a[i].x,&a[i].y);
	a[n]=a[0];
}

double sector_area(const point &a,const point &b)
{
	double theta=atan2(a.y,a.x)-atan2(b.y,b.x);
	while(theta<=0) theta+=2*pi;
	while(theta>2*pi) theta-=2*pi;
	theta=min(theta,2*pi-theta);
	return R*R*theta/2;
}

inline double calc(const point &a,const point &b)
{
	point p[2];int num=0;
	bool ina=sgn(sqrt(dot(a,a))-R)<0,inb=sgn(sqrt(dot(b,b))-R)<0;
	if(ina)
	{
		if(inb) return fabs(det(a,b))/2;
		else
		{
			circle_cross_line(line(a,b),point(0,0),R,p,num);
            return fabs(sector_area(b,p[0]))+fabs(det(a,p[0]))/2.0;
		}
	}
	else
	{
		if(inb)
		{
			circle_cross_line(line(a,b),point(0,0),R,p,num);
			return sector_area(p[0],a)+fabs(det(p[0],b))/2;
		}
		else
		{
			circle_cross_line(line(a,b),point(0,0),R,p,num);
			if(num==2)
				return sector_area(a,p[0])+sector_area(p[1],b)+fabs(det(p[0],p[1]))/2.0;
			else
				return sector_area(a,b);
		}
	}
}

inline double area()
{
	double res=0;
	for(int i=0;i<n;i++)
	{
		int f=sgn(det(a[i],a[i+1]));
		if(f!=0)
			res+=f*(calc(a[i],a[i+1]));
	}
	return res;
}

inline void mainwork()
{
	ans=fabs(area())+eps;
}

inline void print()
{
	printf("%.2f\n",ans);
}

int main()
{
	while(~scanf("%lf",&R))
	{
		prework();
		mainwork();
		print();
	}	
	return 0;
}

 

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/liufengwei1/article/details/87821200

智能推荐

css的mediaquery查询几种不生效的原因_media- query 不生效-程序员宅基地

文章浏览阅读6.3k次。今天上午搞了半个小时,在webstrom中修改了css样式,用chrom浏览器浏览不生效,在网上查了半天,看是不是webstrom的有写保护,或者chorm浏览器不支持,最后发现少and后面少了一个空格导致css mediaquery查询不生效,即使在chrom中的network中css样式文件已经被加载了。@media screen and 不起作用原因汇总。首先确认是不是css本身的问题,而不..._media- query 不生效

Arduino 语法、函数使用、最常用函数、arduino教程、arduino笔记、参考手册_arduino常用函数手册-程序员宅基地

文章浏览阅读1w次,点赞38次,收藏194次。函数部分数字 I/OpinMode()描述将指定的引脚配置成输出或输入。详情请见digital pins。语法pinMode(pin, mode)参数pin:要设置模式的引脚mode:INPUT或OUTPUT返回无例子ledPin = 13 // LED连接到数字脚13 void setup() { pinMode(ledPin,OUTPUT);//设置数字脚为输..._arduino常用函数手册

MySQL数据库概览_后,我们首先就来看一下,mysql数据库中有多少个数据库(你可以理解为,mysql数据库是-程序员宅基地

文章浏览阅读257次。SQL语句分类:DCL: 数据控制语言 GRANT、DENY、REVOKEDDD: 数据定义语言 CREATE、ALTER、DROP、TRUNCATEDML: 数据操纵语言 SELECT、DELETE、UPDATE、SELECTMysql默认数据库有4个(version 5.7及8)mysql: mysql的核心数据库,类似于sql server中的master表,主要负责存储数据库的用户、权限设置、关键字等mysql自己需要使用的控制和管理信息。(常用的,在mysql.user表中_后,我们首先就来看一下,mysql数据库中有多少个数据库(你可以理解为,mysql数据库是

H3C S5500-52C-EI SSH 服务器发送了断开连接数据包_服务器发送了断开连接的数据包-程序员宅基地

文章浏览阅读2.9k次。H3C S5500-52C-EI SSH 服务器发送了断开连接数据包The connection is closed by SSH Server.(code:2)_服务器发送了断开连接的数据包

iOS-底层原理 21:Method-Swizzling 方法交换_ios方法交换 父类 死循环-程序员宅基地

文章浏览阅读1.3k次。iOS 底层原理 文章汇总method-swizzling 是什么?method-swizzling的含义是方法交换,其主要作用是在运行时将一个方法的实现替换成另一个方法的实现,这就是我们常说的iOS黑魔法,在OC中就是利用method-swizzling实现AOP,其中AOP(Aspect Oriented Programming,面向切面编程)是一种编程的思想,区别于OOP(面向对象编程)OOP和AOP都是一种编程的思想ios_lowLevelOOP编程思想更加倾向于对业务模块的._ios方法交换 父类 死循环

前端移动端适配 - 媒体查询适配方案_pc端写的固定宽度,媒体查询适配移动端,width写100%没效果-程序员宅基地

文章浏览阅读3.4k次,点赞6次,收藏18次。工作中难免会有写静态页面的需求,有时候移动端适配真的是做的心累,如果自己新做一个页面倒还好,整体布局会按照自己习惯来,但有时候不得不修改别人的代码,尤其是别人没适配好的代码,找样式以及命名规范等问题够折磨一整天了。_pc端写的固定宽度,媒体查询适配移动端,width写100%没效果

随便推点

PyTorch深度学习遥感影像地物分类与目标检测、分割及遥感影像问题深度学习优化技术_深度学习框架遥感-程序员宅基地

文章浏览阅读930次,点赞2次,收藏8次。深度卷积网络采用“端对端”的特征学习,通过多层处理机制揭示隐藏于数据中的非线性特征,能够从大量训练集中自动学习全局特征(这种特征被称为“学习特征”),是其在遥感影像自动目标识别取得成功的重要原因,也标志特征模型从手工特征向学习特征转变。为使广大学者能理解卷积神经网络背后的数学模型和计算机算法,掌握利用PyTorch为基础的遥感影像地物分类,遥感图像目标检测,以及遥感图像目标分割等应用。1.现有几个优秀模型结构的演变原理,包括AlexNet,VGG,googleNet,ResNet,DenseNet等模型。_深度学习框架遥感

Oauth2.0 ( RFC-6749 ) 中文译文_rfc6749中文-程序员宅基地

文章浏览阅读1w次,点赞14次,收藏26次。本文为RFC6749(OAuth2.0)的中文翻译,本文在不影响原文语义的情况下尽可能地采用更符合中文习惯的方式进行表述,如有翻译不妥当的地方请在评论中指出。PS:有些过于简单,或者不适合翻译的内容直接以原文呈现,如11节往后,基本上能读懂英语单词就不会有阅读障碍,所以就不翻译了。PPS:有些明明能简单的表述出来的东西非要绕十八个弯写出来,搞学术的都是鬼才。。。本文由spawpaw@ho..._rfc6749中文

Java音乐播放器,窗体程序 完整源码_java音乐播放器源码-程序员宅基地

文章浏览阅读1.6k次,点赞2次,收藏20次。java 编写音乐播放器,窗体程序, 完整源码游戏,可以直接用来做课程设计或者毕业设计;代码功能完善,下载后可以直接运行!!_java音乐播放器源码

csgo准星设置代码_【老许教你玩CSGO】如何调整合适的准星参数-程序员宅基地

文章浏览阅读6.6k次。想要找到合适的准星,首先我们得知道调整准星的参数有哪些,这是最基本的准星调整方法,参数在网上一搜便可全部知晓。但是大家在面对一大串的准星参数时,完全是看花了眼,不知道以什么样的姿势来输入,这里我就来给大家整合一下,帮助更多的玩家能找到属于自己的那个准星。其实真正所需要的只有几个参数,只需要这几个便可完成准星的设置1 cl_crosshairsize 这个参数是用来设置准星的大小2 cl_cross..._csgo准星代码怎么用 csdn

docker curl: (56) Recv failure: 连接被对方重设_偶然出现curl pecv failure:连接被对方重设-程序员宅基地

文章浏览阅读7k次。现象:docker服务正常起来了,但是在外部curl的时候一直报错误解决方案:服务可能使用了配置文件或者在代码里面写死了。把原有127.0.0.1地址改成0.0.0.0。重新运行即可_偶然出现curl pecv failure:连接被对方重设

怎么去除pdf中的水印-程序员宅基地

文章浏览阅读112次。有时候下载一个pdf文件会发里面有水印,想要使用里面内容的时候很不方便,那么怎么去除pdf中的水印呢?想知道话就来看看下面的文章吧。 其实想要去除pdf中的水印是有两种方法的,首先给大家介绍的是网页版的 要借助的软件:迅捷pdf在线转换器 软件功能介绍:是一个不用下载软件就可以完成各种文..._ab pdf 下面水印

推荐文章

热门文章

相关标签