CF427D

SA的奇技淫巧,其实就是板子。

题意:

给定两个字符串,求最短的满足各只出现一次的连续公共字串

解析:

一般情况下,SA都是用来求最长公共前缀的,好像和这道题所求的最短公共子串没有任何关系。

但我们依然可以通过类比思路得出:

想一想为什么要寻找zz最大的元素?

因为如果小于最大值,那么最大值就会包含这个序列。

所以答案就是,一个元素,没有z值比这个元素大的,自然就是要选z的最大值

从上述思路,寻找如何让答案尽量小

z的次大值自然是不行的,但是发现次大值+1是满足条件的

一方面,它比最大值小,所以被唯一的最大值包含;另一方面,它比次大值大,所以仅被最大值包含

所以可证,次大值+1也是唯一的。

所以,按如上方法,对 $ S_1 $ 的每一个后缀求最大值和次大值,再用次大值更新答案

然后我们就能愉快的用SA解决这个问题了。

CODE:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring> using namespace std; #define LL long long
#define N 10010 string s1,s2,str;
int SA[N],rk[N],tp[N],cnt[N];
int len,tot,m,a[N],height[N]; void qsort() {
for(int i = 1 ; i <= m ; i++) cnt[i] = 0;
for(int i = 1 ; i <= tot; i++) cnt[rk[i]]++;
for(int i = 1 ; i <= m ; i++) cnt[i] += cnt[i - 1];
for(int i = tot ; i >= 1 ; i--) SA[cnt[rk[tp[i]]]--] = tp[i];
}
inline bool cmp(int *f,int x,int y,int w) {
return f[x] == f[y] && f[x + w] == f[y + w];
}
void build_SA() {
m = 127;
for(int i = 1 ; i <= tot ; i++) {
rk[i] = a[i];
tp[i] = i;
}
qsort();
for(int w = 1 , p = 0 ; p < tot ; w += w,m = p) {
p = 0;
for(int i = tot - w + 1 ; i <= tot ; i++) tp[++p] = i;
for(int i = 1 ; i <= tot ; i++) {
if(SA[i] > w) tp[++p] = SA[i] - w;
}
qsort();
swap(rk,tp);
rk[SA[1]] = p = 1;
for(int i = 2 ; i <= tot ; i++)
rk[SA[i]] = cmp(tp,SA[i],SA[i - 1],w) ? p : ++p;
}
int j = 0, k = 0;
for(int i = 1 ; i <= tot ; height[rk[i++]] = k) {
for(k = k ? k - 1 : k, j = SA[rk[i] - 1] ; a[i + k] == a[j + k] ; k++);
}
}
inline bool check(int k,int div) {
int cnt1 = 0,cnt2 = 0;
for(int i = 1 ; i <= tot ; i++) {
if(height[i] < k) {
if(cnt1 == 1 && cnt2 == 1) return true;
cnt1 = cnt2 = 0;
if(SA[i] <= div) cnt1++;
else if(SA[i] >= div) cnt2++;
continue;
}
if(SA[i] <= div) cnt1++;
else if(SA[i] >= div) cnt2++;
}
return cnt1 == 1 && cnt2 == 1;
} int main() {
cin>>s1>>s2;
len = s1.length();
str = s1 + '#' + s2;//加入'#'表示两个字符串的分界点。
tot = len + s2.length() + 1;
for(int i = 1 ; i <= tot ; i++) a[i] = str[i - 1];
build_SA();
int ans = -1;
for(int i = 1 ; i <= len ; i++) {
if(check(i,len)) {
ans = i;
break;
}
}
printf("%d \n",ans);
//system("pause");
return 0;
}

CF427D的更多相关文章

随机推荐

  1. python 的常见排序算法实现

    python 的常见排序算法实现 参考以下链接:https://www.cnblogs.com/shiluoliming/p/6740585.html 算法(Algorithm)是指解题方案的准确而完 ...

  2. VBA开发项目分享-1

    这个项目的目的是使用VBA制作一个股票筛选器,股票的指标数据存放在多个工作表,输入多个指标的查询条件,可以从相应的工作表里查询出符合条件的股票数据并返回.项目涉及的VBA知识结构有字典.数组.OLED ...

  3. mysql 添加grant权限

    GRANT USAGE ON *.* TO 'xxxx'@'x.%.%.%' WITH GRANT OPTION;

  4. 基于PROMISE解决回调地狱问题

    回调地狱问题: 在使用JavaScript时,为了实现某些逻辑经常会写出层层嵌套的回调函数,如果嵌套过多,会极大影响代码可读性和逻辑,这种情况也被成为回调地狱.比如说你要把一个函数 A 作为回调函数, ...

  5. java代码实现mock数据

    废话不多说,直接上代码. /** * 发get请求,获取文本 * * @param getUrl * @return 网页context */ public static String sendGet ...

  6. AS中Git操作和Git工作流

    查阅并学习安卓中如何使用Git,并扩展学了下Gitflow工作流程,以下是我觉得有用,系统,详细的相关文章. AS中Git操作系列 https://www.jianshu.com/p/0cc226ee ...

  7. 9.EL表达式 和 JSTL核心标签库

    EL表达式 1./*获取数据*/ (某个web域中的对象,访问javabean的属性.访问List集合.访问Map集合.访问数组) <html> <head> <titl ...

  8. gulp die('click').live('click' composer

    gulp  die('click').live('click' composer packagist.org https://getcomposer.org/ 下载后 php composer.pha ...

  9. 关于百度Tongji Api的文档补充

    百度统计的Tongji Api好像没有人维护了,文档缺胳膊少腿也没人理. 今天在这里指出其中一点,因为这一点花时间也没有傻思考的乐趣的. 引用自百度Tongji API文档 这个文档缺了很多东西,其中 ...

  10. MyBatis 插入失败后爆出 500 ,如何捕获异常?

    我们在使用 Mybatis 的时候,会出现以下场景 数据表里有一些字段被设置为了 不可为 null 但是我们的用户在提交表单的时候没有提交所需的 字段数据 然后 Mybatis 在数据库做操作的时候就 ...