总结:
离散点对应回稀疏矩阵时可以分块处理,根据限制条件筛选掉不必处理的块,缩短计算时间。
一、题目要求
题目背景
暑假要到了。可惜由于种种原因,小 P 原本的出游计划取消。失望的小 P 只能留在西西艾弗岛上度过一个略显单调的假期……直到……
某天,小 P 获得了一张神秘的藏宝图。
问题描述
西西艾弗岛上种有 n 棵树,这些树的具体位置记录在一张绿化图上。
简单地说,西西艾弗岛绿化图可以视作一个大小为 (L+1)×(L+1) 的 01 矩阵 A,
地图左下角(坐标 (0,0))和右上角(坐标 (L,L))分别对应 A[0][0] 和 A[L][L]。
其中 A[i][j]=1 表示坐标 (i,j) 处种有一棵树,A[i][j]=0 则表示坐标 (i,j) 处没有树。
换言之,矩阵 A 中有且仅有的 n 个 1 展示了西西艾弗岛上 n 棵树的具体位置。
传说,大冒险家顿顿的宝藏就埋藏在某棵树下。
并且,顿顿还从西西艾弗岛的绿化图上剪下了一小块,制作成藏宝图指示其位置。
具体来说,藏宝图可以看作一个大小为 (S+1)×(S+1) 的 01 矩阵 B(S 远小于 L),对应着 A 中的某一部分。
理论上,绿化图 A 中存在着一处坐标 (x,y)(0≤x,y≤L−S)与藏宝图 B 左下角 (0,0) 相对应,即满足:
对 B 上任意一处坐标 (i,j)(0≤i,j≤S),都有 A[x+i][y+j]=B[i][j]。
当上述条件满足时,我们就认为藏宝图 B 对应着绿化图 A 中左下角为 (x,y)、右上角为 (x+S,y+S) 的区域。
实际上,考虑到藏宝图仅描绘了很小的一个范围,满足上述条件的坐标 (x,y) 很可能存在多个。
请结合西西艾弗岛绿化图中 n 棵树的位置,以及小 P 手中的藏宝图,判断绿化图中有多少处坐标满足条件。
特别地,藏宝图左下角位置一定是一棵树,即 A[x][y]=B[0][0]=1,表示了宝藏埋藏的位置。
输入格式
从标准输入读入数据。
输入的第一行包含空格分隔的三个正整数 n、L 和 S,分别表示西西艾弗岛上树的棵数、绿化图和藏宝图的大小。
由于绿化图尺寸过大,输入数据中仅包含 n 棵树的坐标而非完整的地图;即接下来 n 行每行包含空格分隔的两个整数 x 和 y,表示一棵树的坐标,满足 0≤x,y≤L 且同一坐标不会重复出现。
最后 (S+1) 行输入小 P 手中完整的藏宝图,其中第 i 行(0≤i≤S)包含空格分隔的 (S+1) 个 0 和 1,表示 B[S−i][0]⋯B[S−i][S]。
需要注意,最先输入的是 B[S][0]⋯B[S][S] 一行,B[0][0]⋯B[0][S] 一行最后输入。
输出格式
输出到标准输出。
输出一个整数,表示绿化图中有多少处坐标可以与藏宝图左下角对应,即可能埋藏着顿顿的宝藏。
样例 1 输入
5 100 2
0 0
1 1
2 2
3 3
4 4
0 0 1
0 1 0
1 0 0
Data
样例 1 输出
3
二、我的解法(0)
看到题目一开始往图论方向想,但是很快发现并不是图论问题。本题的大地图只给出了离散点,而笔者浅薄的知识里只知道怎么把矩阵离散化,不知道怎么根据离散点反过来对应稀疏矩阵。由于没见过,想当然地把问题想复杂了,最终也找不到合适的解法。吸取教训,在实际考试时哪怕找不到思路也要先暴力求解,能拿一部分分总好过零分。
在csdn上搜到的一些解法中有根据前缀和先确定树的总数相等再将矩阵逐个元素比较的(我没看懂具体实现,只懂了大致思想,欢迎看懂的朋友评论区里赐教),有调用map的函数直接进行比较的(感觉还是自己实现比较能学到东西,没有细看),最后找到了根据树构建局部地图再比较的方法,简单易懂,豁然开朗,下面给出该方法的实现。
三、满分解法
#include<iostream>
using namespace std;
int a[1001][2],b[51][51];
int temp[51][51];
int main(){
int n,l,s;
cin>>n>>l>>s;
for(int i=1;i<=n;i++){
cin>>a[i][0]>>a[i][1];
}
for(int i=s;i>=0;i--){//注意行是逆序
for(int j=0;j<=s;j++){
cin>>b[i][j];
}
}
int res=0;
for(int i=1;i<=n;i++){//遍历选取基准树点
if(a[i][0]+s>l||a[i][1]+s>l) continue;//该树对应的局部地图超出范围,不予考虑
for(int x=0;x<=s;x++)//把局部地图初始化
for(int y=0;y<=s;y++)
temp[x][y]=0;
for(int j=1;j<=n;j++){//找出存在于局部地图里的树
int dx=a[j][0]-a[i][0];
int dy=a[j][1]-a[i][1];
if(dx<=s&&dy<=s&&dx>=0&&dy>=0)
temp[dx][dy]=1;
}
bool success=true;
for(int x=0;x<=s;x++){
for(int y=0;y<=s;y++){
if(temp[x][y]!=b[x][y]){//比较局部地图与藏宝图
success=false;
break;
}
}
if(!success) break;
}
if(success) res++;//比较成功,藏宝点增加
}
cout<<res;
return 0;
}
分析:
题目超时的原因在于存储大地图的稀疏矩阵太大,遍历比较耗时长。前缀和解法减少了遍历的次数,本解法利用局部地图减少了遍历的范围。这带给我一个启发,由离散点对应回稀疏矩阵时,也可以分块地对应,逐块的处理效果同整个稀疏矩阵一次性处理相同,却能根据限制条件筛选掉很多不必处理的块,大大减少了处理时间。