题目大意
给定一个长度为
n
的序列和
需要维护区间加,区间除法(下取整)以及查询区间和和区间最小值。
Data Constraint
n,m≤105
题解
如果一个区间的最大值和最小值做区间除后的增量相同,就直接打一个区间减法标记。
卡卡常。
时间复杂度: O(nlog2n)
SRC
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std ;
#define N 100000 + 10
#define Min( a , b ) ( a < b ? a : b )
#define Max( a , b ) ( a > b ? a : b )
typedef long long ll ;
struct Tree {
ll del ;
ll Minv , Maxv , Sum ;
int len ;
} T[4*N] ;
int a[N] ;
int n , m ;
ll ret ;
inline int Read() {
int ret = 0 , flag = 1 ;
char ch = getchar() ;
while ( ch < '0' || ch > '9' ) {
if ( ch == '-' ) flag = -1 ;
ch = getchar() ;
}
while ( ch >= '0' && ch <= '9' ) {
ret = ret * 10 + ch - '0' ;
ch = getchar() ;
}
return ret * flag ;
}
inline void Build( int v , int l , int r ) {
if ( l == r ) {
T[v].Minv = T[v].Maxv = T[v].Sum = a[l] ;
T[v].len = 1 ;
return ;
}
int ls = v << 1 , rs = (v << 1) + 1 ;
T[v].len = r - l + 1 ;
int mid = (l + r) >> 1 ;
Build( ls , l , mid ) ;
Build( rs , mid + 1 , r ) ;
T[v].Minv = Min( T[ls].Minv , T[rs].Minv ) ;
T[v].Maxv = Max( T[ls].Maxv , T[rs].Maxv ) ;
T[v].Sum = T[ls].Sum + T[rs].Sum ;
}
inline void ADD( int v , int l , int r , int x , int y , int val ) {
if ( l == x && r == y ) {
T[v].del += val ;
T[v].Minv += val ;
T[v].Maxv += val ;
T[v].Sum += T[v].len * val ;
return ;
}
int ls = v << 1 , rs = (v << 1) + 1 ;
if ( T[v].del ) {
T[ls].Minv += T[v].del , T[ls].Maxv += T[v].del ;
T[rs].Minv += T[v].del , T[rs].Maxv += T[v].del ;
T[ls].Sum += T[v].del * T[ls].len ;
T[rs].Sum += T[v].del * T[rs].len ;
T[ls].del += T[v].del , T[rs].del += T[v].del ;
T[v].del = 0 ;
}
int mid = (l + r) >> 1 ;
if ( y <= mid ) ADD( ls , l , mid , x , y , val ) ;
else if ( x > mid ) ADD( rs , mid + 1 , r , x , y , val ) ;
else {
ADD( ls , l , mid , x , mid , val ) ;
ADD( rs , mid + 1 , r , mid + 1 , y , val ) ;
}
T[v].Minv = Min( T[ls].Minv , T[rs].Minv ) ;
T[v].Maxv = Max( T[ls].Maxv , T[rs].Maxv ) ;
T[v].Sum = T[ls].Sum + T[rs].Sum ;
}
inline ll Calc( ll x , ll d ) {
if ( x * d > 0 ) return x / d ;
return x / d - ( x % d != 0 ) ;
}
inline void DIV( int v , int l , int r , int x , int y , int div ) {
int mid = (l + r) >> 1 ;
int ls = v << 1 , rs = (v << 1) + 1 ;
if ( l == x && r == y ) {
if ( T[v].Minv - Calc( T[v].Minv , div ) == T[v].Maxv - Calc( T[v].Maxv , div ) ) {
ll del = T[v].Minv - Calc( T[v].Minv , div ) ;
T[v].Minv -= del ;
T[v].Maxv -= del ;
T[v].Sum -= del * T[v].len ;
T[v].del -= del ;
} else {
if ( T[v].del ) {
T[ls].Minv += T[v].del , T[ls].Maxv += T[v].del ;
T[rs].Minv += T[v].del , T[rs].Maxv += T[v].del ;
T[ls].Sum += T[v].del * T[ls].len ;
T[rs].Sum += T[v].del * T[rs].len ;
T[ls].del += T[v].del , T[rs].del += T[v].del ;
T[v].del = 0 ;
}
DIV( ls , l , mid , x , mid , div ) ;
DIV( rs , mid + 1 , r , mid + 1 , y , div ) ;
T[v].Minv = Min( T[ls].Minv , T[rs].Minv ) ;
T[v].Maxv = Max( T[ls].Maxv , T[rs].Maxv ) ;
T[v].Sum = T[ls].Sum + T[rs].Sum ;
}
return ;
}
if ( T[v].del ) {
T[ls].Minv += T[v].del , T[ls].Maxv += T[v].del ;
T[rs].Minv += T[v].del , T[rs].Maxv += T[v].del ;
T[ls].Sum += T[v].del * T[ls].len ;
T[rs].Sum += T[v].del * T[rs].len ;
T[ls].del += T[v].del , T[rs].del += T[v].del ;
T[v].del = 0 ;
}
if ( y <= mid ) DIV( ls , l , mid , x , y , div ) ;
else if ( x > mid ) DIV( rs , mid + 1 , r , x , y , div ) ;
else {
DIV( ls , l , mid , x , mid , div ) ;
DIV( rs , mid + 1 , r , mid + 1 , y , div ) ;
}
T[v].Minv = Min( T[ls].Minv , T[rs].Minv ) ;
T[v].Maxv = Max( T[ls].Maxv , T[rs].Maxv ) ;
T[v].Sum = T[ls].Sum + T[rs].Sum ;
}
inline void FindMin( int v , int l , int r , int x , int y ) {
if ( l == x && r == y ) {
ret = Min( ret , T[v].Minv ) ;
return ;
}
int ls = v << 1 , rs = (v << 1) + 1 ;
if ( T[v].del ) {
T[ls].Minv += T[v].del , T[ls].Maxv += T[v].del ;
T[rs].Minv += T[v].del , T[rs].Maxv += T[v].del ;
T[ls].Sum += T[v].del * T[ls].len ;
T[rs].Sum += T[v].del * T[rs].len ;
T[ls].del += T[v].del , T[rs].del += T[v].del ;
T[v].del = 0 ;
}
int mid = (l + r) >> 1 ;
if ( y <= mid ) FindMin( v + v , l , mid , x , y ) ;
else if ( x > mid ) FindMin( v + v + 1 , mid + 1 , r , x , y ) ;
else {
FindMin( v + v , l , mid , x , mid ) ;
FindMin( v + v + 1 , mid + 1 , r , mid + 1 , y ) ;
}
T[v].Minv = Min( T[ls].Minv , T[rs].Minv ) ;
T[v].Maxv = Max( T[ls].Maxv , T[rs].Maxv ) ;
T[v].Sum = T[ls].Sum + T[rs].Sum ;
}
inline void FindSum( int v , int l , int r , int x , int y ) {
if ( l == x && r == y ) {
ret += T[v].Sum ;
return ;
}
int ls = v << 1 , rs = (v << 1) + 1 ;
if ( T[v].del ) {
T[ls].Minv += T[v].del , T[ls].Maxv += T[v].del ;
T[rs].Minv += T[v].del , T[rs].Maxv += T[v].del ;
T[ls].Sum += T[v].del * T[ls].len ;
T[rs].Sum += T[v].del * T[rs].len ;
T[ls].del += T[v].del , T[rs].del += T[v].del ;
T[v].del = 0 ;
}
int mid = (l + r) >> 1 ;
if ( y <= mid ) FindSum( v + v , l , mid , x , y ) ;
else if ( x > mid ) FindSum( v + v + 1 , mid + 1 , r , x , y ) ;
else {
FindSum( v + v , l , mid , x , mid ) ;
FindSum( v + v + 1 , mid + 1 , r , mid + 1 , y ) ;
}
T[v].Minv = Min( T[ls].Minv , T[rs].Minv ) ;
T[v].Maxv = Max( T[ls].Maxv , T[rs].Maxv ) ;
T[v].Sum = T[ls].Sum + T[rs].Sum ;
}
int main() {
n = Read() , m = Read() ;
for (int i = 1 ; i <= n ; i ++ ) a[i] = Read() ;
Build( 1 , 1 , n ) ;
for (int i = 1 ; i <= m ; i ++ ) {
int op = Read() ;
if ( op == 1 ) {
int l = Read() + 1 , r = Read() + 1 , c = Read() ;
ADD( 1 , 1 , n , l , r , c ) ;
} else if ( op == 2 ) {
int l = Read() + 1 , r = Read() + 1 , d = Read() ;
DIV( 1 , 1 , n , l , r , d ) ;
} else if ( op == 3 ) {
int l = Read() + 1 , r = Read() + 1 ;
ret = 1e18 ;
FindMin( 1 , 1 , n , l , r ) ;
printf( "%lld\n" , ret ) ;
} else {
int l = Read() + 1 , r = Read() + 1 ;
ret = 0 ;
FindSum( 1 , 1 , n , l , r ) ;
printf( "%lld\n" , ret ) ;
}
}
return 0 ;
}
以上.