ZR#317.【18 提高 2】A(计算几何 二分)

题意

Sol

非常好的一道题,幸亏这场比赛我没打,不然我估计要死在这个题上qwq

到不是说有多难,关键是细节太多了,我和wcz口胡了一下我的思路,然后他写了一晚上没调出来qwq

解法挺套路的,先提出一个$x$

然后维护一堆直线对应的上凸壳

在凸壳上二分即可。

由于这题的$x$很小,直接处理出答案就行了

/*

*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<map>
#include<vector>
#include<set>
#include<queue>
#include<cmath>
//#include<ext/pb_ds/assoc_container.hpp>
//#include<ext/pb_ds/hash_policy.hpp>
#define Pair pair<int, int>
#define MP(x, y) make_pair(x, y)
#define fi first
#define se second
#define int long long 
#define LL long long 
#define ull unsigned long long 
#define rg register 
#define pt(x) printf("%d ", x);
//#define getchar() (p1 == p2 && (p2 = (p1 = buf) + fread(buf, 1, 1<<22, stdin), p1 == p2) ? EOF : *p1++)
//char buf[(1 << 22)], *p1 = buf, *p2 = buf;
//char obuf[1<<24], *O = obuf;
//void print(int x) {if(x > 9) print(x / 10); *O++ = x % 10 + '0';}
//#define OS  *O++ = ' ';
using namespace std;
//using namespace __gnu_pbds;
const int MAXN = 1e6 + 10, INF = 1e9 + 10, mod = 1e9 + 7, BB = 32323;
const double eps = 1e-9;
inline int read() {
    char c = getchar(); int x = 0, f = 1;
    while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * f;
}
int N;
struct Node {
    double a, b;
    bool operator < (const Node &rhs) const {
        return a == rhs.a ? b < rhs.b : a < rhs.a;
    }
}P[MAXN], s1[MAXN], s2[MAXN];
int t1 = 0, t2 = 0, ans[MAXN];
double cross(Node x, Node y) {
//    printf("%lf\n", 1.0 * (y.b - x.b) / (x.a - x.b));
    return 1.0 * (y.b - x.b) / (x.a - y.a);
}
void Get() {
    sort(P + 1, P + N + 1);
    s1[++t1] = P[1];
    for(int i = 2; i <= N; i++) {
        if(t1 && P[i].a == s1[t1].a) t1--;
        while(t1 > 1 && cross(P[i], s1[t1]) <= cross(s1[t1], s1[t1 - 1])) t1--;
        s1[++t1] = P[i];
    }
    for(int i = 1; i <= N; i++) P[i].b = -P[i].b;
    sort(P + 1, P + N + 1);
    s2[++t2] = P[1];
    for(int i = 2; i <= N; i++) {
        if(t1 && P[i].a == s2[t2].a) t2--;
        while(t2 > 1 && cross(P[i], s2[t2]) <= cross(s2[t2], s2[t2 - 1])) t2--;
        s2[++t2] = P[i];
    }
}
int Query(Node p, int x) {
    return p.a * x * x + p.b * x;
}
void MakeAns() {
    for(int i = 1, c = 1; i <= BB; i++) {
        //printf("%lf\n", cross(s1[c], s1[c + 1]));
        while(c < t1 && cross(s1[c], s1[c + 1]) <= (double)i) c++;
        ans[i + BB] = Query(s1[c], i);
    }
    for(int i = 1, c = 1; i <= BB; i++) {
        while(c < t2 && cross(s2[c], s2[c + 1]) <= (double)i) c++;
        ans[BB - i] = Query(s2[c], i);
    }
}
main() {
//    freopen("a.in", "r", stdin);
    N = read(); int Q = read();
    for(int i = 1; i <= N; i++) P[i].a = read(), P[i].b = read();
    Get();
    MakeAns();
    while(Q--) {
        int x = read();
        printf("%lld\n", ans[x + BB]);
    }
    return 0;
}
/*
2 2 1
1 1
2 1 1
*/

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏ThoughtWorks

TW洞见〡3D打印的各种问题及解决方案

文章作者来自ThoughtWorks:贺思聪 ,图片来自网络。 3D打印机已经买回来几个月了,基本上每天都要打印一些东西,期间遇到了很多的问题积累了很多的经验...

408120
来自专栏HansBug's Lab

Tyvj P1813 [JSOI2008]海战训练

P1813 [JSOI2008]海战训练 时间: 1000ms / 空间: 131072KiB / Java类名: Main 描述 为了准备高层峰会,元首命令武...

376120
来自专栏窗户

scratch写的图灵机

15930
来自专栏数据小魔方

用R语言复盘美国总统大选结果~

这两天各种社交媒体都被美国大选的消息刷屏,各种段子满天飞,把平时不怎么关注政治的小编都吸引了。 美国大选的投票数据,给小编的写作提供了非常宝贵的案例数据,毕竟四...

82250
来自专栏逍遥剑客的游戏开发

Bullet中创建中心不在原点的碰撞体

20850
来自专栏编程微刊

【前端图表】echarts散点图鼠标划过散点显示信息

21630
来自专栏数据结构与算法

逆元的三种解法(附详细证明)

友情提示: Latex加载稍慢,请耐心等待 什么是逆元? 若x满足 我们称x是a在 意义下的逆元 逆元的基本解法 https://loj.ac/pr...

1.1K80
来自专栏Crossin的编程教室

【每周一坑】蜥蜴流感与贝叶斯定理

春季是流感的高发季节。不要觉得只是小小的“感冒”,严重起来甚至也会危及生命,而且还没有特效药。因此,身体不适请及时到医院检查。

17530
来自专栏深度学习之tensorflow实战篇

python jieba分词(结巴分词)、提取词,加载词,修改词频,定义词库

“结巴”中文分词:做最好的 Python 中文分词组件,分词模块jieba,它是python比较好用的分词模块, 支持中文简体,繁体分词,还支持自定义词库。 ...

5.4K70
来自专栏AI科技评论

开发 | 如何利用 TVM 优化深度学习GPU op?教你用几十行Python代码实现2-3倍提升

数天前,陈天奇团队宣布推出 TVM,在微博上表示,「我们今天发布了 TVM,和 NNVM 一起组成深度学习到各种硬件的完整优化工具链,支持手机,cuda, op...

46180

扫码关注云+社区

领取腾讯云代金券