hdu3473Minimum Sum【划分树sum维护】
一个上午唯一的成果T^T还有8天就过年了
说说这个题，比基础的多了一个sum数组，用来记录当前区间在这个数之前进入左子树的元素值之和，
/***************
hdu3473
2016.1.29
561MS 30344K 2956 B C++
***************/
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int MAXN = 1000010;
int sorted[MAXN]; //对原集合中元素排序后的值
int val[30][MAXN]; //val记录第k层当前位置的值
int toleft[30][MAXN]; //记录元素所在区间当前位置前的元素进入到左子树的个数
long long sum[30][MAXN]; //记录元素所在区间当前位置前的元素进入到左子树的元素和
int n,t,cnt,m,x,y;
long long lsum,rsum;
int lnum,rnum;
void build(int l, int r, int d) {
if (l == r) return ;
int mid = (l + r) >> 1;
int same = mid - l + 1;
for (int i=l; i<=r; i++)
if (val[d][i] < sorted[mid])
same--;
int lp = l, rp = mid+1;
for (int i=l; i<=r; i++) {
if (i == l) sum[d][i]=toleft[d][i] = 0;
else toleft[d][i] = toleft[d][i-1],sum[d][i]=sum[d][i-1];
if (val[d][i] < sorted[mid]) {
toleft[d][i]++;
val[d+1][lp++] = val[d][i];
sum[d][i]+=val[d][i];
}
else if (val[d][i] > sorted[mid])
{
val[d+1][rp++] = val[d][i];
}
else {
if (same) {
same--;
toleft[d][i]++;
val[d+1][lp++] = val[d][i];
sum[d][i]+=val[d][i];
}
else
{
val[d+1][rp++] = val[d][i];
}
}
}
build(l, mid, d+1);
build(mid+1, r, d+1);
}
int query(int a, int b, int k, int l, int r, int d) {
if (a == b) return val[d][a];
int mid = (l + r) >> 1;
int s, ss;
long long sss;
if (a == l) {
s = toleft[d][b];
ss = 0;
sss=sum[d][b];
}
else {
s = toleft[d][b] - toleft[d][a-1];
ss = toleft[d][a-1];
sss=sum[d][b]-sum[d][a-1];//为啥没管右区间的事？因为sum数组储存的是”当前区间此数之前进入左子树的元素和“
//和右区间没关系啊。
}
if (s >= k) {
a = l + ss;
b = l + ss + s - 1;
return query(a, b, k, l, mid, d+1);
}
else {
lnum+=s;//既然进入这个else了 说明前面进入左子树的个数小于k 那么k-th左侧的个数肯定要+[a,b]中进入左孩子的个数，即+s
lsum+=sss;//同理，k-th左侧数之和必然要加上（区间[a,b]中小于第k大元素值之和）即sss
a = mid+1 + a - l - ss;
b = mid+1 + b - l - toleft[d][b];
return query(a, b, k-s, mid+1, r, d+1);
}
}
long long s[MAXN];
int main()
{
//freopen("cin.txt","r",stdin);
cnt=1;
scanf("%d",&t);
while(t--)
{
printf("Case #%d:\n",cnt++);
scanf("%d",&n);
s[0]=0;
for(int i=1;i<=n;i++)
{
scanf("%d",&sorted[i]);
val[0][i]=sorted[i];
s[i]=s[i-1]+sorted[i];
}
sort(sorted+1,sorted+n+1);
build(1,n,0);
scanf("%d",&m);
while(m--)
{
scanf("%d%d",&x,&y);
x++;y++;
lsum=lnum=0;
int ave=query(x,y,(y-x+2)/2,1,n,0);
rnum=y-x+1-lnum;//这个区间里除了进入左边的就是进入右边的
rsum=s[y]-s[x-1]-lsum;//s数组储存的是截止到这个数为止的所有元素和
printf("%I64d\n",rsum-lsum+(long long)(lnum-rnum)*ave);//
}
puts("");
}
return 0;
}
本来想照着邝斌的思路改来着，结果一直RE
/***************
hdu3473
2016.1.29
***************/
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int MAXN = 200010;
int sorted[MAXN]; //对原集合中元素排序后的值
int val[30][MAXN]; //val记录第k层当前位置的值
int toleft[30][MAXN]; //记录元素所在区间当前位置前的元素进入到左子树的个数
long long sum[30][MAXN]; //记录比当前元素小的元素的和
int n,t,cnt,m,x,y;
long long lsum,rsum,ans;
int lnum,rnum;
void build(int l, int r, int d) {
if (l == r)
{
sum[d][l]=val[d][l];
return ;
}
int mid = (l + r) >> 1;
int same = mid - l + 1;
for (int i=l; i<=r; i++)
{
if (val[d][i] < sorted[mid])
same--;
sum[d][i]=val[d][i];
if(i>l) sum[d][i]+=sum[d][i-1];
}
int lp = l, rp = mid+1;
for (int i=l; i<=r; i++) {
if (i == l) toleft[d][i] = 0;
else toleft[d][i] = toleft[d][i-1];
if (val[d][i] < sorted[mid]) {
toleft[d][i]++;
val[d+1][lp++] = val[d][i];
// sum[d][i]+=val[d][i];
}
else if (val[d][i] > sorted[mid])
{
val[d+1][rp++] = val[d][i];
}
else {
if (same) {
same--;
toleft[d][i]++;
val[d+1][lp++] = val[d][i];
// sum[d][i]+=val[d][i];
}
else
{
val[d+1][rp++] = val[d][i];
}
}
}
build(l, mid, d+1);
build(mid+1, r, d+1);
}
int query(int a, int b, int k, int l, int r, int d) {
if (a == b) return val[d][a];
int mid = (l + r) >> 1;
int cnt=toleft[d][b]-toleft[d][a-1];
int ss=toleft[d][a-1]-toleft[d][l-1];
int ee=a-l-ss;
int s=toleft[d][b]-toleft[d][a-1];
int e=b-a+1-s;
if(cnt>=k)
{
if(e>0)
{
if(ee>0) ans+=sum[d+1][mid+e+ee]-sum[d+1][mid+ee];
else ans+=sum[d+1][mid+e];
}
int newl=l+toleft[d][a-1]-toleft[d][l-1];
int newr=newl+cnt-1;
return query(newl,newr,k,l,mid,d+1);
}
else
{
if(s>0)
{
if(ss>0) ans-=sum[d+1][l+ss+s-1]-sum[d+1][l+ss-1];
else ans-=sum[d+1][l+s-1];
}
int newr=b+toleft[d][r]-toleft[d][b];
int newl=newr-(b-a-cnt);
return query(newl,newr,k-cnt,mid+1,r,d+1);
}
}
long long s[MAXN];
int main()
{
//freopen("cin.txt","r",stdin);
cnt=1;
scanf("%d",&t);
while(t--)
{
printf("Case #%d:\n",cnt++);
scanf("%d",&n);
s[0]=0;
for(int i=1;i<=n;i++)
{
scanf("%d",&sorted[i]);
val[0][i]=sorted[i];
s[i]=s[i-1]+sorted[i];
}
sort(sorted+1,sorted+n+1);
build(1,n,0);
scanf("%d",&m);
while(m--)
{
scanf("%d%d",&x,&y);
x++;y++;
lsum=lnum=0;
//int ave=query(x,y,(y-x+2)/2,1,n,0);
rnum=y-x+1-lnum;
rsum=s[y]-s[x-1]-lsum;
ans=0;
int tt=query(x,y,(y-x)/2+1,1,n,0);
if((y-x+1)%2==0) ans-=tt;
printf("%I64d\n",ans);
}
puts("");
}
return 0;
}
最近状态不是一般的不好，虽然说每天也在学习，进度明显赶不上计划，可是再这样下去就学不完了啊啊啊啊啊，早睡早起