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

智能推荐

基于SSM和MySQL的企业人事管理系统的设计与实现_基于ssm的企业人事管理系统的设计与实现参考文献-程序员宅基地

文章浏览阅读484次,点赞2次,收藏2次。管理员进入主界面,软件开始运行,提供用户登录功能,不同的用户登录操作的功 能不同,非管理员用户登录只能查看一些公告信息等,而管理员登录后,可以进行用户 管理丶部门管理、职位管理、员工管理、公告管理等功能。基于其他企业人事管理软件的不足,要求能够制作一个可以方便、快捷的对员工信 息进行添加、修改、删除的操作,为了能够更好的存储职工的信息,可以将职工的信息添 加到 Word 文档,这样,不但便于保存,还可以通过 Word 文档进行打印。员工信息的管理:维护员工的基本信息,用户可以进行员工档案信息的录入及更改,_基于ssm的企业人事管理系统的设计与实现参考文献

【C语言】手撕二叉树

【C语言】手撕二叉树

Postgresql源码(127)投影ExecProject的表达式执行分析

无论是投影还是别的计算,表达式执行的入口和计算逻辑都是统一的,这里已投影为分析表达式执行的流程。

hive启动beeline报错

出现上面的问题执行以下代码。

【Hadoop】-Hive客户端:HiveServer2 & Beeline 与DataGrip & DBeaver[14]

DataGrip是由JetBrains公司推出的数据库管理软件,DataGrip支持几乎所有主流的关系数据库产品,如DB2、Derby、MySQL、Oracle、SQL Server等,也支持几乎所有主流的大数据生态圈SQL软件,并且提供了简单易用的界面,开发者上手几乎不会遇到任何困难。3、连接成功,在里面我们可以看到我们前面章节所创建的表,这样子就可以在里面操作我们的sql语句的。5、连接成功,在里面我们可以看到我们前面章节所创建的表,这样子就可以在里面操作我们的sql语句的。

java lambda无法使用_java – 为什么不允许lambda函数?-程序员宅基地

文章浏览阅读1.2k次。我一直在Vaadin的GUI中工作,有一些来自我的IT主管的课程.这一切都很棒,但是,今天,我遇到过我不能在addListener方法类型中使用lambda表达式.此方法是自定义的,作为使用它的对象.这是实施:public class ResetButtonForTextField extends AbstractExtension {private final List listeners= n..._java: -source 1.5 中不支持 lambda 表达式

随便推点

FRTC8563时钟芯片的主要特性和应用场景

一款实时时钟(RTC)芯片,它采用SOP-8封装,这种封装形式使得芯片具有较小的体积和良好的引脚连接稳定性,便于集成到各种电子设备中。:该芯片采用低功耗技术,使得在待机状态下功耗极低,有助于延长电池寿命,特别适合用于便携式设备或长时间运行的系统。:FRTC8563基于稳定的晶振工作,能够提供准确的时钟和日历信息,包括年、月、日、星期、小时、分钟和秒等。:芯片支持较宽的电压范围,使其能够适应不同设备的电源需求。提供准确的时间戳和日历信息,支持设备的时间同步和事件记录。:在便携式仪器仪表中,由于其低功耗特性,

wetrtc简介

wetrtc简介

单片机基于ST25DV动态标签的无线通信_st25dv能量采集-程序员宅基地

文章浏览阅读651次。利用I2C有线链路,任何NFC智能手机或NFC/RFID HF专业读卡器以及MCU均可以访问存储在这些标签中的数据,并且支持掉电保存。这些标签的集成度和性价比极高,可提供丰富的功能集,适用于各种应用。此外,开发人员可从评估板、软件工具、移动应用和其他资源构成的完整生态系统中受益,从而加快应用开发速度。在物联网产品的开发过程中,物联网设备非接触式向外界提供可变的交换信息非常有必要,比如路由器向手机提供可配置的WiFi信息,巡检点向手机提供动态的传感器数据等等。_st25dv能量采集

自己搭建 Linux 服务器踩坑记录_建立服务器踩过的坑-程序员宅基地

文章浏览阅读149次。前言妈蛋,自己搭建一个Linux服务器居然能遇到这么多坑。特此整理下,方便下次遇到同样的错误时能够回过头来快速定位问题并解决问题Number 1,服务器重启之后,Xshell 连接不上注:在服务器重启之前,我只安装了 jdk ,配置了 /etc/profile 环境变量,我一直以为是这个原因,后面把jdk 配置注释掉也没用正确的方向应该是先查看 ssh 服务有没有启动键入命令systemctl status sshd.service如果你的显示跟红框一样 【dead..._建立服务器踩过的坑

MT4606-VB_MOSFET产品应用与参数解析-程序员宅基地

文章浏览阅读187次。通过控制20Vgs (±V)的门源电压,可以实现开关管的导通和截止,实现对电流的控制和开关状态的转换。MT4606详细参数说明 - 极性 N+P沟道- 额定电压 ±30V- 额定电流 9A (N沟道), -6A (P沟道)- 导通电阻 15mΩ @ 10V (N沟道), 42mΩ @ 10V (P沟道), 19mΩ @ 4.5V (N沟道), 50mΩ @ 4.5V (P沟道)- 门源电压 20Vgs (±V)- 阈值电压 ±1.65Vth (V)- 封装类型 SOP8。_mt4606

达梦启云平台中,部署使用HIVE笔记_达梦sql中hiveing-程序员宅基地

文章浏览阅读637次。启云平台部署hive_达梦sql中hiveing