这是一道探究性质的题目,但是思路很清晰明了
性质:对于一个满足要求的子矩形,他的折叠方式的顺序是可以任意改变的
证明:
其他的情况也可以类推了
因此我们只需要考虑这个大矩形横着划分的方式*竖着划分的方式就是答案了
然后对于每一个划分的位置,都是可以合并的
如果翻折时越过了除大矩形边界外的边界是没有影响的
可以证明此时一定可以先自行翻折成较短的长度
剩下的就是简单的MANACHER了
#include<bits/stdc++.h>
using namespace std;
const int mod=1000000007;
char s[2000005];
int n,m,col[2000005],row[2000005];
int b[2000005],p[2000005];
void manacher(int *a,int n)
{
for(int i=1;i<=n;i++)
b[i*2]=a[i];
for(int i=0;i<=n;i++)
b[i*2+1]=-1;
b[0]=-2;b[n*2+2]=-3;
int id=0,maxid=0;
for(int i=1;i<=n*2+1;i++)
{
p[i]=maxid>i?min(maxid-i,p[id+id-i]):1;
while(b[i+p[i]]==b[i-p[i]]) p[i]++;
if(i+p[i]>maxid) maxid=i+p[i],id=i;
}
}
int l[2000005],r[2000005];
long long sl[2000005],sr[2000005];
long long calc(int *a,int n)
{
manacher(a,n);
for(int i=0;i<=n;i++)
{
p[i]=p[i*2+1]>>1;
}
l[0]=1;sl[0]=1;
for(int i=1;i<=n;i++)
{
l[i]=(sl[i-1]-(i-p[i]-1>=0?sl[i-p[i]-1]:0))>0;
sl[i]=sl[i-1]+l[i];
}
r[n]=1;sr[n]=1;
for(int i=n-1;i>=0;i--)
{
r[i]=(sr[i+1]-(i+p[i]+1<=n?sr[i+p[i]+1]:0))>0;
sr[i]=sr[i+1]+r[i];
}
long long ret=0;
for(int i=1;i<=n;i++)
{
ret+=sl[i-1]*r[i];
}
return ret;
}
int main()
{
cin>>n>>m;
for(int i=1;i<=n;i++)
{
scanf("%s",s+1);
for(int j=1;j<=m;j++)
{
row[i]=(29ll*row[i]+s[j])%mod;
col[j]=(29ll*col[j]+s[j])%mod;
}
}
printf("%lld\n",calc(row,n)*calc(col,m));
return 0;
}