2019 ACM-ICPC 南京 现场赛 K. Triangle (二分)
2024-10-07 17:52:03
题意
给定一个三角形和一个点 \(p\)，如果该点不在三角形边上直接输出 \(-1\)，否则在三角形上找一点 \(q\)，使得线段 \(pq\) 平分三角形面积。
思路
看完题想都没想直接二分了。
就是比赛时复制粘贴的时候改错了，贡献了三发罚时。
板子来自 kuangbin 的模板。
#include <bits/stdc++.h>
using namespace std;
typedef double db;
const db eps = 1e-8;
const db inf = 1e20;
const db pi = acos(-1.0);
int sgn(db x){
if(fabs(x) < eps)return 0;
if(x < 0)return -1;
else return 1;
}
struct Point{
db x, y;
Point(){}
Point(db _x, db _y){
x = _x;
y = _y;
}
void input(){
scanf("%lf%lf", &x, &y);
}
bool operator == (Point b)const{
return sgn(x-b.x) == 0 && sgn(y-b.y) == 0;
}
Point operator -(const Point &b)const{
return Point(x-b.x, y-b.y);
}
//叉积
db operator ^(const Point &b)const{
return x*b.y - y*b.x;
}
//点积
db operator *(const Point &b)const{
return x*b.x + y*b.y;
}
//返回两点的距离
db dis(Point p){
return hypot(x-p.x, y-p.y);
}
Point operator +(const Point &b)const{
return Point(x+b.x, y+b.y);
}
Point operator *(const db &k)const{
return Point(x*k, y*k);
}
Point operator /(const db &k)const{
return Point(x/k, y/k);
}
};
struct Line{
Point s,e;
Line(){}
Line(Point _s,Point _e){
s = _s;
e = _e;
}
// 点在线段上的判断
bool pointonseg(Point p){
return sgn((p-s)^(e-s)) == 0 && sgn((p-s) * (p-e)) <= 0;
}
// 求两直线的交点
Point crosspoint(Line v){
double a1 = (v.e-v.s)^(s-v.s);
double a2 = (v.e-v.s)^(e-v.s);
return Point((s.x*a2-e.x*a1)/(a2-a1),(s.y*a2-e.y*a1)/(a2-a1));
}
};
// 求点a和点b的中点
Point get_mid(Point a, Point b) {
return (a + b) * 0.5;
}
int main() {
int T;
scanf("%d", &T);
while(T--) {
Point a, b, c, p;
Line ab, bc, ac;
a.input(), b.input(), c.input(), p.input();
ab = Line(a, b); bc = Line(b, c); ac = Line(a, c);
db area = fabs((b - a) ^ (c - a) * 0.5); // 三角形面积
if(ab.pointonseg(p)) { // 点p在线段ab上
if(a.dis(p) < b.dis(p)) { // 点p靠近点a，则另一点一定在线段bc上
Point l = c, r = b;
Point mid = get_mid(l, r);
int times = 1000;
while(times--) {
mid = get_mid(l, r);
db s = fabs((mid - p) ^ (b - p));
if(sgn(s - area) > 0) {
l = mid;
} else if(sgn(s - area) == 0) {
break;
} else {
r = mid;
}
}
printf("%.10lf %.10lf\n", mid.x, mid.y);
} else { // 点p靠近点b，则另一点一定在线段ac上
Point l = c, r = a;
Point mid = get_mid(l, r);
int times = 1000;
while(times--) {
mid = get_mid(l, r);
db s = fabs((mid - p) ^ (a - p));
if(sgn(s - area) > 0) {
l = mid;
} else if(sgn(s - area) == 0) {
break;
} else {
r = mid;
}
}
printf("%.10lf %.10lf\n", mid.x, mid.y);
}
} else if(bc.pointonseg(p)) { // 另外两种情况复制粘贴就行了
if(b.dis(p) < c.dis(p)) {
Point l = a, r = c;
Point mid = get_mid(l, r);
int times = 1000;
while(times--) {
mid = get_mid(l, r);
db s = fabs((mid - p) ^ (c - p));
if(sgn(s - area) > 0) {
l = mid;
} else if(sgn(s - area) == 0) {
break;
} else {
r = mid;
}
}
printf("%.10lf %.10lf\n", mid.x, mid.y);
} else {
Point l = a, r = b;
Point mid = get_mid(l, r);
int times = 1000;
while(times--) {
mid = get_mid(l, r);
db s = fabs((mid - p) ^ (b - p));
if(sgn(s - area) > 0) {
l = mid;
} else if(sgn(s - area) == 0) {
break;
} else {
r = mid;
}
}
printf("%.10lf %.10lf\n", mid.x, mid.y);
}
} else if(ac.pointonseg(p)) {
if(a.dis(p) < c.dis(p)) {
Point l = b, r = c;
Point mid = get_mid(l, r);
int times = 1000;
while(times--) {
mid = get_mid(l, r);
db s = fabs((mid - p) ^ (c - p));
if(sgn(s - area) > 0) {
l = mid;
} else if(sgn(s - area) == 0) {
break;
} else {
r = mid;
}
}
printf("%.10lf %.10lf\n", mid.x, mid.y);
} else {
Point l = b, r = a;
Point mid = get_mid(l, r);
int times = 1000;
while(times--) {
mid = get_mid(l, r);
db s = fabs((mid - p) ^ (a - p));
if(sgn(s - area) > 0) {
l = mid;
} else if(sgn(s - area) == 0) {
break;
} else {
r = mid;
}
}
printf("%.10lf %.10lf\n", mid.x, mid.y);
}
} else {
printf("-1\n");
}
}
return 0;
}
赛后发现直接推也可以，也不用写这么长的代码了 (虽然大部分是复制粘贴)。感觉在演队友。
最新文章
- A Simple OpenGL Shader Example II
- 【URLDecoder】java.lang.IllegalArgumentException: URLDecoder: Illegal hex characters in es
- 解决Unity3D导出apk失败:Failed to re-package resources
- 扩展django的User的部分方法
- js运动框架之一条乱跑的虫子
- 《DOM启蒙》 随笔
- 又发现一个msdn的坑
- JS 获取服务器时间
- **bootstrap常见常用样式总结
- Exception与相关
- Sublime Text 2&3中输入法不跟随光标移动的问题的解决方法
- 图片延迟加载插件jquery.lazyload.js的使用方法
- 超大批量删除redis中无用key+配置
- 如何去除本地文件与svn服务器的关联
- plus、max、Pro、Edge
- python之创建文件写入内容
- 不应该使用String.valueOf的场景
- LOCK TABLES 和 UNLOCK TABLES
- 013 JstlView
- Oracle分割字符串 REGEXP_SUBSTR用法
热门文章
- centos7运行yum报如下提示：Run "yum repolist all" to see the repos you have
- 测开之路三十四：html常用标签
- day 52协程
- Appium初始化设置：手写代码连接手机、appium-desktop连接手机
- HTML创建文本框的3种方式
- 微信小程序实战篇-分类页面制作
- /usr/bin/python^M: bad interpreter: No such file or dir解决办法
- 42.Flatten Binary Tree to Linked List
- c# WPF DataGrid 获取选中单元格信息
- xampp环境下，配置Zend Studio调试php(XDebug) 转摘：http://www.cnblogs.com/tuyithief/archive/2011/06/02/2068431.html