2018牛客网暑期ACM多校训练营(第二场)G Transform(二分)
题意
在一个数轴上有n个集装箱,第 i 个集装箱的位置为x[i],且在集装箱内装有a[i]件货物,现在将这些集装箱内的货物进行移动(将一件货物从第 i 个集装箱移动到第 j 个集装箱的花费就为2*abs(x[i]-x[j]) ),求在总花费不超过T的情况下,最多能将多少货物移动到同一个集装箱内。
分析
既然要使得花费在不超过T的情况尽可能多的移动货物,那么我们肯定是将一个区间内的所有货物移到坐标中位的集装箱上。那么我们就可以对答案进行二分,然后枚举所要移动的区间的左端点,再找到中位点和右端点,然后判断这个区间移动的花费是否小于T。
可以预处理一下前缀和以及前缀花费、距离等。依次来计算区间的花费。
二分判断的过程中,由于最终的答案不一定是将整个区间内所有的货物都移动到一个集装箱,所以我们还要判断比需求量多的那部分是从左端点移过来的还剩从右端点移过来的,然后再根据花费的大小情况进行二分就可以得到最终答案了。
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cmath>
#include<string>
#include<cstring>
using namespace std;
#define ll long long
const int N = 5e5 + ;
ll dis[N]; ///dis[i]表示第i个集装箱距离x=0的距离
ll num[N]; ///num[i]表示第i个集装箱有的产品个数
ll sum[N]; ///sum[i]表示第1~i个集装箱有的产品数总和
ll cost[N]; ///cost[i]表示第1~i个集装箱的所有产品移到x=0处所"获得"的费用
ll tot,q;
int n, L, R, u; ///u为产品汇总的最优点
ll cost_l(int i) {///费用=将L+1~i中所有产品移动到i处所需费用+(将i~R中所有的产品移动到i处所需的费用 - 区间L+1~R比x多的产品数从R处移动到i处的费用)
return ((sum[i] - sum[L])*dis[i] - (cost[i] - cost[L])) + ((cost[R] - cost[i]) - (sum[R] - sum[i])*dis[i] - (sum[R] - sum[L] - q)*(dis[R] - dis[i]));
}
ll cost_r(int i) {///同上,只是规则是优先从右边界开始取产品
return -((sum[R] - sum[i])*dis[i] - (cost[R] - cost[i])) - ((cost[i-] - cost[L]) - (sum[i-] - sum[L])*dis[i] + (sum[R] - sum[L] - q)*(dis[i] - dis[L+]));
}
bool check(ll &x) { ///判断需求x是否能在所给的费用t内达到
q = x;
L = , R = , u = ;
while () { ///从左边界开始向右移动区间,优先取区间左边的产品
while (R < n&&sum[R] - sum[L] < x) R++;
if (sum[R] - sum[L] < x)break; ///若是当前的L~n无法满足x,那么L++也不可能满足
while (u < L)u++;
while (u < R&&cost_l(u)>cost_l(u + ))u++;
if (cost_l(u) <= tot)return true;
L++;
}
L = n - , R = n, u = n;
while () { ///从右边界开始向左移动区间,优先取区间右边的产品
while (L > && sum[R] - sum[L] < x) L--;
if (sum[R] - sum[L] < x)break; ///若是当前的L~R无法满足x,那么R--也不可能满足
while (u > R)u--;
while (u > L && cost_r(u) > cost_r(u - ))u--;
if (cost_r(u) <= tot)return true;
R--;
}
return false;
}
int main()
{
scanf("%d%lld", &n, &tot);
tot /= ;
for (int i = ; i <= n; i++)
scanf("%lld", &dis[i]);
sum[]=cost[]=;
for (int i = ; i <= n; i++) {
scanf("%lld", &num[i]);
sum[i] = sum[i - ] + num[i];
cost[i] = cost[i - ] + num[i] * dis[i];
}
ll l = ,r = sum[n] + ;
while (l + < r) {
ll mid = (l + r) >> ;
if (check(mid)) l = mid;
else r = mid;
}
printf("%lld\n", l);
return ;
}
2018牛客网暑期ACM多校训练营(第二场)G Transform(二分)的更多相关文章
- 2018牛客网暑期ACM多校训练营(第二场)I- car ( 思维)
2018牛客网暑期ACM多校训练营(第二场)I- car 链接:https://ac.nowcoder.com/acm/contest/140/I来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 ...
- 2018牛客网暑期ACM多校训练营(第一场)D图同构,J
链接:https://www.nowcoder.com/acm/contest/139/D来源:牛客网 同构图:假设G=(V,E)和G1=(V1,E1)是两个图,如果存在一个双射m:V→V1,使得对所 ...
- 2018 牛客网暑期ACM多校训练营(第一场) E Removal (DP)
Removal 链接:https://ac.nowcoder.com/acm/contest/139/E来源:牛客网 题目描述 Bobo has a sequence of integers s1, ...
- 2018牛客网暑期ACM多校训练营(第十场)A Rikka with Lowbit (树状数组)
链接:https://ac.nowcoder.com/acm/contest/148/A 来源:牛客网 Rikka with Lowbit 时间限制:C/C++ 5秒,其他语言10秒 空间限制:C/C ...
- 2018牛客网暑期ACM多校训练营(第十场)J Rikka with Nickname(二分,字符串)
链接:https://ac.nowcoder.com/acm/contest/148/J?&headNav=acm 来源:牛客网 Rikka with Nickname 时间限制:C/C++ ...
- 2018牛客网暑期ACM多校训练营(第二场)J Farm(树状数组)
题意 n*m的农场有若干种不同种类作物,如果作物接受了不同种类的肥料就会枯萎.现在进行t次施肥,每次对一个矩形区域施某种类的肥料.问最后枯萎的作物是多少. 分析 作者:xseventh链接:https ...
- 2018牛客网暑期ACM多校训练营(第一场)B Symmetric Matrix(思维+数列递推)
题意 给出一个矩阵,矩阵每行的和必须为2,且是一个主对称矩阵.问你大小为n的这样的合法矩阵有多少个. 分析 作者:美食不可负064链接:https://www.nowcoder.com/discuss ...
- 2018牛客网暑期ACM多校训练营(第三场) A - PACM Team - [四维01背包][四约束01背包]
题目链接:https://www.nowcoder.com/acm/contest/141/A 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言524288K ...
- 2018牛客网暑期ACM多校训练营(第五场) F - take - [数学期望][树状数组]
题目链接:https://www.nowcoder.com/acm/contest/143/F 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言524288K ...
- 2018牛客网暑期ACM多校训练营(第五场) E - room - [最小费用最大流模板题]
题目链接:https://www.nowcoder.com/acm/contest/143/E 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言524288K ...
随机推荐
- ISAP算法
为什么叫ISAP ISAP(Improved Shortest Augment Path):改进的最短增广路,属于增广路算法 算法 Dinic算法中,我们每次都需要BFS出层次图,而在ISAP中,我们 ...
- bzoj 3123 [Sdoi2013]森林(主席树+启发式合并+LCA)
Description Input 第一行包含一个正整数testcase,表示当前测试数据的测试点编号.保证1≤testcase≤20. 第二行包含三个整数N,M,T,分别表示节点数.初始边数.操作数 ...
- LG P2473 [SCOI2008]奖励关
题目链接:P2473 [SCOI2008]奖励关 题意:有n个宝物 每次等概率抛出其中之一一共抛出k次每个宝物有一个价值 和一个前提集合只有集齐了集合中的所有宝物 才可以领取这个宝物 范围:1 < ...
- 【BZOJ2829】[SHOI2012]信用卡凸包(凸包)
[BZOJ2829][SHOI2012]信用卡凸包(凸包) 题面 BZOJ 洛谷 题解 既然圆角的半径都是一样的,而凸包的内角和恰好为\(360°\),所以只需要把圆角的圆心弄下来跑一个凸包,再额外加 ...
- js处理异步的几种方式
Javascript语言的执行环境是"单线程"(single thread,就是指一次只能完成一件任务.如果有多个任务,就必须排队,前面一个任务完成,再执行后面一个任务,以此类推) ...
- 【php】php分隔字符串为数组
工作中会经常分隔字符串为数组,我们可以用php内置函数str_split(),可是有时候字符串中包含中文,切割后会乱码,比如 print_r(str_split('dw氛围fesf',3)); 输出 ...
- hdu 1686 Oulipo (kmp)
Problem Description The French author Georges Perec (1936–1982) once wrote a book, La disparition, w ...
- C# DateTimePicker控件获取他的年,月,日,时,分,秒
CustomFormat属性设置为: yyyy-MM-dd HH:mm:ss 记住还要修改一个属性值,DateFormat属性 可选项改为Custom,默认是Long
- 用lemon测交互题
题目类型:传统. 答案比较类型:逐行比较类型(忽略多余空格和制表符). 配置:交互. 编译器参数: -o %s %s.* ..\..\data\%s\judge.cpp -Wl,--stack= ju ...
- python面向对象编程 -- 基本概念
面向对象的编程简要概括就是将要处理的问题抽象为数据和操作的集合,用类对其进行封装.其中数据和操作都称为类的属性,它们是一般是不变的. 对类进行实例化生成我们所说的对象,对象有自己的属性.对象的属性一般 ...