【BZOJ4566】找相同字符（后缀数组）

题面

BZOJ

题解

后缀数组的做法，应该不是很难想
首先看到两个不同的串，当然是接在一起求 SA,height SA,height
那么，考虑一下暴力
在两个串各枚举一个后缀，他们的 lcp lcp就是对答案产生的贡献

现在优化一下，按照 SA SA的顺序枚举来处理 lcp lcp
利用一个单调栈维护一下，每次记录一下前面有多少个的贡献和当前答案一样就好啦
只是有点难写。。。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define MAX 440000
#define ll long long
int n,gr[MAX];
int SA[MAX],hg[MAX],rk[MAX];
int a[MAX],t[MAX],x[MAX],y[MAX];
bool cmp(int i,int j,int k){return y[i]==y[j]&&y[i+k]==y[j+k];}
void GetSA()
{int m=50;for(int i=1;i<=n;++i)t[x[i]=a[i]]++;for(int i=1;i<=m;++i)t[i]+=t[i-1];for(int i=n;i>=1;--i)SA[t[x[i]]--]=i;for(int k=1;k<=n;k<<=1){int p=0;for(int i=1;i<=n;++i)y[i]=0;for(int i=n-k+1;i<=n;++i)y[++p]=i;for(int i=1;i<=n;++i)if(SA[i]>k)y[++p]=SA[i]-k;for(int i=0;i<=m;++i)t[i]=0;for(int i=1;i<=n;++i)t[x[y[i]]]++;for(int i=1;i<=m;++i)t[i]+=t[i-1];for(int i=n;i>=1;--i)SA[t[x[y[i]]]--]=y[i];swap(x,y);x[SA[1]]=p=1;for(int i=2;i<=n;++i)x[SA[i]]=cmp(SA[i],SA[i-1],k)?p:++p;if(p>=n)break;m=p;}for(int i=1;i<=n;++i)rk[SA[i]]=i;for(int i=1,j=0;i<=n;++i){if(j)--j;while(a[i+j]==a[SA[rk[i]-1]+j])++j;hg[rk[i]]=j;}
}
int Q1[MAX],top1,S1[MAX];
int Q2[MAX],top2,S2[MAX];
ll ans,sum1,sum2;
char ch[MAX];
int main()
{scanf("%s",ch+1);for(int i=1,l=strlen(ch+1);i<=l;++i)a[++n]=ch[i]-96,gr[n]=1;a[++n]=27;scanf("%s",ch+1);for(int i=1,l=strlen(ch+1);i<=l;++i)a[++n]=ch[i]-96,gr[n]=2;GetSA();for(int i=1,tot=0;i<n;++i){ans+=(gr[SA[i]]==1)?sum2:sum1;ll tmp1=0;while(top1&&Q1[top1]>=hg[i+1]){tmp1+=S1[top1];sum1-=1ll*(Q1[top1]-hg[i+1])*S1[top1];--top1;}++top1;S1[top1]=tmp1;Q1[top1]=hg[i+1];ll tmp2=0;while(top2&&Q2[top2]>=hg[i+1]){tmp2+=S2[top2];sum2-=1ll*(Q2[top2]-hg[i+1])*S2[top2];--top2;}++top2;S2[top2]=tmp2;Q2[top2]=hg[i+1];if(gr[SA[i]]==1)++S1[top1],sum1+=hg[i+1];else ++S2[top2],sum2+=hg[i+1];}printf("%lld\n",ans);return 0;
}

【BZOJ4566】找相同字符（后缀数组）相关推荐

[bzoj4566][HAOI2016]找相同字符(后缀数组)
题目传送门题解这里:把两个串用一个很大的字符连接起来,求一个后缀数组. 考虑怎样暴力的算答案. 在 rank r a n k rank数组中从前往后枚举起点,对于每个枚举的起点,都暴力的往后扫 ...
【HAOI2016/BZOJ4566】找相同字符后缀数组+单调栈
原题走这里鉴于我实在不是很懂单调栈和单调队列这一系列东西,所以我决定稍微具体讲一下单调栈. 恩,本题实质上就是求两个字符串的公共子串数,其中只要出现位置不同,就算是不同的子串. 处理多个字符串的经典 ...
BZOJ4566: [Haoi2016]找相同字符(后缀自动机)
题意题目链接 Sol 直接在SAM上乱搞枚举前缀,用SAM统计可以匹配的后缀,具体在匹配的时候维护和当前节点能匹配的最大值然后再把parent树上的点的贡献也统计上,这部分可以爆跳parent树 ...
【bzoj4566】找相同字符后缀自动机
AC通道:http://www.lydsy.com/JudgeOnline/problem.php?id=4566 [题解] 我们还是先把A串建成SAM,然后让B串在SAM上跑因为相同子串数=不同长 ...
[BZOJ4566][HAOI2016]找相同字符后缀自动机
题目要求的就是B的每个字串在A中的出现次数之和. 我们考虑先建出A串的SAM,每个点所代表的串的个数就是 |Righti|∗(Maxi−Maxfai) |Right_i|*(Max_i-Max_{fa ...
【BZOJ4566】找相同字符，后缀数组
传送门思路: 偶然翻到的一个题苦思冥想算法之时-- 旁边不(jing)会(tong)后缀数组的聪爷爷:这不是后缀数组吗? 赶紧来练一下遗忘的后缀数组(然后手打板子又错了,只能回到博客上重新翻一波以 ...
BZOJ4566: [Haoi2016]找相同字符
BZOJ4566: [Haoi2016]找相同字符 Description 给定两个字符串,求出在两个字符串中各取出一个子串使得这两个子串相同的方案数. 两个方案不同当且仅当这两个子串中有一个位置不同 ...
【BZOJ1031】[JSOI2007]字符加密Cipher 后缀数组
[BZOJ1031][JSOI2007]字符加密Cipher Description 喜欢钻研问题的JS同学,最近又迷上了对加密方法的思考.一天,他突然想出了一种他认为是终极的加密办法 :把需要加密的 ...
长字符串匹配（BWT编码、后缀数组、倍增算法、FM索引）
用 O(m) 时间复杂度找出一个长度为 m 的短字符串在一个长度为 n 的长字符串中的精确匹配(n>>m),限制长短字符串仅由 A.C.G.T 这四种字符组成. 输入:长短字符串输出:短 ...

【BZOJ4566】找相同字符（后缀数组）

题面

题解

【BZOJ4566】找相同字符（后缀数组）相关推荐

最新文章

热门文章