Problem B: 专家系统

Description

一个专家系统是指,你雇佣了\(n\)个专家,他们每个人会做出一个结果,然后你从中选取较多的专家的结果组合而成最终的结果。专家系统广泛应用于传统机器学习领域、决策领域以及老师找学生做出一套试题的答案等等。

现在我们要在平面上找到一个点的坐标\((x, y)\),于是我们雇佣了\(n\)个专家,每个人都会给出一个坐标。为了得出一个较为准确的结果,我们选取其中的\(k\)名专家,他们得出的最小\(x\)坐标记为\(x_{\min}\),最小y坐标记为\(y_{\min}\),最大\(x\)坐标记为\(x_{\max}\),最大\(y\)坐标记为\(y_{\max}\),那么我们有理由相信,最终的这个点坐标一定落在\(x_{\min}\le x\le x_{\max}\), \(y_{\min}\le y\le y_{\max}\)之内。

但是选取不同的\(k\)名专家可能会导致不同的结果。为此,我们定义一个不确定度\(c=\max(x_{\max}-x_{\min},y_{\max}-y_{\min})\)

我们希望选出的\(k\)名专家使得\(c\)值尽可能的小,这样才能比较精确地获取这个点的坐标范围。现在,这个任务交给了你。

简述

现在有\(n\)个坐标\((x_i, y_i)\),你要从中选出\(k\)个。

假设你选出的全部\(k\)个坐标中,\(x\)坐标最小值为\(x_{\min}\),\(x\)坐标最大值为\(x_{\max}\),\(y\)坐标最小值为\(y_{\min}\),\(y\)坐标最大值为\(y_{\max}\)。

那么你得出的不确定度\(c=\max(x_{\max}-x_{\min},y_{\max}-y_{\min})\)

你的目的是让\(c\)尽可能小。

由于可能有很多种选法能够达成这一目的,你只需要输出这个最小的\(c\)即可。

Input

第一行两个正整数\(n,k\), 空格分隔。

接下来\(n\)行,每行\(2\)个整数\(x_i, y_i\),空格分隔,表示一个坐标。

Output

仅一行,一个数,表示最小的\(c\)。

HINT

对于\(20\%\)的数据, \(n\le 50\)。

对于\(40\%\)的数据,\(n\le 300\)。

对于\(60\%\)的数据,\(n\le 2000\)。

对于\(80\%\)的数据,\(n\le 20000\)。

对于\(100\%\)的数据,\(1\le n\le 100000\)。

对于\(100\%\)的数据,\(1\le k\le n\),\(x_i\)和\(y_i\)的绝对值不超过\(2000000000\)


首先吐槽,开2e9的范围..

考试的时候写了个看起来很套路的60分,二分答案然后拿两个指针按\(x\)扫描一下,然后朴素维护一下\(y\),结果因为不明原因爆0,现在都不晓得为啥。

事实上这个题非常的noip,考虑扫描\(x\)区间的时候如何维护\(y\)

每个\(y\)都被选上的时候可以钦定Ta为权值的右端点,设这个权值为\(d\),然后当前二分的值为\(c\),那么所有\(y\)值在\([d-c,d]\)的点都可以被选上。那么我们就可以考虑每个点对其他值域的贡献,其实就是\([d,d+c]\),贡献是个数\(1\)

然后我们可以离散化\(y\)坐标,把这些东西放在线段树里面,就是维护一个区间加和最大值就可以了。


Code:

#include <cstdio>
#include <algorithm>
const int N=1e5+10;
#define ll long long
const ll inf=(1ll<<40);
struct node
{
ll x,y;
bool friend operator <(node a,node b){return a.x<b.x;}
}dx[N];
int mx[N<<2],tag[N<<2];
using std::max;
#define ls id<<1
#define rs id<<1|1
void build(int id,int l,int r)
{
mx[id]=tag[id]=0;
int mid=l+r>>1;
if(l^r) build(ls,l,mid),build(rs,mid+1,r);
}
void pushdown(int id)
{
if(tag[id])
{
tag[ls]+=tag[id],tag[rs]+=tag[id];
mx[ls]+=tag[id],mx[rs]+=tag[id];
tag[id]=0;
}
}
void change(int id,int L,int R,int l,int r,int d)
{
if(l==L&&r==R)
{
tag[id]+=d,mx[id]+=d;
return;
}
pushdown(id);
int Mid=L+R>>1;
if(r<=Mid) change(ls,L,Mid,l,r,d);
else if(l>Mid) change(rs,Mid+1,R,l,r,d);
else change(ls,L,Mid,l,Mid,d),change(rs,Mid+1,R,Mid+1,r,d);
mx[id]=max(mx[ls],mx[rs]);
}
int n,m,k;ll dy[N];
bool check(ll d)
{
build(1,1,m);
int l=1,r=1,p;
p=std::upper_bound(dy+1,dy+1+m,dy[dx[l].y]+d)-dy-1;
change(1,1,m,dx[l].y,p,1);
while(r<n&&dx[r+1].x-dx[l].x<=d)
{
++r;
p=std::upper_bound(dy+1,dy+1+m,dy[dx[r].y]+d)-dy-1;
change(1,1,m,dx[r].y,p,1);
}
if(mx[1]>=k) return true;
while(r<n)
{
++r;
while(dx[r].x-dx[l].x>d)
{
p=std::upper_bound(dy+1,dy+1+m,dy[dx[l].y]+d)-dy-1;
change(1,1,m,dx[l++].y,p,-1);
}
p=std::upper_bound(dy+1,dy+1+m,dy[dx[r].y]+d)-dy-1;
change(1,1,m,dx[r].y,p,1);
if(mx[1]>=k) return true;
}
return false;
}
int main()
{
scanf("%d%d",&n,&k);
for(int i=1;i<=n;i++) scanf("%lld%lld",&dx[i].x,&dx[i].y),dy[i]=dx[i].y;
std::sort(dx+1,dx+1+n);
std::sort(dy+1,dy+1+n);
m=std::unique(dy+1,dy+1+n)-dy-1;
for(int i=1;i<=n;i++) dx[i].y=std::lower_bound(dy+1,dy+1+m,dx[i].y)-dy;
ll l=1,r=4000000000ll;
while(l<r)
{
ll mid=l+r>>1;
if(check(mid)) r=mid;
else l=mid+1;
}
printf("%lld\n",l);
return 0;
}

2018.12.31

Problem B: 专家系统 解题报告的更多相关文章

  1. ZOJ Problem Set - 1025解题报告

    ZOJ Problem Set - 1025 题目分类:基础题 原题地址:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=10 ...

  2. ACM: A Simple Problem with Integers 解题报告-线段树

    A Simple Problem with Integers Time Limit:5000MS Memory Limit:131072KB 64bit IO Format:%lld & %l ...

  3. BestCoder18 1002.Math Problem(hdu 5105) 解题报告

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5105 题目意思:给出一个6个实数:a, b, c, d, l, r.通过在[l, r]中取数 x,使得 ...

  4. Problem - 433C - Codeforces解题报告

    对于这题本人刚开始的时候的想法是:先把最大两数差的位置找到然后merge计算一个值再与一连串相同的数做merge后计算一个值比较取最大值输出:可提交后发现不对,于是本人就搜了一下正解发现原来这题的正确 ...

  5. Problem A: 选举 解题报告

    Problem A: 选举 题意 给出一个投票过程.有\(n\)个选民和\(m\)个候选人,每个选民\(i\)有个不重且有序的可投集合\(\{a_i\}\). 对于第一轮投票,选民\(i\)会投给\( ...

  6. Problem A: 种树 解题报告

    Problem A: 种树 Description 很久很久以前,一个蒟蒻种了一棵会提问的树,树有\(n\)个节点,每个节点有一个权值,现在树给出\(m\)组询问,每次询问两个值:树上一组点对\((x ...

  7. Problem C: 多线程 解题报告

    Problem C: 多线程 Description 多线程是一种常见的加速手段,利用多个线程同时处理不同的任务可以一定程度上减少总耗时,达到提高效率的目的.然而,多个线程间的执行顺序是完全不可控的, ...

  8. Problem A: 踢罐子 解题报告

    Problem A: 踢罐子 Description 平面上有\(n\)个点,其中任意2点不重合,任意3点不共线. 我们等概率地选取一个点A,再在剩下的\(n-1\)个点中等概率地选取一个点B,再在剩 ...

  9. Problem C Dist 解题报告

    Problem C Dist Description 有一个\(n\)个点带边权的连通无向图,边集用\(k\)个集合\(s_1,s_2,\dots,s_k\)和\(k\)个整数\(w_1,w_2,\d ...

随机推荐

  1. linux下安装redis安装使用

    1.下载redis 下载地址:http://redis.io/download,下载最新稳定版本 2.解压redis 1)  cd redis-x.x.x 2) make 3.启动redis 1) c ...

  2. WPF Get jiayuan outbox list(send mail box)

      Request URL: http://www.jiayuan.com/msg/outbox/list.php Request Method: POST form data: type=all&a ...

  3. MVVM Light Toolkit使用指南

    原文:MVVM Light Toolkit使用指南 原文地址:  https://blog.csdn.net/ldld1717/article/details/77040077 概述 MVVM Lig ...

  4. JavaEE笔记(十四)

    #SSH配置文件整合笔记实例 spring-BaseBean.xml <?xml version="1.0" encoding="UTF-8"?> ...

  5. 奔跑吧Linux

    刚拿到MiZ702,就被他的"外貌"深深的吸引,核心板加底板的形式让她看上去,强大而神秘-- 华丽的外表之下是否有着与之相当的内含呢,我们拿Linux将其检验一番! 板载的TF卡里 ...

  6. mfc 类型间的强制转换

    一. static_cast运算符 用法:static_cast < type-id > ( expression ) 该运算符把expression 转换为type-id类型,但没有运行 ...

  7. 01-BAT算法特训班

  8. HTML基础之JS

    HTML中的三把利器的JS 又称为JavaScript,看着好像和Java有点联系,实际上他和java半毛钱关系都没有,JavaScript和我们学习的Python.Go.Java.C++等,都是一种 ...

  9. 阿里云rds 备份和还原

    阿里云rds 备份和还原 转发:https://www.cnblogs.com/lin1/p/8617764.html 转发:https://help.aliyun.com/knowledge_det ...

  10. Linux/Mac 挂载远程服务器目录到本地

    1. 安装 sudo apt-get installsshfs 2. 创建SSHFS 挂载目录 sudo mkdir/mnt/siyuan 3.使用SSHFS 挂载远程的文件系统 sudo sshfs ...