[NOI2005]月下柠檬树 (自适应辛普森)

P4207 [NOI2005]月下柠檬树

在这里插入图片描述

如图,我们要求的面积就是这些圆形跟梯形的组合,由于投射到地面上,显然有h′=htanθh' = \frac{h}{tan \theta}h=tanθh,由此我们就可以开始推导这个f(x)f(x)f(x)函数了。

所以转换为我们要推导出直线a,ba, ba,b的函数表达式了。

有公切线ababab,矩形abgcabgcabgc,所以显然有eg=rc−reeg = r_c - r_eeg=rcre,由三角形acd∼ceg∼befacd \sim ceg \sim befacdcegbef,由此我们可以得到cd,efcd, efcd,ef从而得到ad,bfad, bfad,bf,知道这里对圆公切线问题就已经解决了。

所以接下来我们知道带到函数中利用自适应辛普森求解即可,详细细节见代码描述。

/*Author : lifehappy
*/
#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include <bits/stdc++.h>using namespace std;typedef long long ll;const int inf = 0x3f3f3f3f;
const double eps = 1e-6;const int N = 510;int n;double angle;struct Point {double x, y;
}circle[N], line[N << 1];void get_point() {//求得切线的左右两个坐标。for(int i = 1; i < n; i++) {double x1 = circle[i].x, x2 = circle[i + 1].x, r1 = circle[i].y, r2 = circle[i + 1].y;double ce = r1 * (r1 - r2) / (x2 - x1), xa = x1 + ce, ya = sqrt(r1 * r1 - ce * ce);double gd = r2 * (r1 - r2) / (x2 - x1), xb = x2 + gd, yb = sqrt(r2 * r2 - gd * gd);line[2 * i - 1] = {xa, ya}, line[2 * i] = {xb, yb};}
}double f(double x) {double ans = 0.0;for(int i = 1; i <= n; i++) {//在圆里面,if(x < circle[i].x + circle[i].y && x > circle[i].x - circle[i].y) {ans = max(ans, sqrt(circle[i].y * circle[i].y - (x - circle[i].x) * (x - circle[i].x)));}}for(int i = 1; i <= 2 * (n - 1); i += 2) {//被切线包围,if (x >= line[i].x && x <= line[i + 1].x) {ans = max(ans, (line[i + 1].y - line[i].y) / (line[i + 1].x - line[i].x) * (x - line[i].x) + line[i].y);}}//这些值一定是取最大嘛,因为它可能同时符合圆,切线的要求,所以我们得取覆盖面积最大的。return ans;
}double sim(double l, double r) {return (r - l) * (f(l) + f(r) + f((l + r) / 2.0) * 4.0) / 6.0;
}double asr(double l, double r, double eps, double ans) {double mid = (l + r) / 2.0;double ansl = sim(l, mid), ansr = sim(mid, r);if(fabs(ansl + ansr - ans) < 15.0 * eps) return ansl + ansr + (ansl + ansr - ans) / 15.0;return asr(l, mid, eps / 2.0, ansl) + asr(mid, r, eps / 2.0, ansr);
}int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);scanf("%d %lf", &n, &angle);angle = tan(angle);n++;for(int i = 1; i <= n; i++) {   scanf("%lf", &circle[i].x);circle[i].x /= angle;//改变x的值,求个前缀和得到对应的准确的坐标。circle[i].x += circle[i - 1].x;}for(int i = 1; i < n; i++) {scanf("%lf", &circle[i].y);}circle[n].y = 0;get_point();double l = circle[1].x - circle[1].y, r = circle[n].x;//确定积分区间。for(int i = 1; i <= n; i++) {l = min(l, circle[i].x - circle[i].y);r = max(r, circle[i].x + circle[i].y);}printf("%.2f", 2.0 * asr(l, r, eps, sim(l, r)));return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/314240.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Kong 1.3发布,原生gRPC代理、上游TLS交叉认证

Kong 1.3 发布了&#xff0c;此版本亮点包括支持原生 gRPC 代理、上游 TLS 交叉认证&#xff0c;以及一系列新功能和性能改进。原生 gRPC 代理越来越多的用户转向微服务架构&#xff0c;并且希望有对原生 gRPC 代理的支持&#xff0c;Kong 1.3 解决了这个问题&#xff0c;为支持…

对Windows桌面应用程序进行UI自动化测试

所谓UI自动化测试&#xff0c;就是模拟一个用户&#xff0c;对应用程序的UI进行操作&#xff0c;以完成特定场景的功能性集成测试。要对Windows桌面应用程序进行UI自动化测试&#xff0c;目前可选的技术主要是两种&#xff1a;VS自带的CodedUI Test和AppiumWinAppDriver。但是&…

项目实战中如何使用抽象类和接口

引子&#xff1a;时常会有这么一个疑惑&#xff0c;抽象类和接口功能好像&#xff0c;真正用起来该如何抉择呢&#xff1f;&#xff1f;好问题。。来看看书上怎么说的&#xff08;C#7.0本质论&#xff09;虽然方法可在基类中声明为抽象成员&#xff0c;但是&#xff01;&#x…

番茄日志发布1.0.3版本-增加Kafka支持

番茄日志&#xff08;TomatoLog&#xff09;能做什么可能你是第一次听说TomatoLog&#xff0c;没关系&#xff0c;我可以从头告诉你&#xff0c;通过了解番茄日志&#xff0c;希望能帮助有需要的朋友&#xff0c;番茄日志处理将大大降低你采集、分析、处理日志的过程。介绍Toma…

ArangoDB 3.5发布:流事务API、蒙面数据、搜索性能大幅提升、最短路径功能

ArangoDB 3.5 发布了。ArangoDB 是一个分布式原生的多模型数据库&#xff0c;具有灵活的文档、图形和键值数据模型。使用方便的 SQL 查询语言或 JavaScript 扩展构建高性能应用程序。此版本亮点包括&#xff1a;期待已久的 Streaming Transactions API&#xff0c;可以直接使用…

ASP.NET Core on K8S深入学习(7)Dashboard知多少

本篇已加入《.NET Core on K8S学习实践系列文章索引》&#xff0c;可以点击查看更多容器化技术相关系列文章。在第二篇《部署过程解析与Dashboard》中介绍了如何部署Dashboard&#xff0c;但是没有更多地介绍如何使用Dashboard&#xff0c;本文就来对Dashboard的使用进行补充。…

【学习笔记】Docker - 02. 在容器中运行软件(上)

2.1 控制容器: 构建一个网站监视器 需求: 客户想让你做一个网站, 这个网站需要被紧密的监视, 如果服务器宕机了, 那么它们的团队会收到相关的邮件. 这里用到了3个容器. 第一个运行NGINX; 第二个运行一个叫做mailer的程序. 这两个容器都是detached的. Detached 表示容器将在后台…

lintcode 有效的括号序列

心血来潮&#xff0c;半夜做了一道Lintcode的题目&#xff0c;调试完睡觉&#xff0c;欢迎大家批评指正。 public boolean isValidParentheses(String s) {// Write your code hereStack stacknew Stack();if(s.length()0){return true;}else if(s.length()1){return false;}fo…

2019 年 8 月编程语言排行榜,C#重回增长之路

至于火热的 Python&#xff0c;其占有率还是保持着一如既往的增长势头&#xff0c;由 7 月的 9.260% 上升到现在的 10.020%。C 的占有率则稍有下降&#xff0c;由 6.705% 下降到了 6.057%&#xff0c;C# 和 Visual Basic .NET 分别呈现出增长和下降的趋势。可能是本期榜单没什么…

asp.net core系列 71 Web架构分层指南

一.概述本章Web架构分层指南&#xff0c;参考了“Microsoft应用程序体系结构指南”&#xff08;该书是在2009年出版的&#xff0c;当时出版是为了帮助开发人员和架构师更快速&#xff0c;更低风险地使用Microsoft平台和.NET Framework设计和构建有效&#xff0c;高质量的应用程…

C#高级语法之泛型、泛型约束,类型安全、逆变和协变(思想原理)

一、为什么使用泛型&#xff1f;泛型其实就是一个不确定的类型&#xff0c;可以用在类和方法上&#xff0c;泛型在声明期间没有明确的定义类型&#xff0c;编译完成之后会生成一个占位符&#xff0c;只有在调用者调用时&#xff0c;传入指定的类型&#xff0c;才会用确切的类型…

谈自由,ASP.NET Core才是未来?

首先我要说一下自己对自由的理解&#xff1a;自由是我可以选择不干什么&#xff0c;但我要保留我可以干什么的可能性。比如说我现在只有一个码农的角色&#xff0c;但我仍然要保留可以扮演其他角色的可能&#xff0c;比如成为一个作者&#xff0c;当我写下文章的时候已经是了&a…

今天,全网曝光这几个公众号

有人统计过&#xff0c;我们平均每天花在看内容上的时间是5-6小时与其每天被各种看过就忘的内容占据时间不如看点真正对你有价值的信息下面小编为你推荐几个高价值的公众号&#xff0c;这些公众号都是专注.NET技术它们提供的信息能真正提高你生活的质量当你迷茫的时候刷刷这些大…

.NET Core 小程序开发零基础系列(1)——开发者启用并牵手成功

最近几个月本人与团队一直与小程序打交道&#xff0c;对小程序的实战开发算比较熟悉&#xff0c;也因一些朋友经常问我各种小程序问题&#xff0c;无不能一一回答&#xff0c;想了很久&#xff0c;决定还是空余时间来写写文章吧&#xff0c;偶尔发现一个人安静的时候写文章特爽…

学习笔记之12个月提升计划

Java世界博大精深&#xff0c;有太多的东西要学。如果一头扎进去&#xff0c;很可能会淹没在Java技术的海洋里。于是&#xff0c;最近一直在思考列一个提纲&#xff0c;作为高级工程师到资深、再到架构之路的路标。 学习笔记一栏&#xff0c;即为本计划的博客记录。将自己的计划…

从“梁漱溟:思考问题有八层境界”所联想到的

最近一段时间以来写的文章比较少了&#xff0c;这固然是有一些客观原因&#xff0c;但确实有我不可说的一些自我反省和认识等主观因素。记得8月初有一次友人聚餐&#xff0c;席间有朋友聊到公众号的运营心得体会&#xff0c;其中有一条是&#xff1a;避免粉丝减少的黄金法则之一…

程序员过关斩将--cookie和session的关系其实很简单

喜欢就点关注吧!月高风下&#xff0c;下班路上....菜菜哥&#xff0c;告诉你一个秘密&#xff0c;但是不允许告诉任何人这么秘密&#xff0c;你有男票了&#xff1f;~不是&#xff0c;昨天我偷偷去面试了&#xff0c;结果挂了这不是好事吗&#xff0c;上天让公司留住你.....好吧…

.NET导出Excel的四种方法及评测

前言导出Excel是.NET的常见需求&#xff0c;开源社区、市场上&#xff0c;都提供了不少各式各样的Excel操作相关包。本文&#xff0c;我将使用NPOI、EPPlus、OpenXML、Aspose.Cells四个市面上常见的库&#xff0c;各完成一个导出Excel示例。然后对其代码风格和性能做一个横向比…

[2021.1.27多校省选模拟10]染色(min-max容斥/二项式反演)

[2021.1.27多校省选模拟10]染色 突然发现我对概率期望的理解不是很好。。。 部分分1&#xff1a;可以直接进行状压dp&#xff0c;然后按照题意模拟即可。 部分分2&#xff1a;首先可以发现这个问题是min_max容斥形式&#xff0c;然后对于min(T)的问题&#xff0c;我们将问题转…