【题目链接】 http://www.spoj.pl/problems/PHRASES/

【题目大意】

  求在每个字符串中出现至少两次的最长的子串

【题解】

  注意到这么几个关键点:最长,至少两次,每个字符串。
  首先对于最长这个条件,我们可以想到二分答案,
  然后利用后缀数组所求得的三个数组判断是否满足条件。
  其次是出现两次,每次出现这个条件的时候,
  我们就应该要想到这是最大值最小值可以处理的,
  将出现在同一个字符串中的每个相同字符串的起始位置保存下来,
  如果最小值和最大值的差距超过二分长度L,则表明在这个字符串中这个条件是可行的。
  将所有的字符串通过拼接符连接在一起,做一遍后缀数组,
  现在我们根据h数组将大于二分长度的前后后缀分为一组,
  每当存在分组中的后缀数量大于2*n,
  就说明这个字符串有可能是我们需要的答案,那么对它进行检验,
  检验所有可能合法的字符串就可以完成对一个长度的判断了。

【代码】

#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
const int N=1000010;
int n,m,rank[N],sa[N],h[N],tmp[N],cnt[N],ans,a[N],s[N]; char str[N];
void suffixarray(int n,int m){int i,j,k;n++;for(i=0;i<2*n+5;i++)rank[i]=sa[i]=h[i]=tmp[i]=0;for(i=0;i<m;i++)cnt[i]=0;for(i=0;i<n;i++)cnt[rank[i]=s[i]]++;for(i=1;i<m;i++)cnt[i]+=cnt[i-1];for(i=0;i<n;i++)sa[--cnt[rank[i]]]=i;for(k=1;k<=n;k<<=1){for(i=0;i<n;i++){j=sa[i]-k;if(j<0)j+=n;tmp[cnt[rank[j]]++]=j;}sa[tmp[cnt[0]=0]]=j=0;for(i=1;i<n;i++){if(rank[tmp[i]]!=rank[tmp[i-1]]||rank[tmp[i]+k]!=rank[tmp[i-1]+k])cnt[++j]=i;sa[tmp[i]]=j;}memcpy(rank,sa,n*sizeof(int));memcpy(sa,tmp,n*sizeof(int));if(j>=n-1)break;}for(j=rank[h[i=k=0]=0];i<n-1;i++,k++)while(~k&&s[i]!=s[sa[j-1]+k])h[j]=k--,j=rank[sa[j]+1];
}int first=0,len[N],u,K;
vector<int> S[N];
int Min[15],Max[15];
bool check(int L){int cur=-1;for(int i=1;i<=u;i++){if(h[i]<L)S[++cur].clear();S[cur].push_back(i);}for(int i=0;i<=cur;i++){if(S[i].size()>=2*n){memset(Min,-1,sizeof(Min));memset(Max,-1,sizeof(Max));for(int j=0;j<S[i].size();j++){int k=S[i][j];int x=upper_bound(a,a+n+1,sa[k])-a-1;Min[x]=Min[x]==-1?sa[k]:min(Min[x],sa[k]);Max[x]=Max[x]==-1?sa[k]:max(Max[x],sa[k]);}bool flag=1;for(int i=0;i<n;i++){if(Min[i]==-1||Max[i]-Min[i]<L){flag=0;break;}}if(flag)return 1;}}return 0;
}int T;
int main(){scanf("%d",&T);while(T--){scanf("%d",&n);int tmp=200; u=0;for(int i=0;i<n;i++){scanf("%s",str);len[i]=strlen(str);for(int j=0;j<len[i];j++)s[u++]=(int)str[j];s[u++]=tmp++;}tmp=0; s[u]=0; //注意处理完的字符串最后封零for(int i=0;i<=n;i++){a[i]=tmp;if(i<n)tmp=tmp+(i==0?len[i]:len[i]+1);}suffixarray(u,310);int l=1,r=10000,ans=0;while(l<=r){int mid=(l+r)>>1;if(check(mid))ans=mid,l=mid+1;else r=mid-1;}printf("%d\n",ans);}return 0;
}

  

转载于:https://www.cnblogs.com/forever97/p/PHRASES.html

SPOJ 220 Relevant Phrases of Annihilation(后缀数组+二分答案)相关推荐

  1. SPOJ - PHRASES Relevant Phrases of Annihilation(后缀数组+二分)

    题目链接:点击查看 题目大意:给出n个字符串,求出n个字符串中的最长公共子串,且该字串必须在每个字符串中出现至少两次且不能重叠 题目分析:结合之前整理过的: 一个字符串中两个不重叠的重复子串的最大长度 ...

  2. SPOJ - PHRASES Relevant Phrases of Annihilation —— 后缀数组 出现于所有字符串中两次且不重叠的最长公共子串...

    题目链接:https://vjudge.net/problem/SPOJ-PHRASES PHRASES - Relevant Phrases of Annihilation no tags  You ...

  3. SPOJ220 Relevant Phrases of Annihilation(后缀数组)

    引用罗穗骞论文中的话: 先将n 个字符串连起来,中间用不相同的且没有出现在字符串中的字符隔开,求后缀数组.然后二分答案,再将后缀分组.判断的时候,要看是否有一组后缀在每个原来的字符串中至少出现两次,并 ...

  4. bzoj3277 串 (后缀数组+二分答案+ST表)

    常见操作:先把所有串都连到一起,但中间加上一个特殊的符号(不能在原串中/出现过)作为分割 由于全部的子串就等于所有后缀的所有前缀,那我们对于每一个后缀,去求一个最长的前缀,来满足这个前缀在至少K个原串 ...

  5. [BZOJ4310]跳蚤-后缀数组-二分答案

    跳蚤 Description 很久很久以前,森林里住着一群跳蚤.一天,跳蚤国王得到了一个神秘的字符串,它想进行研究.首先,他会把串分成不超过 k 个子串,然后对于每个子串 S,他会从S的所有子串中选择 ...

  6. BZOJ3998 TJOI2015弦论(后缀数组+二分答案)

    先看t=1的情况.显然得求出SA(因为我不会SAM).我们一位位地确定答案.设填到了第len位,二分这一位填什么之后,在已经确定的答案所在的范围(SA上的某段区间)内二分,找到最后一个小于当前串的后缀 ...

  7. HDU - 5030 Rabbit's String(后缀数组+二分)

    题目链接:点击查看 题目大意:给出一个字符串,现在要求将其分为不大于k个连续的子串,对于每个子串求出字典序最大的子串,现在要求所有子串的最大子串的最大值最小,输出这个最大子串 题目分析:最大值最小,标 ...

  8. HDU - 5008 Boring String Problem(后缀数组+二分)

    题目链接:点击查看 题目大意:给出一个字符串,接下来给出 q 个询问,每次询问字符串中第 k 大的子串,要求输出该字串的左右端点,如果有多个答案,输出左端点最小的一个 题目分析:因为在求出后缀数组后, ...

  9. 洛谷 P4094 [HEOI2016/TJOI2016]字符串 后缀数组+二分+主席树

    题目链接 后缀数组 题目分析: sa[i] – 第i小的后缀的编号 rank[i] --编号为i的后缀排第几: height[i] – 第i和第i-1的最长lcp最长公共前缀: 1.二分答案,答案肯定 ...

最新文章

  1. 线程池的介绍及简单实现
  2. Touch Event
  3. Weblogic 集群部署说明 --转
  4. 关于HTML下overflow-y:auto无效、清除HTML默认边距、解决去除手机访问网页时的左右多余空白的三个方法汇总
  5. hibernate的映射之二(一对多双向关联)
  6. 图说苹果工作站-MAC PRO
  7. 介绍一个办法查TREX search的log
  8. 杭州高职单招计算机试题,河北省高职单招职业技能练习题20
  9. JPA 2.1:不同步的持久性上下文
  10. webstorm 内存溢出怎么弄_webstrom 内存溢出,软件崩溃卡死解决的方法
  11. poj 3268 SliverrCowParty Dijkstra最短路(单向路转向)
  12. 搭建Discuz论坛的两种方式
  13. javascript焦点图自动播放
  14. zabbix安装详解
  15. Adobe 全家桶系列卸载工具
  16. Make sure ‘SystemCfg‘ is registered using qRegisterMetaType
  17. Android 10 低内存应用白名单和应用保活
  18. 篮球记分牌c语言程序和报告,基于单片机的篮球比赛电子记分牌最终版(备份存档)...
  19. 绘真妙笔千山怎么在电脑上玩 绘真妙笔千山安卓模拟器教程
  20. 网名闲话之“茶乡浪子”

热门文章

  1. 解决虚拟机 正在决定eht0 的ip信息失败 无链接-- 添加虚拟网卡
  2. UITableView,UICollectionView,UIScrollView快速返回顶部
  3. Lecture 2 Asymptotic Notation
  4. (马世龙)Linux下CACTI完全搭建技术文档二
  5. 《游戏编程模式》一7.8 并发状态机
  6. 致敬2016,拥抱2017
  7. Objective-c 程序结构
  8. Statspack安装心得
  9. GitHub动作简介
  10. leetcode257. 二叉树的所有路径(回溯算法)