poj1279

板子题,求多边形内核面积。

话说jls的板子返回的是边,然后我就在冥思苦想怎么根据割边求面积啊。。

然后发现自己果然是个傻逼,求一下交点存起来就好了。。。

//板子题到此为止了

  1 #include <iostream>
  2 #include <cstdio>
  3 #include <algorithm>
  4 #include <vector>
  5 #include <cmath>
  6 #include <deque>
  7 
  8 using namespace std;
  9 typedef double db;
 10 const db eps=1e-6;
 11 const db pi=acos(-1);
 12 int sign(db k){
 13     if (k>eps) return 1; else if (k<-eps) return -1; return 0;
 14 }
 15 int cmp(db k1,db k2){return sign(k1-k2);}
 16 struct point{
 17     db x,y;
 18     point operator + (const point &k1) const{return (point){k1.x+x,k1.y+y};}
 19     point operator - (const point &k1) const{return (point){x-k1.x,y-k1.y};}
 20     point operator * (db k1) const{return (point){x*k1,y*k1};}
 21     point operator / (db k1) const{return (point){x/k1,y/k1};}
 22     db getP()const { return sign(y)==1||(sign(y)==0&&sign(x)==-1);}
 23 };
 24 db cross(point k1,point k2){ return k1.x*k2.y-k1.y*k2.x;}
 25 db dot(point k1,point k2){ return k1.x*k2.x+k1.y*k2.y;}
 26 db rad(point k1,point k2){ return atan2(cross(k1,k2),dot(k1,k2));}
 27 int compareangle(point k1,point k2){
 28     return k1.getP()<k2.getP()||(k1.getP()==k2.getP()&&sign(cross(k1,k2))>0);
 29 }
 30 point getLL(point k1,point k2,point k3,point k4){
 31     db w1=cross(k1-k3,k4-k3),w2=cross(k4-k3,k2-k3);
 32     return (k1*w2+k2*w1)/(w1+w2);
 33 }
 34 struct line{
 35     point p[2];
 36     line(point k1,point k2){p[0]=k1;p[1]=k2;}
 37     point &operator[](int k){ return p[k];}
 38     int include(point k){ return sign(cross(p[1]-p[0],k-p[0])>0);}
 39     point dir(){ return p[1]-p[0];}
 40 };
 41 point getLL(line k1,line k2){
 42     return getLL(k1[0],k1[1],k2[0],k2[1]);
 43 }
 44 int parallel(line k1,line k2){ return sign(cross(k1.dir(),k2.dir()))==0;}
 45 int sameDir(line k1,line k2){
 46     return parallel(k1,k2)&&sign(dot(k1.dir(),k2.dir()))==1;
 47 }
 48 int operator <(line k1,line k2){
 49     if(sameDir(k1,k2))return k2.include(k1[0]);
 50     return compareangle(k1.dir(),k2.dir());
 51 }
 52 int checkpos(line k1,line k2,line k3){ return k3.include(getLL(k1,k2));}
 53 vector<line> getHL(vector<line> &L){
 54     sort(L.begin(),L.end());deque<line> q;
 55     for(int i=0;i<L.size();i++){
 56         if(i&&sameDir(L[i],L[i-1]))continue;
 57         while (q.size()>1&&!checkpos(q[q.size()-2],q[q.size()-1],L[i]))q.pop_back();
 58         while (q.size()>1&&!checkpos(q[1],q[0],L[i]))q.pop_front();
 59         q.push_back(L[i]);
 60     }
 61     while (q.size()>2&&!checkpos(q[q.size()-2],q[q.size()-1],q[0]))q.pop_back();
 62     while (q.size()>2&&!checkpos(q[1],q[0],q[q.size()-1]))q.pop_front();
 63     vector<line> ans;for(int i=0;i<q.size();i++)ans.push_back(q[i]);
 64     return ans;
 65 }
 66 int t,n;
 67 point p[1551];
 68 bool cw(){//时针
 69     db s=0;
 70     for(int i=1;i<n-1;i++){
 71         s+=cross(p[i]-p[0],p[i+1]-p[0]);
 72     }
 73     return s>0;
 74 }
 75 vector<line>l;
 76 int main(){
 77     scanf("%d",&t);
 78     while (t--){
 79         scanf("%d",&n);
 80         for(int i=0;i<n;i++){
 81             scanf("%lf%lf",&p[i].x,&p[i].y);
 82         }
 83         if(!cw())reverse(p,p+n);
 84         for(int i=0;i<n;i++){
 85             l.push_back(line(p[i],p[(i+1)%n]));
 86         }
 87         l=getHL(l);
 88         if(l.size()<3){
 89             printf("0.00\n");
 90         } else{
 91             vector<point> a;
 92             for(int i=0;i<l.size();i++){
 93                 a.push_back(getLL(l[i],l[(i+1)%l.size()]));
 94             }
 95             db ans = 0;
 96             for(int i=1;i<a.size()-1;i++)
 97                 ans+=cross(a[i]-a[0],a[i+1]-a[0]);
 98             ans/=2;
 99             printf("%.2f\n",ans);
100         }
101         l.clear();
102     }
103 }
View Code

 

转载于:https://www.cnblogs.com/MXang/p/10453657.html

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

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

相关文章

蓝桥杯第八届省赛JAVA真题----日期问题

标题&#xff1a;日期问题小明正在整理一批历史文献。这些历史文献中出现了很多日期。小明知道这些日期都在1960年1月1日至2059年12月31日。令小明头疼的是&#xff0c;这些日期采用的格式非常不统一&#xff0c;有采用年/月/日的&#xff0c;有采用月/日/年的&#xff0c;还有…

蓝桥杯第八届省赛JAVA真题----拉马车

标题&#xff1a;拉马车小的时候&#xff0c;你玩过纸牌游戏吗&#xff1f;有一种叫做“拉马车”的游戏&#xff0c;规则很简单&#xff0c;却很吸引小朋友。其规则简述如下&#xff1a;假设参加游戏的小朋友是A和B&#xff0c;游戏开始的时候&#xff0c;他们得到的随机的纸牌…

算法提高 合并石子【动态规划】

问题描述在一条直线上有n堆石子&#xff0c;每堆有一定的数量&#xff0c;每次可以将两堆相邻的石子合并&#xff0c;合并后放在两堆的中间位置&#xff0c;合并的费用为两堆石子的总数。求把所有石子合并成一堆的最小花费。输入格式输入第一行包含一个整数n&#xff0c;表示石…

回溯法在解决八皇后问题中的应用

回溯法&#xff1a;有这样一类题目&#xff0c;它们要求在相对问题的输入规模按照指数速度增长&#xff08;或者更快&#xff09;的域中&#xff0c;找出一个具有指定特性的元素。例如&#xff1a;在图顶点的所有排列中求一个哈密顿回路&#xff0c;在背包问题的一个实例中求其…

python数据结构与算法之问题求解

懂得计算机的童鞋应该都知道&#xff0c;一条计算机程序由数据结构跟算法两大部分组成。所以&#xff0c;其实不管你使用哪种计算机语言编写程序&#xff0c;最终这两部分才是一个程序设计的核心。所以&#xff0c;一个不懂得数据结构与算法的程序员不是一个好工程师。因此&…

蓝桥杯算法提高----2n皇后

2n皇后问题描述给定一个n*n的棋盘&#xff0c;棋盘中有一些位置不能放皇后。现在要向棋盘中放入n个黑皇后和n个白皇后&#xff0c;使任意的两个黑皇后都不在同一行、同一列或同一条对角线上&#xff0c;任意的两个白皇后都不在同一行、同一列或同一条对角线上。问总共有多少种放…

运用tp5上传图片,并生成缩略图

最近想做个相册&#xff0c;需要用到上传图像&#xff0c;并且考虑到性能问题&#xff0c;还要生成缩略图&#xff0c;就学习下。在网上看了很多大神写的文章&#xff0c;经过各种调试总算出来了&#xff0c;分享下。不好之处&#xff0c;多多指教 ​ ​ ps&#xff1a;运用tp5…

求解最长回文子串----Manacher 算法

最长回文子串问题&#xff1a;给定一个字符串&#xff0c;求它的最长回文子串长度。 如果一个字符串正着读和反着读是一样的&#xff0c;那么我们称之为回文串。例如&#xff1a;abba、aaaa、abvcba、123321等 暴力法&#xff1a;遍历字符串的所有子串&#xff0c;对每个字串进…

Peter's smokes -poj 2509

题意&#xff1a;彼得有n支雪茄&#xff0c;每k个烟头可以换一支新雪茄&#xff0c;问彼得最多可以吸多少支雪茄 &#xff1f; 当时自己做时&#xff0c;错在了直接在while循环开始前&#xff0c;便将雪茄的初始数量给加上了&#xff0c;然而应该是先处理后再加上最终剩余的雪茄…

模式匹配算法----KMP算法以及next数组的解法

KMP算法&#xff1a;求字符串匹配&#xff08;也叫模式匹配&#xff09;的算法&#xff0c;即给定一个字符串&#xff0c;求其某一子串在其中出现的位置。 普通模式匹配 例如&#xff1a;给定字符串为abcabaaabaabcac&#xff0c;求其子串abaabcac在其中出现的位置。 结果为…

Spring Boot使用layui的字体图标时无法正常显示 解决办法

在html文件使用字体图标并且预览时正常&#xff0c;但是启动工程后显示不正常&#xff0c;浏览器调试界面显示字体文件无法decode&#xff1a; Failed to decode downloaded font: xxxxx 如图所示&#xff1a; 显示结果&#xff1a; 原因&#xff1a;经过maven的filter&#xf…

蓝桥杯第七届国赛JAVA真题----机器人塔

机器人塔X星球的机器人表演拉拉队有两种服装&#xff0c;A和B。 他们这次表演的是搭机器人塔。 类似&#xff1a;队内的组塔规则是&#xff1a; A 只能站在 AA 或 BB 的肩上。B 只能站在 AB 或 BA 的肩上。你的任务是帮助拉拉队计算一下&#xff0c;在给定A与B的人数时&…

致自己人生的一句话

千锤万凿出深山 烈火焚烧若等闲 转载于:https://www.cnblogs.com/wuhefeng/p/10467737.html

单链表常见操作

最近在重新学习数据结构和算法的知识&#xff0c;数据结构和算法问题本身比较枯燥和乏味&#xff0c;而且比较难不容易掌握&#xff0c;但是属于程序员内功的一部分&#xff0c;学习起来容易上瘾。 1. 单链表定义 package algorithm.datastructors;/*** 单向链表* author i3247…

蓝桥杯第五届省赛JAVA真题----最长公共子序列

一个串的子串是指该串的一个连续的局部。如果不要求连续&#xff0c;则可称为它的子序列。 比如对串&#xff1a; "abcdefg" 而言&#xff0c;"ab","abd","bdef" 等都是它的子序列。 特别地&#xff0c;一个串本身&#xff0c;以及空串…

Python divmod() 函数

Python divmod() 函数 Python 内置函数 python divmod() 函数把除数和余数运算结果结合起来&#xff0c;返回一个包含商和余数的元组(a // b, a % b)。 在 python 2.3 版本之前不允许处理复数。 函数语法 divmod(a, b)参数说明&#xff1a; a: 数字b: 数字实例 >>>div…

蓝桥杯第五届省赛JAVA真题----单位分数

形如&#xff1a;1/a 的分数称为单位分数。 可以把1分解为若干个互不相同的单位分数之和。 例如&#xff1a; 1 1/2 1/3 1/9 1/18 1 1/2 1/3 1/10 1/15 1 1/3 1/5 1/7 1/9 1/11 1/15 1/35 1/45 1/231 等等&#xff0c;类似这样的分解无穷无尽。我们增加一个约…

.NET面试题总结

1、c#垃圾回收机制 从以下方面入手展开&#xff1a; 1、压缩合并算法 2、代的机制 3、GC调用终结器 2、委托和事件 先说它的定义&#xff1a;委托的本质是类&#xff0c;类型安全的指针&#xff0c;然后从用途上考虑&#xff0c;事件是包装的委托&#xff0c;但事件不是委托…

蓝桥杯第六届省赛JAVA真题----垒骰子

垒骰子 赌圣atm晚年迷恋上了垒骰子&#xff0c;就是把骰子一个垒在另一个上边&#xff0c;不能歪歪扭扭&#xff0c;要垒成方柱体。 经过长期观察&#xff0c;atm 发现了稳定骰子的奥秘&#xff1a;有些数字的面贴着会互相排斥&#xff01; 我们先来规范一下骰子&#xff1a;1 …

JAVA应试技巧----大数开方

大数开方--牛顿迭代法&#xff1a; 如果一个数的位数为偶数个&#xff0c;那么这个数开方之后就有n/2位&#xff1b; 如果一个数的位数为奇数个&#xff0c;那么这个数开方之后就有n/21位&#xff1b; 我们以num1000为例&#xff0c;位数为4&#xff0c;偶数&#xff0c;那么开…