2018.06.29 NOIP模拟 繁星(前缀和)
繁星
【问题描述】
要过六一了,大川正在绞尽脑汁想送给小伙伴什么礼物呢。突然想起以前拍过一张夜空中的繁星的照片,这张照片已经被处理成黑白的,也就是说,每个像素只可能是两个颜色之一,白或黑。像素(x,y)处是一颗星星,当且仅当,像素(xxx,yyy),(x−1x-1x−1,yyy),(x+1x+1x+1,yyy),(xxx,y−1y-1y−1),(xxx,y+1y+1y+1)都是白色的。因此一个白色像素有可能属于多个星星,也有可能有的白色像素不属于任何一颗星星。但是这张照片具有研究价值,所以大川不想把整张照片都送给小伙伴,而只准备从中裁下一小块长方形照片送给他。但为了保证效果,大川认为,这一小块相片中至少应该有 kkk 颗星星。
现在大川想知道,到底有多少种方法裁下这一小块长方形相片呢?
【输入格式】
从文件 star.instar.instar.in 中输入数据。
输入的第一行包含三个正整数 nnn,mmm,kkk,意义见题目所示。
接下来 nnn 行,每行一个长度为 mmm 的字符串,字符串仅由’.‘和’‘构成,’.‘表示这个像素为黑色,’'表示这个像素为白色。
【输出格式】
输出到文件 star.outstar.outstar.out 中。
输出的第一行包含一个整数,表示大川有多少种满足题意的裁剪方法。
【样例输入】
555 666 333
∗∗∗...***...∗∗∗...
∗∗∗∗..****..∗∗∗∗..
.**.*.
∗∗∗∗∗∗******∗∗∗∗∗∗
.∗.∗∗∗.*.***.∗.∗∗∗
【样例输出】
333
【样例说明】
图中共有 444 颗星星,分别位于第 222 行第 222 列、第 222 行第 333 列、第 444 行第 222 列、第444 行第 555 列。
有 333 种符合题意的选择方法(以左上角行列 - 右下角行列方式给出): (1,11,11,1)-(5,45,45,4)
(1,11,11,1)-(5,55,55,5) (1,11,11,1)-(5,65,65,6)
【数据规模与约定】
对于 202020%的数据,满足 NNN,MMM<=202020. 对于 404040%的数据,满足 N,M<=100N,M<=100N,M<=100. 对于 707070%的数据,满足 N,MN,MN,M<=200200200. 对于 100100100%的数据,满足 N,MN,MN,M<=500500500,0<k<N∗M0<k<N*M0<k<N∗M.
202020%做法:
枚举左上角、右下角,暴力统计其中的星星个数。O(N6N^6N6)
404040%做法:
首先发现查询矩形内星星个数可以用部分和在O(N2N^2N2)预处理后O(111)回答。
枚举左上角、右下角,用部分和计算其中星星个数。O(N4N^4N4)
707070%做法:
枚举子矩阵的左、右边界。然后发现对于某个给定的下边界,可行的上边界必然是连续一段。
因此枚举左、右边界后用部分和维护这一条里的星星,枚举下边界,二分查找上边界。O(N3lgNN^3lgNN3lgN)
100100100%做法:
发现左右边界确定后,随着下边界的向下移动,上边界也只会向下移动或不动。
于是枚举左、右边界后用部分和维护星星,然后用一个滑窗维护可行的上边界。O(N3N^3N3)
此题十分垃圾,考试时二维前缀和加双指针水过,下面来一发考试代码。
#include<cstdio>
#include<cstdlib>
#include<string>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<ctime>
#include<iostream>
#include<queue>
#include<stack>
#include<vector>
#include<set>
using namespace std;
#define N 505
inline long long read(){
long long ans=0,w=1;
char ch=getchar();
while(!isdigit(ch)){
if(ch=='-')w=-1;
ch=getchar();
}
while(isdigit(ch)){
ans=(ans<<3)+(ans<<1)+ch-'0';
ch=getchar();
}
return ans*w;
}
inline void write(long long x){
if(x<0){
x=-x;
putchar('-');
}
if(x>9)write(x/10);
putchar(x%10+'0');
}
int n,m,t,b[N],sum[N][N];
long long cnt=0;
bool map[N][N],ok[N];
char s[N];
inline int total(int x1,int y1,int x2,int y2){return sum[x2][y2]-sum[x1-1][y2]-sum[x2][y1-1]+sum[x1-1][y1-1];}
inline void solve1(){
for(int len=1;len<=n;++len)
for(int i=len+1;i<n;++i){
for(int j=1;j<=m;++j)b[j]=sum[i][j]-sum[i-len][j];
if(b[m]<t)continue;
int l=2,r=2;
while(r<=m-1){
if(b[r]-b[l-1]>=t){
cnt+=m-r;
++l;
}
else ++r;
}
}
}
int main(){
freopen("star.in","r",stdin);
freopen("star.out","w",stdout);
n=read();
m=read();
t=read();
for(int i=1;i<=n;++i){
scanf("%s",s);
for(int j=0;j<m;++j){
if(s[j]=='*')map[i][j+1]=true;
else map[i][j+1]=false;
}
}
for(int i=2;i<=n;++i)
for(int j=2;j<=m;++j){
if(map[i][j]&&map[i-1][j]&&map[i+1][j]&&map[i][j-1]&&map[i][j+1])sum[i][j]=1;
sum[i][j]=sum[i-1][j]+sum[i][j-1]+sum[i][j]-sum[i-1][j-1];
}
if(sum[n][m]<t){
printf("0");
return 0;
}
solve1();
write(cnt);
return 0;
}
2018.06.29 NOIP模拟 繁星(前缀和)的更多相关文章
- 2018.06.29 NOIP模拟 区间(前缀和差量)
区间(interval.cpp) 时限:2000ms 空间限制:512MB [问题描述] 给出一个长度为 n 的序列 a[1]-a[n]. 给出 q 组询问,每组询问形如<x,y>< ...
- 2018.06.29 NOIP模拟 Gcd(容斥原理)
Gcd 题目背景 SOURCE:NOIP2015-SHY-2 题目描述 给出n个正整数,放入数组 a 里. 问有多少组方案,使得我从 n 个数里取出一个子集,这个子集的 gcd 不为 1 ,然后我再从 ...
- 2018.06.29 NOIP模拟 旅馆(线段树)
旅馆 [问题描述] OIEROIEROIER 们最近的旅游计划,是到长春净月潭,享受那里的湖光山色,以及明 媚的阳光.你作为整个旅游的策划者和负责人,选择在潭边的一家著名的旅馆住 宿.这个巨大的旅馆一 ...
- 2018.06.29 NOIP模拟 排列(线段树)
排列(premu.cpp) [题目描述] 对于一个 1 到 n 的排列,逆序数的定义为:排列中第 i 位 ai的逆序数就是 a1-ai-1中比 ai大的数的个数.另外用 pi表示 a1,-,ai的逆序 ...
- 2018.06.29 NOIP模拟 Minimum(最小生成树)
Minimum 题目背景 SOURCE:NOIP2015-SHY-2 题目描述 给出一幅由 n 个点 m 条边构成的无向带权图. 其中有些点是黑点,另外点是白点. 现在每个白点都要与他距离最近的所有黑 ...
- 2018.06.29 NOIP模拟 1807(简单递推)
1807 题目背景 SOURCE:NOIP2015-SHY-2 题目描述 给出一个由数字('0'-'9')构成的字符串.我们说一个子序列是好的,如果他的每一位都是 1.8.0.7 ,并且这四个数字按照 ...
- 2018.06.29 NOIP模拟 边的处理(分治+dp)
边的处理(side.cpp) [问题描述] 有一个 n 个点的无向图,给出 m 条边,每条边的信息形如<x,y,c,r><x,y,c,r><x,y,c,r>. 给出 ...
- 2018.08.29 NOIP模拟 pmatrix(线性筛)
[问题描述] 根据哥德巴赫猜想(每个不小于 6 的偶数都可以表示为两个奇素数之和),定义 哥德巴赫矩阵 A 如下:对于正整数对(i,j),若 i+j 为偶数且 i,j 均为奇素数,则 Ai,j = 1 ...
- 2018.06.27 NOIP模拟 节目(支配树+可持久化线段树)
题目背景 SOURCE:NOIP2015-GDZSJNZX(难) 题目描述 学校一年一度的学生艺术节开始啦!在这次的艺术节上总共有 N 个节目,并且总共也有 N 个舞台供大家表演.其中第 i 个节目的 ...
随机推荐
- linux c++下载http文件并显示进度<转>
#include <stdio.h> #include <sys/types.h> #include <sys/socket.h> #include <net ...
- 4 python 类的专有方法介绍
1.__init__ : 构造函数,在生成对象时调用 该方法是在对象产生之后才会执行,只用来为对象进行初始化操作,可以有任意代码,但不一定有返回值. 所谓初始化构造函数就是在构造对象的同时被对象自动 ...
- Eclipse安装Svn无法连接问题
在Eclipse市场上安装完SVN插件后连接SVN时出现以下错误: SVN: '0x00400006: Validate Repository Location' operation finished ...
- [转]Gulp思维 —— Gulp高级技巧
感受过gulp.js带来的兴奋过后,你需要的不仅仅是它的光鲜,而是切切实实的实例.这篇文章讨论了一些使用gulp.js时常踩的坑,以及一些更加高级和定制化的插件和流的使用技巧. 基本任务 gulp的基 ...
- WEB性能测试工具
做Web开发,难免要对自己开发的页面进行性能检测,自己写工具检测,工作量太大.网上有几款比较成熟的检测工具,以下就介绍一下,与大家分享. 互联网现有工具 基于网页分析工具: 1. 阿里测 ...
- AUC理解
https://www.zhihu.com/question/39840928 机器学习和统计里面的auc怎么理解?
- python引入模块时import与from ... import的区别(转)
import datetime是引入整个datetime包,如果使用datetime包中的datetime类,需要加上模块名的限定. 1 import datetime 2 3 print datet ...
- CentOS7.0重置Root的密码 (忘记密码)
首先进入开启菜单,按下e键进入编辑现有的内核,如下图所示 然后滚动列表,找到ro,将它替换成rw,并加上init=/sysroot/bin/sh,最终变为如下图 然后按CTRL+X进入到单用户模式,在 ...
- ie6浏览器的安装
试过各种方法都不行,最后用ie8卸载工具,不仅卸载了ie8还自动安装了ie6,非常棒!!!我的空间有那款卸载工具.
- Island Transport
Island Transport http://acm.hdu.edu.cn/showproblem.php?pid=4280 Time Limit: 20000/10000 MS (Java/Oth ...