12年杭州 D 二维树状数组+离散化 HDU 4456

12年杭州 D

赛中的时候感觉数组开不下，并没有什么很好的思路
实际上，离散化每个需要操作的点之后，这就变成了一个简单的二维树状数组题目。
关于离散化的操作，由于服务器内存卡的比较死，试了几种姿势后都不行，还是采用网上的线性膜的方法。按理说这种方法可以构造出极端数据卡掉的，但是过了

失分点：
1. 没有想到二维线树状数组离散化
2. 离散化的离散方法，离散化的数组大小。
学习点：
1. 离散化的方式总结以下几种：
1）map
2）set+lower_bound
3) 将所有查询存入数组后，用unique函数去重再进行二分查找
4）康拓展开（只针对特定题目）
5）线性膜
2. 二维线段树的写法。

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <string>
#include <algorithm>
#include <iostream>
#include <vector>
using namespace std;
#define mp make_pair
#define fi first
#define se second
#define pb push_back
//const int MAXN = 80000 + 5;
const int MAXN = 2222222;
const int mod = 2222222;
int n, m, up;
int sum[MAXN], id[MAXN];
int cnt;
int lowbit(int i){return i & (-i);}
int Hash(int u)
{
    int temp = u % mod;
    while(1) {
        if(id[temp] == -1) {
            id[temp] = u;
            return temp;
        }
        else if(id[temp] > 0 && id[temp] == u) return temp;
        else if(id[temp] > 0 && id[temp] != u) temp++;
        if(temp >= mod) temp = 0;
    }
}
int getHash(int u)
{
    int temp = u % mod;
    while(1) {
        if(id[temp] == -1) return -1;
        else if(id[temp] > 0 && id[temp] == u) return temp;
        temp++;
        if(temp >= mod) temp = 0;
    }
}
void addid(int x, int y, int w)
{
    for(int i = x ; i <= up ; i += lowbit(i)) {
        for(int j = y ; j <= up ; j += lowbit(j)) {
            sum[Hash(i * up + j)] += w;
        }
    }
}
void add(int x, int y, int w)
{
    for(int i = x ; i <= up ; i += lowbit(i)) {
        for(int j = y ; j <= up ; j += lowbit(j)) {
            sum[Hash(i * up + j)] += w;
        }
    }
}
int query(int x, int y)
{
    x = min(x, up); y = min(y, up);
    x = max(x, 0); y = max(y, 0);
    int res = 0;
    for(int i = x ; i ; i -= lowbit(i)) {
        for(int j = y ; j ; j -= lowbit(j)) {
            int mark = getHash(i * up + j);
            if(mark == -1) continue;
            res += sum[mark];
        }
    }
    return res;
}
void trans(int &x, int &y)
{
    int tx = y - x + n;
    int ty = y + x;
    x = tx, y = ty;
}
int main()
{
    while(scanf("%d", &n) != EOF && n) {
        scanf("%d", &m);
        up = n << 1;
        memset(sum, 0, sizeof sum);
        memset(id, -1, sizeof id);
        cnt = 0;
        for(int i = 0 ; i < m ; i++) {
            int op, x, y, z;
            scanf("%d%d%d%d", &op, &x, &y, &z);
            trans(x, y);
            if(op == 1) add(x, y, z);
            else {
                int t1 = query(x+z,y+z) + query(x-z-1,y-z-1);
                int t2 = query(x-z-1,y+z) + query(x+z,y-z-1);
                printf("%d\n", t1 - t2);
            }
        }

    }
    return 0;
}