cf13B Letter A(分类+简单计算几何,,)
题意:
给三个线段(每个线段的两个端点的坐标),问这三个线段能否组成字母A。
组成字母A的条件:
1.两个线段有公共端点。
2.这两个线段夹角小于等于90度。
3.第三个线段的两个端点分别在这两个线段上,且各自分割的大小比率不超过4:1
思路:
直接。。。。。
应该开始积累计算几何的代码了,,,,
代码:
struct Point{
double x,y;
};
struct segment{
Point a,b;
}
S[5];
bool samePoint(Point a,Point b){
if(fabs(a.x-b.x)<eps && fabs(a.y-b.y)<eps){
return true;
}
return false;
}
bool PointOnSegment(Point a,segment L){
if((a.x-L.a.x)*(L.b.y-L.a.y)==(a.y-L.a.y)*(L.b.x-L.a.x)
&& min(L.a.x,L.b.x)<=a.x && a.x<=max(L.a.x,L.b.x)
&& min(L.a.y,L.b.y)<=a.y && a.y<=max(L.a.y,L.b.y)){
return true;
}
return false;
}
double rates(Point a,segment L){
double l1=(a.x-L.a.x)*(a.x-L.a.x)+(a.y-L.a.y)*(a.y-L.a.y);
double l2=(a.x-L.b.x)*(a.x-L.b.x)+(a.y-L.b.y)*(a.y-L.b.y);
if(l1<eps || l2<eps){
return 999999.0;
}
if(l1<l2){
return l2/l1;
}
else{
return l1/l2;
}
}
ll Cos(double x1,double y1,double x2,double y2){
if(x1*x2+y1*y2>=0){
return true;
}
return false;
}
bool solve(){
//1 and 2
bool t1=samePoint(S[1].a,S[2].a);
bool t2=samePoint(S[1].a,S[2].b);
bool t3=samePoint(S[1].b,S[2].a);
bool t4=samePoint(S[1].b,S[2].b);
if(t1||t2||t3||t4){
bool flag=false;
if(t1){
double x1,y1,x2,y2;
x1=S[1].b.x-S[1].a.x;
y1=S[1].b.y-S[1].a.y;
x2=S[2].b.x-S[2].a.x;
y2=S[2].b.y-S[2].a.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(t2){
double x1,y1,x2,y2;
x1=S[1].b.x-S[1].a.x;
y1=S[1].b.y-S[1].a.y;
x2=S[2].a.x-S[2].b.x;
y2=S[2].a.y-S[2].b.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(t3){
double x1,y1,x2,y2;
x1=S[1].a.x-S[1].b.x;
y1=S[1].a.y-S[1].b.y;
x2=S[2].b.x-S[2].a.x;
y2=S[2].b.y-S[2].a.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(t4){
double x1,y1,x2,y2;
x1=S[1].a.x-S[1].b.x;
y1=S[1].a.y-S[1].b.y;
x2=S[2].a.x-S[2].b.x;
y2=S[2].a.y-S[2].b.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(flag){
if(PointOnSegment(S[3].a,S[1]) && PointOnSegment(S[3].b,S[2])){
if(rates(S[3].a,S[1])<=16.0 && rates(S[3].b,S[2])<=16.0){
return true;
}
}
if(PointOnSegment(S[3].b,S[1]) && PointOnSegment(S[3].a,S[2])){
if(rates(S[3].b,S[1])<=16.0 && rates(S[3].a,S[2])<=16.0){
return true;
}
}
}
}
//1 and 3
t1=samePoint(S[1].a,S[3].a);
t2=samePoint(S[1].a,S[3].b);
t3=samePoint(S[1].b,S[3].a);
t4=samePoint(S[1].b,S[3].b);
if(t1||t2||t3||t4){
bool flag=false;
if(t1){
double x1,y1,x2,y2;
x1=S[1].b.x-S[1].a.x;
y1=S[1].b.y-S[1].a.y;
x2=S[3].b.x-S[3].a.x;
y2=S[3].b.y-S[3].a.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(t2){
double x1,y1,x2,y2;
x1=S[1].b.x-S[1].a.x;
y1=S[1].b.y-S[1].a.y;
x2=S[3].a.x-S[3].b.x;
y2=S[3].a.y-S[3].b.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(t3){
double x1,y1,x2,y2;
x1=S[1].a.x-S[1].b.x;
y1=S[1].a.y-S[1].b.y;
x2=S[3].b.x-S[3].a.x;
y2=S[3].b.y-S[3].a.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(t4){
double x1,y1,x2,y2;
x1=S[1].a.x-S[1].b.x;
y1=S[1].a.y-S[1].b.y;
x2=S[3].a.x-S[3].b.x;
y2=S[3].a.y-S[3].b.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(flag){
if(PointOnSegment(S[2].a,S[1]) && PointOnSegment(S[2].b,S[3])){
if(rates(S[2].a,S[1])<=16.0 && rates(S[2].b,S[3])<=16.0){
return true;
}
}
if(PointOnSegment(S[2].b,S[1]) && PointOnSegment(S[2].a,S[3])){
if(rates(S[2].b,S[1])<=16.0 && rates(S[2].a,S[3])<=16.0){
return true;
}
}
}
}
//2 and 3
t1=samePoint(S[2].a,S[3].a);
t2=samePoint(S[2].a,S[3].b);
t3=samePoint(S[2].b,S[3].a);
t4=samePoint(S[2].b,S[3].b);
if(t1||t2||t3||t4){
bool flag=false;
if(t1){
double x1,y1,x2,y2;
x1=S[2].b.x-S[2].a.x;
y1=S[2].b.y-S[2].a.y;
x2=S[3].b.x-S[3].a.x;
y2=S[3].b.y-S[3].a.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(t2){
double x1,y1,x2,y2;
x1=S[2].b.x-S[2].a.x;
y1=S[2].b.y-S[2].a.y;
x2=S[3].a.x-S[3].b.x;
y2=S[3].a.y-S[3].b.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(t3){
double x1,y1,x2,y2;
x1=S[2].a.x-S[2].b.x;
y1=S[2].a.y-S[2].b.y;
x2=S[3].b.x-S[3].a.x;
y2=S[3].b.y-S[3].a.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(t4){
double x1,y1,x2,y2;
x1=S[2].a.x-S[2].b.x;
y1=S[2].a.y-S[2].b.y;
x2=S[3].a.x-S[3].b.x;
y2=S[3].a.y-S[3].b.y;
if(Cos(x1,y1,x2,y2)){
flag=true;
}
}
if(flag){
if(PointOnSegment(S[1].a,S[2]) && PointOnSegment(S[1].b,S[3])){
if(rates(S[1].a,S[2])<=16.0 && rates(S[1].b,S[3])<=16.0){
return true;
}
}
if(PointOnSegment(S[1].b,S[2]) && PointOnSegment(S[1].a,S[3])){
if(rates(S[1].b,S[2])<=16.0 && rates(S[1].a,S[3])<=16.0){
return true;
}
}
}
}
return false;
}
int T;
int main(){
cin>>T;
while(T--){
scanf("%lf%lf%lf%lf",&S[1].a.x,&S[1].a.y,&S[1].b.x,&S[1].b.y);
scanf("%lf%lf%lf%lf",&S[2].a.x,&S[2].a.y,&S[2].b.x,&S[2].b.y);
scanf("%lf%lf%lf%lf",&S[3].a.x,&S[3].a.y,&S[3].b.x,&S[3].b.y);
if(solve()){
puts("YES");
}else{
puts("NO");
}
}
return 0;
}
cf13B Letter A(分类+简单计算几何,,)的更多相关文章
- CF13B Letter A
CF13B Letter A 洛谷传送门 题目描述 Little Petya learns how to write. The teacher gave pupils the task to writ ...
- ●POJ 1556 The Doors(简单计算几何+最短路)
●赘述题目 10*10的房间内,有竖着的一些墙(不超过18个).问从点(0,5)到(10,5)的最短路. 按照输入样例,输入的连续5个数,x,y1,y2,y3,y4,表示(x,0--y1),(x,y2 ...
- 2018.07.04 POJ 2398 Toy Storage(二分+简单计算几何)
Toy Storage Time Limit: 1000MS Memory Limit: 65536K Description Mom and dad have a problem: their ch ...
- 2018.07.03 BZOJ 1007: [HNOI2008]水平可见直线(简单计算几何)
1007: [HNOI2008]水平可见直线 Time Limit: 1 Sec Memory Limit: 162 MB Description 在xoy直角坐标平面上有n条直线L1,L2,-Ln, ...
- python - 实现文本分类[简单使用第三方库完成]
第三方库 pandas sklearn 数据集 来自于达观杯 训练:train.txt 测试:test.txt 概述 TF-IDF 模型提取特征值建立逻辑回归模型 代码 # _*_ coding:ut ...
- HDU2948Geometry Darts(简单计算几何)
题目大意就是说两个人掷飞镖,飞镖在所给定的图形内就记一分,现在给定N个图形(圆.三角形和矩形),问每一次比赛(没人分别掷三次)谁赢. #include <map> #include < ...
- Bzoj4558:分类讨论 计算几何 组合数学
国际惯例的题面: 这题让我爆肝啦......这种计数显然容斥,正好不含任何坏点的我们不会算,但是我们能算至少含零个坏点的,至少含一个坏点的,至少含两个坏点的......所以最终的答案就是(至少含零个坏 ...
- 2018.07.04 POJ 1654 Area(简单计算几何)
Area Time Limit: 1000MS Memory Limit: 10000K Description You are going to compute the area of a spec ...
- 2018.07.04 POJ 3304 Segments(简单计算几何)
Segments Time Limit: 1000MS Memory Limit: 65536K Description Given n segments in the two dimensional ...
随机推荐
- 一文彻底搞懂Hive的数据存储与压缩
目录 行存储与列存储 行存储的特点 列存储的特点 常见的数据格式 TextFile SequenceFile RCfile ORCfile 格式 数据访问 Parquet 测试 准备测试数据 存储空间 ...
- Java基础系列(30)- 命令行传参
命令行传参 有时候你希望运行一个程序的时候再传递给它消息.这就要靠传递命令行参数main()函数实现 package method; public class CommandLine { public ...
- 腾讯云与 Grafana Labs 达成深度合作, 推出全新 Grafana 托管服务
9 月 23 日,腾讯云宣布与业界领先的开源数据可视化公司 Grafana Labs 达成深度合作协议,共同开发和验证全新的 Grafana 托管服务,通过 Grafana Labs 开源软件与腾讯云 ...
- vue 改变路由参数
更改路由传递的参数: const query = JSON.parse(JSON.stringify(this.$route.query)) // 获取路由参数信息 query.id = Name / ...
- 分析 ajax 请求并抓取 “今日头条的街拍图”
今日头条抓取页面: 分析街拍页面的 ajax 请求: 通过在 XHR 中查看内容,获取 url 链接,params 参数信息,将两者进行拼接后取得完整 url 地址.data 中的 article_u ...
- P3175-[HAOI2015]按位或【min-max容斥,FWT】
正题 题目链接:https://www.luogu.com.cn/problem/P3175 题目大意 开始有一个\(n\)位二进制数\(s=0\),每次有\(p_i\)概率选取数字\(i\)让\(s ...
- 解决报错:The import javax.servlet.annotation cannot be resolved
maven项目,引入javax.servlet.annotation.WebServlet的jar包,使用@WebServlet注解来实现对传统web.xml中servlet和url的映射 报错:Th ...
- Springboot中使用Redisson实现分布式锁
1. 概述 老话说的好:便宜没好货,有价值的商品,即使再贵,也有人会买. 言归正传,今天继续讨论有关"锁"的话题,synchronized 和 ReentrantLock 大家应该 ...
- 基于Hyperledger Fabric实现ERC721
介绍 超级账本(Hyperledger)项目是首个面向企业应用场景的开源分布式账本平台.由linux基金会牵头,包括 IBM 等 30家初始企业成员共同成立的. 区块链网络主要有三种类型:公共区块链. ...
- 一次OutOfMemoryError: GC overhead limit exceeded
现象: 由于需要将mysql表中的过期数据在凌晨定时读取出过滤后转入到MongoDB,一个转换SQL达到百行,而且有几十个,集中运行后程序反馈异常: Handler dispatch failed; ...