HDU-4027: Can you answer these queries?（区间）

线段树

Description

A lot of battleships of evil are arranged in a line before the battle. Our commander decides to use our secret weapon to eliminate the battleships. Each of the battleships can be marked a value of endurance. For every attack of our secret weapon, it could decrease the endurance of a consecutive part of battleships by make their endurance to the square root of it original value of endurance. During the series of attack of our secret weapon, the commander wants to evaluate the effect of the weapon, so he asks you for help.
You are asked to answer the queries that the sum of the endurance of a consecutive part of the battleship line.

Notice that the square root operation should be rounded down to integer.

Input

The input contains several test cases, terminated by EOF.
For each test case, the first line contains a single integer N, denoting there are N battleships of evil in a line. (1 <= N <= 100000)
The second line contains N integers Ei, indicating the endurance value of each battleship from the beginning of the line to the end. You can assume that the sum of all endurance value is less than $2^{63}$.
The next line contains an integer M, denoting the number of actions and queries. (1 <= M <= 100000)
For the following M lines, each line contains three integers T, X and Y. The T=0 denoting the action of the secret weapon, which will decrease the endurance value of the battleships between the X-th and Y-th battleship, inclusive. The T=1 denoting the query of the commander which ask for the sum of the endurance value of the battleship between X-th and Y-th, inclusive.

Output

For each test case, print the case number at the first line. Then print one line for each query. And remember follow a blank line after each test case.

Sample Input

10
1 2 3 4 5 6 7 8 9 10
5
0 1 10
1 1 10
1 1 5
0 5 8
1 4 8 

Sample Output

Case #1:
19
7
6 

#include<stdio.h>
#include<algorithm>
#include<math.h>
#define LL long long
using namespace std;
LL a[100005];
struct node
{
    int l,r,lazy;
    LL sum;
}s[100005*4];
void build(int id,int l,int r)
{
    s[id].l=l;
    s[id].r=r;
    if(l==r)
    {
        s[id].sum=a[l];
        if(s[id].sum==1)
            s[id].lazy=0;
        else
            s[id].lazy=1;
    }
    else
    {
        int mid=(l+r)/2;
        build(id*2,l,mid);
        build(id*2+1,mid+1,r);
        s[id].sum=s[id*2].sum+s[id*2+1].sum;
        s[id].lazy=s[id*2].lazy+s[id*2+1].lazy;
    }
}
void rep(int id,int l,int r)
{
    if(s[id].lazy==0) //都是0或者1，不用开方了！
        return;
    if(s[id].r<l||s[id].l>r)
        return;
    if(s[id].l==s[id].r)
    {
        s[id].sum=(LL)sqrt(s[id].sum);
        if(s[id].sum==1)
            s[id].lazy=0;
        return ;
    }
    else
    {
        int mid=(s[id].l+s[id].r)/2;
        if(r<=mid)
            rep(id*2,l,r);
        else if(l>mid)
            rep(id*2+1,l,r);
        else
        {
            rep(id*2,l,r);
            rep(id*2+1,l,r);
        }
    }
    s[id].sum=s[id*2].sum+s[id*2+1].sum;
    s[id].lazy=s[id*2].lazy+s[id*2+1].lazy;
}
LL sum(int id,int l,int r)
{
    if(s[id].r<l||s[id].l>r)
        return 0;
    if(l<=s[id].l&&s[id].r<=r)
        return s[id].sum;
    int mid=(s[id].l+s[id].r)/2;
    if(r<=mid)
        return sum(id*2,l,r);
    else if(l>mid)
        return sum(id*2+1,l,r);
    else
        return sum(id*2,l,r)+sum(id*2+1,l,r);
}
int main()
{
    int n,i,m,t,x,y,tt=0;
    while(scanf("%d",&n)!=EOF)
    {
        tt++;
        for(i=1;i<=n;i++)
            scanf("%lld",&a[i]);
        build(1,1,n);
        scanf("%d",&m);
        printf("Case #%d:\n",tt);
        while(m--)
        {
            scanf("%d%d%d",&t,&x,&y);
            if(x>y) swap(x,y);
            if(t==0)
                rep(1,x,y);
            if(t==1)
                printf("%lld\n",sum(1,x,y));
        }
        printf("\n");
    }
    return 0;
}