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);
    }

    

}