Palindrome

简要题意: 

我们有一个字符串S,字符串的长度不超过500000。

求满足S[i]=S[2n−i]=S[2n+i−2](1≤i≤n)(n≥2)的子串个数。 

分析:

我们能通过简单的数学知识,得到:

该子串是两个回文串拼在一起的,例如abcbabc中,前5项为一个回文串,后5项有一个回文串。

第n项以及第2*n-1项为回文串的中心。

我们可以用Manacher 求得以每个$i$为中心点的回文串半径len[i]。

求得len[i]后,写一个方程:

令i<j,
1. i+len[i]>=j;
2. j-len[j]<=i;
可以推得:
1. i+1<=j<=i+len[i];
2. j-len[j]<=i

所以,这道题便转化成:

已知一个len[i]数组,求满足上列条件的个数! 

怎么求一个区间内小于一个数x的个数呢?

主席树! 

主席树中记录的是每个数字出现过的个数,区间[1,5]即记录数字[1,5]出现的总次数,

我们主席树类似前缀和的思想,依次往后合并。

最终我们求得的答案为:∑query(root[i+len[i]],x)-query(root[i],x);

#include<bits/stdc++.h>
using namespace std;
#define re register int
#define LL long long
const int N=1000005;
int seg,root[N];
struct segment{int v,ls,rs;}tr[N<<4];
int build(int x,int y)
{
int p=++seg;tr[p].v=0;
if(x<y)
{
int mid=(x+y)>>1;
tr[p].ls=build(x,mid);
tr[p].rs=build(mid+1,y);
}
return p;
}
int update(int pre,int l,int r,int x)
{
int p=++seg;
tr[p]=(segment){tr[pre].v+1,tr[pre].ls,tr[pre].rs};
if(l<r)
{
int mid=(l+r)>>1;
if(x<=mid)tr[p].ls=update(tr[pre].ls,l,mid,x);
else tr[p].rs=update(tr[pre].rs,mid+1,r,x);
}
return p;
}
int query(int p1,int p2,int l,int r,int x)
{
if(x<l)return 0;
if(r<=x)return tr[p1].v-tr[p2].v;
int mid=(l+r)>>1;
int res=query(tr[p1].ls,tr[p2].ls,l,mid,x);
if(mid+1<=x)res+=query(tr[p1].rs,tr[p2].rs,mid+1,r,x);
return res;
}
char S[N],s[N<<1];
int val[N<<1],pos[N];
inline void work()
{
scanf("%s",S);int n=strlen(S);
int k=0,id=0,mx=0;
for(re i=0;i<n;++i)s[++k]='#',s[++k]=S[i];s[0]='!';s[++k]='#';
for(re i=1;i<=k;++i)
{
if(i<mx)val[i]=min(val[id*2-i],mx-i);
else val[i]=1;
while(s[i-val[i]]==s[i+val[i]])val[i]++;
if(i+val[i]>mx)mx=i+val[i],id=i;
}
int MAX=0;
for(re i=2;i<=k;i+=2)
{
int id=(i>>1);
pos[id]=(val[i]>>1)-1;
MAX=max(MAX,id-pos[id]);
} seg=0;root[0]=build(1,MAX);
for(re i=1;i<=n;++i)root[i]=update(root[i-1],1,MAX,i-pos[i]);
LL ans=0;
for(re i=1;i<=n;++i)
if(pos[i])ans+=1ll*query(root[i+pos[i]],root[i],1,MAX,i);
printf("%lld\n",ans);
}
signed main()
{
int T;scanf("%d",&T);
while(T--)work();
}

“向前跑,迎着冷眼和嘲笑”

“失败后郁郁寡欢,那是懦夫的表现”----《追梦赤子心》

热身训练3 Palindrome的更多相关文章

  1. 数位dp & 热身训练7

    数位dp 数位dp是一种计数用的dp,一般就是要统计一段区间$[L,R]$内,满足一定条件的数的个数,或者各个数位的个数. 数位dp使得暴力枚举变为满足一定状态的记忆化,更加优秀. 数位dp常常会考虑 ...

  2. 热身训练4 Article

    Article 在这个学期即将结束时,DRD开始写他的最后一篇文章. DRD使用著名的Macrohard的软件World来写他的文章. 不幸的是,这个软件相当不稳定,它总是崩溃. DRD需要在他的文章 ...

  3. 热身训练4 Eighty seven

    Eighty seven 简要题意: n个卡片,其中第i个卡片的数值为$a[i]$.一共q次询问,每次询问将删去其中3个卡片(可能删除若干相同的卡片)后,问能否选出10个卡片,数值之和等于87. n≤ ...

  4. 热身训练2 The All-purpose Zero

    The All-purpose Zero 简要题意:  长度为n的数组,每个数字为S[i],$0$是一种很神奇的数字,你想要的,它都可以变! 问这个序列的最长上升子序列长度为多少? 分析: 我们将除了 ...

  5. 热身训练2 GCD

    题目描述 简要题意:  n个数字,a1,a2,...,an m次询问(l,r),每次询问需回答 1.gcd(al,al+1,al+2,...,ar);2.gcd(ax,ax+1,ax+2,...,ay ...

  6. 热身训练2 Another Meaning

    题目来源 简要题意: 众所周知,在许多情况下,一个词语有两种意思.比如"hehe",不仅意味着"hehe",还意味着"excuse me". ...

  7. 热身训练1 Calculator

    题目出处:Calculator 简要题意: 你有一个确定的函数,f(x)=+...*...^...,其中共有n个操作,从左到右依次计算. 共有m次询问,我们每次询问,1.会修改f(x)中的操作:2.输 ...

  8. 热身训练1 ping ping ping

    点此进入 题意: 一棵树,n+1 个节点,以0号节点为根,给出端点(a,b),节点a到节点b的路径上,至少有一个点是"坏掉的",求"坏掉的点"最少 分析: St ...

  9. 热身训练1 Problem B. Harvest of Apples

    http://acm.hdu.edu.cn/showproblem.php?pid=6333 题意: 求 C(0,n)+C(1,n)+...+C(m,n) 分析: 这道题,我们令s(m,n) = C( ...

随机推荐

  1. CSS导航菜单(二级菜单)

    index.html <div class="nav"> <ul> <li> <a href="#">Java& ...

  2. 计算机网络-TCP篇

    TCP篇 之前的总结文章:TCP简单版本介绍-三次握手等 基本认识 TCP 是⾯向连接的(⼀定是「⼀对⼀」才能连接).可靠的.基于字节流的传输层通信协议. RFC 793 是如何定义「连接」的:⽤于保 ...

  3. go build 与go install

    相同点都能生成可执行文件 不同点go build 不能生成包文件, go install 可以生成包文件go build 生成可执行文件在当前目录下, go install 生成可执行文件在bin目录 ...

  4. Spring Cloud Eureka 之服务端自我注册

    Eureka服务端实现了一种自我注册机制,涉及配置项: eureka.client.register-with-eureka spring.application.name Eureka Server ...

  5. C++快速读入

    使用C++的标准cin进行读入速度比较慢,尤其是在大数据的情况下,所以我们需要使用一种方法,按照字符读入,最后再"组装"成整数.由于字符读入比数字要快,所以这样做可以提高读入速度. ...

  6. 猪齿鱼 SaaS 版效能平台发布

    ​日前,猪齿鱼Choerodon全场景效能平台Saas版发布,提供体系化方法论和协作.测试.DevOps及容器工具,帮助企业拉通需求.设计.开发.部署.测试和运营流程,一站式提高管理效率和质量.从团队 ...

  7. Selenium多浏览器处理 (Chrome/Firefox/IE)

    测试用例文件:test_selenium/test_hogwarts.py 使用pytest框架 定义一个变量,通过外部传入变量,确定使用哪个浏览器 browser = os.getenv(" ...

  8. python读取文件编码转换问题

    encode(编码)   decode(解码)   encoding(编码格式) #-*- coding:utf-8 -*- import chardet #用于查看编码 with open(&quo ...

  9. CF1251F-Red-White Fence【NTT】

    前言 刚开始看错题推了半天的生成函数 正题 题目链接:https://www.luogu.com.cn/problem/CF1251F 题目大意 $n$个白色木板,$k$个红色木板,给出这些木板的高度 ...

  10. P5363-[SDOI2019]移动金币【阶梯博弈,dp,组合数学】

    正题 题目链接:https://www.luogu.com.cn/problem/P5363 题目大意 \(1\times n\)的网格上有\(m\)个硬币,两个人轮流向前移动一个硬币但是不能超过前一 ...