BZOJ3922 Karin的弹幕 [线段树]

F.A.Qs	Home	Discuss	ProblemSet	Status	Ranklist	1 Contest	入门OJ	ModifyUser   Sakura_Lemon	Logout	捐赠本站

Notice:九月份月赛将于本周五晚六点到十一点进行，欢迎大家参加，鸣谢Claris主持，欢迎进群633587999交流！

Problem 3922. – Karin的弹幕

3922: Karin的弹幕

Time Limit: 10 Sec   Memory Limit: 512 MB
Submit: 607   Solved: 231
[ Submit][ Status][ Discuss]

Description

Karin在战斗之余的闲暇时光里喜欢上B站看鬼畜视频，尤其喜欢发弹幕。她这天对一个视频的弹幕产生了兴趣，她记录了每个时间点的弹幕数量，并且可能对一段呈等差数列的时间的弹幕数量求最大值；她还可能修改某个时间点的弹幕数量。

为了不在Yuuna面前暴露出她弱爆了的数学能力，保持她傲娇的属性，你需要帮助她。

精简题意：给定一个序列，支持以下操作：①对一段下标是等差数列的子序列进行求最大值操作(参见输入格式)；②单点修改。

Input

第一行是一个整数n，

第二行是一个长度为n的整数序列a1…an，

第三行是一个整数m，

接下来m行，每行首先有一个整数op，

然后，若op=0，则之后有两个整数p，v，代表将a[p]的值加上v，

若op=1，则之后有两个整数x0，d，代表询问max{a[x0],a[x0+d],a[x0+2d],…,a[x0+kd]}(x0+kd<=n,x0+(k+1)d>n)。

数据中可能有多余空格。

Output

对每个op=1，单独输出一行，代表该等差子序列的最大值。

Sample Input

【输入样例1】

10

1 6 1 4 9 4 8 2 8 5

10

1 3 3

0 5 4

0 3 8

1 2 5

1 4 8

1 7 5

1 3 6

0 1 2

1 5 3

1 4 9



【输入样例2】

10

-9 -6 2 -10 -2 -6 10 6 -4 -2

10

1 2 3

1 6 3

0 7 8

0 4 -6

0 10 -5

1 10 4

0 3 -8

1 2 4

0 10 -5

1 1 2

Sample Output

【输出样例1】

8

8

4

8

9

13

4



【输出样例2】

6

-4

-7

-6

18

HINT

【数据范围】

1<=n<=70000，

1<=m<=70000，

保证任何时刻abs(a[i])(1<=i<=n)<=2147483647，

0<=op<=1，

1<=p<=n，

abs(v)<=2147483647

1<=x0<=n,

1<=d<=n，

保证涉及的所有数在C++的int内。

2015.4.2新加四组数据

Source

本题题意就是对每个等差数列维护一棵线段树。
网上的大佬说差值为5好，那就5吧。。。。。
第一次把线段树写入struct中。。。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define ls rt<<1
#define rs rt<<1|1
const int INF = 0x3f3f3f3f,N=3e5,M=5;
using namespace std;
int src[N],n,m;
struct Segment_Tree{
    int sum[N],st[N],en[N],c[N],t[N];
    void up(int rt){
        sum[rt]=max(sum[ls],sum[rs]);
    }
    void build(int rt=1,int l=1,int r=n){
        if(l==r){
            sum[rt]=t[l],c[l]=rt;return ;
        }
        int mid=l+r>>1;
        build(ls,l,mid),build(rs,mid+1,r),up(rt);
    }
    void init(int d){
        int p=0;
        for(register int i=1;i<=d;i++){
            for(register int k=i;k<=n;k+=d)
                st[k]=++p,t[p]=src[k];
            for(int k=i;k<=n;k+=d)en[k]=p;
        }
        build();
    }
    void add(int x,int w){
        sum[x=c[st[x]]]+=w;
        for(x>>=1;x;x>>=1)up(x);
    }
    int query(int l,int r,int rt=1,int L=1,int R=n){
        if(L==l&&r==R)return sum[rt];
        register int mid=(L+R)>>1;
        if(r<=mid)return query(l,r,ls,L,mid);
        else if(l>mid)return query(l,r,rs,mid+1,R);
        else return max(query(l,mid,ls,L,mid),query(mid+1,r,rs,mid+1,R));
    }
    int ask(int x){return query(st[x],en[x]);}
}a[M+1];
template<class T>inline void read(T &res){
    static char ch;T flag=1;
    while((ch=getchar())<'0'||ch>'9')if(ch=='-')flag=-1;res=ch-48;
    while((ch=getchar())>='0'&&ch<='9')res=res*10+ch-48;res*=flag;
}
int main(){
    read(n);
    for(register int i=1;i<=n;i++)read(src[i]);
    for(register int i=1;i<=M;i++)a[i].init(i);
    for(read(m);m--;){
        int x,y,z;
        read(x),read(y),read(z);
        if(x==0){
            src[y]+=z;
            for(register int i=1;i<=M;i++)a[i].add(y,z);
        }else{
            if(z==0)printf("%d\n",src[y]);
            else if(z>M){
                int ans=-INF;
                while(y<=n)ans=max(ans,src[y]),y+=z;
                printf("%d\n",ans);
            }else printf("%d\n",a[z].ask(y));
        }
    }
    return 0;
}