[神坛]pta
*相邻两条边围成的三角形面积会是最小的
极角排序+叉积计算三角形面积
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N = 2e5 + 10;
int n, m;
double ans = -1;
struct node{
int x, y;
int pos;
}a[N], b[N];
int jd(node s){
if(s.x > 0 && s.y > 0) return 1;
if(s.x < 0 && s.y > 0) return 2;
if(s.x < 0 && s.y < 0) return 3;
if(s.x > 0 && s.y < 0) return 4;
return 0;//不加的话,vscode会报错
}
//计算几何:极角排序
bool cmp(node s1, node s2){
if(s1.pos != s2.pos) return s1.pos < s2.pos;
return s1.x * s2.y - s1.y * s2.x > 0;
}
signed main(){
// ios::sync_with_stdio(false);
// cin.tie(0); cout.tie(0);
cin >> n;
for(int i = 1; i <= n; i++) cin >> a[i].x >> a[i].y;
for(int i = 1; i <= n; i++){
int cnt = 0;
for(int j = 1; j <= n; j++){//确定第二个点
if(i == j) continue;
b[++cnt].x = a[i].x - a[j].x;
b[cnt].y = a[i].y - a[j].y;
b[cnt].pos = jd(b[cnt]);
}
sort(b + 1, b + n, cmp);
for(int j = 1; j < n - 1; j ++){
if(ans == -1 || abs(b[j].x * b[j + 1].y - b[j].y * b[j + 1].x ) * 0.5 < ans){
//叉乘计算三角形面积
ans = abs(b[j].x * b[j + 1].y - b[j].y * b[j + 1].x ) * 0.5;
}
}
}
printf("%.3lf", ans);
return 0;
}
原文地址:http://www.cnblogs.com/N-lim/p/16906963.html
1. 本站所有资源来源于用户上传和网络,如有侵权请邮件联系站长!
2. 分享目的仅供大家学习和交流,请务用于商业用途!
3. 如果你也有好源码或者教程,可以到用户中心发布,分享有积分奖励和额外收入!
4. 本站提供的源码、模板、插件等等其他资源,都不包含技术服务请大家谅解!
5. 如有链接无法下载、失效或广告,请联系管理员处理!
6. 本站资源售价只是赞助,收取费用仅维持本站的日常运营所需!
7. 如遇到加密压缩包,默认解压密码为"gltf",如遇到无法解压的请联系管理员!
8. 因为资源和程序源码均为可复制品,所以不支持任何理由的退款兑现,请斟酌后支付下载
声明:如果标题没有注明"已测试"或者"测试可用"等字样的资源源码均未经过站长测试.特别注意没有标注的源码不保证任何可用性