2020牛客寒假算法基础集训营3 J题（线段树）

有n个k值和n个b值，m次询问+修改，设f[i][x]=k[i]*x+b[i]
操作1:修改k[i],b[i]
操作2:给一段区间[l,r]，输出f[r][f[r-1][…f[l]]…]的值
人生第一颗自建的线段树，第一次写博客，挖土挖了2个小时，对拍1小时，debug了3个小时，最后发现是ask函数写错了
题解:瞎jb乱推可以得到一个这样的式子

然后开始瞎搞

在这里#include<bits/stdc++.h>
#include<ctime>
#define ll long long
#define ull unsigned long long
#define ld long double
#define rep1(i,n) for(ll i=1;i<=n;i++)
#define rep0(i,n) for(ll i=0;i<n;i++)
#define rep(i,l,r) for(ll i=l;i<=r;i++)
using namespace std;
const ll mod=1e9+7;
const ll maxn=2e5+7;
ll b[maxn],k[maxn];
struct SegmentTree
{
    ll l,r;
    ll pk,pb;
}t[maxn*4];

void build(ll p,ll l,ll r)
{
    t[p].l=l,t[p].r=r;
    if(l==r){t[p].pk=k[l],t[p].pb=b[l];return;}
    ll mid=(l+r)/2;
    build(p*2,l,mid);
    build(p*2+1,mid+1,r);
    t[p].pk=1ll*t[p*2].pk*t[p*2+1].pk%mod;
    t[p].pb=(1ll*t[p*2].pb*t[p*2+1].pk%mod+t[p*2+1].pb)%mod;
}

void change(ll p,ll x,ll k1,ll b1)
{
    if(t[p].l==t[p].r){t[p].pk=k1,t[p].pb=b1;return;}
    ll mid=(t[p].l+t[p].r)/2;
    if(x<=mid)change(p*2,x,k1,b1);
    else change(p*2+1,x,k1,b1);
    t[p].pk=1ll*t[p*2].pk*t[p*2+1].pk%mod;
    t[p].pb=(1ll*t[p*2].pb*t[p*2+1].pk%mod+t[p*2+1].pb)%mod;
}

pair<ll,ll> Merge(pair<ll,ll> a,pair<ll,ll> b)
{
    return {a.first*1ll*b.first%mod,(1ll*a.second*b.first%mod+b.second)%mod};
}

pair<ll,ll> ask(ll p,ll l,ll r)
{
    if(l<=t[p].l&&r>=t[p].r)return {t[p].pk,t[p].pb};
    ll mid=(t[p].l+t[p].r)/2;
    pair<ll,ll>ans1=make_pair(1,0),ans2=make_pair(1,0);
    //这里的2个if我一值用了return，所以一直错
    if(l<=mid)ans1=ask(p*2,l,r);
    if(r>mid)ans2=ask(p*2+1,l,r);
    return Merge(ans1,ans2);
}

int main()
{
//    freopen("in.txt","r",stdin);
    ll n,m;scanf("%lld%lld",&n,&m);
    ll op,x,y,z;
    for(ll i=1;i<=n;i++)scanf("%lld",&k[i]);
    rep1(i,n)scanf("%lld",&b[i]);
    build(1,1,n);
    while(m--)
    {
        scanf("%lld%lld%lld",&op,&x,&y);
        if(op==2)
        {
            pair<ll,ll> ans=ask(1,x,y);
            printf("%lld\n",(ans.first+ans.second)%mod);
        }
        else scanf("%lld",&z),change(1,x,y,z);
    }
    return 0;
}