HDU 3473
分析可知,最优的x应该在区间中的数排列后最中间的地方选择。由于区间的数个数有奇偶之分,于是当奇数个时,就是中位数了。偶数个时,就是排列后中间两数区间的任意一个。
可以应用划分树求得前半部分,树状数组统计。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL __int64
#define lowbit(x) (x&(-x))
using namespace std; const int N=100050;
int sorted[N],n,m;
struct node{
int val[N];
int num[N];
LL sum[N];
}tp[20];
LL ANS;
LL arrayTree[N]; void build(int lft,int rht,int p){
if(lft==rht) return ;
int i,mid=(lft+rht)/2;
int isame=mid-lft+1,same=0;
for(i=lft;i<=rht;i++){
if(tp[p].val[i]<sorted[mid]) isame--;
}
int ln=lft,rn=mid+1;
for(i = lft; i <= rht; i++) {
if(i == lft) {
tp[p].num[i] = 0;
tp[p].sum[i]=0;
} else {
tp[p].num[i] = tp[p].num[i-1];
tp[p].sum[i]=tp[p].sum[i-1];
}
if(tp[p].val[i] < sorted[mid]) {
tp[p].num[i]++;
tp[p].sum[i]+=tp[p].val[i];
tp[p+1].val[ln++] = tp[p].val[i];
}else if(tp[p].val[i] > sorted[mid]) {
tp[p+1].val[rn++] = tp[p].val[i];
}else {
if(same < isame) {
same++;
tp[p].num[i]++;
tp[p].sum[i]+=tp[p].val[i];
tp[p+1].val[ln++] = tp[p].val[i];
}else {
tp[p+1].val[rn++] = tp[p].val[i];
}
}
}
build(lft, mid, p+1);
build(mid+1, rht, p+1);
} int query(int a, int b, int k, int p, int lft, int rht) {
if(lft == rht) {
return tp[p].val[a];
}
LL sss=0;
int s, ss, b2, bb, mid = (lft+rht)/2;
if(a == lft) {
s = tp[p].num[b];
ss = 0;
sss=tp[p].sum[b];
} else {
s = tp[p].num[b] - tp[p].num[a-1];
ss = tp[p].num[a-1];
sss=tp[p].sum[b]-tp[p].sum[a-1];
}
if(s >= k) {
a = lft + ss;
b = lft + ss + s - 1;
return query(a, b, k, p+1, lft, mid);
} else {
bb = a - lft - ss;
b2 = b - a + 1 - s;
a = mid + bb + 1;
b = mid + bb + b2;
ANS+=sss;
return query(a, b, k-s,p+1, mid+1, rht);
}
} void update(int i,int w){
for(int k=i;k<=n;k+=lowbit(k)){
arrayTree[k]+=w;
}
} LL arraySum(int index){
if(index==0) return 0;
LL ret=0;
for(int i=index;i>=1;i-=lowbit(i))
ret+=arrayTree[i];
return ret;
} LL getsum(int l,int r){
// cout<<arraySum(r)<<endl;
return arraySum(r)-arraySum(l-1);
} int main(){
int T,Case=0;
scanf("%d",&T);
while(T--){
scanf("%d",&n);
memset(arrayTree,0,sizeof(arrayTree));
for(int i=1;i<=n;i++){
scanf("%d",&tp[0].val[i]);
sorted[i]=tp[0].val[i];
update(i,sorted[i]);
}
ANS=0;
sort(sorted+1,sorted+1+n);
build(1,n,0);
scanf("%d",&m);
int l,r;
printf("Case #%d:\n",++Case);
for(int i=1;i<=m;i++){
scanf("%d%d",&l,&r);
l++,r++;
LL spre=getsum(l,r);
ANS=0;
if((r-l+1)%2==0){
int mid=query(l,r,(r-l+1)/2+1,0,1,n);
printf("%I64d\n",spre-ANS-ANS);
}
else{
int mid=query(l,r,(r-l+1)/2+1,0,1,n);
printf("%I64d\n",spre-ANS-ANS-mid);
}
}
puts("");
}
return 0;
}
HDU 3473的更多相关文章
- hdu 3473 划分树
Minimum Sum Time Limit: 16000/8000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) Tot ...
- HDU 3473 Minimum Sum 划分树,数据结构 难度:1
http://acm.hdu.edu.cn/showproblem.php?pid=3473 划分树模板题目,需要注意的是划分树的k是由1开始的 划分树: 参考:http://blog.csdn.ne ...
- hdu 3473 Minimum Sum
传送门 之前看挑战的时候看到一道分桶法的题目,其实我不是很明白分桶法应该怎么写.看到poj后面的讨论版上写着划分树裸题,而我以前就听说过了划分树,就干脆拿来学习一下.在写这篇博客的时候,其实我还是对这 ...
- hdu 3473 区间条件极值 - 区间 差的绝对值 之和的最小
题目传送门//res tp hdu 目的 对长度为n的区间,给定q个子区间,求一x,使得区间内所有元素与x的差的绝对值之和最小. 多测. n 1e5 q 1e5 ai [1,1e9] (i∈[1,n] ...
- hdu 3473 裸的划分树
思路: 用Sum[dep][i]记录从tree[po].l到i中进入左子树的和. #include<iostream> #include<algorithm> #include ...
- HDU 3473 Minimum Sum (划分树)
题意:给定一个数组,有Q次的询问,每次询问的格式为(l,r),表示求区间中一个数x,使得sum = sigma|x - xi|最小(i在[l,r]之间),输出最小的sum. 思路:本题一定是要O(nl ...
- HDU 3473 Minimum Sum(划分树)
Minimum Sum Time Limit: 16000/8000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Tota ...
- hdu 3473 (划分树)2
Minimum Sum Time Limit: 16000/8000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)Tota ...
- HDU 3473 Minimum Sum (划分树求区间第k大带求和)(转)
题意:在区间中找一个数,求出该区间每个数与这个数距离的总和,使其最小 找的数字是中位数(若是偶数个,则中间随便哪个都可)接着找到该区间比此数大的数的总和 区间中位数可以使用划分树,然后在其中记录:每层 ...
- HDU 3473 Minimum Sum 划分树
题意: 给出一个长度为\(n(1 \leq n \leq 10^5)\)的序列\(a\) 有若干次查询l r:找到一个\(x\)使得\(\sum \limits_{l \leq i \leq r} \ ...
随机推荐
- Windows(7/8/10)搭建kibana 6.x版本(elasticsearch的可视化服务)
在搭建kibana之前,我们先了解下什么是kibana Kibana 是一款开源的数据分析和可视化平台,它是 Elastic Stack 成员之一,设计用于和 Elasticsearch 协作.您可以 ...
- 安卓中Canvas实现清屏效果
可以在代码里面添加: paint.setXfermode(new PorterDuffXfermode(PorterDuff.Mode.CLEAR)); canvas.drawPaint(paint) ...
- 15年用canvas画的
请恕我当年的工作太轻松,用canvas手打了一个图,技术含量并没有什么,现在看看,甚是怀念_(¦3」∠)_ <!DOCTYPE html> <html> <head&g ...
- Laravel5.1学习笔记13 系统架构5 Contract
Contract 简介 为什么要用 Contract? Contract 参考 如何使用 Contract 简介 Laravel 中的 Contract 是一组定义了框架核心服务的接口.例如,Illu ...
- Oracle12C用户创建、授权、登录
Oracle12C用户创建.授权.登录 1.以系统用户登录 C:\Users\LEI>sqlplus sys/dwh as sysdba; SQL*Plus: Release 12.1.0.2. ...
- JS高级——面向对象方式解决tab栏切换问题
注意事项 1.给li元素注册事件,函数里面的this指的li元素,那么我们可以在注册事件之前将Tab对象用that=this进行保存 2.使用沙箱模式,所以暴露给外面的变量使用的是window.tab ...
- 【技术累积】【点】【java】【28】Map遍历
Map遍历 map的遍历一般有几种吧 遍历entrySet for(Map.Entry<String,String> entry : map.entrySet()){ } Iterator ...
- js常用方法和技巧
随着AJAX的流行,js又得到了很多人的重视,js最大的优势就是它能够对html上的所有元素进行操作,包括创建标签元素,更改元素属性等,这样就使得我们能够利用js来实现很多的动态效果,来提供给用户更强 ...
- Windows系统设置与北京Internet时间同步
找到提供Ntp服务器的网址:http://support.ntp.org/我从中找到东北大学的: ntp.neu.edu.cn ( 202.118.1.46 ) 长话短说,下面开始设置: 修改后的值 ...
- 用shell编写dhcp自动获取脚本
#!/bin/bash#net=$(ifconfig ens33 | awk -F'[ .]+' '/inet\>/{print $3"."$4"."$5 ...