質問の意味：
そこ\（N- \）の長さ\（a_iを\）竹の、ある\（m個\）操作は、各操作は次の4つの数字、\（Lは、R、X、Y \） 。今、あなたは間隔よ\（[L、R] \ ） 竹のカット\（Y \）の各カットのために、あなたは高さを選択する必要があり、ナイフを\（H \） 、あなたは間隔が必要になります\を（[L、R] \）以上の高さに\（H \）竹のカット、最終的にあなたは、すべてが同じ高さをknifedことを確認する必要があります。今すぐあなたの最初尋ねる\（X \）ナイフが必要な高さをカット。
分析：
質問の意味のための1つの\（クエリ（Lは、R、X、Y）\） 、各ナイフは、必要な高さに固定されているカット\（\ FRAC {\ sum_ { i = L} ^ {R} Y}} {a_iを\）。我々は最初にあるときに、\（X \）動作中に、それが切断されていなければならない\（\ FRAC {\ sum_ { i = L} ^ {R} a_iを* X} {Y} \） の高さ。そして、我々は、操作の高さを切るたびに単調に増加するので、の我々は可能性として半分なければならないので、\（X \）回必要なカットの高さ\（ヘイ\） 。
現在決定するために、\（丙を\）正当である、我々は現在の高さに応じてのみ必要（\開平）\全高が遮断された計算する前に、\（sumhei \）を比較することにより、最終的に、および\（sumheiを\）と\ （\ FRAC {\ sum_ {I = L} ^ {R} a_iを* X} {Y} \） の間のサイズ。
そして必要な\（sumhei \） 、我々は間隔知る必要が\（[L、R]を\ ） に等しい高さよりも大きい、\（特開平\）ツリーの数\（サイズ\）と重量\（ヴァル\）、とlet \（ヴァル・サイズ*の丙の\ ） の答えです。数や間隔のサイズよりも大きい数の要件は、我々は唯一のツリーのメンテナンスの会長に目を向ける必要があります。
大統領バイナリツリークエリセットので、全体の時間に複雑である\（\ mathcal {O}（ nlogn ^ 2）\）
コード：
#include <bits/stdc++.h>
#define maxn 200005
using namespace std;
typedef long long ll;
struct ST{
int l,r,sz;
ll val;
}tr[maxn*30];
int T[maxn],tot,n,m,a[maxn];
ll sum[maxn];
int update(int l,int r,int pre,int pos,int val){
int rt=++tot;
tr[rt]=tr[pre];
tr[rt].val+=val;
tr[rt].sz++;
if(l==r) return rt;
int mid=(l+r)>>1;
if(pos<=mid) tr[rt].l=update(l,mid,tr[rt].l,pos,val);
else tr[rt].r=update(mid+1,r,tr[rt].r,pos,val);
return rt;
}
void query(int l,int r,int rt,int pre,int K,int &ans1,ll &ans2){
if(l==r){
ans1+=tr[rt].sz-tr[pre].sz;
ans2+=tr[rt].val-tr[pre].val;
return ;
}
int mid=(l+r)>>1;
if(K<=mid){
query(l,mid,tr[rt].l,tr[pre].l,K,ans1,ans2);
ans1+=tr[tr[rt].r].sz-tr[tr[pre].r].sz;
ans2+=tr[tr[rt].r].val-tr[tr[pre].r].val;
}else query(mid+1,r,tr[rt].r,tr[pre].r,K,ans1,ans2);
return;
}
int main()
{
scanf("%d%d",&n,&m);
for(int i=1;i<=n;i++){
scanf("%d",&a[i]);
sum[i]=sum[i-1]+a[i];
T[i]=update(1,100005,T[i-1],a[i],a[i]);
}
for(int i=1;i<=m;i++){
int l,r,x,y;
scanf("%d%d%d%d",&l,&r,&x,&y);
double tmp=1.0*(sum[r]-sum[l-1])/y*x;
double ll=0,rr=100005;
for(int j=0;j<100;j++){
double mid=(ll+rr)/2;
int tt=ceil(mid),res1=0;
long long res2=0;
query(1,100005,T[r],T[l-1],tt,res1,res2);
double ans=res2-res1*mid;
if(ans<tmp) rr=mid;
else ll=mid;
}
printf("%.8f\n",ll);
}
return 0;
}