Problem G Ray in the tube(思维)

最新推荐文章于 2019-08-24 01:03:15 发布

Link_Ray

最新推荐文章于 2019-08-24 01:03:15 发布

阅读量406

点赞数

CC 4.0 BY-SA版权

分类专栏：思维

本文链接：https://blog.youkuaiyun.com/Link_Ray/article/details/89404052

思维专栏收录该内容

9 篇文章

订阅专栏

博客探讨了二维坐标系中，从传感器A到B发射的激光遇到直线反射的情况，求解最多能经过多少个传感器。通过分析射线与传感器的关系，提出了一种优化的解决方案，即在每个传感器处枚举特定数量的射线，达到O((n+m)log(1e9)log(n))的时间复杂度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

https://codeforces.com/gym/101911/problem/L

题意

在二维坐标系中给出两个直线 y = a, y = b。在这两条直线上分别有n个和m个传感器，可以任意选两个传感器A,B，发射一束A->B的激光，这个激光遇到直线就会反射，问最多会经过多少个传感器。

1<=n,m<=3e5
在这里插入图片描述

题解

先解决一个问题，已知一条射线，如何快速求出经过的传感器个数？

对于上下直线分开讨论。

假设如图所示，从下往上的射线在下直线的跨度是 $c = (x_b-x_a)*2$ ，对于每个传感器令其坐标为 $x_u$ ，只需满足 $x_u \equiv x_a \ mod\ c$ 即在这个射线上。
对于上直线，同理，只需满足 $x_u+c/2 \equiv x_a \ mod\ c$

这两条直线间隔的距离是无所谓的，这个只决定射线的斜率而已。
一个很直观的做法就是枚举上下所有的两个点，然后取最大值，但这个复杂度不允许 $O(n^2m^2)$ 。

画画图可以发现，一些c比较大的射线能用比它小的射线代替，例如c = 6可以被 c = 2代替，但是c = 4 不能被 c = 2 代替。

多尝试一下发现只要 $c = 2^l$ 的射线都不可被替代， $c = 2^l*奇数$ 都可以被 $c = 2^l$ 的所代替。

所以我们只需在每个传感器处枚举 $l o g (1 e 9)$ 根射线即可。时间复杂度为 $O ((n + m) l o g (1 e 9) l o g (n))$

代码

#include <bits/stdc++.h>
using namespace std;
#define FOR0(a,b) for(int i = a; i < b; ++i)
#define FORE(a,b) for(int i = a; i <= b; ++i)
typedef long long ll;
typedef pair<int,int> pii;

const int maxn = 1e5+5;

map<ll,ll> mp1,mp2;
ll n,m;
ll a[maxn], b[maxn];;

int main() {
    int y;
    scanf("%lld%lld", &n,&y);
    for(int i = 0; i < n; ++i)
        scanf("%lld", &a[i]);
    scanf("%lld%lld", &m, &y);
    for(int i = 0; i < m; ++i)
        scanf("%lld", &b[i]);

    ll c = 2;
    ll ans = 2;
    for(; c <= 1e9; c <<= 1) {
        mp1.clear();mp2.clear();
        for(int i = 0; i < n; ++i)
            mp1[a[i]%c]++;
        for(int i = 0; i < m; ++i)
            mp2[b[i]%c]++;

        for(int i = 0; i < n; ++i) 
            ans = max(ans, mp1[a[i]%c]+mp2[(a[i]+c/2)%c]);
        for(int i = 0; i < m; ++i)
            ans = max(ans, mp2[b[i]%c]+mp1[(b[i]+c/2)%c]);
    }
    printf("%lld\n", ans);
    return 0;
}