●BZOJ 3796 Mushroom追妹纸
题链:
http://www.lydsy.com/JudgeOnline/problem.php?id=3796
题解:
题意:
给出三个串 A,B,C
找出一个最长串 S,
使得 S是A,B 的子串,但是 C不是S的子串。
首先,对于第一二个限制,只需要把 A,B串用一个分隔符连接在一起。
求出sa[],rank[],height[]数组,
那么在排好序的后缀中,相邻的两个后缀如果一个属于A串(令在A串p位置),另一个属于B串的话,
则可能贡献答案。 但应该贡献多少呢?
由于有 C 的限制,即在LCP中选出的最长前缀里不能出现 C串。
所以考虑如下做法。
定义 appear[i] 表示 A串的 i位置后面第一次出现 C串的开始位置。
求出这个以后,那么就可以贡献答案了:
ANS=max(ANS,min(hei[i],appear[i]-p+LenC-1))
以下是求pre数组的做法:
把 A串作为文本串,C串作为匹配串,用 KMP去匹配。
并在 A串中每一个出现 C串的开始位置打上标记
然后反向跑一边,得出appear[]数组。
总的复杂度 (Nlog2N+N)
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#define MAXN 100500
#define INF 0x3f3f3f3f
#define filein(x) freopen(#x".in","r",stdin);
#define fileout(x) freopen(#x".out","w",stdout);
using namespace std;
char S[MAXN],A[MAXN],B[MAXN],C[MAXN];
int sa[MAXN],rak[MAXN],hei[MAXN],app[MAXN];
int L1,L2,L3,ANS;
void build(int N,int M){
static int cc[MAXN],ta[MAXN],tb[MAXN],*x,*y,h,p;
x=ta; y=tb; h=0;
for(int i=0;i<M;i++) cc[i]=0;
for(int i=0;i<N;i++) cc[x[i]=S[i]]++;
for(int i=1;i<M;i++) cc[i]+=cc[i-1];
for(int i=N-1;i>=0;i--) sa[--cc[x[i]]]=i;
for(int k=1;p=0,k<N;k<<=1){
for(int i=N-k;i<N;i++) y[p++]=i;
for(int i=0;i<N;i++) if(sa[i]>=k) y[p++]=sa[i]-k;
for(int i=0;i<M;i++) cc[i]=0;
for(int i=0;i<N;i++) cc[x[y[i]]]++;
for(int i=1;i<M;i++) cc[i]+=cc[i-1];
for(int i=N-1;i>=0;i--) sa[--cc[x[y[i]]]]=y[i];
swap(x,y); y[N]=-1; x[sa[0]]=0; M=1;
for(int i=1;i<N;i++)
x[sa[i]]=y[sa[i]]==y[sa[i-1]]&&y[sa[i]+k]==y[sa[i-1]+k]?M-1:M++;
if(M>=N) break;
}
for(int i=0;i<N;i++) rak[sa[i]]=i;
for(int i=0,j;i<N;i++){
if(h) h--;
if(rak[i]){
j=sa[rak[i]-1];
while(S[i+h]==S[j+h]) h++;
}
hei[rak[i]]=h;
}
}
void KMP(char *T,char *S){
static int nxt[MAXN],i,j,k,lT,lS;
memset(app,0x3f,sizeof(app));
nxt[0]=-1; lT=strlen(T); lS=strlen(S);
j=0; k=-1;
while(j<lS){
if(k==-1||S[j]==S[k]) j++,k++,nxt[j]=k;
else k=nxt[k];
}
i=0; j=0;
while(i<lT){
if(j==-1||T[i]==S[j]){
i++; j++;
if(!S[j]) app[i-lS]=i-lS;
}
else j=nxt[j];
}
for(int i=lT-1;i>=0;i--) app[i]=min(app[i],app[i+1]);
}
bool bel(int p){
if(p<L1) return 0;
else return 1;
}
int main()
{
int N;
scanf(" %s %s %s",A,B,C);
L1=strlen(A); L2=strlen(B); L3=strlen(C); N=L1+L2+1;
for(int i=0;i<L1;i++) S[i]=A[i]; S[L1]='&';
for(int i=0;i<L2;i++) S[L1+i+1]=B[i]; S[N]=0;
build(N,300);
KMP(A,C);
for(int i=1,p,q;i<N;i++){
if(hei[i]<=ANS||!(bel(sa[i-1])^bel(sa[i]))) continue;
if(!bel(sa[i-1])) p=sa[i-1]; else p=sa[i]; q=app[p];
ANS=max(ANS,min(hei[i],q-p+L3-1));
}
printf("%d",ANS);
return 0;
}
●BZOJ 3796 Mushroom追妹纸的更多相关文章
- [BZOJ 3796]Mushroom追妹纸
[BZOJ 3796]Mushroom追妹纸 题目 Mushroom最近看上了一个漂亮妹纸.他选择一种非常经典的手段来表达自己的心意——写情书.考虑到自己的表达能力,Mushroom决定不手写情书.他 ...
- bzoj 3796 Mushroom追妹纸——后缀数组
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3796 长度一般都是 1e5 ,看这个是 5e4 ,一看就是把两个串接起来做. 自己本来想的是 ...
- bzoj 3796 Mushroom追妹纸 —— 后缀数组
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3796 先把三个串拼在一起,KMP 求 s1 , s2 中每个位置和 s3 的匹配情况: 注意 ...
- bzoj 3796: Mushroom追妹纸 AC自动机+后缀自动机+dp
题目大意: 给定三个字符串s1,s2,s3,求一个字符串w满足: w是s1的子串 w是s2的子串 s3不是w的子串 w的长度应尽可能大 题解: 首先我们可以用AC自动机找出s3在s1,s2中出现的位置 ...
- BZOJ 3796 Mushroom追妹纸 哈希+二分(+KMP)
先把两个串能匹配模式串的位置找出来,然后标记为$1$(标记在开头或末尾都行),然后对标记数组求一个前缀和,这样可以快速查到区间内是否有完整的一个模式串. 然后二分子串(答案)的长度,每次把长度为$md ...
- bzoj 3796: Mushroom追妹纸【二分+后缀数组+st表】
把三个串加上ASCII大于z的分隔符连起来,然后求SA 显然每个相同子串都是一个后缀的前缀,所以枚举s1的每个后缀的最长和s2相同的前缀串(直接在排序后的数组里挨个找,最近的两个分别属于s1和s2的后 ...
- 【BZOJ3796】Mushroom追妹纸 二分+hash
[BZOJ3796]Mushroom追妹纸 Description Mushroom最近看上了一个漂亮妹纸.他选择一种非常经典的手段来表达自己的心意——写情书.考虑到自己的表达能力,Mushroom决 ...
- 【bzoj3796】Mushroom追妹纸 hash/sa+kmp+二分
Description Mushroom最近看上了一个漂亮妹纸.他选择一种非常经典的手段来表达自己的心意--写情书.考虑到自己的表达能力,Mushroom决定不手写情书.他从网上找到了两篇极佳的情书, ...
- BZOJ3796 : Mushroom追妹纸
将S1与S2用#号拼接在一起形成S串 将S3与S串跑KMP求出S3在S串中每次出现的位置l[i] 对于S串每个后缀i,求出f[i]表示该串不包含S3串的最长前缀 然后求出S串的后缀数组 先从小到大扫描 ...
随机推荐
- c++第0次作业
1.你认为大学的学习生活.同学关系.师生应该是怎样? 随着大学生活的慢慢到来,我开始领悟到大学并不是自由的天堂,相反,我们更加的走进社会这个牢笼.在这个牢笼中有着从前的我们并不需要在意和考虑的规则与问 ...
- 第十二条:考虑实现Comparable接口
与前面讨论的方法不同,compareTo()方法并没有在Object类中定义.相反,它是Comparable接口中唯一的方法. 一个类的实例对象要想是可以比较大小的,那么这个类需要实现Comparab ...
- Spring MVC Restful Put方法无法获取参数值
Spring MVC Restful 无法通过@ReqeustParam获取参数值 原因是Tomcat只支持POST/GET获取参数值,对于PUT这些方法需要通过HttpPutFormContentF ...
- System.Reflection名称空间下的程序集类Assembly应用.
利用反射中的程序集类(Assembly--抽象类)动态加载类库(.dll)或者可执行程序(.exe). 优点:①.可以消除if条件的逻辑判断.②.减少内存资源.③.有利于程序扩展. 缺点... 使用静 ...
- vue-入门
数据绑定 <!--步骤1:创建html文件--> <!DOCTYPE html> <html lang="en"> <head> ...
- OAuth2.0学习(1-10)新浪开放平台微博认证-手机应用授权和refresh_token刷新access_token
1.当你是使用微博官方移动SDK的移动应用时,授权返回access_token的同时,还会多返回一个refresh_token: JSON 1 2 3 4 5 6 { "access ...
- Archaius 原理
Archaius 原理 Archaius是什么? Archaius提供了动态修改配置的值的功能,在修改配置后,不需要重启应用服务.其核心思想就是轮询配置源,每一次迭代,检测配置是否更改,有更改重新更新 ...
- Mysql 测试题
一. 表结构和数据 作业要求 /* Navicat Premium Data Transfer Source Server : localhost Source Server Type : MySQL ...
- POJ-3069 Saruman's Army---区间选点
题目链接: https://vjudge.net/problem/POJ-3069 题目大意: 在一条直线上,有n个点.从这n个点中选择若干个,给他们加上标记.对于每一个点,其距离为R以内的区域里必须 ...
- hdu1060 Leftmost Digit---求N的N次方的首位(对数)
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=1060 题目描述:求N的N次方的第一位. 思路: 第一次做这种类型的题目,学到了如何运用对数. 首先推 ...