Noip前的大抱佛脚----字符串

时间:2024-12-12 15:37:32

字符串

Tags:Noip前的大抱佛脚

经验

用FFT求解字符串匹配问题

  • 一一对应

把其中一个\(Reverse\)后,对于每个字符跑一遍FFT,打上\(Tag\)

如果在某个位置上有串长个\(Tag\)那便是匹配上了一处

  • 模糊匹配

\(Fuzzy Search\) 在跑\(FFT\)前把模糊门限值的区间内全部置为1,然后同样的操作

两(多)串DP时状态合并

插入AC自动机,老套路了

最长公共子序列转LIS

求两个串的最长公共子序列,把第二个串的每个值映射到第一个串上 该值的位置

然后对第二个串求LIS即可(最长公共子串就是第二个串的最长连续段)

位运算最大值

异或最大值:建Trie贪心

与最大值:建Trie贪心(1的儿子siz大于2则走)/高维前缀和逐位贪心

或最大值:高维前缀和逐位贪心

但是求(x+B)^A的最大值呢(SCOI2016美味)(当然&操作也是一样的,这种题通常值域很小)

同样贪心地做

开值域线段树,贪心到某一位,需要该位为0或者1,则对应地可以算出x的范围,查询是否有在这个范围之内的x即可

挂链哈希

期望\(O(1)\),当然支持查询多个关键字

const int Mo=YYB;
struct HashTable
{
struct Line{int next,val;}a[Mo];
int head[Mo],cnt;
void reset() {memset(head,0,sizeof(head));cnt=0;}
void Add(int p,int val) {a[++cnt]=(Line){head[p],val};head[p]=cnt;}
int Query(int x)
{
for(int i=head[x%Mo];i;i=a[i].next)
if(a[i].val==x) return 1;return 0;
}
}Hash;

哈希处理回文串

正反哈希前缀和即可求出区间哈希值,然后查询起点到回文中心的正哈希和回文中心到终点的反哈希即可

树哈希

例:一棵无根树的本质不同的独立集个数(\(k\)棵相同子树方案为\(x\)则乘一个可重组合)

字符串模板库

KMP

【模板】KMP字符串匹配

用于两串匹配问题,做法是对子串求next后匹配母串,复杂度\(O(n+m)\)

const int N=1e6+10;
char s1[N],s2[N];
int nxt[N];
int main()
{
scanf("%s%s",s1+1,s2+1);
int l1=strlen(s1+1),l2=strlen(s2+1);
for(int i=2;i<=l2;i++)
{
int j=nxt[i-1];
while(j&&s2[j+1]!=s2[i]) j=nxt[j];
nxt[i]=(s2[j+1]==s2[i])+j;
}
for(int i=1,j=0;i<=l1;i++)
{
while(j&&s2[j+1]!=s1[i]) j=nxt[j];
if(s2[j+1]==s1[i]) j++;
if(j==l2) printf("%d\n",i-l2+1),j=nxt[j];
}
for(int i=1;i<=l2;i++) printf("%d ",nxt[i]);
return puts(""),0;
}

最小循环表示

工艺

\(O(n)\)求一个环从某点断开按一定方向形成的字典序最小的链

int i,j=2,k,l,p,s[610000];
int main()
{
cin>>l;for(i=1;i<=l;i++) cin>>s[i],s[i+l]=s[i];
for(i=1;j<=l&&i<=l&&k<=l;)
{
if(s[i+k]==s[j+k]) {k++;continue;}
s[i+k]<s[j+k]?j+=k+1:i+=k+1;
if(i==j) i++;k=0;
}
for(;p<l;p++) cout<<s[min(i,j)+p]<<" ";
}

Mancher

【模板】manacher算法

求出以每个位置为中心的最长回文串,复杂度\(O(n)\),证明:每次操作要么不动\(while\),要么给两个单调的指针至少\(+1\)


const int N=3e7+10;
char s[N],t[N];
int l,p[N],R,C,Ans;
int main()
{
scanf("%s",t+1);
for(int i=1,len=strlen(t+1);i<=len;i++)
s[++l]='#',s[++l]=t[i];s[++l]='#';
for(int i=1;i<=l;i++)
{
p[i]=i<=R?min(p[C*2-i],R-i):1;
while(s[i+p[i]]==s[i-p[i]]&&i+p[i]<=l&&i-p[i]>=1) p[i]++;
if(i+p[i]-1>R) R=i+p[i]-1,C=i;
Ans=max(Ans,p[i]-1);
}
cout<<Ans<<endl;
}

AC自动机

【模板】AC自动机(加强版)

用于处理多串匹配单串等多串问题,复杂度\(O(n*26)\)

方式是先建\(Trie\),求出\(fail\)指针后建成\(Trie\)图

int n,node,fail[N],ch[N][26];
queue<int> Q;
void Insert(char *s,int id)
{
int x=0,l=strlen(s+1);
for(int i=1;i<=l;i++)
{
int &p=ch[x][s[i]-'a'];
if(!p) p=++node;x=p;
}
}
void Get_fail()
{
for(int i=0;i<26;i++) if(ch[0][i]) Q.push(ch[0][i]);
while(!Q.empty())
{
int x=Q.front();Q.pop();
for(int i=0;i<26;i++)
if(ch[x][i]) fail[ch[x][i]]=ch[fail[x]][i],Q.push(ch[x][i]);
else ch[x][i]=ch[fail[x]][i];
}
}

后缀数组

【模板】后缀排序

用于处理字符串后缀的东西(不过这东西Noip不会考,省选题倒是经常出现)

复杂度\(O(nlogn)\),基于一种倍增和桶排的思想对后缀排序


const int N=1e6+10;
int m=300,t[N],x[N],y[N],rk[N],h[N],SA[N],l;char s[N];
int cmp(int i,int j,int k) {return y[i]==y[j]&&y[i+k]==y[j+k];}
void Sort()
{
for(int i=1;i<=m;i++) t[i]=0;
for(int i=1;i<=l;i++) t[x[i]]++;
for(int i=1;i<=m;i++) t[i]+=t[i-1];
for(int i=l;i>=1;i--) SA[t[x[y[i]]]--]=y[i];
}
void Get_SA()
{
for(int i=1;i<=l;i++) x[i]=s[i],y[i]=i;Sort();
for(int k=1,p=0;k<=l;k<<=1)
{
for(int i=l-k+1;i<=l;i++) y[++p]=i;
for(int i=1;i<=l;i++) if(SA[i]>k) y[++p]=SA[i]-k;
Sort();swap(x,y);x[SA[1]]=p=1;
for(int i=2;i<=l;i++) x[SA[i]]=cmp(SA[i],SA[i-1],k)?p:++p;
if(p==l) break;m=p;p=0;
}
for(int i=1;i<=l;i++) rk[SA[i]]=i;
for(int i=1,j=0;i<=l;i++)
{
while(s[i+j]==s[SA[rk[i]-1]+j]) j++;
h[rk[i]]=j;if(j) j--;
}
}
int main()
{
scanf("%s",s+1);l=strlen(s+1);Get_SA();
for(int i=1;i<=l;i++) printf("%d ",SA[i]);
}

后缀自动机

【模板】后缀自动机

用于处理子串的问题。这家伙比SA好写复杂度还优秀适用范围还广些

不过Noip还是不会考,复杂度\(O(n)\)


const int N=2e6+10;
int l,lst=1,node=1,t[N],A[N],len[N],fa[N];
int ch[N][26],siz[N];char s[N];
void Extend(int c)
{
int f=lst,p=++node;lst=p;
len[p]=len[f]+1;siz[p]=1;
while(f&&!ch[f][c]) ch[f][c]=p,f=fa[f];
if(!f) {fa[p]=1;return;}
int x=ch[f][c],y=++node;
if(len[f]+1==len[x]) {fa[p]=x;node--;return;}
len[y]=len[f]+1;fa[y]=fa[x];fa[x]=fa[p]=y;
memcpy(ch[y],ch[x],sizeof(ch[y]));
while(f&&ch[f][c]==x) ch[f][c]=y,f=fa[f];
}
int main()
{
scanf("%s",s+1);l=strlen(s+1);
for(int i=1;i<=l;i++) Extend(s[i]-'a');
}