1414: [ZJOI2009]对称的正方形

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 331  Solved: 149
[Submit][Status]

Description

Orez很喜欢搜集一些神秘的数据，并经常把它们排成一个矩阵进行研究。最近，Orez又得到了一些数据，并已经把它们排成了一个n行m列的矩阵。通过观察，Orez发现这些数据蕴涵了一个奇特的数，就是矩阵中上下对称且左右对称的正方形子矩阵的个数。 Orez自然很想知道这个数是多少，可是矩阵太大，无法去数。只能请你编个程序来计算出这个数。

Input

文件的第一行为两个整数n和m。接下来n行每行包含m个正整数，表示Orez得到的矩阵。

Output

文件中仅包含一个整数answer，表示矩阵中有answer个上下左右对称的正方形子矩阵。

Sample Input

5 5
4 2 4 4 4
3 1 4 4 3
3 5 3 3 3
3 1 5 3 3
4 2 1 2 4

Sample Output

27

数据范围
对于30%的数据 n，m≤100
对于100%的数据 n，m≤1000 ，矩阵中的数的大小≤109

 
 
　　現在我確實是遇到難一點的題套單調隊列就搞不清楚，這道題看了以後很容易想到計算每個點向兩個方向延展的迴文串最長長度，然而如何快速統計卻是個問題，邊長2a的目標正方形滿足其行上連續a的範圍內的格子列方向延展長度大於a，且a小於中心橫向迴文串延展長度。（說的我都暈了），解決這個問題，我用的是單調隊列上二分，單調隊列主要是用來快速解決區間最值問題。
　　編的時候出現了三個問題：
  1. manacher算法中忘記更行mx，id變量，導致其退化成O（n^2)，在這個問題調試的時候我確信了manacher算法複雜度是嚴格O(n)的。
  2. 單調隊列中兩相鄰元素所夾區間應屬於後面一個元素。
  3. 這樣得題數組大小一定要準確，如果所有下標範圍都開大一倍一定會MLE
#include<iostream>

#include<cstdio>

#include<cstdlib>

#include<cstring>

#include<ctime>

#include<cmath>

#include<algorithm>

#include<set>

#include<map>

#include<vector>

#include<string>

#include<queue>

using namespace std;

#ifdef WIN32

#define LL "%I64d"

#else

#define LL "%lld"

#endif

#define MAXN 2100

#define MAXV MAXN*2

#define MAXE MAXV*2

#define INF 0x3f3f3f3f

#define INFL 0x3f3f3f3f3f3f3f3fLL

typedef long long qword;

inline int nextInt()

{

        char ch;

        int x=;

        bool flag=false;

        do

                ch=(char)getchar(),flag=(ch=='-')?true:flag;

        while(ch<''||ch>'');

        do x=x*+ch-'';

        while (ch=(char)getchar(),ch<='' && ch>='');

        return x*(flag?-:);

}

int n,m;

typedef int arr_t[MAXN];

typedef int arr2[MAXN*];

typedef arr_t map_t[MAXN];

map_t a,b;

map_t sa,sb,tsb,tsa;

void manacher(arr2 seq,int n,arr2 &res)

{

        int id,mx;

        int i,j;

        id=mx=-;

        int cnt=;

        for (i=;i<n;i++)

        {

                if (i<mx)

                {

                        res[i]=min(mx-i,res[id*-i]);

                }else res[i]=;

                while (i-res[i]->= && i+res[i]+<n && seq[i+res[i]+]==seq[i-res[i]-])

                {

                        res[i]++;

                        //cnt++;

                }

                if (i+res[i]>mx)//Do not forget

                {

                        mx=i+res[i];

                        id=i;

                }

        }

        return ;

}

void init_p(map_t mp,int n,int m,map_t &sl)

{

        int i,j;

        arr2 seq,res;

        for (i=;i<n;i++)

        {

                for (j=;j<m;j++)

                {

                        seq[j*+]=mp[i][j];

                        seq[j*]=-INF;

                }

                seq[*m]=-INF;

                manacher(seq,*m+,res);

                for (j=;j<m;j++)

                {

                        sl[i][j]=res[j*+];

                }

        }

}

int pseq[MAXN];

int head,tail;

map_t fa,fb,fc;

void pm(map_t &a)

{

        int i,j;

        for (i=;i<n;i++)

        {

                for (j=;j<m;j++)

                {

                        printf("%d ",a[i][j]);

                }

                printf("\n");

        }

        printf("\n");

}

void work(map_t &sa,map_t &sb,map_t &fa,int n,int m)

{

        int i,j;

        for (i=;i<n;i++)

                for (j=;j<m/;j++)

                {

                        swap(sb[i][j],sb[i][m-j-]);

                        swap(sa[i][j],sa[i][m-j-]);

                }

        for (i=;i<n;i++)

        {

                head=,tail=;

                pseq[]=-;

                pseq[]=;

                fa[i][m-]=;

                for (j=;j<m;j++)

                {

                        while (tail>=head && sb[i][j]<sb[i][pseq[tail]])tail--;

                        pseq[++tail]=j;

                        int l,r,mid;

                        l=-,r=tail;

                        while (l+<r)

                        {

                                int mid=(l+r)>>;

                                if (sb[i][pseq[mid]]>=j-pseq[mid])

                                        r=mid;

                                else

                                        l=mid;

                        }

                        fa[i][m-j-]=min(j-(pseq[l]+),min(sa[i][j],sb[i][pseq[r]]));

                }

        }

        for (i=;i<n;i++)

                for (j=;j<m/;j++)

                {

                        swap(sb[i][j],sb[i][m-j-]);

                        swap(sa[i][j],sa[i][m-j-]);

                }

        for (i=;i<n;i++)

        {

                head=,tail=;

                pseq[]=-;

                pseq[]=;

                fa[i][]=;

                for (j=;j<m;j++)

                {

                        while (tail>=head && sb[i][j]<sb[i][pseq[tail]])tail--;

                        pseq[++tail]=j;

                        int l,r,mid;

                        l=-,r=tail;

                        while (l+<r)

                        {

                                int mid=(l+r)>>;

                                if (sb[i][pseq[mid]]>=j-pseq[mid])

                                        r=mid;

                                else

                                        l=mid;

                        }

                        fa[i][j]=min(fa[i][j],min(j-(pseq[l]+),min(sa[i][j],sb[i][pseq[r]])));

                }

        }

}

int main()

{

        //freopen("input.txt","r",stdin);

        //freopen("output.txt","w",stdout);

        int i,j,k;

        int x,y,z;

        int ans=;

        scanf("%d%d",&n,&m);

        memset(a,INF,sizeof(a));

        for (i=;i<n;i++)

        {

                for (j=;j<m;j++)

                {

                        scanf("%d",&x);

                        a[i*+][j*+]=x;

                        b[j*+][i*+]=x;

                }

        }

        n*=;m*=;n++;m++;

        init_p(a,n,m,sa);

        init_p(b,m,n,tsb);

        for (i=;i<n;i++)

                for (j=;j<m;j++)

                        sb[i][j]=tsb[j][i];

        for (i=;i<n;i++)

        {

                for (j=;j<m;j++)

                {

                        sa[i][j]/=;

                        sb[i][j]/=;

                }

        }    

        work(sa,sb,fa,n,m);

        for (i=;i<n;i++)

                for (j=;j<m;j++)

                        tsa[i][j]=a[i][j];

        for (i=;i<n;i++)

                for (j=;j<m;j++)

                        b[i][m-j-]=a[m-j-][i]=tsa[i][j];

        swap(n,m);

        init_p(a,n,m,sa);

        init_p(b,m,n,tsb);

        for (i=;i<n;i++)

                for (j=;j<m;j++)

                        sb[i][j]=tsb[j][i];

        for (i=;i<n;i++)

        {

                for (j=;j<m;j++)

                {

                        sa[i][j]/=;

                        sb[i][j]/=;

                }

        }

        work(sa,sb,fb,n,m);

        for (i=;i<n;i++)

                for (j=;j<m;j++)

                        tsb[i][j]=fb[i][j];

        swap(n,m);

        for (i=;i<n;i++)

                for (j=;j<m;j++)

                        fb[i][j]=tsb[m-j-][i];

        //pm(fa);

        //pm(fb);

        for (i=;i<n;i++)

        {

                for (j=;j<m;j++)

                {

                        fc[i][j]=max(,min(fa[i][j]-(i%==),fb[i][j]-(j%==)));

                        fc[i][j]=(fc[i][j]+)/;

                }

        }

        //pm(fc);

        for (i=;i<n;i++)

        {

                for (j=;j<m;j++)

                {

                        if (i%+j%==)continue;

                        ans+=fc[i][j];

                }

        }

        cout<<ans<<endl;

}

bzoj 1414: [ZJOI2009]对称的正方形 manacher算法+單調隊列的更多相关文章

  1. bzoj 1414: [ZJOI2009]对称的正方形

    Description Orez很喜欢搜集一些神秘的数据,并经常把它们排成一个矩阵进行研究.最近,Orez又得到了一些数据,并已经把它们排成了一个n行m列的矩阵.通过观察,Orez发现这些数据蕴涵了一 ...

  2. 【bzoj 1414】对称的正方形 单调队列+manacher

    Description Orez很喜欢搜集一些神秘的数据,并经常把它们排成一个矩阵进行研究.最近,Orez又得到了一些数据,并已经把它们排成了一个n行m列的矩阵.通过观察,Orez发现这些数据蕴涵了一 ...

  3. 【BZOJ1414】[ZJOI2009]对称的正方形（哈希）

    [BZOJ1414][ZJOI2009]对称的正方形(哈希) 题面 BZOJ 洛谷 题解 深思熟虑一波,发现一个矩阵如果左右对称的话,那么它每行都是一个回文串,同理,如果上下对称的话,那么每列都是一个 ...

  4. 【BZOJ1414/3705】[ZJOI2009]对称的正方形 二分+hash

    [BZOJ1414/3705][ZJOI2009]对称的正方形 Description Orez很喜欢搜集一些神秘的数据,并经常把它们排成一个矩阵进行研究.最近,Orez又得到了一些数据,并已经把它们 ...

  5. BZOJ 2342 双倍回文（manacher算法）

    题目链接:http://61.187.179.132/JudgeOnline/problem.php?id=2342 题意:定义双倍回文串为:串的长度为4的倍数且串的前一半.后一半.串本身均是回文的. ...

  6. [luoguP2601] [ZJOI2009]对称的正方形（二维Hash + 二分 || Manacher）

    传送门 很蒙蔽,不知道怎么搞. 网上看题解有说可以哈希+二分搞,也有的人说用Manacher搞,Manacher是什么鬼?以后再学. 对于这个题,可以从矩阵4个角hash一遍,然后枚举矩阵中的点,再二 ...

  7. 题解-------[ZJOI2009]对称的正方形

    传送门 题目大意 找到所有的上下左右都相同的正方形. 思路:二分+二维Hash 这道题我们首先想到不能暴力判断一个正方形是否合法. 然后我们发现当一个正方形合法时,以这个正方形为中心且比它小的正方形也 ...

  8. manacher算法学习(求最长回文子串长度)

    Manacher总结 我的代码 学习:yyb luogu题目模板 xzy的模板 #include<iostream> #include<cstdlib> #include< ...

  9. 【学习笔记】Manacher算法

    本文部分图片来源 代码来源(代码是学姐哒.. 一.引入 Manacher算法是用来求最长回文子串的算法,时间复杂度O(n). 回文子串指的是''aacaa'',''noon'',这种正着反着读都一样的 ...

随机推荐

  1. OpenStack_Swift源代码分析——创建Ring及加入�设备源代码算法具体分析

    1 创建Ring 代码具体分析 在OpenStack_Swift--Ring组织架构中我们具体分析了Ring的具体工作过程,以下就Ring中添加�设备,删除设备,已经又一次平衡的实现过程作具体的介绍. ...

  2. python学习笔记--Django入门四 管理站点

    上一节  Django入门三 Django 与数据库的交互:数据建模 "管理员界面"是基础功能中的重要部分. django.contrib 包 Django自动管理工具是djang ...

  3. angularjs 根据变量改变 动态加载模板

    directive return { restrict: 'E', replace: true, templateUrl: 'app/view/order.html', link: function ...

  4. 在后台CS文件里面,隐藏和显示Repeater里面控件

    <asp:Repeater ID="Repeater1" runat="server"><ItemTemplate><asp:Pa ...

  5. PHP负载均衡

    待更新 版权声明:本文为博主原创文章,未经博主允许不得转载.

  6. 原生js日期时间插件鼠标点击文本框弹出日期时间表格选择日期时间

    原文出处 (这是我从互联网上搜来的,感觉能满足各方面的需求.个人感觉挺不错的,所以后期修改了一下向大家推荐!) 效果图: html代码: <!DOCTYPE html PUBLIC " ...

  7. Bootstrap-全局css样式之按钮

    这里所说的按钮只是Bootstrap设计的能使标签或元素呈现按钮样式的属性,所以为 <a>.<button> 或 <input> 元素添加按钮类(button cl ...

  8. MongoDB的安装和基本操作

    一.使用前的准备(windows下的安装)  1.下载 目前MongoDB的官网不知道问什么不能进行下载了,但是可以在MongoDB中文论坛进行下载, 地址如下:http://www.mongoing ...

  9. js - get-the-value-from-the-url-parameter(可以在非模态对话框中使用)

    ref: http://stackoverflow.com/questions/979975/how-to-get-the-value-from-the-url-parameter 函数: funct ...

  10. 关于for循环中的闭包问题

    还是昨天的那个简单的小项目,已经花了一天的时间了 - - .从&&的用法,到CSStext,到今天马上要谈的闭包(closure),通过一个小东西,真真发现了自己的各方面不足.昨天发完 ...