2019CCPC A . Angle Beats(斜率哈希)

最新推荐文章于 2021-11-19 10:46:13 发布
最新推荐文章于 2021-11-19 10:46:13 发布
阅读量407 收藏
点赞数
文章标签： 哈希算法 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/demaci/article/details/121071762
版权

题目大意：给定平面内的 n 个点，m次询问，每次询问给出一个点，询问这个点与给定的n个点中的某两个点构成直角三角形的方案数。

n,m<=2000，4000ms

设询问的点为 q 点

分为两种情况，q 点为直角顶点和 q 为锐角顶点。

第一种情况使用一个 map 存下所有其它点与 q 点的斜率，然后在map中找乘积为 -1 的斜率的对数，第二种情况先预处理 n 个点中每个点的每种斜率的点的个数，然后询问的时候遍历一遍算出对应斜率并求查询负倒数斜率的个数即可。

然鹅这两种操作使用map会超时：整体复杂度是O(n^2logn+3nmlogn)

#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define pii pair<int, int>
const int maxn = 2010;
pii a[maxn];
map<pii, int> num[maxn];
map<pii, int> g;
int main()
{
    int n, m;
    scanf("%d%d", &n, &m);
    for (int i = 1; i <= n; i++)
    {
        scanf("%d%d", &a[i].fi, &a[i].se);
    }
    for (int i = 1; i <= n; i++) //预处理
    {
        for (int j = 1; j <= n; j++)
        {
            if (i == j || a[i] == a[j]) continue;
            int tx = a[i].fi - a[j].fi;
            int ty = a[i].se - a[j].se;
            if (tx && ty)
            {
                int t = __gcd(abs(tx), abs(ty));
                tx /= t, ty /= t; //约分
                //规范斜率使分母非负
                if ((tx > 0 && ty < 0) || (tx < 0 && ty < 0)) tx = -tx, ty = -ty;
            }
            else
            {
                if (tx) tx = 1;
                if (ty) ty = 1;
            }
            num[i][{ tx, ty }]++;
        }
    }
    while (m--)
    {
        g.clear();
        int x, y;
        int ans = 0;
        scanf("%d%d", &x, &y);
        pii p = { x, y };
        for (int i = 1; i <= n; i++)
        {
            if (a[i] == p) continue;
            int tx = a[i].fi - x;
            int ty = a[i].se - y;
            if (tx && ty)
            {
                int t = __gcd(abs(tx), abs(ty));
                tx /= t, ty /= t;
            }
            else
            {
                if (tx) tx = 1;
                if (ty) ty = 1;
            }
            if ((tx > 0 && ty < 0) || (tx < 0 && ty < 0)) tx = -tx, ty = -ty;
            g[{ tx, ty }]++;
        }
        for (auto gg : g)
        {
            if (gg.fi.fi == 0 && gg.fi.se == 1 && g.count({ 1, 0 })) // 特判垂直和水平
            {
                ans += gg.se * g[{ 1, 0 }];
            }
            else if (gg.fi.fi == 1 && gg.fi.se == 0 && g.count({ 0, 1 }))
            {
                ans += gg.se * g[{ 0, 1 }];
            }
            else
            {
                int tx = -gg.fi.se;
                int ty = gg.fi.fi;
                if ((tx > 0 && ty < 0) || (tx < 0 && ty < 0)) tx = -tx, ty = -ty;
                if (g.count({ tx, ty }))
                    ans += gg.se * g[{ tx, ty }];
            }
        }
        ans /= 2;      //每种情况会被计算两次
        for (int i = 1; i <= n; i++)
        {
            if (p == a[i]) continue;
            int tx = a[i].fi - x;
            int ty = a[i].se - y;
            if (tx && ty)
            {
                int t = __gcd(abs(tx), abs(ty));
                tx /= t, ty /= t;
            }
            else
            {
                if (tx) tx = 1;
                if (ty) ty = 1;
            }
            swap(tx, ty);
            if (tx && ty)
            {
                tx = -tx;
                if ((tx > 0 && ty < 0) || (tx < 0 && ty < 0)) tx = -tx, ty = -ty;
            }
            if (num[i].count({ tx, ty }))
                ans += num[i][{ tx, ty }];
        }
        printf("%d\n", ans);
    }
    return 0;
}

使用unordered_map存储斜率的话就能少一个log的复杂度，就可以过了，但是使用存储坐标哈希的方法会莫名的wa掉，然后改成质数的哈希并更改一下遍历的方式就可以过了： 

// Problem: A. Angle Beats
// URL: https://codeforces.com/gym/102361/problem/A
//
// Powered by CP Editor (https://cpeditor.org)

#include <bits/stdc++.h>
using namespace std;
#define fi first
#define se second
#define pii pair<int, int>
#define ull unsigned long long
const int maxn = 2010;
pii a[maxn];
unordered_map<ull, int> num[maxn];
unordered_map<ull, int> g;
int gcd(int a, int b)
{
    return b ? gcd(b, a % b) : a;
}
inline ull has(int x, int y)
{
    return 233ull * x + 2333ull * y + 998244353;
}
int main()
{
    int n, m;
    scanf("%d%d", &n, &m);
    for (int i = 1; i <= n; i++)
    {
        scanf("%d%d", &a[i].fi, &a[i].se);
    }
    for (int i = 1; i <= n; i++)
    {
        for (int j = 1; j <= n; j++)
        {
            if (i == j || a[i] == a[j]) continue;
            int tx = a[i].fi - a[j].fi;
            int ty = a[i].se - a[j].se;
            if (tx && ty)
            {
                int t = gcd(tx, ty);
                tx /= t, ty /= t;
            }
            else
            {
                if (tx) tx = 1;
                if (ty) ty = 1;
            }
            num[i][has(tx, ty)]++;
        }
    }
    while (m--)
    {
        g.clear();
        int x, y;
        int ans = 0;
        scanf("%d%d", &x, &y);
        pii p = { x, y };
        for (int i = 1; i <= n; i++)
        {
            if (a[i] == p) continue;
            int tx = a[i].fi - x;
            int ty = a[i].se - y;
            if (tx && ty)
            {
                int t = gcd(abs(tx), abs(ty));
                tx /= t, ty /= t;
            }
            else
            {
                if (tx) tx = 1;
                if (ty) ty = 1;
            }
            if ((tx > 0 && ty < 0) || (tx < 0 && ty < 0)) tx = -tx, ty = -ty;
            g[has(tx, ty)]++;
        }
        for (int i = 1; i <= n; i++)
        {
            if (a[i] == p) continue;
            int tx = y - a[i].se;
            int ty = a[i].fi - x;
            if (tx && ty)
            {
                int t = gcd(abs(tx), abs(ty));
                tx /= t, ty /= t;
                ull t1 = has(tx, ty), t2 = has(-tx, -ty);
                if (g.count(t1)) ans += g[t1];
                if (g.count(t2)) ans += g[t2];
            }
            else
            {
                if (tx) tx = 1;
                if (ty) ty = 1;
                ull t = has(tx, ty);
                if (g.count(t)) ans += g[t];
            }
        }
        ans /= 2;
        for (int i = 1; i <= n; i++)
        {
            if (p == a[i]) continue;
            int tx = y - a[i].se;
            int ty = a[i].fi - x;
            if (tx && ty)
            {
                int t = gcd(tx, ty);
                tx /= t, ty /= t;
                ull t1 = has(tx, ty), t2 = has(-tx, -ty);
                if (num[i].count(t1)) ans += num[i][t1];
                if (num[i].count(t2)) ans += num[i][t2];
            }
            else
            {
                if (tx)
                    tx = 1;
                else if (ty)
                    ty = 1;
                ull t = has(tx, ty);
                if (num[i].count(t)) ans += num[i][t];
            }
        }
        printf("%d\n", ans);
    }
    return 0;
}

另外官解是极角排序，还有优化排序方式来省掉约分的gcd复杂度的map的版本可以过，但是擦边3.5s过。

义往尘沙
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
2021CCPC华为云挑战赛.pdf
08-21
2021CCPC华为云挑战赛.pdf
2019第五届CCPC河南省赛
10-28
2019第五届CCPC省赛 word 文档，找了很久一直没找到，这是HUEL—ACM实验室20级杨宇航同学自己扒的，感谢学弟
2019CCPC秦皇岛站 A.Angle Beats(hdu 6731)(极角排序)
krypton12138的博客
10-01 714
链接 Angle Beats 题意 平面上给2000个固有点，点坐标范围-1e9~+1e9，给2000个询问，每个询问给个点A，问和给定的2000个固有点能组成多少个直角三角形？ 思路 对于给定点A，如果点A当作直角顶点，直接极角排序，对每个固有点判断旋转90°的位置有多少点，累加即可。 如果A不是直角顶点，应以固有的20...
Gym - 102361A Angle Beats(几何)
Falcon的博客
01-29 846
题目链接：点击查看 题目大意：给出n个点，再给出m个询问，每次询问给出一个点 x，我们需要回答包括点 x 的直角三角形有多少个 题目分析：题目比较直接，数据也比较小，支持n*n的算法处理，首先我们必须知道，如果要包括点 x 所组成三角形，那么点 x 可以在直角点上，也可以在非直角点上，比较容易想到的方法是离线处理，我们需要分类讨论，也就是上述两个情况，对于每个询问而言，以点 x 为相对原点，对...
2019CCPCQHD Angle Beats
Sofan He
09-23 227
原题链接 题目来源：CCPC2019 秦皇岛站 A题 问题重述 给定二维平面上nnn个点，然后给出q次询问，每次询问给出一个点P（n+q个点都不同），询问P和之前的n个点能构成多少个直角三角形。 其中0≤n,q≤2000;∣x∣,∣y∣≤1e90 \leq n,q \leq 2000; |x|,|y|\leq 1e90≤n,q≤2000;∣x∣,∣y∣≤1e9 问题分析 假设存在两个点A,B跟这次的这个P构成了直角三角形, 那么可以得到以下的式子: {Py−AyPx−Ax⋅Py−ByPx−Bx=−1,P是顶
2019 CCPC秦皇岛站 A - Angle Beats （计算几何）
qq_41645482的博客
10-02 742
题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6731 题目大意：给n个点，查询q个点，问包含这个点的直角三角形有多少个。 题目思路： 对于一个查询点，要么他是直角三角形的直角顶点，要么是锐角顶点，先把询问离线，然后假设该点做直角顶点，遍历n个点，处理出一个向量，之后查询就是查询与当前向量之积为-1的向量有多少个，对于n个点分别做...
2019 CCPC 秦皇岛 A. Angle Beats
qq_43624815的博客
10-27 160
题面链接 题目大意: 在二维平面上 给你n个点 有q次查询 每次查询会有一个坐标Q,询问 有多少点对 和 Q点 能组成直角三角形; 思路:分两个情况: 一:Q点是直角顶点 二:Q点不是直角顶点; 第一个情况:只需要看有多少个和Q点组成的斜率互乘为-1, 第二个情况:就是看Q点和任意一个点 组成的直线的斜率 和 这个点 和 点集合中的点组成的直线有没有存在互乘为-1的 注意点: 会有卡时间需要注意不要用map 建议用哈希查找表 或者 hash_map 或者 unordered_map维护 时间.
2020CCPC网络赛题解.pdf
07-12
2020CCPC网络赛题解
2019CCPC哈尔滨站题解（F、I、J、K）（csdn）————程序.pdf
12-01
2019CCPC哈尔滨站题解（F、I、J、K）（csdn）————程序
【2019CCPC秦皇岛：A】Angle Beats（离线+斜率Hash+分类讨论）
Cassie_zkq的博客
10-05 613
题目地址：https://codeforces.com/gym/102361/problem/A 题目： n个给定点，q个询问点，每次询问给出一个坐标A,问从n中选定两个点B,C,有多少种方案使得ABC是个直角三角形。 解题思路： 标称的思路是从极角考虑的，训练的时候我以为是从极角考虑这个问题的，但是atan2，再加上极角+-，误差太大，我第二个样例就挂了（技艺不精QAQ）。 离线...
2019秦皇岛ccpc A题：Angle Beats[计算几何：统计符合直角三角形的个数]+[向量hash+3hash]
weixin_45630722的博客
10-09 312
题目链接 题目大意：就是已经知道一些点集合{p}，然后q次询问每次给出一个点A，问你这个点和点集合里面多少对点可以构成三角形？ 解题思路：很无奈这道题卡常了我们常规思路就是：我们分类讨论A这个点，因为A这个有两种情况就是A是直角点，A不是直角点 1.A是直角点：那么我们就用A这个点然后和点集合里面的每个点去连一条边然后对于每一条边我们存一下斜率然后用map去hash存一下 2.由于斜率有精度误差所有对于一个边我们在直接存x和y值，但是求得是x/y相等的那么统一除以gcd就可以了 3.但是普通的map+g
A Angle Beats(三角形计数)
liexss的博客
11-15 389
传送门:QAQ 题意：给你n个点，然后有q次询问一个点，询问能和这n个点组成多少个直角三角形。 2≤n≤2000,1≤q≤2000 思路：首先回想到分类讨论，询问点作为直角点和不为直角点，然后就会想到存斜率，利用直角斜率的性质去计数。 果然第一发map+pair tle了，常数太大，外加自己处理太粗糙。 上网查题解时，发现存在两种做法，一种就是用unorder_map或者自己has...
Angle Beats（计算几何 极角排序）
jkchen's Haven
09-29 1211
original link - http://acm.hdu.edu.cn/showproblem.php?pid=6731 题意： 给出nnn个点BBB，mmm次询问，每次给出一个点AAA，问有多少个AB1B2AB_1B_2AB1​B2​为直角三角形。 AB1B2,AB2B1AB_1B_2,AB_2B_1AB1​B2​,AB2​B1​视为一种，保证n+mn+mn+m个点不重复。 解析： 分别考虑...
[2019CCPC秦皇岛] A Angle Beats 极角排序
qq_40859782的博客
10-08 312
给出n≤2e3n\leq2e3n≤2e3个点，以及q≤2e3q\leq2e3q≤2e3个点，保证所有的点都互不相同，对于每次询问的这个点，求问有多少对(u,v)(u,v)(u,v)，和当前这个点构成一个直角三角形。 分成两种情形： 1.询问的点是直角的顶点：对n个点以询问点为中心极角排序。然后枚举每个点，二分出逆时针90°的点有多少个，只考虑一个方向就不会重复计算。 2.询问的点是不是直角的顶点：...
2019秦皇岛A - Angle Beats Gym - 102361A（极角排序，多少个直角三角形）
tomjobs的博客
10-21 821
Given n points P1,P2,⋯,Pn on 2D plane and q queries. In i-th query, a point Ai is given, and you should determine the number of tuples (u,v) that 1≤u<v≤n and Ai,Pu,Pv form a non-degenerate right-an...
「2019 CCPC秦皇岛现场赛A」Angle Beats【哈希 或 二分】
wzw
10-11 341
题目连接 https://codeforces.com/gym/102361/problem/A 题意 就是给你n(n≤2000)n(n\leq 2000)n(n≤2000)个整数点，q(q≤2000)q(q\leq 2000)q(q≤2000)个询问，每次询问给你一个点p(x,y)p(x,y)p(x,y)，求nnn个点中有多少点对(u,v)(u,v)(u,v)，使得u,v,pu,v,pu,...
秦皇岛站2019CCPC A.Angel Beats
weixin_30892037的博客
09-28 694
题意：平面内给定n个点，q次询问，给次给定一个点P，问这个点与平面内n个点可以组成多少直角三角形，其中（n+q）个点互不相等 思路： 分别考虑P点作直角顶点和非直角顶点。这个题思路很简单，就是看如何实现简单而且不会tle！！！ 对于直角顶点和非直角顶点代码都比较简单，求后者有点离线的思想。 这里想说的就是map的用法，自定义小于运算符，使得在map中查找的时候，统一斜率的向量都会...
2019CCPC秦皇岛 A. Angle Beats
RunningBeef的博客
11-19 288
#include<bits/stdc++.h> #define IOS std::ios_base::sync_with_stdio(false), std::cin.tie(0), std::cout.tie(0);// 快读 using namespace std; typedef long long ll; const ll inf = 4e9 + 1; const int N = 2E3 + 10; struct ty{ ll x,y;// x = (xa - x2) y =
2017-CCPC-FINAL
最新发布
10-04
引用:2017 CCPC Final Contest Onsite Round, Sunday, December 3rd, 2017 Problem A. Dogs and Cages Jerry likes dogs. He has N dogs numbered 0, 1, . . . , N − 1. He also has N cages numbered 0, 1, . . . ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值