BZOJ1218: [HNOI2003]激光炸弹

1218: [HNOI2003]激光炸弹

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 1041 Solved: 500
[Submit][Status]

Description

一种新型的激光炸弹，可以摧毁一个边长为R的正方形内的所有的目标。现在地图上有n(N<=10000)个目标，用整数Xi,Yi(其值在[0,5000])表示目标在地图上的位置，每个目标都有一个价值。激光炸弹的投放是通过卫星定位的，但其有一个缺点，就是其爆破范围，即那个边长为R的正方形的边必须和x，y轴平行。若目标位于爆破正方形的边上，该目标将不会被摧毁。 0

Input

输入文件的第一行为正整数n和正整数R，接下来的n行每行有3个正整数，分别表示

Output

输出文件仅有一个正整数，表示一颗炸弹最多能炸掉地图上总价值为多少的目标（结果不会超过32767）。

Sample Input

2 1
0 0 1
1 1 1

Sample Output

HINT

Source

题解：

真是不能再暴力。

树状数组的思想：

f[i][j]表示1..i,1..j这个矩阵的所有数和，这样暴力枚举每个正方形的右下角的点，正方形内的和可以在O(1)时间内算出来

WA了一次，因为有可能出现有数的只在几行，<r

代码：

 1 #include<cstdio>
 2 #include<cstdlib>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<algorithm>
 6 #include<iostream>
 7 #include<vector>
 8 #include<map>
 9 #include<set>
10 #include<queue>
11 #define inf 1000000000
12 #define maxn 6100
13 #define maxm 500+100
14 #define eps 1e-10
15 #define ll long long
16 #define pa pair<int,int>
17 using namespace std;
18 inline int read()
19 {
20     int x=0,f=1;char ch=getchar();
21     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
22     while(ch>='0'&&ch<='9'){x=10*x+ch-'0';ch=getchar();}
23     return x*f;
24 }
25 int n,r,mx,my,ans,f[maxn][maxn];
26 int main()
27 {
28     freopen("input.txt","r",stdin);
29     freopen("output.txt","w",stdout);
30     n=read();r=read();mx=0;my=0;
31     for(int i=1;i<=n;i++)
32     {
33         int x=read(),y=read(),z=read();x++;y++;
34         f[x][y]+=z;
35         mx=max(mx,x);my=max(my,y);
36     }
37     //cout<<mx<<' '<<my<<endl;
38     for(int i=1;i<=min(mx+r,maxn);i++)
39      for(int j=1;j<=min(my+r,maxn);j++)
40       {
41           f[i][j]+=f[i][j-1]+f[i-1][j]-f[i-1][j-1];
42           if(i>=r&&j>=r)
43           ans=max(ans,f[i][j]-f[i-r][j]-f[i][j-r]+f[i-r][j-r]);
44           //cout<<i<<' '<<j<<' '<<f[i][j]<<endl;
45       }  
46     printf("%d
",ans);
47     return 0;
48 }

View Code