最大子矩阵和、Radar Installation、激光炸弹

1.最大子矩阵和

一个M*N的矩阵，找到此矩阵的一个子矩阵，并且这个子矩阵的元素的和是最大的，输出这个最大的值。

例如：3*3的矩阵：

-1 3 -1
2 -1 3
-3 1 2

和最大的子矩阵是：

3 -1
-1 3
1 2

Input

第1行：M和N，中间用空格隔开（2 <= M,N <= 500)。 第2 ~ N+1行：矩阵中的元素，每行M个数，中间用空格隔开。(-10^9 <= M[i] <= 10^9)

Output

输出和的最大值。如果所有数都是负数，就输出0。

Sample

Input

3 3
-1 3 -1
2 -1 3
-3 1 2

Output

7

 题解：

#include<iostream>
#include<cstring>
using namespace std;
int m,n,a[503][503],s[503][503],dp[503],mx;

void solve(int j){   //求最大子序列和
	memset(dp,0,sizeof(dp));
	for(int k=1;k<=m;k++){
		dp[k]=max(dp[k-1]+s[j][k],s[j][k]);
		if(dp[k]>mx) mx=dp[k];
	}
}

int main(){
	cin>>m>>n;
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++) cin>>a[i][j];
	}
	for(int i=1;i<=n;i++){
		memset(s,0,sizeof(s));
		for(int j=i;j<=n;j++){
			for(int k=1;k<=m;k++){
				s[j][k]=s[j-1][k]+a[j][k];
			}
			solve(j);
		}
	}
	cout<<mx;
	return 0;
}

---

2.Radar Installation

如今我们位于沿海地区，需要安装大炮，使得火力可以覆盖整个区域。海岸线可以视为是无限长的直线。陆地位于海岸线的一侧，海洋位于另一侧。海洋里有若干个岛屿，每个小岛可以视为海洋中的一个点。我们需要在海岸线上安装大炮，每个大炮智能覆盖距离d，因此海洋中的小岛被大炮安装所覆盖的条件是两者间的距离不超过 d 。 我们将海岸线视为 x 轴。海洋的一侧位于 x 轴上方，陆地的一侧位于下方。给定海洋中每个小岛的位置（以xy坐标给出），并给定大炮的覆盖距离，你需要写程序，找出安装大炮的最少数量，使得所有的小岛都被覆盖。

Input

输入由多个测试用例组成。每个测试用例的第一行，包含了两个整数 n (1<=n<=1000) 和 d，其中 n 是海洋中小岛的数目，d 是大炮可以覆盖的距离。接下来是 n 行，每行包含了两个整数x,y，表示每个小岛的坐标。每组测试用例之间，以一个空行间隔。（0<=x,y,d<=1e9)

输入终止于包含两个 0 的一行。

Output

对于每个测试用例，输出一行，包括测试用例的编号，以及安装大炮所需的最小数量。如果不能全部覆盖，则输出"-1"。

Sample Input

3 2
1 2
-3 1
2 1

1 2
0 2

0 0

Sample Output

Case 1: 2
Case 2: 1

 题解：

分析：把问题转换为区间放雷问题。 

#include<iostream> //题解1
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
int n,d;
struct map{
	double one,two;
}edg[1003];

bool cmp(map a,map b){
	return a.one<b.one;
}

int main(){
	for(int k=1; ;k++){
		scanf("%d%d",&n,&d);
		if(!n&&!d) break;
		int flag=0;
		for(int x,y,i=1;i<=n;i++){
			scanf("%d%d",&x,&y);
			if(y>d) flag=1;
			else{
				edg[i].one=x-sqrt(1.0*d*d-1.0*y*y);
				edg[i].two=x+sqrt(1.0*d*d-1.0*y*y);
			}
		}
		if(flag){
			printf("Case %d: -1\n",k);
			continue;
		}
		sort(edg+1,edg+n+1,cmp);
		int sum=1;
		double l=edg[1].one,r=edg[1].two;
		for(int i=2;i<=n;i++){
			if(edg[i].one>=l&&edg[i].two<=r){
				l=edg[i].one;
				r=edg[i].two;
			}
			else if(edg[i].one<=r&&edg[i].two>=r){
				l=edg[i].one;
			}
			else{
				l=edg[i].one;
				r=edg[i].two;
				sum++;
			}
		}
		printf("Case %d: %d\n",k,sum);
	}
	return 0;
} 

#include<iostream>  //题解2
#include<cstdio>
#include<cmath>
#include<algorithm>
using namespace std;
int n,d;
struct map{
	double one,two;
}edg[1003];

bool cmp(map a,map b){
	if(a.one==b.one) a.one>b.one;
	return a.two<b.two;
}

int main(){
	for(int k=1; ;k++){
		scanf("%d%d",&n,&d);
		if(!n&&!d) break;
		int flag=0;
		for(int x,y,i=1;i<=n;i++){
			scanf("%d%d",&x,&y);
			if(y>d) flag=1;
			else{
				edg[i].one=x-sqrt(1.0*d*d-1.0*y*y);
				edg[i].two=x+sqrt(1.0*d*d-1.0*y*y);
			}
		}
		if(flag){
			printf("Case %d: -1\n",k);
			continue;
		}
		sort(edg+1,edg+n+1,cmp);
		int sum=1;
		double t=edg[1].two;
		for(int i=2;i<=n;i++){
			if(edg[i].one>t){
				t=edg[i].two;
				sum++;
			}
		}
		printf("Case %d: %d\n",k,sum);
	}
	return 0;
} 

---

3.激光炸弹

Description

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

Input

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

Output

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

Sample

 Input

2 1
0 0 1
1 1 1

Output

1

 题解：

 分析：这里用到了二维前缀和的知识。（前缀和与差分知识传送阵）

#include<iostream>
#include<cstdio>
using namespace std;
int n,r,a[5003][5003];
int nn,mm,mx;

int main(){
	cin>>n>>r;
	nn=r,mm=r;
	for(int x,y,v,i=1;i<=n;i++){
		scanf("%d%d%d",&x,&y,&v);
		a[++x][++y]+=v;
		nn=max(nn,x);
		mm=max(mm,y);
	}
	for(int i=1;i<=nn;i++){
		for(int j=1;j<=mm;j++){
			a[i][j]+=a[i-1][j]+a[i][j-1]-a[i-1][j-1];
		}
	}
	for(int i=r;i<=nn;i++){
		for(int j=r;j<=mm;j++){
			mx=max(mx,a[i][j]-a[i-r][j]-a[i][j-r]+a[i-r][j-r]);
		}
	}
	cout<<mx;
	return 0;
}