題意:有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);
}
}