ABC 207 D. Congruence Points 计算几何

2023-10-20 17:59

本文主要是介绍ABC 207 D. Congruence Points 计算几何,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

题目链接

在这里插入图片描述
题解:

这道题,看起来既要平移又要旋转,挺麻烦。可以抓住一个性质,如果这两个点集可以通过旋转匹配,那么,这两个点集中的点相对重心的位置分布是相同的

因此,我们可以先求出这两个点集中的点相对其重心的相对坐标,这样我们就能够忽略掉平移这个操作所带来的影响。

下面就要考虑能否通过旋转使两个点集匹配。

从点集 S S S中选一个 x x x坐标不为0的点 p p p,这样可以保证有 t a n tan tan值。由于这题范围很小,因此直接暴力枚举点集 T T T中的所有点转到点 p p p所需的角度 θ \theta θ,然后,判断点集 S S S的剩余点转 θ \theta θ能否与点集 T T T中的某一个点匹配即可。

时间复杂度 O ( n 3 ) O(n^3) O(n3)

实现细节见代码:

#include <bits/stdc++.h>
using namespace std;
#define endl '\n'
#define int long long
typedef long long ll;
const int inf = 0x3f3f3f3f;
const int MAXN = 1e6 + 10;
const int mod = 1e9 + 7;
const double eps = 1e-6;
int a[110], b[110], c[110], d[110];
signed main() 
{ios::sync_with_stdio(false);cin.tie(0), cout.tie(0);int n;int sumx = 0, sumy = 0;cin >> n;sumx = sumy = 0;for (int i = 0; i < n; i++) {cin >> a[i] >> b[i];sumx += a[i];sumy += b[i];a[i] *= n;b[i] *= n;}for (int i = 0; i < n; i++) {a[i] -= sumx;b[i] -= sumy;}sumx = sumy = 0;for (int i = 0; i < n; i++) {cin >> c[i] >> d[i];sumx += c[i];sumy += d[i];c[i] *= n;d[i] *= n;}for (int i = 0; i < n; i++) {c[i] -= sumx;d[i] -= sumy;}for (int i = 0; i < n; i++) { // 找到一个x坐标不为0的点if (a[i]) {swap(a[i], a[0]);swap(b[i], b[0]);break;}}for (int i = 0; i < n; i++) {double angle = atan2(d[i], c[i]) - atan2(b[0], a[0]);bool flag = true;for (int j = 0; j < n; j++) {double newa = a[j] * cos(angle) - b[j] * sin(angle);double newb = a[j] * sin(angle) + b[j] * cos(angle);bool flag1 = false;for (int k = 0; k < n; k++) {if (fabs(newa - c[k]) <= eps && fabs(newb - d[k]) <= eps) {flag1 = true;break;}}flag &= flag1;}if (flag) {cout << "Yes" << endl;return 0;}}cout << "No" << endl;return 0;
}

这篇关于ABC 207 D. Congruence Points 计算几何的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/248748

相关文章

如何用Java结合经纬度位置计算目标点的日出日落时间详解

《如何用Java结合经纬度位置计算目标点的日出日落时间详解》这篇文章主详细讲解了如何基于目标点的经纬度计算日出日落时间,提供了在线API和Java库两种计算方法,并通过实际案例展示了其应用,需要的朋友... 目录前言一、应用示例1、天安门升旗时间2、湖南省日出日落信息二、Java日出日落计算1、在线API2

uva 10387 Billiard(简单几何)

题意是一个球从矩形的中点出发,告诉你小球与矩形两条边的碰撞次数与小球回到原点的时间,求小球出发时的角度和小球的速度。 简单的几何问题,小球每与竖边碰撞一次,向右扩展一个相同的矩形;每与横边碰撞一次,向上扩展一个相同的矩形。 可以发现,扩展矩形的路径和在当前矩形中的每一段路径相同,当小球回到出发点时,一条直线的路径刚好经过最后一个扩展矩形的中心点。 最后扩展的路径和横边竖边恰好组成一个直

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

poj 3304 几何

题目大意:给出n条线段两个端点的坐标,问所有线段投影到一条直线上,如果这些所有投影至少相交于一点就输出Yes!,否则输出No!。 解题思路:如果存在这样的直线,过投影相交点(或投影相交区域中的点)作直线的垂线,该垂线(也是直线)必定与每条线段相交,问题转化为问是否存在一条直线和所有线段相交。 若存在一条直线与所有线段相交,此时该直线必定经过这些线段的某两个端点,所以枚举任意两个端点即可。

POJ 2318 几何 POJ 2398

给出0 , 1 , 2 ... n 个盒子, 和m个点, 统计每个盒子里面的点的个数。 const double eps = 1e-10 ;double add(double x , double y){if(fabs(x+y) < eps*(fabs(x) + fabs(y))) return 0 ;return x + y ;}struct Point{double x , y

poj 2653 几何

按顺序给一系列的线段,问最终哪些线段处在顶端(俯视图是完整的)。 const double eps = 1e-10 ;double add(double x , double y){if(fabs(x+y) < eps*(fabs(x) + fabs(y))) return 0 ;return x + y ;}struct Point{double x , y ;Point(){}Po

音视频入门基础:WAV专题(10)——FFmpeg源码中计算WAV音频文件每个packet的pts、dts的实现

一、引言 从文章《音视频入门基础:WAV专题(6)——通过FFprobe显示WAV音频文件每个数据包的信息》中我们可以知道,通过FFprobe命令可以打印WAV音频文件每个packet(也称为数据包或多媒体包)的信息,这些信息包含该packet的pts、dts: 打印出来的“pts”实际是AVPacket结构体中的成员变量pts,是以AVStream->time_base为单位的显