[Luogu] 金字塔神话
这是出给pj的题ccccc
#include <bits/stdc++.h> #define ll long long
#define INF 2147483647
#define y1 uihaoissa
#define y2 oaioiuwio using namespace std;
const int wx[]= {,,,-},wy[]= {,,-,}; int T,n,m,k,maxh,ltt;
int a[][],sum[][],cal[];
char tow[][][];
bool f[][][];
int cnt[][][];
short pre1[][][],pre2[][][],fl[][][];
bool debug; inline int read() {
int x=;
char ch=getchar();
while(ch<''||ch>'')ch=getchar();
while(ch>=''&&ch<='') {
x=x*+ch-'';
ch=getchar();
}
return x;
} inline int sqr(int x) {
return x*x;
} inline int toedge(int x,int y) {
return min(min(x,y),min(n-x+,m-y+));
} inline int islarge1(int x1,int y1,int y2,int z) {
return pre1[x1][y2][z]-pre1[x1][y1-][z];
} inline int islarge2(int x1,int x2,int y1,int z) {
return pre2[x2][y1][z]-pre2[x1-][y1][z];
} inline int getsum(int x1,int y1,int x2,int y2) {
return sum[x2][y2]-sum[x1-][y2]-sum[x2][y1-]+sum[x1-][y1-];
} int check(int x,int y,int z) {
if (x>) {
int h;
if (cnt[x-][y][z]==-)
switch (tow[x-][y][z]) {
case :
if (fl[x-][y][z]!=z-) {
h=fl[x-][y][z]+;
if (islarge1(x-h,y-h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
}
break;
case :
h=fl[x-][y][z];
if (islarge2(x-h,x+h,y-h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
break;
case :
h=fl[x-][y][z];
if (islarge2(x-h,x+h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
break;
case :
h=fl[x-][y][z]-;
if (islarge1(x+h,y-h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
break;
}
}
if (y>) {
int h;
if (cnt[x][y-][z]==-) {
if (fl[x][y][z]==) tow[x][y][z]=;
switch (tow[x][y-][z]) {
case :
h=fl[x][y-][z];
if (islarge1(x-h,y-h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
break;
case :
if (fl[x][y-][z]!=z-) {
h=fl[x][y-][z]+;
if (islarge2(x-h,x+h,y-h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
}
break;
case :
h=fl[x][y-][z]-;
if (islarge2(x-h,x+h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
break;
case :
h=fl[x][y-][z];
if (islarge1(x+h,y-h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
break;
}
}
}
if (z>) {
if (cnt[x][y][z-]==-) {
int h=fl[x][y][z-];
switch (tow[x][y][z-]) {
case :
if (islarge1(x-h,y-h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
break;
case :
if (islarge2(x-h,x+h,y-h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
break;
case :
if (islarge2(x-h,x+h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
break;
case :
if (islarge1(x+h,y-h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
break;
}
}
}
for (int h=z-; h>=; h--) { //h为中心到该圈的差数
if (cnt[x][y][h+]!=-) break;
if (islarge1(x-h,y-h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
if (islarge2(x-h,x+h,y-h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
if (islarge2(x-h,x+h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
if (islarge1(x+h,y-h,y+h,z-h)>) return fl[x][y][z]=h,tow[x][y][z]=,cnt[x][y][z]=-,;
}
int tot=cal[z]-getsum(x-z+,y-z+,x+z-,y+z-);
return tot>k?:(cnt[x][y][z]=tot,);
} int main() {
for (int i=; i<=; i++)
cal[i]=cal[i-]+(*i-)*(*i-);
T=read();
while (T--) {
maxh=;
memset(cnt,-,sizeof(cnt));
n=read(),m=read(),k=read();
for (int i=; i<=n; i++)
for (int j=; j<=m; j++) {
a[i][j]=read();
sum[i][j]=sum[i-][j]+sum[i][j-]-sum[i-][j-]+a[i][j];
for (int k=,to=toedge(i,j); k<=to; k++) {
f[i][j][k]=(a[i][j]>k);
pre1[i][j][k]=pre1[i][j-][k]+f[i][j][k];
pre2[i][j][k]=pre2[i-][j][k]+f[i][j][k];
}
}
for (int i=; i<=n; i++)
for (int j=; j<=m; j++) {
cnt[i][j][]=-a[i][j];
if (a[i][j]<=) maxh=max(maxh,);
for (int l=max(max(a[i][j],),maxh+),to=toedge(i,j); l<=to;) {
int tmp=check(i,j,l);
if (!tmp) maxh=max(maxh,l),tmp++;
l+=tmp;
}
}
printf("%d\n",maxh);
}
return ;
}
[Luogu] 金字塔神话的更多相关文章
- luogu P1216 [IOI1994][USACO1.5]数字三角形 Number Triangles (递推)
链接:https://www.luogu.org/problemnew/show/P1216 题面: 题目描述 观察下面的数字金字塔. 写一个程序来查找从最高点到底部任意处结束的路径,使路径经过数字的 ...
- Atitit.软件开发的三层结构isv金字塔模型
Atitit.软件开发的三层结构isv金字塔模型 第一层,Implements 层,着重与功能的实现.. 第二次,spec层,理论层,设计规范,接口,等.流程.方法论 顶层,val层,价值观层,原则, ...
- 在Excel中制作金字塔条形图
使用场景:一项市场调查研究中,男性和女性.赞同和反对.满意和不满意的两方面的消费者,他们在某些项目上的指标分布特性一项产品组合决策中,乐观场景和悲观场景下各产品的获利情况一个产品试销活动中,不同门店渠 ...
- Luogu 魔法学院杯-第二弹(萌新的第一法blog)
虽然有点久远 还是放一下吧. 传送门:https://www.luogu.org/contest/show?tid=754 第一题 沉迷游戏,伤感情 #include <queue> ...
- luogu p1268 树的重量——构造,真正考验编程能力
题目链接:http://www.luogu.org/problem/show?pid=1268#sub -------- 这道题费了我不少心思= =其实思路和标称毫无差别,但是由于不习惯ACM风格的题 ...
- SIFT中的尺度空间和传统图像金字塔
SIFT中的尺度空间和传统图像金字塔 http://www.zhizhihu.com/html/y2010/2146.html 最近自己混淆了好多概念,一边弄明白的同时,也做了一些记录,分享一下.最近 ...
- Atitit 图像金字塔原理与概率 attilax的理解总结qb23
Atitit 图像金字塔原理与概率 attilax的理解总结qb23 1.1. 高斯金字塔 ( Gaussianpyramid): 拉普拉斯金字塔 (Laplacianpyramid):1 1.2 ...
- 软件项目发展历史<人月神话>这本书好
几乎是计算机软件开发的发展历史 人月神话,增加人手并不一定能提高开发速度. 原因在于,有些任务是无法分解的,存在先后顺序.无法同步进行. 增加人手,增加的是沟通成本,相互牵制.可以分解的任务就 ...
- OpenCV图像金字塔:高斯金字塔、拉普拉斯金字塔与图片尺寸缩放
这篇已经写得很好,真心给作者点个赞.题目都是直接转过来的,直接去看吧. Reference Link : http://blog.csdn.net/poem_qianmo/article/detail ...
随机推荐
- 关于UBOOT,LINUX内核编译,根文件系统的15个小问题
(1)内核默认运行地址和加载地址在哪里设置? 由 arch/arm/kernel/vmlinux.lds.S 生成的 arch/armkernel/vmlinux.lds决定 (2)从FLASH什 ...
- vi 使用系统剪贴板(clipboard)
ref : https://www.jianshu.com/p/771b95e34293 http://www.bubuko.com/infodetail-469867.html 在vi中,如果编译时 ...
- LOJ2557 CTSC2018组合数问题(提交答案)
直接利用simulator退火应该可以得到大量分数. op=1:1,4,5,6,10 即构造序列{ai},最小化Σti,ai+rai,aj. 1:暴搜/退火. 4:观察到图大致成一条链(注意其中有两个 ...
- C# LoadXml System.Xml.XmlException: Data at the root level is invalid. Line 1, position 1.
去掉BOM头 writer = new XmlTextWriter(stream, new UnicodeEncoding(false,false)); 如果是UTF8 writer = new Xm ...
- ggalluvial|TCGA临床数据绘制桑基图(Sankey)
本文首发于”生信补给站“,https://mp.weixin.qq.com/s/yhMgkST-rVD6SaQS7R-eoA 桑基图(Sankey diagram),是一种特定类型的流程图,图中延伸的 ...
- 奇妙的算法【7】-贪婪算法-dp
问题1描述:[贪婪算法,Dijistra算法] ①有一只兔子要从一个N*N的二维矩阵方格中从上跳到下面: ②每次只能向左或向下,越过一个方格,跳到下一个方格中: ③被越过的方格中的数值,表示该兔子越过 ...
- 通过gpio控制一个进程开启或关闭
目标: 板子上有个进程需要通过读取gpio的值, 当gpio值为1 时, 开启指定的进程,当gpio为0时, 杀掉这个指定的进程. #include <stdio.h> int main( ...
- centOS 7单机安装 kong
kong 网关 单机部署 环境:centOS 7:依赖:jdk1.8 安装内容:postgresql数据库, kong 网关,nodeJs和npm,kong Dashboard (可视化管理界面) 版 ...
- nginx分割日志
实现nginx日志按照时间分割存储 backups_log.sh #设置log日志的存储地址 LOG_PATH=/home/soft/nginx/logs #设置历史日志的存储地址 HISTORY_L ...
- springboot2.0(二)
三. Web开发 3.1.静态资源访问 在我们开发Web应用的时候,需要引用大量的js.css.图片等静态资源. 默认配置 Spring Boot默认提供静态资源目录位置需置于classpath下,目 ...