2019牛客暑期多校训练营（第九场）H Cutting Bamboos —— 区间更新主席树+查询第一个比x小（大）的值

This way

题意：

给你一些竹子，q个询问，问你从第l到第r个竹子，如果你要用y次砍完它，并且每次砍下来的长度是相同的，问你第x次砍在哪。

题解：

其实我是在写题解的时候才看的题目的。。
刚做完j队友就将两个式子丢给我说是这道题的做法，然后由于自己没看题目不知道输出的是位置，结束了才A掉。。
一个式子是$val\frac{sum}{y}\ast (y-x)$也就是前面砍了竹子的总长度。
主席树按照从左到右建root，竹子的高度为基础，区间更新。
需要求l-r区间中第一个小于等于它的位置，和竹子的总长，如果这个长度就是我们上面式子算出来的答案，那就输出位置，否则我们需要找到第一个大于等于它的位置，以及高度，因为这个高度一定是小于它位置的高度+1（理解一下就懂了），那么大-小也就是这个高度上竹子的数量，那么答案就是$p_{小}+\frac{val-小}{大-小}$
主席树由与不能push_down，那么我们在查询的时候就将flag参数化，也就是把父亲的flag用传参的方法传到儿子即可。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define lb long double
const int N=2e5+5;
int ls[N*40],rs[N*40],rt[N];
ll sum[N*40],flag[N*40],a[N];
int n,m,tot;
void update(int l,int r,int root,int last,int ql,int qr)
{
    ls[root]=ls[last];
    rs[root]=rs[last];
    sum[root]=sum[last]+(ll)(min(qr,r)-max(ql,l)+1);
    flag[root]=flag[last];
    if(ql<=l&&qr>=r)
    {
        flag[root]++;
        return ;
    }
    int mid=l+r>>1;
    if(mid>=ql)
        update(l,mid,ls[root]=++tot,ls[last],ql,qr);
    if(mid<qr)
        update(mid+1,r,rs[root]=++tot,rs[last],ql,qr);
}
ll pp;
ll q_s(int l,int r,int root,int last,lb val,ll f_r,ll f_l)
{
    if(sum[root]+f_r*(ll)(r-l+1)-sum[last]-f_l*(ll)(r-l+1)<=val)
    {
        pp=r;
        return sum[root]+f_r*(ll)(r-l+1)-sum[last]-f_l*(ll)(r-l+1);
    }
    if(l==r)
    {
        pp=l-1;
        return 0ll;
    }
    ll ans;
    int mid=l+r>>1;
    ll l_sum_root=sum[ls[root]]+(f_r+flag[root])*(ll)(mid-l+1);
    ll l_sum_last=sum[ls[last]]+(f_l+flag[last])*(ll)(mid-l+1);
    if(l_sum_root-l_sum_last>=val)
        return q_s(l,mid,ls[root],ls[last],val,f_r+flag[root],f_l+flag[last]);
    return l_sum_root-l_sum_last+q_s(mid+1,r,rs[root],rs[last],val-(l_sum_root-l_sum_last),f_r+flag[root],f_l+flag[last]);
}
ll q_b(int l,int r,int root,int last,lb val,ll f_r,ll f_l)
{
    if(l==r)
    {
        //printf("%lld\n",1ll*(sum[root]+flag[root]*(r-l+1)-sum[last]-flag[last]*(r-l+1)));
        return sum[root]+f_r-sum[last]-f_l;
    }
    int mid=l+r>>1;
    ll l_sum_root=sum[ls[root]]+(f_r+flag[root])*(ll)(mid-l+1);
    ll l_sum_last=sum[ls[last]]+(f_l+flag[last])*(ll)(mid-l+1);
    ll ans=0;
    if(l_sum_root-l_sum_last<=val)
        ans=l_sum_root-l_sum_last+q_b(mid+1,r,rs[root],rs[last],val-(l_sum_root-l_sum_last),f_r+flag[root],f_l+flag[last]);
    else
        ans=q_b(l,mid,ls[root],ls[last],val,f_r+flag[root],f_l+flag[last]);
    //printf("%lld\n",ans);
    return ans;
}
ll pre[N];
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
    {
        scanf("%lld",&a[i]),pre[i]=pre[i-1]+a[i];
        update(1,1e5,rt[i]=++tot,rt[i-1],1,a[i]);
    }
    while(m--)
    {
        int l,r;
        ll x,y;
        scanf("%d%d%lld%lld",&l,&r,&x,&y);
        ll sum=pre[r]-pre[l-1];
        double val=1.0*sum/y*(y-x);
        ll s=q_s(1,1e5,rt[r],rt[l-1],val,0,0);

        if(s==val)
            printf("%.10f\n",1.0*pp);
        else
        {
            ll b=q_b(1,1e5,rt[r],rt[l-1],val,0,0);
            printf("%.10f\n",1.0*(pp+(val-s)/(b-s)));
        }
    }
    return 0;
}
/*
10 4
2 8 5 1 10 5 9 9 3 5
8 9 19 76

*/