Gym 100971D Laying Cables 二分 || 单调栈
要求找出每个a[i],找到离他最近而且权值比它大的点,若距离相同,输出权利最大的那个
我的做法有点复杂,时间也要500+ms,因为只要时间花在了map上。
具体思路是模拟一颗树的建立过程,对于权值最大的那个,必须是-1,次大的那个,必须是pos_peo[mx];就是最大人口的节点id、
然后维护一个单调的序列,记录当前有多少个位置加入了树。用个set保证单调性。小到大
把结构体按人口排序大到小,枚举没个城市,这样保证加入后,后面加入的直接找位置最短即可,人口最对的bigger than now的。
二分一个位置,> now.pos的,枚举它左右,选择即可。注意就是当距离相同的时候,还要再判断一次。
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define inf (0x3f3f3f3f)
typedef long long int LL; #include <iostream>
#include <sstream>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <string>
const int maxn = + ;
map<int,int>pos_peo;
map<int,int>pos_pos;
struct data {
int pos,peo;
} a[maxn],b[maxn];
int ans[maxn]; struct cmp1 {
bool operator()(int a,int b) {
return a < b; //
}
};
bool cmp2 (data a,data b)
{
return a.peo > b.peo;
}
set<int,cmp1>SS;
void work ()
{
int n;
scanf ("%d",&n);
int mx = -inf;
for (int i=; i<=n; ++i) {
scanf ("%d%d",&a[i].pos,&a[i].peo);
b[i].pos=a[i].pos;
b[i].peo=a[i].peo;
pos_peo[a[i].peo] = i; //id
pos_pos[a[i].pos] = i;
mx = max(mx,a[i].peo);
}
if (n==) {
printf ("-1\n");
return ;
}
ans[pos_peo[mx]] = -;
int sec = -inf;
for (int i=; i<=n; ++i) {
if (sec < a[i].peo && a[i].peo != mx)
sec = a[i].peo;
}
ans[pos_peo[sec]] = pos_peo[mx]; SS.insert(a[pos_peo[mx]].pos);
SS.insert(a[pos_peo[sec]].pos); sort (a+,a++n,cmp2); // peo up to low set<int>::iterator it;
for (int i=; i<=n; ++i) {
int val = a[i].pos;
int ppeo = a[i].peo;
it = SS.lower_bound(val);
int t1 = inf,t2 = inf;
if (it == SS.begin()) { //在开头
ans[pos_peo[a[i].peo]] = pos_pos[*it];
} else if (it == SS.end()) { //在末尾
it --;
ans[pos_peo[a[i].peo]] = pos_pos[*it];
} else {
int tt1 = *it;
t1 = abs(val - (*it));
it--;
int tt2 = *it;
t2 = abs(val - (*it));
if (t1 < t2) {
ans[pos_peo[a[i].peo]] = pos_pos[tt1];
} else if (t1 > t2) {
ans[pos_peo[a[i].peo]] = pos_pos[tt2];
} else { //xiangdeng
int id2 = pos_pos[tt1];
int id1 = pos_pos[tt2];
int cut2 = abs(b[id2].peo - ppeo);
int cut1 = abs(b[id1].peo - ppeo);
if (cut2 > cut1) {
ans[pos_peo[a[i].peo]] = pos_pos[tt1];
} else {
ans[pos_peo[a[i].peo]] = pos_pos[tt2];
}
}
}
SS.insert(a[i].pos);
}
for (int i=; i<=n; ++i) {
printf ("%d ",ans[i]);
}
printf ("\n");
}
int main()
{
#ifdef LOCAL
freopen("data.txt","r",stdin);
#endif
work ();
return ;
}
其实这个可以用单调栈O(n)解决
首先,对于任何一个city,只有两种可能,选择在它左边的第一个城市,或者选择在它右边的第一个城市,当然这些城市都是要合法的。就是要满足人口数 > 当前城市。
所以首先对pos进行排序。这样压栈的时候就能知道相对位置了。用单调栈预处理ri[i]表示右边第一个合法城市。le[i]同理。比较即可。
为什么是第一个合法城市呢?第二个合法城市不行?这是因为距离要最短,要先满足距离。
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <algorithm>
using namespace std;
#define inf (0x3f3f3f3f)
typedef long long int LL; #include <iostream>
#include <sstream>
#include <vector>
#include <set>
#include <map>
#include <queue>
#include <string>
const int maxn = + ;
struct node {
int pos, val, id;
bool operator < (const node &rhs) const {
return pos < rhs.pos;
}
}a[maxn], ri[maxn], le[maxn];
int stack[maxn], ans[maxn];
void work ()
{
int n;
scanf ("%d", &n);
for (int i = ; i <= n; ++i) {
scanf ("%d%d", &a[i].pos, &a[i].val);
a[i].id = i;
}
sort (a + , a + + n);
int top = ;
for (int i = ; i <= n; ++i) {
while (top >= && a[i].val > a[stack[top]].val) --top;
++top;
stack[top] = i;
if (top != ) {
le[i] = a[stack[top - ]];
} else {
le[i].id = -inf;
}
}
top = ;
for (int i = n; i >= ; --i) {
while (top >= && a[i].val > a[stack[top]].val) --top;
++top;
stack[top] = i;
if (top != ) {
ri[i] = a[stack[top - ]];
} else {
ri[i].id = -inf;
}
}
for (int i = ; i <= n; ++i) {
int toans;
if (le[i].id == -inf) {
toans = ri[i].id == -inf ? - : ri[i].id;
} else if (ri[i].id == -inf) {
toans = le[i].id == -inf ? - : le[i].id;
} else {
int posL = abs(le[i].pos - a[i].pos);
int posR = abs(ri[i].pos - a[i].pos);
if (posL > posR) {
toans = ri[i].id;
} else if (posL == posR) {
if (le[i].val > ri[i].val) {
toans = le[i].id;
} else {
toans = ri[i].id;
}
} else {
toans = le[i].id;
}
}
ans[a[i].id] = toans;
}
for (int i = ; i <= n; ++i) {
printf ("%d ", ans[i]);
}
} int main ()
{
#ifdef local
freopen("data.txt","r",stdin);
#endif
work ();
return ;
}
单调栈
Gym 100971D Laying Cables 二分 || 单调栈的更多相关文章
- Code Forces Gym 100971D Laying Cables(单调栈)
D - Laying Cables Time Limit:2000MS Memory Limit:262144KB 64bit IO Format:%I64d & %I64u ...
- Gym 100971D Laying Cables 单调栈
Description One-dimensional country has n cities, the i-th of which is located at the point xi and h ...
- BZOJ 1012--[JSOI2008]最大数maxnumber(二分&单调栈)
1012: [JSOI2008]最大数maxnumber Time Limit: 3 Sec Memory Limit: 162 MBSubmit: 14142 Solved: 6049[Subm ...
- SPOJ MINSUB - Largest Submatrix(二分+单调栈)
http://www.spoj.com/problems/MINSUB/en/ 题意:给出一个n*m的矩阵M,和一个面积k,要使得M的子矩阵M'的最小元素最大并且面积大于等于k,问子矩阵M'的最小元素 ...
- 51Nod 1279:扔盘子(二分||单调栈)
1279 扔盘子 1.0 秒 131,072.0 KB 5 分 1级题 有一口井,井的高度为N,每隔1个单位它的宽度有变化.现在从井口往下面扔圆盘,如果圆盘的宽度大于井在某个高度的宽度,则圆盘被卡住( ...
- D - Laying Cables Gym - 100971D (单调栈)
题目链接:https://cn.vjudge.net/problem/Gym-100971D 题目大意:给你n个城市的信息,每一个城市的信息包括坐标和人数,然后让你找每一个城市的父亲,作为一个城市的父 ...
- Codeforces gym 100971 D. Laying Cables 单调栈
D. Laying Cables time limit per test 2 seconds memory limit per test 256 megabytes input standard in ...
- Gym 100971D 单调栈
D - Laying Cables Time Limit:2000MS Memory Limit:262144KB 64bit IO Format:%I64d & %I64u ...
- BZOJ1767/Gym207383I CEOI2009 Harbingers 斜率优化、可持久化单调栈、二分
传送门--BZOJCH 传送门--VJ 注:本题在BZOJ上是权限题,在Gym里面也不能直接看,所以只能在VJ上交了-- 不难考虑到这是一个\(dp\). 设\(dep_x\)表示\(x\)在树上的带 ...
随机推荐
- NSDictionary 用法
//Dictionary //不可变 //字典里面:是以键值对存放的 //字典存放顺序是无序的 //字典里面:value可以重复,但key不能重复 //字典里面:key一般用字符串表示,value可用 ...
- MATLAB模糊逻辑工具箱函数
说明:本文档中所列出的函数适用于Matlab5.3以上版本,为了简明起见,只列出了函数名,若需要进一步的说明,请参阅MATLAB的帮助文档. 1. GUI工具 Anfisedit 打开ANF ...
- WPF 导出Excel(合并单元格)
WPF 导出Excel(合并单元格) DataTable 导出Excel(导出想要的列,不想要的去掉) ,B1,B2,B3,B4,B5} MisroSoft.Office.Interop.Excel. ...
- BluetoothSetServiceState 函数
DWORD BluetoothSetServiceState( HANDLE hRadio, BLUETOOTH_DEVICE_INFO* pbtdi, GUID* pGuidService, DWO ...
- Mac系统的launchd、守护进程daemon(2013笔记整理)
1. launchd Mac系统下通用的进程管理器,是Mac系统下非常重要的一个进程,一般来说该进程不允许直接以命令行的形式调用.只能通过其控制管理界面,launchctl来进行控制. launchd ...
- Redis简介及基础知识
一.Redis简介 Redis是什么?两句话可以做下概括: 1. 是一个完全开源免费的key-value内存数据库 2. 通常被认为是一个数据结构服务器,主要是因为其有着丰富的数据结构 string( ...
- C语言入门题
1. 如下语句通过算术运算和逻辑运算之后 i 和 j 的结果是() int i=0; int j=0; if((++i>0)||(++j>0 ...
- 正则表达式需要匹配的内容本身就自带了html转义字符,需要转义,否则无法匹配
1.正则表达式需要匹配的内容本身就自带了html转义字符,需要转义,否则无法匹配 例如,对于"requestNo\":\"3b89957436eaacd8311535e0 ...
- 一个web.Config或app.Config自定义段configSections的示例--转
直接用配置文件app.Config或web.Config配置应用系统的运行参数,比自己做一个xml配置文件,简洁方便得多.这两个配置文件不仅有常见的connectionStrings和appSetti ...
- jexus处理静态文件(处理后缀)
AspNet_Exts=txt就能把你指定的扩展名交给asp.net处理.同理,可以写很多个,AspNet_Exts=txt,htm,html