whu 1464 Deal whth Numbers 线段树维护区间除法

题意

　　N个数，M个操作，操作有三类：

　　1.  D,A,B,C  区间[a,b]中 大于0 的数 除以C，计算机整除

　　2.  M,A,B,C  区间[a,b]所有数减去C

　　3.  S,A,B     区间[a,b]求和

解法

　　区间更新和求和，线段树上延迟更新就可以了。主要是除法如何维护。

　   这里的数C <= 5*1E5 , 因为大于0才除，所以每个数最多被除 log(Ai) 次， N个数则最多被除 N*log(Ai) 次了。

　   所以做除法可以直接更新到底，还有就是使用一个标记，区间中都不可被除时，就没必要更新到底了。

　　另外，除数是1的时候，不能除，貌似后台数据好多1，╮(╯▽╰)╭~~~

　　

View Code

 1 #include<cstdio>

 2 #include<cstdlib>

 3 #include<cstdlib>

 4 #include<algorithm>

 5 using namespace std;

 6 

 7 typedef long long LL;

 8 

 9 const int N = 5e5+10;

10 //segment_Tree

11 #define lch rt<<1,l,m

12 #define rch rt<<1|1,m+1,r

13 LL sum[N<<2], add[N<<2];

14 bool flag[N<<2];

15 

16 void push_up(int rt){

17     sum[rt] = sum[rt<<1]+sum[rt<<1|1];

18     flag[rt] = flag[rt<<1] | flag[rt<<1|1];    

19 }

20 void push_down(int rt,int L){

21     if( add[rt] ){

22         add[rt<<1] += add[rt];

23         add[rt<<1|1] += add[rt];

24         sum[rt<<1] += (L-(L>>1))*add[rt];

25         sum[rt<<1|1] +=    (L>>1)*add[rt];

26         add[rt] = 0;

27     }    

28 }

29 void build(int rt,int l,int r){

30     sum[rt] = add[rt] = 0;

31     if(l == r){ scanf("%lld",&sum[rt]); flag[rt] = (sum[rt]>0); return; }

32     int m = (l+r)>>1;

33     build(lch), build(rch);

34     push_up(rt);    

35 }

36 void Minus(int rt,int l,int r,int a,int b,int c){

37     if(a <= l && r <= b){

38         add[rt] += c; sum[rt] += (r-l+1)*c; return;    

39     }

40     push_down(rt,r-l+1);

41     int m = (l+r)>>1;

42     if( a <= m ) Minus( lch,a,b,c );

43     if( m <  b ) Minus( rch,a,b,c );

44     push_up(rt);

45 }

46 void Division(int rt,int l,int r,int a,int b,int c){

47     if( flag[rt] == 0 ) return;

48     if(l == r){

49         if(sum[rt]>0) sum[rt] /= (LL)c; 

50         flag[rt] = (sum[rt] > 0);     return;

51     }    

52     push_down(rt, r-l+1);

53     int m = (l+r)>>1;

54     if( a <= m ) Division( lch,a,b,c );

55     if( m <  b ) Division( rch,a,b,c );

56     push_up(rt);

57 }

58 LL Sum(int rt,int l,int r,int a,int b){

59     if(a <= l && r <= b) return sum[rt];

60     push_down(rt,r-l+1);

61     int m = (l+r)>>1;

62     LL res = 0;

63     if( a <= m ) res += Sum(lch,a,b);

64     if( m <  b ) res += Sum(rch,a,b);

65     return res;    

66 }

67 

68 int main(){

69 //    freopen("D.in","r", stdin);

70 //    freopen("test.out","w",stdout);

71     int T;

72     scanf("%d", &T);

73     for(int Case = 1; Case <= T; Case++){

74         int n, m, a, b, c;

75         char op[10];

76         scanf("%d%d", &n,&m);

77         printf("Case %d:\n",Case);

78         build( 1, 1, n );

79         for(int i = 0; i < m; i++){

80             scanf("%s", op);

81             if( op[0] == 'D' ){

82                 scanf("%d%d%d",&a,&b,&c);

83                 if( c == 1 ) continue;

84                 Division(1,1,n,a,b,c);     

85             }

86             else if( op[0] == 'M' ){

87                 scanf("%d%d%d",&a,&b,&c);

88                 Minus(1,1,n,a,b,-c);    

89             }

90             else{

91                 scanf("%d%d",&a,&b);

92                 LL res = Sum(1,1,n,a,b);

93                 printf("%lld\n",res);

94             }

95         }    

96         puts("");

97     }

98     return 0;    

99 }