【玲珑杯 round#18 A】计算几何你瞎暴力

【Link】:http://www.ifrog.cc/acm/problem/1143?contest=1020&no=0

【Description】

【Solution】



因为每个点的(xi,yi)里面,

0<=xi,yi<=10

则,开一个d[11][11][11]的数组统计每个点出现的次数;

最多只会有11^3个点;

则暴力枚举第一个教室的点(x0,y0,z0);

枚举第二个教室的点(x1,y1,z1);

设x为这两个点的距离

是相同的点的话

dis[x]+=A2d[x0][y0][z0]

是不同的点的话

dis[x]+=d[x0][y0][z0]∗d[x1][y1][z1];

把dis数组搞一个前缀和

我每个点都是从0..11枚举;

所以肯定会重一次的

最后对于每个询问r输出dis[min(r,max[x])]



【NumberOf WA】



1



【Reviw】



又是忘记初始化。



【Code】

#include <bits/stdc++.h>
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define LL long long
#define rep1(i,a,b) for (int i = a;i <= b;i++)
#define rep2(i,a,b) for (int i = a;i >= b;i--)
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define ms(x,y) memset(x,y,sizeof x)
#define Open() freopen("F:\\rush.txt","r",stdin)
#define Close() ios::sync_with_stdio(0)

typedef pair<int,int> pii;
typedef pair<LL,LL> pll;

const int dx[9] = {0,1,-1,0,0,-1,-1,1,1};
const int dy[9] = {0,0,0,-1,1,-1,1,-1,1};
const double pi = acos(-1.0);
const int N = 12;

int n,q;
LL d[N][N][N];
LL dis[10000];

int main(){
    //freopen("F:\\rush.txt","r",stdin);
    int T;
    scanf("%d",&T);
    while (T--){
        rep1(i,0,10)
            rep1(j,0,10)
                rep1(k,0,10)
                    d[i][j][k] = 0;
        rep1(i,0,9999) dis[i] = 0;
        scanf("%d%d",&n,&q);
        rep1(i,1,n){
            int x,y,z;
            scanf("%d%d%d",&x,&y,&z);
            d[x][y][z]++;
        }
        int ma = 0;
        rep1(x0,0,10)
            rep1(y0,0,10)
                rep1(z0,0,10)
                    rep1(x1,0,10)
                        rep1(y1,0,10)
                            rep1(z1,0,10){
                                int t = abs(x0-x1)+abs(y0-y1)+abs(z0-z1);
                                if (x0==x1 && y0==y1 && z0==z1) {
                                        LL temp = d[x0][y0][z0];
                                        dis[t] += temp*(temp-1);
                                        continue;
                                }
                                ma = max(ma,t);
                                dis[t]+=d[x0][y0][z0]*d[x1][y1][z1];
                            }
        rep1(i,1,ma)
            dis[i] = dis[i-1] + dis[i];
        rep1(i,1,q){
            int x;
            cin >> x;
            x = min(x,ma);
            printf("%lld\n",dis[x]/2);
        }
    }
    //Open();
    //Close();
    return 0;
}