K - Krystalova‘s Trivial Problem Gym - 103940K

题目链接
题意:
1:u l r x 表示在[l,r]中加入x
2:q l r 表示询问[l,r]的x^2的和
题解:很明显的线段树
1:(x+k) ^ 2 = x ^ 2 + k ^ 2 + 2 * k * x
所以我们维护一个区间和和一个区间平方和来不断更新数据即可,就是这里的mod有点搞。
下面是WA21的代码:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
#define int long long
const int N=1e5+10,mod=1e9+7;
int a[N];
int n,q;
struct node
{
    int l,r;
    int sum;
    int copsum;
    int lazy;
}tr[N*4];
void pushup(int u)
{
    tr[u].sum=((tr[u<<1].sum+tr[u<<1|1].sum)%mod+mod)%mod;
    tr[u].copsum=((tr[u<<1].copsum+tr[u<<1|1].copsum)%mod+mod)%mod;
}
void pushdown(int u)
{
    node &root=tr[u],&left=tr[u<<1],&right=tr[u<<1|1];
    if(root.lazy)
    {
        left.sum=((left.sum+2*root.lazy%mod*left.copsum%mod+(left.r-left.l+1)*root.lazy%mod*root.lazy%mod)%mod+mod)%mod;
        left.copsum=((left.copsum+root.lazy%mod*(left.r-left.l+1)%mod)%mod+mod)%mod;
        left.lazy+=root.lazy;
        right.sum=((right.sum+2*root.lazy%mod*right.copsum%mod+(right.r-right.l+1)*root.lazy%mod*root.lazy%mod)%mod+mod)%mod;
        right.copsum=((right.copsum+root.lazy%mod*(right.r-right.l+1)%mod)%mod+mod)%mod;
        right.lazy+=root.lazy;
        root.lazy=0;
    }
} 
void build(int u,int l,int r)
{
    //printf("%lld %lld %lld\n",u,l,r);
    tr[u]={l,r,0,0,0};
    if(l==r)
    {
        tr[u].sum=a[l]*a[l]%mod;
        tr[u].copsum=a[l];
        tr[u].lazy=0;
        return ;
    }
    int mid=(l+r)>>1;
    build(u<<1,l,mid);
    build(u<<1|1,mid+1,r);
    pushup(u);
}
int query(int u,int l,int r)
{
    if(tr[u].l>=l&&tr[u].r<=r) return tr[u].sum;
    pushdown(u);
    int mid=tr[u].l+tr[u].r>>1;
    int sum=0;
    if(l<=mid) sum+=query(u<<1,l,r);
    if(r>mid) sum+=query(u<<1|1,l,r);
    return sum;
}
int query1(int u,int l,int r)
{
    if(tr[u].l>=l&&tr[u].r<=r) return tr[u].copsum;
    pushdown(u);
    int mid=tr[u].l+tr[u].r>>1;
    int sum=0;
    if(l<=mid) sum+=query1(u<<1,l,r);
    if(r>mid) sum+=query1(u<<1|1,l,r);
    return (sum+mod)%mod;
}
void modify(int u,int l,int r,int k)
{
    //printf("%lld %lld %lld\n",u,l,r);
    if(tr[u].l>=l&&tr[u].r<=r)
    {
        tr[u].lazy+=k;
        (tr[u].sum=(tr[u].sum%mod+tr[u].copsum*k%mod*2%mod+(tr[u].r-tr[u].l+1)*k%mod*k%mod)%mod+mod)%mod;
        tr[u].copsum=((tr[u].copsum%mod+k*(tr[u].r-tr[u].l+1)%mod)%mod+mod)%mod;
        return;
    }
    pushdown(u);
    int mid=tr[u].l+tr[u].r>>1;
    if(l<=mid) modify(u<<1,l,r,k);
    if(r>mid) modify(u<<1|1,l,r,k);
    pushup(u);
    return;
}
signed main()
{
    cin>>n>>q;
    for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
    build(1,1,n);
    // for(int j=1;j<=n;j++)
    // cout<<query(1,j,j)<<"==>"<<endl;
    for(int i=1;i<=q;i++)
    {
        char op[2];
        scanf("%s",op);
        if(*op=='u')
        {
            int l,r,k;
            scanf("%lld%lld%lld",&l,&r,&k);
            modify(1,l,r,k);
        }
        else
        {
            int l,r;
            scanf("%lld%lld",&l,&r);
            int ans=query(1,l,r);
            printf("%lld\n",(ans+mod)%mod);
        }
    }
    return 0;
}

下面是AC代码:

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
#define int long long
const int N=1e5+10,mod=1e9+7;
int a[N];
int n,q;
struct node
{
    int l,r;
    int sum;
    int copsum;
    int lazy;
}tr[N*4];
void pushup(int u)
{
    tr[u].sum=tr[u<<1].sum+tr[u<<1|1].sum;
    tr[u].sum%=mod;
    tr[u].copsum=tr[u<<1].copsum+tr[u<<1|1].copsum;
    tr[u].copsum%=mod;
}

void pushdown(int u) 
{
    node &root=tr[u],&left=tr[u<<1],&right=tr[u<<1|1];
    if(root.lazy)
    {   
        left.sum+=2*root.lazy%mod*left.copsum%mod;
        left.sum%=mod;
        left.sum+=(left.r-left.l+1)*root.lazy%mod*root.lazy%mod;
        left.sum%=mod;
        left.sum=(left.sum+mod)%mod;
        left.copsum+=root.lazy%mod*(left.r-left.l+1)%mod;
        left.copsum=(left.copsum+mod)%mod;
        left.lazy+=root.lazy;
        left.lazy%=mod;
        
        right.sum+=2*root.lazy%mod*right.copsum%mod;
        right.sum%=mod;
        right.sum+=(right.r-right.l+1)*root.lazy%mod*root.lazy%mod;
        right.sum%=mod;
        right.sum=(right.sum+mod)%mod;
        right.copsum+=root.lazy%mod*(right.r-right.l+1)%mod;
        
        right.copsum=(right.copsum+mod)%mod;
        right.lazy+=root.lazy;
        right.lazy%=mod;
        root.lazy=0;
    }
} 
void build(int u,int l,int r)
{
    //printf("%lld %lld %lld\n",u,l,r);
    tr[u]={l,r,0,0,0};
    if(l==r)
    {
        tr[u].copsum=a[l];
        tr[u].sum=a[l]*a[l]%mod;
        tr[u].lazy=0;
        return ;
    }
    int mid=(l+r)>>1;
    build(u<<1,l,mid);
    build(u<<1|1,mid+1,r);
    pushup(u);
}
int query(int u,int l,int r)
{
    if(tr[u].l>=l&&tr[u].r<=r) return tr[u].sum%mod;
    pushdown(u);
    int mid=tr[u].l+tr[u].r>>1;
    int sum=0;
    if(l<=mid) sum+=query(u<<1,l,r),sum%=mod;
    if(r>mid) sum+=query(u<<1|1,l,r),sum%=mod;
    return sum%mod;
}
void modify(int u,int l,int r,int k)
{
    //printf("%lld %lld %lld\n",u,l,r);
    if(tr[u].l>=l&&tr[u].r<=r)
    {
        tr[u].lazy+=k;
        tr[u].lazy%=mod;
        tr[u].sum+=tr[u].copsum*k%mod*2%mod;
        tr[u].sum%=mod;
        tr[u].sum+=(tr[u].r-tr[u].l+1)%mod*k%mod*k%mod;
        tr[u].sum%=mod;
        tr[u].sum=(tr[u].sum+mod)%mod;
        tr[u].copsum+=k*(tr[u].r-tr[u].l+1)%mod;
        tr[u].copsum%=mod;
        return;
    }
    pushdown(u);
    int mid=tr[u].l+tr[u].r>>1;
    if(l<=mid) modify(u<<1,l,r,k);
    if(r>mid) modify(u<<1|1,l,r,k);
    pushup(u);
    return;
}
signed main()
{
    cin>>n>>q;
    for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
    build(1,1,n);
    for(int i=1;i<=q;i++)
    {
        char op[2];
        scanf("%s",op);
        if(*op=='u')
        {
            int l,r,k;
            scanf("%lld%lld%lld",&l,&r,&k);
            modify(1,l,r,k);
        }
        else
        {
            int l,r;
            scanf("%lld%lld",&l,&r);
            int ans=query(1,l,r);
            printf("%lld\n",(ans+mod)%mod);
        }
    }
    return 0;
}

这里我们可以总结出就是那个我们在mod的过程中可以去每加一次,就mod一次,这样可能不容易超。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值