HDU1166(分块)

 敌兵布阵

Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u

Submit Status

Description

C国的死对头A国这段时间正在进行军事演习，所以C国间谍头子Derek和他手下Tidy又开始忙乎了。A国在海岸线沿直线布置了N个工兵营地,Derek和Tidy的任务就是要监视这些工兵营地的活动情况。由于采取了某种先进的监测手段，所以每个工兵营地的人数C国都掌握的一清二楚,每个工兵营地的人数都有可能发生变动，可能增加或减少若干人手,但这些都逃不过C国的监视。 
中央情报局要研究敌人究竟演习什么战术,所以Tidy要随时向Derek汇报某一段连续的工兵营地一共有多少人,例如Derek问:“Tidy,马上汇报第3个营地到第10个营地共有多少人!”Tidy就要马上开始计算这一段的总人数并汇报。但敌兵营地的人数经常变动，而Derek每次询问的段都不一样，所以Tidy不得不每次都一个一个营地的去数，很快就精疲力尽了，Derek对Tidy的计算速度越来越不满:"你个死肥仔，算得这么慢，我炒你鱿鱼!”Tidy想：“你自己来算算看，这可真是一项累人的工作!我恨不得你炒我鱿鱼呢!”无奈之下，Tidy只好打电话向计算机专家Windbreaker求救,Windbreaker说：“死肥仔，叫你平时做多点acm题和看多点算法书，现在尝到苦果了吧!”Tidy说："我知错了。。。"但Windbreaker已经挂掉电话了。Tidy很苦恼，这么算他真的会崩溃的，聪明的读者，你能写个程序帮他完成这项工作吗？不过如果你的程序效率不够高的话，Tidy还是会受到Derek的责骂的. 

Input

第一行一个整数T，表示有T组数据。 
每组数据第一行一个正整数N（N<=50000）,表示敌人有N个工兵营地，接下来有N个正整数,第i个正整数ai代表第i个工兵营地里开始时有ai个人（1<=ai<=50）。 
接下来每行有一条命令，命令有4种形式： 
(1) Add i j,i和j为正整数,表示第i个营地增加j个人（j不超过30） 
(2)Sub i j ,i和j为正整数,表示第i个营地减少j个人（j不超过30）; 
(3)Query i j ,i和j为正整数,i<=j，表示询问第i到第j个营地的总人数; 
(4)End 表示结束，这条命令在每组数据最后出现; 
每组数据最多有40000条命令 

Output

对第i组数据,首先输出“Case i:”和回车, 
对于每个Query询问，输出一个整数并回车,表示询问的段中的总人数,这个数保持在int以内。 

Sample Input

1
10
1 2 3 4 5 6 7 8 9 10
Query 1 3
Add 3 6
Query 2 7
Sub 10 2
Add 6 3
Query 3 10
End 

Sample Output

Case 1:
6
33
59

这题在之前是用线段树做的，今天艾神讲分块，把艾神的分块做法贴一下

//2016.8.12
#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#define kuai 230
using namespace std;
int arr[50005],n;
int sum[305];
void add(int x,int y)
{
    int belong=x/kuai;
    sum[belong]+=y;
    arr[x]+=y;
    return;
}
int query(int x,int y)
{
    int l=0,r=kuai-1,ans=0;
    for (int i=0;i<kuai;i++,l+=kuai,r+=kuai)
    {
        int L,R;
        L=max(l,x);
        R=min(r,y);
        if (L>R) continue;
        if (L==l && R==r)
            ans+=sum[i];
        else
        {
            for (int j=L;j<=R;j++)
                ans+=arr[j];
        }
    }
    return ans;
}
int main()
{
    int t;
    cin>>t;
    for (int cas=1;cas<=t;cas++)
    {
        cin>>n;
        for (int i=1;i<=n;i++)
            scanf("%d",&arr[i]);
        memset(sum,0,sizeof(sum));
        for (int i=1;i<=n;i++)
        {
            int belong=i/230;
            sum[belong]+=arr[i];
        }
        char ch[10];
        printf("Case %d:
",cas);
        scanf("%s",ch);
        while (ch[0]!='E')
        {
            int x,y;
            scanf("%d%d",&x,&y);
            if (ch[0]=='A')
                add(x,y);
            else
            if (ch[0]=='S')
                add(x,-y);
            else
            printf("%d
",query(x,y));
            scanf("%s",ch);
        }
    }
}

线段树做法。

//2016-05-13
#include<iostream>
#include<string>
#include<cstdio>  
#include<cstring>  
#include<algorithm>
using namespace std;

struct Node{
    int l, r, sum;
}A[500050];

void Pushup(int root)
{
    A[root].sum = A[root<<1].sum + A[root<<1|1].sum;
    return ;
}

void Build(int l, int r, int root)
{
    A[root].l = l;
    A[root].r = r;
    if(l == r)
    {
        scanf("%d", &A[root].sum);
        return;
    }
    int mid = (l + r)>>1;
    Build(l, mid, root<<1);
    Build(mid+1, r, root<<1|1);
    Pushup(root);
    return ;
}

void update(int l, int r, int root , int k)
{
    if(l == A[root].l && r == A[root].r)
    {
        A[root].sum += k;
        return ;
    }
    int mid = (A[root].l + A[root].r) >> 1;
    if(r <= mid)
        update(l, r, root<<1, k);
    else if(l > mid)
        update(l, r, root<<1|1, k);
    else
    {
        update(l, mid, root<<1, k);
        update(mid+1, r, root<<1|1, k);
    }
    Pushup(root);
    return ;    
}

int Query(int l, int r, int root)
{
    if(l == A[root].l && r == A[root].r)
        return A[root].sum;
    int mid = (A[root].l + A[root].r) >> 1;
    int ans = 0;
    if(r <= mid)
        ans += Query(l, r, root<<1);
    else if(l > mid)
        ans += Query(l, r, root<<1|1);
    else
        ans += Query(l, mid, root<<1) + Query(mid+1, r, root<<1|1);
    return ans;
}

int main ()
{
    int T, cas, N;
    int l, k;
    string OP, A = "Add", S = "Sub", Q = "Query", E = "End";
    cas = 0;
    scanf("%d", &T);
    while(T--)
    {
        scanf("%d", &N);
        Build(1, N, 1);
        printf("Case %d:
", ++cas);
        while(cin>>OP)
        {
            if(OP == A)
            {
                scanf("%d %d", &l, &k);
                update(l, l, 1, k);
            }
            else if(OP == S)
            {
                scanf("%d %d", &l, &k);
                update(l, l, 1, -k);
            }
            else if(OP == Q)
            {
                scanf("%d %d", &l, &k);
                printf("%d
", Query(l, k, 1));
            }
            else if(OP == E)
                break;
        }
    }
    return 0;
}

树状数组做法。

//2017-05-17
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

const int N = 50005;
int e[N], n;

int lowbit(int x)
{
    return x&-x;
}

void add(int pos, int w)
{
    while(pos <= n){
        e[pos] += w;
        pos += lowbit(pos);
    }
}

int sum(int pos)
{
    int ans = 0;
    while(pos > 0){
        ans += e[pos];
        pos -= lowbit(pos);
    }
    return ans;
}

int main()
{
    int T, a;
    scanf("%d", &T);
    for(int kase = 1; kase <= T; kase++){
        scanf("%d", &n);
        memset(e, 0, sizeof(e));
        for(int i = 1; i <= n; i++)
        {
            scanf("%d", &a);
            add(i, a);
        }          
        char op[10];
        int x, y;
        printf("Case %d:
", kase);
        while(scanf("%s", op)){
            if(op[0] == 'Q'){
                scanf("%d%d", &x, &y);
                printf("%d
", sum(y)-sum(x-1));
            }else if(op[0] == 'A'){
                scanf("%d%d", &x, &y);
                add(x, y);
            }else if(op[0] == 'S'){
                scanf("%d%d", &x, &y);
                add(x, -y);
            }else if(op[0] == 'E')
                  break;
        }
    }

    return 0;
}