Escape from the Hell

[JAG Asia 2016]

容易证明优先选择差值大的更优

对于最后一瓶我们可以枚举

枚举最后一瓶,然后在树状数组上消去它的影响,然后线段树check是否出现被追上的情况,即查询区间最小值。

需要用到两个线段树,因为当二分找到的位置在最后一瓶后面,需要在线段树上消去最后一瓶的影响。

特别注意当差值为负数的时候前缀和就没有单调性了,所以二分要在单调递增区间二分。

#include <bits/stdc++.h>

#define ll long long
using namespace std;
const int maxn = 1e5 + 7;
const ll inf = 0x3f3f3f3f3f3f3f3f;
int n;
int C[maxn];
struct node {
int a, b;
} s[maxn]; bool cmp(node x, node y) {
return x.a - x.b > y.a - y.b;
} ll c[maxn];
ll sum[maxn];
struct tree {
int l, r;
ll min1, min2;
} t[maxn << 2]; int lowbit(int x) {
return x & (-x);
} ll getsum(int i) {
ll res = 0;
while (i > 0) {
res += c[i];
i -= lowbit(i);
}
return res;
} void update(int i, ll val) {
while (i < maxn) {
c[i] += val;
i += lowbit(i);
}
} void build(int p, int l, int r) {
t[p].l = l, t[p].r = r;
if (l == r) {
t[p].min1 = getsum(l) - sum[l];
t[p].min2 = getsum(l) - sum[l - 1];
return;
}
int mid = (l + r) >> 1;
build(p << 1, l, mid);
build(p << 1 | 1, mid + 1, r);
t[p].min1 = min(t[p << 1].min1, t[p << 1 | 1].min1);
t[p].min2 = min(t[p << 1].min2, t[p << 1 | 1].min2);
} ll ask1(int p, int l, int r) {
if (l <= t[p].l && r >= t[p].r) return t[p].min1;
int mid = (t[p].l + t[p].r) >> 1;
ll val = inf;
if (l <= mid) val = min(val, ask1(p << 1, l, r));
if (r > mid) val = min(val, ask1(p << 1 | 1, l, r));
return val;
} ll ask2(int p, int l, int r) {
if (l <= t[p].l && r >= t[p].r) return t[p].min2;
int mid = (t[p].l + t[p].r) >> 1;
ll val = inf;
if (l <= mid) val = min(val, ask2(p << 1, l, r));
if (r > mid) val = min(val, ask2(p << 1 | 1, l, r));
return val;
} int erfen(int z, int y, ll x) {
int l = z, r = y;
while (l < r) {
int mid = (l + r) >> 1;
if (getsum(mid) >= x)r = mid;
else l = mid + 1;
}
return l;
} void dubug() {
for (int i = 1; i <= n; ++i) {
cout << ask1(1, i, i) << " " << ask2(1, i, i) << endl;
}
} int main() {
ll L;
scanf("%d%lld", &n, &L);
for (int i = 1; i <= n; ++i) {
scanf("%d%d", &s[i].a, &s[i].b);
}
sort(s + 1, s + 1 + n, cmp);
int k=-1;
for (int i = 1; i <= n; ++i) {
if(s[i].a-s[i].b>=0) update(i, s[i].a - s[i].b);
else if(k==-1){
k=i-1;
}
scanf("%d", &C[i]);
sum[i] = sum[i - 1] + C[i];
}
if(k==-1) k=n;
build(1, 1, n);
int minn = n + 1;
for (int i = 1; i <= n; ++i) {
if(s[i].a-s[i].b>0) update(i, s[i].b - s[i].a);
int pp = erfen(1, k + 1, L - s[i].a);
if (pp != k + 1) {
if (pp < i) {
if (ask1(1, 1, pp) > 0) {
minn = min(minn, pp + 1);
}
} else {
if (i==1||ask1(1, 1, i - 1) > 0) {
if (pp==1||ask2(1, i + 1, pp) - (s[i].a - s[i].b) > 0) {
minn = min(minn, pp);
}
}
}
}
if(s[i].a-s[i].b>0)
update(i, s[i].a - s[i].b);
}
if (minn == n + 1) {
printf("-1\n");
} else {
printf("%d\n", minn);
}
return 0;
}

Escape from the Hell的更多相关文章

  1. ACM: Gym 101047E Escape from Ayutthaya - BFS

    Gym 101047E Escape from Ayutthaya Time Limit:2000MS     Memory Limit:65536KB     64bit IO Format:%I6 ...

  2. 简单明了区分escape、encodeURI和encodeURIComponent

    一.前言 讲这3个方法区别的文章太多了,但是大部分写的都很绕.本文试图从实践角度去讲这3个方法. 二.escape和它们不是同一类 简单来说,escape是对字符串(string)进行编码(而另外两种 ...

  3. c#模拟js escape方法

    public static string Escape(string s) { StringBuilder sb = new StringBuilder(); byte[] ba = System.T ...

  4. 【BZOJ-1340】Escape逃跑问题 最小割

    1340: [Baltic2007]Escape逃跑问题 Time Limit: 5 Sec  Memory Limit: 162 MBSubmit: 264  Solved: 121[Submit] ...

  5. LYDSY热身赛 escape

    Description 给出数字N(1<=N<=10000),X(1<=x<=1000),Y(1<=Y<=1000),代表有N个敌人分布一个X行Y列的矩阵上矩形的行 ...

  6. javascript escape()函数和unescape()函数

    javascript escape()函数和unescape()函数 escape() 函数可对字符串进行编码,这样就可以在所有的计算机上读取该字符串. 语法: escape(string) stri ...

  7. HDU 3605 Escape(状压+最大流)

    Escape Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others) Total Sub ...

  8. escape,encodeURI,encodeURIComponent的区别

    escape是对字符串进行编码而另外两种是对URL. encodeURI方法不会对下列字符编码 ASCII字母 数字 ~!@#$&*()=:/,;?+'encodeURIComponent方法 ...

  9. C#针对js escape解码

    在javascript 中通常用escape与unescape进行编码以方便传输. 在asp.net页面接收到这些数据以后可以使用 Microsoft.JScript.GlobalObject.une ...

  10. escape(), encodeURI()和encodeURIComponent()(转)

      escape(), encodeURI()和encodeURIComponent()是在Javascript中用于编码字符串的三个常用的方法,而他们之间的异同却困扰了很多的Javascript初学 ...

随机推荐

  1. 服务器io资源查看

    资源查看命令安装 yum provides /usr/bin/find #查看命令是哪个软件包安装的 执行 yum provides */netstat 命令就可以看到提供命令的工具包net-tool ...

  2. Redis: Reducing Memory Usage

    High Level Tips for Redis Most of Stream-Framework's users start out with Redis and eventually move ...

  3. Mybatis之配置文件

    Mybatis_config.xml <?xml version="1.0" encoding="UTF-8"?><!DOCTYPE conf ...

  4. Oracle专题

    定时任务:JOB文件夹 监视执行过的SQL 语句 select * from v$sqlarea a where module='PL/SQL Developer' order by a.FIRST_ ...

  5. Vue 项目de一些准备工作

    1.安装node,同时也会自动安装npm,npm是node的一种包安装工具. 2.准备一个git,可以用来管理代码. 3.打开vue官网,可以使用vue-cli脚手架工作. 这里介绍一个element ...

  6. Bootstrap-模态框 modal.js

    参考网址:http://v3.bootcss.com/(能抄不写) 1.大模态框 图片效果图: 代码:(button的属性data-target对应的是具体模态框的class) <!-- Lar ...

  7. pytorch学习问题汇总

    问题六: 问题五:这里是怎么得到的? 问题四:为什么会是如下结果? torch.bernoulli(a)怎么是这个结果? 问题1:torch各个类型数据格式如何转换?数据类型在官方文档torch.Te ...

  8. zabbix监控tcp/nginx/memcache连接数自定义监控shell

    #!/bin/bashtcp_status_fun(){ TCP_STAT=$1 #netstat -n | awk '/^tcp/ {++state[$NF]} END {for(key in st ...

  9. SeetaFaceQt:Qt多线程

    为什么要做多线程,说个最简单的道理就是我们不希望在软件处理数据的时候界面处于无法响应的假死状态.有些处理是灰常花时间的,如果把这样的处理放到主线程中执行,就会导致软件一条路走到底,要等到处理完才能接收 ...

  10. 翻译——3_Gaussian Process Regression

    使用不同的机器学习方法进行预测 续上篇2_Linear Regression and Support Vector Regression 高斯过程回归 %matplotlib inline impor ...