Save the President UVA - 11589

问题

https://vjudge.net/problem/UVA-11589

分析

这题目的逻辑比较明显。
建立数组A[x][y][z][t],找到安全的长方体就可以了。
可以暴力求解，但是时间复杂度是O(10^10)左右的，因为有8个循环。
将每个不安全的地方置为1，求子长方体的元素和，就是求前缀数组，如果前缀和是0，说明安全。
注意遍历时，是从1开始的左闭右闭区间，例如两个顶点是（0，0，0）和（1，1，1）的长方体。
和例题UVA10775很像

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <map>
#include <string>
#include <vector>
#include <algorithm>
#include <queue>
using namespace std;
typedef long long LL;
const int maxn=20;
int n,x,y,z,q,kase=0,a[maxn][maxn][maxn][100];

inline void expand(int x,int &b0,int &b1,int &b2,int &b3){
    b0=(x&1);
    b1=(x>>1)&1;
    b2=(x>>2)&1;
    b3=(x>>3)&1;
}

inline int sign(int b0,int b1,int b2,int b3){
    return ((b0+b1+b2+b3)&1)?1:-1;
}

int main(void){
    while(scanf("%d%d%d%d%d",&n,&x,&y,&z,&q) && n){
        memset(a,0,sizeof(a));
        int x1,y1,z1,x2,y2,z2,hour,minute,t1,t2;
        for(int i=0;i<n;++i){
            scanf("%d%d%d",&x1,&y1,&z1);
            scanf("%d%d%d",&x2,&y2,&z2);
            scanf("%d:%d",&hour,&minute);
            t1=hour*4+minute/15;
            scanf("%d:%d",&hour,&minute);
            t2=hour*4+minute/15;
            for(int j=x1+1;j<=x2;++j){
                for(int k=y1+1;k<=y2;++k){
                    for(int l=z1+1;l<=z2;++l){
                        for(int m=t1+1;m<=t2;++m){
                            a[j][k][l][m]=1;
                        }
                    }
                }
            }
        }
        //求前缀和数组
        int b0,b1,b2,b3;
        for(int i=1;i<=x;++i){
            for(int j=1;j<=y;++j){
                for(int k=1;k<=z;++k){
                    for(int l=1;l<=96;++l){
                        for(int m=1;m<16;++m){
                            expand(m,b0,b1,b2,b3);
                            int s=sign(b0,b1,b2,b3);
                            a[i][j][k][l]+=a[i-b0][j-b1][k-b2][l-b3]*s;
                        }
                    }
                }
            }
        }
        printf("3D World %d:\n",++kase);
        for(int p=0;p<q;++p){
            int ans=0;
            scanf("%d%d%d%d:%d",&x1,&y1,&z1,&hour,&minute);
            t1=hour*4+minute/15;
            for(int i=x1;i<=x;++i){
                for(int j=y1;j<=y;++j){
                    for(int k=z1;k<=z;++k){
                        for(int l=t1;l<=96;++l){
                            int b0,b1,b2,b3,s=0;
                            //子长方体的元素和
                            for(int x=0;x<16;++x){
                                expand(x,b0,b1,b2,b3);
                                s-=a[i-x1*b0][j-y1*b1][k-z1*b2][l-t1*b3]*sign(b0,b1,b2,b3);
                            }
                            if(s==0) ++ans;
                        }
                    }
                }
            }
            if(ans) printf("%d safe place(s) found\n",ans);
            else printf("No safe place(s) found\n");
        }
    }
    return 0;
}

//1 1 1 1 1
//0 0 0 1 1 1 9:30 10:30
//1 1 1 0:30
//0 0 0 0 0