2019牛客暑期多校訓練營(第九場)H-Cutting Bamboos 主席樹+二分

題意:有n個竹子,m個查詢,要求每次切割竹子到一定高度,且每次切割的總長度相同,y次切完,問你第x次切割之後的高度是多少

題解:每次切割掉的長度Σai/y,那麼第x次切割完,切割掉的總長度應該爲Σai*x/y,同時我們可以二分第x次切割完的高度,設二分的高度爲h,則切割掉的總長度=高於h的竹子的高度總和-h*高於h的竹子的數量,這個與Σai*x/y相比就可以了

// #pragma comment(linker, "/STACK:1024000000,1024000000")
 
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<cmath>
#include<vector>
#include<set>
#include<map>
#include<stack>
#include<queue>
#include<time.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
#define ms(s) memset(s, 0, sizeof(s))
const int INF = 0x3f3f3f3f;
const int maxn = 2e5+5;
const double eps=1e-8;
struct node {
	int l,r,sum;
    ll val;
}T[maxn*20];
int rt[maxn],a[maxn];
ll pre[maxn];
int n,m,tot;
void build(int &p,int l,int r){
	p=++tot;
	T[p].sum=0;T[p].val=0;
	if(l==r)return ;
	int mid=(l+r)/2;
	build(T[p].l,l,mid);
	build(T[p].r,mid+1,r);
}
void update(int l,int r,int &now,int last,int k){
	T[++tot]=T[last];
	now=tot;
	T[tot].sum++;T[tot].val+=k;
	if(l==r)return;
	int mid=(l+r)/2;
	if(k<=mid){
		update(l,mid,T[now].l,T[last].l,k);
	}else update(mid+1,r,T[now].r,T[last].r,k);
}
void query(int l,int r,int x,int y,int k,ll &ans1,ll &ans2){
    if(k<=l&&k<=r){
        ans1+=T[y].sum-T[x].sum;
        ans2+=T[y].val-T[x].val;
        return;
    }
    if(k>r)return;
	int mid=(l+r)/2;
	//int cnt=T[T[y].l].sum-T[T[x].l].sum;//因爲這裏查詢的是值,所以直接用mid來判斷左右區間
    query(l,mid,T[x].l,T[y].l,k,ans1,ans2); 
    query(mid+1,r,T[x].r,T[y].r,k,ans1,ans2);
}
int main() 
{
	tot=0;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
        pre[i]=pre[i-1]+a[i];
	}
	build(rt[0],1,100001);
	for(int i=1;i<=n;i++){
		update(1,100001,rt[i],rt[i-1],a[i]);
	}
    for(int i=1;i<=m;i++){
        int l,r,x,y;
        scanf("%d%d%d%d",&l,&r,&x,&y);
        double left=0,right=INF;
        for(int j=1;j<=100;j++){
            double mid=(left+right)/2;int midx=ceil(mid),flag=0;
            ll sum=0,val=0;
            //cout<<233<<endl;
            query(1,100001,rt[l-1],rt[r],midx,sum,val);
            double f1=1.0*val-1.0*sum*mid;
            double sumh=1.0*(pre[r]-pre[l-1])/y*x;
            if(f1<=sumh){
                right=mid;
            }else left=mid;
        }
        printf("%.8f\n",left);
    }

    

}

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章