2017.10.12 礼物(zjoi2011) 失败总结

听说n^4能过，于是愉快地敲了n^4 ，还带了两个最优性剪枝。。

结果事实证明是别人家的n^4.。。（~~wys？~~）

很奇怪玄学剪枝有时2^50次方都能过，为什么这个连1e10都过不了。.。

TLE代码：

#include<iostream>
#include<cstdio>
using namespace std;
int x,y,z,xx[155][155][155],yy[155][155][155],zz[155][155][155],i,j,k,lin[155],a,qi,tu[155][155][155],ans;
char ch;
int main()
{
    scanf("%d%d%d",&y,&x,&z);
    for(i=1;i<=x;i++)
    for(j=1;j<=y;j++)
    for(k=1;k<=z;k++)
    {
        scanf("%c",&ch);
        while(ch!='N'&&ch!='P')scanf("%c",&ch);
        if(ch=='N')tu[i][j][k]=1;
        else tu[i][j][k]=0;
    }
    for(i=1;i<=x;i++)
    for(j=1;j<=y;j++)
    {
    for(k=1;k<=z;k++)
    {
    lin[k]=lin[k-1]+tu[i][j][k];
    xx[i][j][k]=xx[i][j-1][k]+lin[k];
    }}
        for(i=1;i<=y;i++)
    for(j=1;j<=x;j++)
    {
    for(k=1;k<=z;k++)
    {
    lin[k]=lin[k-1]+tu[j][i][k];
    yy[i][j][k]=yy[i][j-1][k]+lin[k];
    }}
            for(i=1;i<=z;i++)
    for(j=1;j<=x;j++)
    {
    for(k=1;k<=y;k++)
    {
    lin[k]=lin[k-1]+tu[j][k][i];
    zz[i][j][k]=zz[i][j-1][k]+lin[k];
    }}
    for(i=1;i<=x;i++)
    for(j=1;j<=y;j++)
for(a=min(x-i,y-j);a>=0;a--)
{   if(4*(a+1)*z<=ans)break;
qi=0;
    for(k=1;k<=z;k++)
    {
if(4*(z-qi+1)*(a+1)<=ans)break;
if(zz[k][i+a][j+a]-zz[k][i-1][j+a]-zz[k][i+a][j-1]+zz[k][i-1][j-1]==(a+1)*(a+1))
{
    ans=max(ans,4*(a+1)*(k-qi));
}else  qi=k;
    }
}
    for(i=1;i<=y;i++)
    for(j=1;j<=z;j++)
for(a=min(y-i,z-j);a>=0;a--)
{   if(4*(a+1)*x<=ans)break;
qi=0;
    for(k=1;k<=x;k++)
    {
if(4*(x-qi+1)*(a+1)<=ans)break;
if(xx[k][i+a][j+a]-xx[k][i-1][j+a]-xx[k][i+a][j-1]+xx[k][i-1][j-1]==(a+1)*(a+1))
{
    ans=max(ans,4*(a+1)*(k-qi));
}else  qi=k;
    }
}   

    for(i=1;i<=x;i++)
    for(j=1;j<=z;j++)
for(a=min(x-i,z-j);a>=0;a--)
{
    if(4*(a+1)*y<=ans)break;
qi=0;
    for(k=1;k<=y;k++)
    {
if(4*(y-qi+1)*(a+1)<=ans)break;
if(xx[k][i+a][j+a]-xx[k][i-1][j+a]-xx[k][i+a][j-1]+xx[k][i-1][j-1]==(a+1)*(a+1))
{
    ans=max(ans,4*(a+1)*(k-qi));
}else  qi=k;
    }
}
    printf("%d",ans);
}

正解是单调栈。。

首先先选一个为底面，把他分成好几层然后再枚举矩形右下角坐标，，这时候暴力的做法就是枚举正方形大小，然而这样做忽略了

整体满足，部分也满足的性质。所以对于一个正方形，若以n为边长合适，则以1~n为边长也一定合适

所以就在这一个线所在点的矩形最大的n都考虑进来，就相当于

AC代码：

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define N 155
int k,i,j,zuo[N],you[N],f[N][N][N],sta[N],top,ans,x,y,z;
char tu[N][N][N],tu2[N][N][N];
void work(int z,int x,int y)
{for(k=1;k<=z;k++){for(i=1;i<=x;i++){for(j=1;j<=y;j++){if(tu[i][j][k]=='N')f[k][i][j]=min(f[k][i-1][j-1],min(f[k][i-1][j],f[k][i][j-1]))+1;      else f[k][i][j]=0;}        }   }for(j=1;j<=x;j++)for(k=1;k<=y;k++){
memset(zuo,0,sizeof(zuo));
top=0;
sta[0]=0;
for(i=1;i<=z;i++)
{while(top&&f[sta[top]][j][k]>=f[i][j][k])--top;
zuo[i]=i-sta[top];
sta[++top]=i;
}       memset(you,0,sizeof(you));
top=0;
sta[0]=z+1;
for(i=z;i>=1;i--)
{while(top&&f[sta[top]][j][k]>=f[i][j][k])--top;
you[i]=sta[top]-i;
sta[++top]=i;
}
for(i=1;i<=z;i++)
if(tu[j][k][i]=='N')ans=max(ans,(zuo[i]+you[i]-1)*4*f[i][j][k]);
}
}
int main()
{scanf("%d%d%d",&y,&x,&z);for(i=1;i<=x;i++)for(j=1;j<=y;j++)for(k=1;k<=z;k++){scanf("%c",&tu[i][j][k]);while(tu[i][j][k]!='N'&&tu[i][j][k]!='P')scanf("%c",&tu[i][j][k]);tu2[i][j][k]=tu[i][j][k];        }
work(z,x,y);for(i=1;i<=x;i++)for(j=1;j<=y;j++)for(k=1;k<=z;k++){tu[j][k][i]=tu2[i][j][k];     }
work(x,y,z);for(i=1;i<=x;i++)for(j=1;j<=y;j++)for(k=1;k<=z;k++){tu[i][k][j]=tu2[i][j][k];     }
work(y,x,z);    printf("%d",ans);
}

2017.10.12 礼物(zjoi2011) 失败总结相关推荐

2017.10.12 记者招待会
2017.10.12 记者招待会记者招待会看到这个题目,你会不会吓一跳?一名普通的员工,哪里会参加记者招待会?是的,不是我参加.你脆弱的小心脏可以简单的修正一下,因为后面还有一个劲爆的消息:学生 ...
2017.10.12 小Q的无敌异或失败总结
第一问是可以O(n)的二进制位是独立的,所以直接分开算即可.. 记一个前缀和,开桶装0.1的个数然后每次新加入一个点就直接更新就可以了第二问就比较难了,,手玩只能发现最低位可以这么搞,对 ...
2017.9.12 人员雇佣失败总结
网络流好难啊,还被卡题意了(不卡题意我也不会) 然后强套昨天的图,其实样子基本是对的,,但是实际是错的这个题的模型就是: 一个点选s有贡献, 两个点同时选s有贡献, 两个点同时不选s无贡献, 两个点 ...
2017.3.12 分割矩阵失败总结
感觉再也做不出题了 QAQ.. 其实,第一眼看上去我依然是懵逼的, 发现只要有左上角的坐标和左下角的坐标是可以10^4 预处理区间和的简单划了一下只有O((a*b)^n)的做法. 然后发现沿数 ...
2017/10/12 表格与表单
1 2 文本框文本域男女 1234 转载于:https://www.cnblogs.com/buzhiheyan/p/7658754.html
2017.10.23 硬币购物失败总结
这题常规套路根本用不上一共就四枚硬币,所以要考虑对硬币的算法然而谁能想到容斥啊... 首先查询是很多的,需要对每个查询做到尽量常数级首先排除dp查询,因为状态肯定不支持O(1)查询然后需要想到 ...
2017.10.22 方格染色失败总结
详细题解网上有好多.. 一开始想用一个并查集维护确定关系一个并查集维护相等关系但会出现一些无解判不出来(原因未知) 所以就只能写加权并查集,记录每个点和他根节点的相等关系:0相等 1不等然后 ...
2017.10.15 旅行comf 失败总结
越是取值要求奇怪的题,范围不大的题,就越需要枚举. 这个比值的要求非常特殊,他没有任何最优值可以取,所以不能考虑dp及其相似算法所以就要考虑枚举一个求另一个,易知枚举了最小的,加的必须是比他大的,所 ...
2017.10.11 Problem c 失败总结
一开始就想偏了.. 首先计数题 30s + n<=300显然是让你n^3 dp的先从手玩开始考虑注意到从大到小出现的数字个数是可以判断有无解的 ,,然后可以画出一些区间的条件如(> ...

2017.10.12 礼物(zjoi2011) 失败总结

2017.10.12 礼物(zjoi2011) 失败总结相关推荐

最新文章

热门文章