BZOJ 3922: Karin的弹幕 线段树 暴力乱搞

3922: Karin的弹幕

Time Limit: 10 Sec  Memory Limit: 512 MB
Submit: 606  Solved: 230
[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新加四组数据

新加数据之后数据范围应该是300000了啊

害我本该AC的代码WA两发 泪

推荐blog

发现我多年未碰数学

连基本不等式都不会了

觉得线段树的维护还是蛮新颖的

#include<cmath>#include<ctime>#include<cstdio>#include<cstring>#include<cstdlib>#include<complex>#include<iostream>#include<algorithm>#include<iomanip>#include<vector>#include<string>#include<bitset>#include<queue>#include<map>#include<set>using namespace std;inline int read(){int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch<='9'&&ch>='0'){x=(x<<1)+(x<<3)+ch-'0';ch=getchar();}return x*f;}void print(int x){if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}const int N=300010,M=11,inf=2147483647;int n,a[N];struct seg_tree{int mx;}str[M][N];int pos[M][N],val[M][N],ed[M][N];inline void pushup(int x,int k){str[x][k].mx=max(str[x][k<<1].mx,str[x][k<<1|1].mx);}void build(int x,int k,int l,int r){if(l==r){str[x][k].mx=a[val[x][l]];return ;}int mid=(l+r)>>1;build(x,k<<1,l,mid);build(x,k<<1|1,mid+1,r);pushup(x,k);}void modify(int x,int k,int l,int r,int aim){if(l==r){str[x][k].mx=a[val[x][l]];return ;}int mid=(l+r)>>1;aim<=mid?modify(x,k<<1,l,mid,aim):modify(x,k<<1|1,mid+1,r,aim);pushup(x,k);}int query(int x,int k,int l,int r,int L,int R){if(l>=L&&r<=R){return str[x][k].mx;}int mid=(l+r)>>1;if(mid<L)return query(x,k<<1|1,mid+1,r,L,R);else if(mid>=R)return query(x,k<<1,l,mid,L,R);else return max(query(x,k<<1|1,mid+1,r,L,R),query(x,k<<1,l,mid,L,R));}void initial(int x){register int i,j,k=0;for(i=1;i<=x;++i){for(j=i;j<=n;j+=x){pos[x][j]=++k;val[x][k]=j;}for(j=i;j<=n;j+=x)ed[x][j]=k;}build(x,1,1,n);}int main(){n=read();register int i,opt,j;for(i=1;i<=n;++i)a[i]=read();int m=min(n/int(log2(n)),10);for(i=1;i<=m;++i)initial(i);int Q=read();while(Q--){opt=read();i=read();j=read();if(opt&&j<=m){print(query(j,1,1,n,pos[j][i],ed[j][i]));puts("");}else if(opt&&j>m){int ans=-inf;for(;i<=n;i+=j)ans=max(ans,a[i]);print(ans);puts("");}else if(!opt){a[i]+=j;for(opt=1;opt<=m;++opt)modify(opt,1,1,n,pos[opt][i]);}}return 0;}/*101 6 1 4 9 4 8 2 8 5101 3 30 5 40 3 81 2 51 4 81 7 51 3 60 1 21 5 31 4 98848913410-9 -6 2 -10 -2 -6 10 6 -4 -2101 2 31 6 30 7 80 4 -60 10 -51 10 40 3 -81 2 40 10 -51 1 26-4-7-618*/