hdu1166 敵兵布陣（CDQ分治）

Problem 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
        

Statistic |  Submit |  Discuss | Note

分析：

這道不算是很典型的三維偏序，實際上是一道二維偏序：

• 時間（一開始就有序了）
• x坐标

這樣的話我們就不用樹狀數組了

struct node{
    int x,y,type,id;
    //x坐标 y值 type操作類型 id詢問編号
};
           

在CDQ分治中，我們隻需要一個sum就可以記錄字首和了

tip

數組要開夠，不然會T的呢

能用歸并的就用歸并，畢竟STL中的sort複雜度是O(nlogn)

//這裡寫代碼片
#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

const int N=;
struct node{
    int x,y,type,id;
};
node po[N<<],q[N<<2];
int n,ans[N],tot,totx;

void CDQ(int L,int R)
{
    if (L==R) return;
    int M=(L+R)>>;
    int sum=;
    CDQ(L,M); CDQ(M+,R);

    int t1=L,t2=M+;
    for (int i=L;i<=R;i++)
    {
        if ((t1<=M&&po[t1].x<=po[t2].x)||t2>R){
            sum+=po[t1].y;
            q[i]=po[t1++];
        }
        else 
        {
            if (po[t2].type==) ans[po[t2].id]-=sum;
            else if(po[t2].type==) ans[po[t2].id]+=sum;
            q[i]=po[t2++];
        }
    }
    for (int i=L;i<=R;i++) po[i]=q[i];
}

int main()
{
    int T;
    scanf("%d",&T);
    for (int cas=;cas<=T;cas++)
    {
        scanf("%d",&n);
        memset(po,,sizeof(po));
        memset(ans,,sizeof(ans)); tot=; totx=;

        for (int i=;i<=n;i++)
        {
            scanf("%d",&po[i].y);
            po[i].x=i; po[i].type=;
        }
        tot=n;
        char s[];
        int x,y;
        while (scanf("%s",s)!=EOF&&s[]!='E')
        {
            scanf("%d%d",&x,&y);
            if (s[]=='A')
            {
                tot++;
                po[tot].x=x; po[tot].y=y; po[tot].type=;
            }
            else if (s[]=='S')
            {
                tot++;
                po[tot].x=x; po[tot].y=-y;  po[tot].type=;
            }
            else
            {
                tot++; totx++;
                po[tot].x=x-; po[tot].type=; po[tot].id=totx;
                tot++;
                po[tot].x=y; po[tot].type=; po[tot].id=totx;
            }
        }

        CDQ(,tot);   //num(時間)已經有序 

        printf("Case %d:\n",cas);
        for (int i=;i<=totx;i++) printf("%d\n",ans[i]);
    }
    return ;
}