因为没有简化题意一直没去做,直到今天讲这道题才口胡出来

要求对称,很明显这样一个“子序列”的对称中心只可能有一个,那么先枚举这个对称中心。

然后我们需要判断两个位置是否相同。看上去好像很困难。

考虑设计哈希函数 \(f(x,y)\),使得 \(f(V_a,V_b)\neq 0\) 且 \(f(V_a,V_a)=f(V_b,V_b)=0\)。

首先等于 \(0\),那么一定有一个 \((x-y)\)。其次,注意到字符集只有 \(2\),为了方便,我们令 \(f(x,y)=-f(y,x)\)。那么就包含了 \((x-y)^2\)。

所以 \(f(x,y)=(V_x-V_y)^2=V_x^2+V_y^2-2V_xV_y\),这就够了。

我们需要计算的是 \(match[i]=\sum_{j=0}f(i-j,i+j)\),这有点像卷积。

拆开得到:

\[\sum_{j=0}(V_{i-j}^2+V_{i+j}^2)-2V_{i+j}V_{i-j}
\]

后面这玩意儿可以用前缀和搞定,前面这玩意儿做一个卷积就好啦。

这样做会算上连续的一段,使用 manacher 计算出来并且减掉就好啦。

复杂度 \(O(n\log n)\)。

#include<cstring>
#include<cstdio>
typedef unsigned ui;
const ui M=1e5+5,G=3,mod=998244353,MOD=1e9+7;
ui n,V[M],S[M],pw2[M],f[M<<2];ui ans;char s[M];
ui p[M<<1];char c[M<<1];
inline ui min(const ui&a,const ui&b){
return a>b?b:a;
}
inline ui pow(ui a,ui b=mod-2){
ui ans(1);for(;b;b>>=1,a=1ull*a*a%mod)if(b&1)ans=1ull*ans*a%mod;return ans;
}
inline ui Add(const ui&a,const ui&b){
return a+b>=mod?a+b-mod:a+b;
}
inline ui Del(const ui&a,const ui&b){
return b>a?a-b+mod:a-b;
}
inline void swap(ui&a,ui&b){
ui c=a;a=b;b=c;
}
inline void DFT(ui*f,const ui&n){
for(ui len=n>>1;len>=1;len>>=1){
const ui w1=pow(G,(mod-1>>1)/len);
for(ui k=0;k<n;k+=len<<1){
for(ui w(1),i=0;i<len;++i){
const ui x=f[i|k],y=f[i|k|len];
f[i|k]=Add(x,y);f[i|k|len]=1ull*w*Del(x,y)%mod;
w=1ull*w*w1%mod;
}
}
}
}
inline void IDFT(ui*f,const ui&n){
for(ui len=1;len<n;len<<=1){
const ui w1=pow(G,(mod-1>>1)/len);
for(ui k=0;k<n;k+=len<<1){
for(ui w(1),i=0;i<len;++i){
const ui x=f[i|k],y=1ull*w*f[i|k|len]%mod;
f[i|k]=Add(x,y);f[i|k|len]=Del(x,y);
w=1ull*w*w1%mod;
}
}
}
const ui inv=pow(n);
for(ui i=0;i<n;++i)f[i]=1ull*f[i]*inv%mod;
for(ui i=1;(i<<1)<n;++i)swap(f[i],f[n-i]);
}
inline ui manacher(char*s){
ui R(0),mid(0),ans(0);c[0]='^';c[1]='#';
for(ui i=1;i<=n;++i)c[i<<1]=s[i],c[i<<1|1]='#';
for(ui i=1;i<=(n<<1);++i){
p[i]=i<=R?min(p[(mid<<1)-i],p[mid]+mid-i):1;
while(c[i-p[i]]==c[i+p[i]])++p[i];
if(i+p[i]-1>R)mid=i,R=i+p[i]-1;
ans=(ans+(p[i]>>1))%MOD;
}
return ans;
}
signed main(){
scanf("%s",s+1);n=strlen(s+1);pw2[0]=1;
for(ui i=1;i<=n;++i){
V[i]=s[i]=='a'?1:2;pw2[i]=pw2[i-1]*2%MOD;
f[i]=V[i];S[i]=S[i-1]+V[i]*V[i];
}
ui len(1);
while(len<n+n+1)len<<=1;
DFT(f,len);
for(ui i=0;i<len;++i)f[i]=1ull*f[i]*f[i]%mod;
IDFT(f,len);
for(ui i=1;i<=n;++i){
ui len=min(n-i,i-1);
ans=(ans+pw2[len+1-((S[i+len]-S[i-len-1])-f[i<<1])]-1)%MOD;
if(len!=n-i)++len;
ans=(ans+pw2[len-((S[i+len]-S[i-len])-f[i<<1|1])]-1)%MOD;
}
printf("%u",(MOD+ans-manacher(s))%MOD);
}

LGP4199题解的更多相关文章

  1. 2016 华南师大ACM校赛 SCNUCPC 非官方题解

    我要举报本次校赛出题人的消极出题!!! 官方题解请戳:http://3.scnuacm2015.sinaapp.com/?p=89(其实就是一堆代码没有题解) A. 树链剖分数据结构板题 题目大意:我 ...

  2. noip2016十连测题解

    以下代码为了阅读方便,省去以下头文件: #include <iostream> #include <stdio.h> #include <math.h> #incl ...

  3. BZOJ-2561-最小生成树 题解(最小割)

    2561: 最小生成树(题解) Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1628  Solved: 786 传送门:http://www.lyd ...

  4. Codeforces Round #353 (Div. 2) ABCDE 题解 python

    Problems     # Name     A Infinite Sequence standard input/output 1 s, 256 MB    x3509 B Restoring P ...

  5. 哈尔滨理工大学ACM全国邀请赛(网络同步赛)题解

    题目链接 提交连接:http://acm-software.hrbust.edu.cn/problemset.php?page=5 1470-1482 只做出来四道比较水的题目,还需要加强中等题的训练 ...

  6. 2016ACM青岛区域赛题解

    A.Relic Discovery_hdu5982 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Jav ...

  7. poj1399 hoj1037 Direct Visibility 题解 (宽搜)

    http://poj.org/problem?id=1399 http://acm.hit.edu.cn/hoj/problem/view?id=1037 题意: 在一个最多200*200的minec ...

  8. 网络流n题 题解

    学会了网络流,就经常闲的没事儿刷网络流--于是乎来一发题解. 1. COGS2093 花园的守护之神 题意:给定一个带权无向图,问至少删除多少条边才能使得s-t最短路的长度变长. 用Dijkstra或 ...

  9. CF100965C题解..

    求方程 \[ \begin{array}\\ \sum_{i=1}^n x_i & \equiv & a_1 \pmod{p} \\ \sum_{i=1}^n x_i^2 & ...

随机推荐

  1. indexOf() 和 lastIndexOf()

    一,定义和用法 indexOf() 方法可返回某个指定的字符串值在字符串中首次出现的位置. lastIndexOf() 方法可返回一个指定的字符串值最后出现的位置,在一个字符串中的指定位置从后向前搜索 ...

  2. jdk1.5出现的新特性---->增强for循环

    import java.util.HashMap; import java.util.HashSet; import java.util.Iterator; import java.util.Map; ...

  3. android中listView下拉刷新

    Android的ListView是应用最广的一个组件,功能强大,扩展性灵活(不局限于ListView本身一个类),前面的文章有介绍分组,拖拽,3D立体,游标,圆角,而今天我们要介绍的是另外一个扩展Li ...

  4. Nginx+Tomcat负载均衡、动静分离群集

    Nginx+Tomcat负载均衡.动静分离群集 目录 Nginx+Tomcat负载均衡.动静分离群集 一.Tomcat 1. Tomcat简介 2. Tomcat重要目录 二.Nginx负载均衡原理 ...

  5. 大话devops

    一.敏捷的局限性的促使devops诞生 敏捷的局限性:敏捷只注重开发阶段的敏捷,未涉及到整个产品生命周期流程其他环节导致采用敏捷开发流程后效果不明显. devops成为企业数字化转型的助推器,扮演基础 ...

  6. Solution -「SV 2020 Round I」「SRM 551 DIV1」「TC 12141」SweetFruits

    \(\mathcal{Description}\)   link.   给定 \(n\) 个水果,每个结点可能有甜度 \(v_i\),或不甜(\(v_i=-1\)).现在把这些水果串成一棵无根树.称一 ...

  7. Solution -「CF 510E」Fox And Dinner

    \(\mathcal{Description}\)   Link.   给定正整数集合 \(\{a_n\}\),求一种把这些数放置在任意多个圆环上的方案,使得每个环的大小大于 \(2\) 且环上相邻两 ...

  8. Process Doppelg&#228;nging

    进程注入:Process Doppelgänging   攻击者可以通过Process Doppelgänging将恶意代码注入到进程中,从而逃避基于进程的防护,并且进行可能的特权提升.Process ...

  9. c++ 子类与父类之间的类型转换

    子类与父类之间的类型转换 先给一段代码 class Base { public: int a = 10; }; class pub_Derv : public Base { Base *getBase ...

  10. 【曹工杂谈】Mysql-Connector-Java时区问题的一点理解--写入数据库的时间总是晚13小时问题

    背景 去年写了一篇"[曹工杂谈]Mysql客户端上,时间为啥和本地差了整整13个小时,就离谱",结果最近还真就用上了. 不是我用上,是组内一位同事,他也是这样:有个服务往数据库in ...