1. 程式人生 > >Matrix 【POJ - 2155】【二維線段樹+永久化標記】

Matrix 【POJ - 2155】【二維線段樹+永久化標記】

題目連結


  挺好的一道題,一開始用lazy標記往下推,總是推不出樣例的正解,然後就去看了相關部落格,發現卻確實如此,在這裡是無法用lazy標記來層層推的,並且還會出現超記憶體的情況,所以,便改用了永久化標記來解這道題。

  還有一件是,關於discuss裡的討論區有一塊這樣的討論:

1
2 10
C 2 1 2 2
Q 2 3
C 2 1 2 1
Q 1 1
C 1 1 2 1
C 1 2 1 2
C 1 1 2 2
Q 1 1
C 1 1 2 1
Q 2 1

相信很多人是過不了上面的資料的吧!——但是你們就沒發現數據本身就錯了嗎。。。它開的N*N的空間不對啊,怎能詢問空間外的點呢?應改為:

正確樣例

1
3 10
C 2 1 2 2
Q 2 3
C 2 1 2 1
Q 1 1
C 1 1 2 1
C 1 2 1 2
C 1 1 2 2
Q 1 1
C 1 1 2 1
Q 2 1

那麼,我相信大家就都能過了吧。


#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#define lowbit(x) ( x&(-x) )
#define efs 1e-6
#define pi 3.141592653589793
#define e 2.718281828459045
using namespace std;
typedef unsigned long long ull;
typedef long long ll;
const int maxN = 1005;
int N, Q;
short int tree[maxN<<2][maxN<<2];
void pushup(int rt, int X0)
{
    tree[X0][rt] = tree[X0][rt<<1] + tree[X0][rt<<1|1];
}
void update_In(int rt, int X0, int l, int r, int ql, int qr)
{
    if(ql == l && qr == r)
    {
        tree[X0][rt]^=1;
        return;
    }
    int mid = (l + r)>>1;
    if(ql>mid) update_In(rt<<1|1, X0, mid+1, r, ql, qr);
    else if(qr<=mid) update_In(rt<<1, X0, l, mid, ql, qr);
    else
    {
        update_In(rt<<1, X0, l, mid, ql, mid);
        update_In(rt<<1|1, X0, mid+1, r, mid+1, qr);
    }
}
void update_Out(int rt, int l, int r, int qlx, int qly, int qrx, int qry)
{
    if(qlx == l && qrx == r)
    {
        update_In(1, rt, 1, N, qly, qry);
        return;
    }
    int mid = (l + r)>>1;
    if(qlx>mid) update_Out(rt<<1|1, mid+1, r, qlx, qly, qrx, qry);
    else if(qrx<=mid) update_Out(rt<<1, l, mid, qlx, qly, qrx, qry);
    else
    {
        update_Out(rt<<1, l, mid, qlx, qly, mid, qry);
        update_Out(rt<<1|1, mid+1, r, mid+1, qly, qrx, qry);
    }
}
int query_In(int rt, int X0, int l, int r, int qy, int sum)
{
    sum^=tree[X0][rt];
    if(l == r) return sum;
    int mid = (l + r)>>1;
    if(qy>mid) return query_In(rt<<1|1, X0, mid+1, r, qy, sum);
    else return query_In(rt<<1, X0, l, mid, qy, sum);
}
int query_Out(int rt, int l, int r, int qx, int qy, int sum)
{
    sum^=query_In(1, rt, 1, N, qy, 0);
    if(l == r) return sum;
    int mid = (l + r)>>1;
    if(qx>mid) return query_Out(rt<<1|1, mid+1, r, qx, qy, sum);
    else return query_Out(rt<<1, l, mid, qx, qy, sum);
}
int main()
{
    int T;  scanf("%d", &T);
    int Cas = 0;
    while(T--)
    {
        if(Cas++ > 0) printf("\n");
        scanf("%d%d", &N, &Q);
        memset(tree, 0, sizeof(tree));
        while(Q--)
        {
            char s[3];
            scanf("%s", s);
            if(s[0] == 'C')
            {
                int lx, ly, rx, ry;
                scanf("%d%d%d%d", &lx, &ly, &rx, &ry);
                update_Out(1, 1, N, lx, ly, rx, ry);
            }
            else
            {
                int x, y;
                scanf("%d%d", &x, &y);
                printf("%d\n", query_Out(1, 1, N, x, y, 0));
            }
        }
    }
    return 0;
}