ccf 202104-2 邻域均值

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

问题描述

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

输入格式

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

 

 

样例输入
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

评测用例规模与约定

 

对每个点遍历，根据r以及Aij确定领域范围，对邻域范围的点再进行遍历

for循环次数过多，超出时间限制 

#include<iostream>
#include<algorithm>
#include<bits/stdc++.h>
using namespace std;
//70分 暴力解法
int main()
{
    int n,L,r,t;
    cin>>n>>L>>r>>t;
    int a[n][n];
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
        {
            cin>>a[i][j];
        }
    }
    int cnt=0;
    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)//遍历每一个点
        {
            int num=0,nm=0;
            for(int u=max(i-r,0);u<=min(r+i,n-1);u++)//确定范围
            {
                for(int v=max(j-r,0);v<=min(j+r,n-1);v++)
                {
                    num+=a[u][v];//灰度值累计
                    nm++;//区域内点数累加
                }
            }

            if(num<=t*nm)//不要用num/nm<=t 出错
            {
                cnt++;
            }
            num=nm=0;
        }
    }
    cout<<cnt<<endl;
    return 0;
}

想办法减少循环次数-----------前缀和 num【i】+=num【i-1】

每一行使用前缀和a【i】+=a【i-1】记录输入信息，

确定邻域范围及其内点的个数

            up=max(i-r,1);
            int down=min(i+r,n);
            int left=max(j-r,1);
            int right=min(j+r,n);
            int nm=(right-left+1)*(down-up+1);

然后对邻域内的每一行灰度值进行计算a【u】【right】-a【u】【left-1】（up<=u<=down)

最后通过比较即得到结果。

​
#include<bits/stdc++.h>
using namespace std;
//100分（每行前缀和）
int a[601][601]={0};
int main()
{
    int n,L,r,t;
    cin>>n>>L>>r>>t;
    int sum=0;
    memset(a,0,sizeof(a));
    memset(x,0,sizeof(x));
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            cin>>a[i][j];
            a[i][j]+=a[i][j-1];//每一行灰度值累加记录
        }
    }
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=n;j++)
        {
            int up=max(i-r,1);
            int down=min(i+r,n);
            int left=max(j-r,1);
            int right=min(j+r,n);
            int nm=(right-left+1)*(down-up+1);
            int t1=0;
          for(int u=up;u<=down;u++)//邻域内灰度值累加
          {
              t1+=a[u][right]-a[u][left-1];//left-1需要注意一下
          }
          if(t1<=nm*t)sum++;
        }
    }
    cout<<sum<<endl;
    return 0;
}

​

虽然只减少了一重循环，但时间复杂度却大大减小。

另外，前缀和在降低时间复杂度方面非常重要，大家可以了解一下哈！

如有错误，欢迎指正！！！