Cube
Problem Description
Given an N*N*N cube A, whose elements are either 0 or 1. A[i, j, k]
means the number in the i-th row , j-th column and k-th layer.
Initially we have A[i, j, k] = 0 (1 <= i, j, k <= N).
We
define two operations, 1: “Not” operation that we change the A[i, j,
k]=!A[i, j, k]. that means we change A[i, j, k] from 0->1,or
1->0. (x1<=i<=x2,y1<=j<=y2,z1<=k<=z2).
0: “Query” operation we want to get the value of A[i, j, k].
Input
Multi-cases.
First line contains N and M, M lines follow indicating the operation below.
Each operation contains an X, the type of operation. 1: “Not” operation and 0: “Query” operation.
If X is 1, following x1, y1, z1, x2, y2, z2.
If X is 0, following x, y, z.
Output
For each query output A[x, y, z] in one line. (1<=n<=100 sum of m <=10000)
Sample Input
2 5
1 1 1 1 1 1 1
0 1 1 1
1 1 1 1 2 2 2
0 1 1 1
0 2 2 2
Sample Output
題意：在一個三維長方體里把0元素改成1或者把1改成0或者詢問某個位置是0還是1.
分析:
非樹狀數組莫屬，其他數據結構都太帥了，用不著。
先說說一維上的樹狀數組，什么問題都得成最簡單的開始，理解后才能推廣到多維或復雜的。
把對某個方體的修改操作看做是對某個方體的操作數+1，那一個位置操作數的奇偶性決定了它的值。
樹狀數組1位置是管1的，2位置是管1到2的,3是管3的，4是管1到4的......所以對一個區間操作，這個區間首先按樹狀數組的分區規則，會被分成若干個不相交的區間，在樹狀數組里對應的就是若干不相交的小樹。
比如對1到3操作，會分成分別對[1,2],[3,3]兩段的操作；對2到4的操作，可以分成[2,2],[3,3],[4,4]，但是對這種情況，按這種方式去修改樹狀數組是麻煩的，可以轉化成分別對[1,4]操作，再對[1,1]操作，這樣[1,1]被修改了兩次，所以相當于沒修改過(配合這題的奇偶性)，也可處理成對[1,4]的操作數+1，再對[1,1]的操作數-1。
知道每段的信息保存操作次數，可以通過疊加就知道每個位置的操作數。
比如對1到3操作了一次，c[2] = 1(對應[1,2]的區間)，c[3] = 1(對應[3,3]區間)；再對1操作，則c[1] = 1（對應[1,1]區間）。則1的操作數是c[1] + c[2]。
到此還不理解的，建議最好看看
樹狀數組的詳細結構，然后畫一畫就應該清楚多了，表達上的不便請原諒，這個算是比較水的，不過對樹狀數組的理解是有點用處的，還有就是在樹狀數組上對方塊的切割，例如對對2到4的操作需要處理成[1,4]操作，再對[1,1]操作 ，將這個操作推廣到3維。
一維的是區間管區間，二維是方塊管理方塊，三維是方體管理方體，本質跟一維一樣。
#include <stdio.h>
#include <stdlib.h>
#include <iostream>
#define maxn 110
using namespace std;
int c[maxn][maxn][maxn];
int n;
void set()
{
int i, j, k;
for (i = 0; i <= n; i++)
{
for (j = 0; j <= n; j++)
{
for (k = 0; k <= n; k++)
{
c[i][j][k] = 0;
}
}
}
}
inline int lowbit(int x)
{
return x & (-x);
}
void update(int x, int y, int z)//分塊更新或操作
{
int i, j, k;
for (i = x; i > 0; i -= lowbit(i))
{
for (j = y; j > 0; j -= lowbit(j))
{
for (k = z; k > 0; k -= lowbit(k))
{
c[i][j][k]++;
}
}
}
}
int sum(int x, int y, int z)//向上疊加每個父節點的操作次數
{
int i, j, k, var = 0;
for (i = x; i <= n; i += lowbit(i))
{
for (j = y; j <= n; j += lowbit(j))
{
for (k = z; k <= n; k += lowbit(k))
{
var += c[i][j][k];
}
}
}
return var;
}
int main()
{
int m, op, x1, y1, z1, x2, y2, z2, ans;
while (scanf("%d%d", &n, &m) - EOF)
{
set();
while (m--)
{
scanf("%d%d%d%d", &op, &x1, &y1, &z1);
if (op == 1)
{
scanf("%d%d%d", &x2, &y2, &z2);
update(x2, y2, z2);
update(x2, y1 - 1, z2);
update(x1 - 1, y2, z2);
update(x1 - 1, y1 - 1, z2);
//--------
update(x2, y2, z1 - 1);
update(x2, y1 - 1, z1 - 1);
update(x1 - 1, y2, z1 - 1);
update(x1 - 1, y1 - 1, z1 - 1);
}
else
{
ans = sum(x1, y1, z1);
printf("%d\n", ans % 2);
}
}
}
return 0;
}
/*
2 100
1 2 2 2 2 2 2
0 1 1 1
1
2 10
1 2 1 2 2
Q 2 2
1 2 1 2 1
Q 1 1
1 1 1 2 1
1 1 2 1 2
1 1 1 2 2
Q 1 1
1 1 1 2 1
Q 2 1
*/