[回文自动机 Manacher] BZOJ4166: 月宫的符卡序列

news/2025/3/14 16:39:58/

hash被卡…
本来以为是回文自动机裸题
发现fail树上一条链的节点表示的回文子串的中点是不一样的…
不过回文树上的链是一样的

那么用建出回文树(我用回文自动机建的,manacher建不知道为什么WA了),然后找到以每个点为中点的最大回文子串,这个用manacher找

在对应节点加上贡献就行了

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>using namespace std;typedef unsigned long long ull;
typedef long long ll;const int N=1000010,P1=19980403,P2=998244353,base1=131,base2=129;int t,n,p,cnt=1,nxt[N][30],fa[N],len[N],fail[N],val[N];
char a[N];ull hs1[N],pw1[N],hs2[N],pw2[N];struct Hst{int G[P1+10],nxt[N<<1],pos[N<<1],cnt;ull val[N<<1];void clear(){ memset(G,0,sizeof(G)); cnt=0; }void insert(ull x,int p){ull u=x>>32,v=x-(u<<32);//for(int i=G[u];i;i=nxt[i])//  if(val[i]==x) return ;++cnt; val[cnt]=v; nxt[cnt]=G[u]; G[u]=cnt; pos[cnt]=p;}int find(ull x){ull u=x>>32,v=x-(u<<32);for(int i=G[u];i;i=nxt[i])if(val[i]==v) return pos[i];}
}H;inline ull Getsub(int l,int r){ull a=(hs1[r]+P1-hs1[l-1]*pw1[r-l+1]%P1)%P1;ull b=(hs2[r]+P2-hs2[l-1]*pw2[r-l+1]%P2)%P2;return a<<32|b;
}inline void extend(int c,int ps){while(a[ps-len[p]-1]!=a[ps]) p=fail[p];if(!nxt[p][c]){int np=++cnt,k=fail[p]; len[np]=len[p]+2;while(a[ps-len[k]-1]!=a[ps]) k=fail[k];fail[np]=nxt[k][c]; fa[np]=p;nxt[p][c]=np;H.insert(Getsub(ps-len[np]+1,ps),np);}p=nxt[p][c];
}char b[N<<1];int r[N<<1];inline void manacher(){int t=0;for(int i=1;i<=n;i++)b[++t]=a[i],b[++t]='%';int mxr=0,p=0;for(int i=1;i<=t;i++){int j=0;if(i<=mxr)j=min(i+r[2*p-i],mxr+1);j=max(j,i+1);while(2*i-j>0 && j<=t && b[2*i-j]==b[j]) j++;r[i]=j-i;int L=i-r[i]+1,R=i+r[i]-1;if(b[L]=='%') L++,R--;if(L<=R) val[H.find(Getsub(L+1>>1,R+1>>1))]^=(i+1>>1)-1;if(i+r[i]-1>mxr) mxr=i+r[i]-1,p=i;}
}inline char nc(){static char buf[100000],*p1=buf,*p2=buf;return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}inline void read(int &x){char c=nc(); x=0;for(;c>'9'||c<'0';c=nc());for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc());
}inline int read(char *a){char c=nc(); int x=0;for(;c>'z'||c<'a';c=nc());for(;c>='a'&&c<='z';a[++x]=c,c=nc()); return x;
}int main(){read(t);while(t--){n=read(a); cnt=1; p=0;memset(nxt,0,sizeof(nxt));memset(val,0,sizeof(val));H.clear();len[1]=-1; fail[0]=fail[1]=1; pw1[0]=pw2[0]=1;for(int i=1;i<=n;i++){hs1[i]=(1LL*hs1[i-1]*base1+a[i])%P1,pw1[i]=1LL*pw1[i-1]*base1%P1;hs2[i]=(1LL*hs2[i-1]*base2+a[i])%P2,pw2[i]=1LL*pw2[i-1]*base2%P2;}for(int i=1;i<=n;i++) extend(a[i]-'a',i);int ans=0;manacher();for(int i=cnt;i;i--) val[fa[i]]^=val[i];for(int i=2;i<=cnt;i++) ans=max(ans,val[i]);printf("%d\n",ans);}return 0;
}

http://www.ppmy.cn/news/243415.html

相关文章

【BZOJ4166】月宫的符卡序列 Manacher+hash

【BZOJ4166】月宫的符卡序列 题解&#xff1a;题倒不难&#xff0c;就是有点恶心。 首先学习回文串的时候一定学到了这样一个结论&#xff1a;一个长度为n的串的本质不同的回文子串数量不超过n个。 那么我们就可以试图将所有回文串的价值都计算出来&#xff0c;这就需要我们先计…

洛谷P4166 [SCOI2007]最大土地面积

将四边形拆成两个三角形。旋转卡壳经典题。 #include <bits/stdc.h> using namespace std; typedef long long LL; typedef int lint; const int maxn 2001; const double eps 1e-12; const double PI acos(-1.0); int sgn(double x) {if(fabs(x) < eps)return 0;…

BZOJ 1069 Luogu P4166 最大土地面积 (凸包)

题目链接: (bzoj)https://www.lydsy.com/JudgeOnline/problem.php?id1069 (luogu)https://www.luogu.org/problemnew/show/P4166 题解: 水题&#xff0c;凸包极角排序之后枚举凸四边形对角线\(i,j\)然后找面积最大的点\(k\)&#xff0c;\(k\)随着\(i,j\)是单调的 但是有个易错…

Android中的Intent(显示隐式)

Android中的Intent(显示&隐式) 显示Intent 显示Intent是明确目标Activity的类名 通过Intent(Context packageContext, Class<?> cls)构造方法Intent intent new Intent(this, SecondActivity.class) startActivity(intent);通过Intent的setComponent()方法Componen…

NLP的idea,看了就能水一篇论文

1.问题 在中文情感分析任务中,已有方法仅从单极、单尺度来考虑情感特征&#xff0c;无法充分挖掘和利用情感特征信息&#xff0c;模型性能不理想。 单级单尺度&#xff1a;只从一个方面学习文本的特征 多级多尺度&#xff1a;应该是分别从不同方面学习文本的特征&#xff0c…

一键提升分辨率,呈现更清晰、更细腻的视觉盛宴

牛学长视频修复工具是一款使用AI技术对你的视频分辨率就行调整放大清晰的软件&#xff0c;最高支持8K超高清。 现如今&#xff0c;视频成为人们记录生活、表达创意的重要方式之一。然而&#xff0c;我们常常遇到一个问题&#xff1a;旧有的视频素材分辨率低&#xff0c;画质模…

priority_queue的模拟实现和仿函数

priority_queue模拟 首先查看源代码&#xff0c;源代码就在queue剩下的部分中 push_heap是STL库中的堆算法&#xff0c;STL库中包装有支持堆的算法&#xff0c;在algorithm.h中&#xff1a; 只要不断用堆的形式插入数据&#xff0c;就会形成堆。 priority_queue模拟——初版 pr…

网络协议——STP协议是什么?是如何实现的?

作者&#xff1a;Insist-- 个人主页&#xff1a;insist--个人主页 作者会持续更新网络知识和python基础知识&#xff0c;期待你的关注 目录 一、STP协议是什么 二、为什么需要STP协议 三、STP的实现过程 ​编辑 1、选举跟桥 2、给非跟桥交换机选举跟端口 3、给每个网段选…