CCF202104-2 邻域均值

题目选自CCF2020104-2

问题分析：按照题意处理即可。需要注意处理时间，如果没有合适处理则有可能得不了满分100分。所以需要使用前缀和，或者使用二维前缀和来实现。

题记：前缀和用于n次不同的数列区间求和十分有效，可以避免大量重复求和计算。二维前缀和用于n次不同的子矩阵求和十分有效，可以避免大量重复求和计算。

试题编号： 202104-2
试题名称：邻域均值
时间限制： 1.0s
内存限制： 512.0MB
问题描述：

试题背景

顿顿在学习了数字图像处理后，想要对手上的一副灰度图像进行降噪处理。不过该图像仅在较暗区域有很多噪点，如果贸然对全图进行降噪，会在抹去噪点的同时也模糊了原有图像。因此顿顿打算先使用邻域均值来判断一个像素是否处于较暗区域，然后仅对处于较暗区域的像素进行降噪处理。

问题描述

待处理的灰度图像长宽皆为n个像素，可以表示为一个n×n大小的矩阵A，其中每个元素是一个[0,L)范围内的整数，表示对应位置像素的灰度值。
对于矩阵中任意一个元素Aij（0≤i,j＜n），其邻域定义为附近若干元素的集和：
N e i g h b o r ( i , j , r ) = { A x y ∣ 0 ≤ x , y ＜ n 　 a n d 　 ∣ x − i ∣ ≤ r 　 a n d 　 ∣ y − j ∣ ≤ r } Neighbor(i,j,r)=\{A_{xy}|0≤x,y＜n　and　|x-i|≤r　and　|y-j|≤r\}
Neighbor(i,j,r)={A xy∣0≤x,y＜n　and　∣x−i∣≤r　and　∣y−j∣≤r}

这里使用了一个额外的参数r来指明Aij附近元素的具体范围。根据定义，易知Neighbor(i,j,r)最多有(2r+1)2个元素。

如果元素Aij邻域中所有元素的平均值小于或等于一个给定的阈值t，我们就认为该元素对应位置的像素处于较暗区域。
下图给出了两个例子，左侧图像的较暗区域在右侧图像中展示为黑色，其余区域展示为白色。

现给定邻域参数r和阈值t，试统计输入灰度图像中有多少像素处于较暗区域。

输入格式

输入共n+1行。

输入的第一行包含四个用空格分隔的正整数n、L、r和t，含义如前文所述。

第二到第n+1行输入矩阵A。
第i+2（0≤i＜n）行包含用空格分隔的n个整数，依次为Ai0，Ai1，…，Ai(n-1)。

输出格式

输出一个整数，表示输入灰度图像中处于较暗区域的像素总数。

样例输入

4 16 1 6
0 1 2 3
4 5 6 7
8 9 10 11
12 13 14 15

样例输出

7

样例输入

11 8 2 2
0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0
0 7 0 0 0 7 0 0 7 7 0
7 0 7 0 7 0 7 0 7 0 7
7 0 0 0 7 0 0 0 7 0 7
7 0 0 0 0 7 0 0 7 7 0
7 0 0 0 0 0 7 0 7 0 0
7 0 7 0 7 0 7 0 7 0 0
0 7 0 0 0 7 0 0 7 0 0
0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0

样例输出

83

评测用例规模与约定

70%的测试数据满足n≤100、r≤10。

全部的测试数据满足0＜n≤600、 0＜r≤100且2≤t＜L≤256。

解题代码：

#include <iostream>
#include <algorithm>

using namespace std;

const int N = 610;

int n, L, r, t;
int s[N][N];

int get_sum(int x1, int y1, int x2, int y2)
{
    return s[x2][y2] - s[x1 - 1][y2] - s[x2][y1 - 1] + s[x1 - 1][y1 - 1];
}

int get_cnt(int x1, int y1, int x2, int y2)
{
    return (x2 - x1 + 1) * (y2 - y1 + 1);
}

int main()
{
    scanf("%d%d%d%d", &n, &L, &r, &t);
    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= n; j ++ )
        {
            int x;
            scanf("%d", &x);
            s[i][j] = x + s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1];
        }

    int res = 0;
    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= n; j ++ )
        {
            int x1 = max(1, i - r), y1 = max(1, j - r);
            int x2 = min(n, i + r), y2 = min(n, j + r);
            if (get_sum(x1, y1, x2, y2) <= t * get_cnt(x1, y1, x2, y2))
                res ++ ;
        }
    printf("%d\n", res);
    return 0;
}