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. Java生成随机数列表

    生成随机数列表 1.Java8以前 (1)Math.random private List<UserEntity> random1() { ArrayList<UserEntity& ...

  2. 集合源码阅读——ArrayList

    ArrayList 关键点: >>扩容每次扩容1.5倍 >>modcount的作用 >>ArrayList的父类AbstractList的成员变量 >> ...

  3. canvas之五角星的绘制

    <html> <head> <meta charset=utf-8> <title>绘制简单图形线及矩形</title> <style ...

  4. Asp.net Core 微信小程序支付

    最近要做一个微信小程序支付的功能 在网上找了一下 .net Core做微信支付的博客 和 demo 几乎没有 自己研究了好几天 参考了 很多 大牛的博客 勉强做出来了  因为参数都没有 比如 opid ...

  5. Xcode8警告⚠️ Empty paragraph passed to '@xxx' command

    问题 Xcode8升级后,之前添加的注释会有很多警告 解决方法 基础知识,就是在编译选项中,添加警告屏蔽 解决步骤 显示警告信息 显示警告信息.png 查看警告类型 查看警告类型.png 屏蔽警告 W ...

  6. DDOS 攻击的防范教程--转载自阮一峰的博客

    一个多月前,我的个人网站遭受 DDOS 攻击,下线了50多个小时.这篇文章就来谈谈,如何应对这种攻击. 需要说明的是,我对 DDOS 并不精通,从没想过自己会成为攻击目标.攻击发生以后,很多素昧平生的 ...

  7. WPF实战案例-在线程内同步集合数据到UI线程

    有这样一个场景,在vm中,我们为了ui的体验,会异步访问后端接口,获取数据集合,如果这个集合绑定到界面,并且在线程内,怎么处理? 有人讲:this.Dispatcher.Invoke,如果在vm内呢? ...

  8. cbv装饰器 中间件 跨站请求伪造

    给cbv下面的函数加装饰器 写一个验证用户登录的程序 前端页面 # 写一个装饰器验证session def login_auth(func): def inner(request,*args,**kw ...

  9. 2.03_01_Python网络爬虫urllib2库

    一:urllib2库的基本使用 所谓网页抓取,就是把URL地址中指定的网络资源从网络流中抓取出来.在Python中有很多库可以用来抓取网页,我们先学习urllib2. urllib2 是 Python ...

  10. vue2.0+webpack+vuerouter+vuex+axios构建项目基础

    前言 本文讲解的是vue2.0+webpack+vuerouter+vuex+axios构建项目基础 步骤 1.全局安装webpack,命令 npm install webpack -g 注意,web ...