bzoj 4032(A的一个最短的子串,它不是B的子串 || A的一个最短的子串,它不是B的子序列 || A的一个最短的子序列,它不是B的子串||A的一个最短的子序列,它不是B的子序列)
在虐各种最长公共子串、子序列的题虐的不耐烦了之后,你决定反其道而行之。
一个串的“子串”指的是它的连续的一段,例如bcd是abcdef的子串,但bde不是。 一个串的“子序列”指的是它的可以不连续的一段,例如bde是abcdef的子串,但bdd不是。 下面,给两个小写字母串A,B,请你计算: (1) A的一个最短的子串,它不是B的子串 (2) A的一个最短的子串,它不是B的子序列 (3) A的一个最短的子序列,它不是B的子串 (4) A的一个最短的子序列,它不是B的子序列
Input
有两行,每行一个小写字母组成的字符串,分别代表A和B。
Output
输出4行,每行一个整数,表示以上4个问题的答案的长度。如果没有符合要求的答案,输出-1.
Sample Input
aabbcc abcabc
Sample Output
2 4 2 4
HINT
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdio>
#include<map>
#define N 20013
#define p1 2000001001
#define ull unsigned long long
using namespace std;
int cnt,n,m,np,p,q,nq,l[N],r[N],pos[N],root;
int ch[N][],fa[N],len[N],last,ch1[N][],f[][];
ull mi[N];
char s[N],s1[N];
map<int,int> mp;
map<ull,int> mp1;
void extend(int i)
{
int c=s1[i]-'a'+;
p=last; np=last=++cnt; l[np]=l[p]+;
for (;!ch[p][c]&&p;p=fa[p]) ch[p][c]=np;
if (!p) fa[np]=root;
else {
q=ch[p][c];
if (l[q]==l[p]+) fa[np]=q;
else {
nq=++cnt; l[nq]=l[p]+;
memcpy(ch[nq],ch[q],sizeof(ch[q]));
fa[nq]=fa[q];
fa[q]=fa[np]=nq;
for (;ch[p][c]==q;p=fa[p]) ch[p][c]=nq;
}
}
}
void solve1()
{
/*mi[0]=1;
for (int i=1;i<=max(n,m)+2;i++) mi[i]=mi[i-1]*p1;
for (int i=1;i<=m;i++){
ull x=(ull)s1[i]*mi[1]; mp1[x]=1;
//cout<<x<<endl;
for (int j=i+1;j<=m;j++) {
x=x*p1+s1[j]*mi[i];
mp1[x]=(j-i+1);
// cout<<x<<endl;
}
}
int mn=m+1;
for (int i=1;i<=n;i++) {
ull x=(ull)s[i]*mi[1];
if (!mp1[x]) mn=min(mn,1);
for (int j=i+1;j<=n;j++){
x=x*p1+s[j]*mi[1];
if (!mp1[x]&&j-i+1<=n) mn=min(mn,j-i+1);
}
}
if (mn>=m+1) PRintf("-1\n");
else printf("%d\n",mn);*/
memset(f,,sizeof(f));
int mn=m+;
for (int i=;i<=n;i++){
int a=,b=;
for (int j=;j<=n;j++) {
if (s[i]==s1[j]) f[i][j]=f[i-][j-]+;
a=max(f[i][j],a);
}
if (a!=i) mn=min(mn,a);
}
if (mn>=m+) printf("-1\n");
else printf("%d\n",mn+);
}
void solve2()//A的最短子串,不是B的子序列
{
memset(f,,sizeof(f));
int mn=m+;
for (int i=;i<=n;i++){
int a=,b=;
for (int j=;j<=m;j++) {
if (s[i]==s1[j]) f[i][j]=b+;
a=max(a,f[i][j]);
b=max(b,f[i-][j]); //f[i][j]表示A中的第i个字符匹配到B中第j个字符,A匹配时必须连续
}
if (a!=i) mn=min(mn,a);//匹配的长度不是i说明,在i+前一字符就是一个合法的子串
}
if (mn>=m+) printf("-1\n");
else printf("%d\n",mn+);
}
void solve3()//A的最短的子序列,不是b的子串。
{
memset(len,/,sizeof(len));
len[]=; int mn=m+; int t;
for (int i=;i<=n;i++)
for (int j=;j<=cnt;j++)
if (!(t=ch[j][s[i]-'a'+])) mn=min(mn,len[j]);//如果到该点匹配不上了,说明在当前基础上加入s[i]就能得到了一个合法的子序列
else len[t]=min(len[t],len[j]+);//len表示的是用a的子序列去匹配后缀自动机中的节点,到节点i能得到的最短长度
if (mn==m+) printf("-1\n");
else printf("%d\n",mn+);
}
void solve4()
{
for (int i=m;i>=;i--){//只要是在当前位置之后出现的字符都可以用来构造子序列,如果我们要选择字符x,那么必然优先选择最靠近当前位置的字符。
for (int j=;j<=;j++)
if (mp[j]) ch1[i][j]=mp[j];
mp[s1[i]-'a'+]=i;
}
memset(len,/,sizeof(len));
len[]=;
int mn=m+; int t;
for (int i=;i<=n;i++)
for (int j=m;j>=;j--)//这里与第三问不同,必须倒序枚举,因为这里相当于是个背包,那么对于一个字符来说,不能连续更新同一层中的位置,否则原串只有一个字符x,正序枚举会出现xx的情况。
if (!(t=ch1[j][s[i]-'a'+])) mn=min(mn,len[j]);
else len[t]=min(len[t],len[j]+);
if (mn==m+) printf("-1\n");
else printf("%d\n",mn+);
}
int main()
{
freopen("sus4.in","r",stdin);
//freopen("sus.out","w",stdout);
scanf("%s",s+);
scanf("%s",s1+);
n=strlen(s+); m=strlen(s1+);
root=last=++cnt;
for (int i=;i<=m;i++)
extend(i);
solve1();
solve2();
solve3();
solve4();
}
先留坑 ,待我变的更厉害好吧
bzoj 4032(A的一个最短的子串,它不是B的子串 || A的一个最短的子串,它不是B的子序列 || A的一个最短的子序列,它不是B的子串||A的一个最短的子序列,它不是B的子序列)的更多相关文章
- 编写一个ComputerAverage抽象类,类中有一个抽象方法求平均分average,可以有参数。定义 Gymnastics 类和 School 类,它们都是 ComputerAverage 的子类。Gymnastics 类中计算选手的平均成绩的方法是去掉一个最低分,去掉一个最高分,然后求平均分;School 中计算平均分的方法是所有科目的分数之和除以总科目数。 要求:定义ComputerAv
题目: 编写一个ComputerAverage抽象类,类中有一个抽象方法求平均分average,可以有参数. 定义 Gymnastics 类和 School 类,它们都是 ComputerAverag ...
- BZOJ 4032: [HEOI2015]最短不公共子串
4032: [HEOI2015]最短不公共子串 Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 446 Solved: 224[Submit][Sta ...
- BZOJ 4032: [HEOI2015]最短不公共子串 后缀自动机 暴力
4032: [HEOI2015]最短不公共子串 题目连接: http://www.lydsy.com/JudgeOnline/problem.php?id=4032 Description 在虐各种最 ...
- bzoj 4032 [HEOI2015]最短不公共子串——后缀自动机
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4032 不是 b 的子串的话就对 b 建后缀自动机,在 a 上枚举从每个位置开始的子串或者找子 ...
- bzoj 4032 [ HEOI 2015 ] 最短不公共子串 —— 后缀自动机+序列自动机
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4032 序列自动机其实就是每个位置记录一下某字母后面第一个出现位置,为了子序列能尽量长. 对字 ...
- BZOJ 4032 Luogu P4112 [HEOI2015]最短不公共子串 (DP、后缀自动机)
这其实是道水题... 题目链接: (bzoj)https://www.lydsy.com/JudgeOnline/problem.php?id=4032 (luogu)https://www.luog ...
- BZOJ.4032.[HEOI2015]最短不公共子串(DP 后缀自动机)
题目链接 1.求A的最短子串,它不是B的子串. 子串是连续的,对B建SAM,枚举起点,在SAM上找到第一个无法匹配点即可.O(n)用SAM能做吗..开始想错了. 2.求A的最短子串,它不是B的子序列. ...
- bzoj 4032: [HEOI2015]最短不公共子串【dp+SAM】
第一.二问: 就是最小的最长公共长度+1,设f[i][j]为a匹配到i,b匹配到j,第一问的转移是f[i][j]=(a[i]==b[j]?f[i-1][j-1]+1:0),第二问的转移是f[i][j] ...
- BZOJ 4032: [HEOI2015]最短不公共子串 (dp*3 + SAM)
转博客大法好 第4个子任务中,为什么只转移最近的一个位置,自己YY吧(多YY有益身体健康). #include <bits/stdc++.h> using namespace std; t ...
随机推荐
- EZOJ #257
传送门 分析 先进行缩点 之后从终点倒着跑 对于一组边如果有一个点不能到达则这组边直接废掉 最后看只用没废掉的边能不能从起点走到终点 代码 #include<iostream> #incl ...
- Selenium运用-漫画批量下载
今天我们要爬去的网站是http://comic.sfacg.com/.漫画网站一般都是通过JavaScript和AJAX来动态加载漫画的,这也就意味着想通过原来爬取静态网站的方式去下载漫画是不可能的, ...
- Java泛型(泛型接口、泛型类、泛型方法)
转载 转载出处:https://www.cnblogs.com/JokerShi/p/8117556.html 泛型接口: 定义一个泛型接口: 通过类去实现这个泛型接口的时候指定泛型T的具体类型. 指 ...
- Java菜鸟之java基础语法,运算符(三)
赋值运算符 (一)JAVA种的赋值运算符 = ,代表代表的等于,一般的形式是 左边变量名称 = 右边的需要赋的指或者表达式,如果左侧的变量类型级别比较高,就把右侧的数据转换成左侧相同的高 ...
- JavaScript - this详解 (一)
侃侃JavaScript中的this this为何如此多变? this总是跟它的执行上下文有关,而在JavaScript总会有开辟新的执行上下文的东西,比如函数,所以,this才如此的变化. 执行环境 ...
- WP8注册文件关联---分享图片
其实这个博友写得差不多了http://www.cnblogs.com/lipan/archive/2013/05/15/3080275.html 里面提供大部分代码,没有提供源代码,小弟就拷贝过来调试 ...
- Ado.NET SqlDataReader详解
ado.net的数据提供程序有三个分别是SqlServer数据提供程序,OLE DB提供程序,ODBC提供程序. 本次记录的是SqlServer提供程序中的一些知识点. ①SqlDataReader必 ...
- WebService 常用的设置
1.修改WebService接收长度 <binding name="IAuthServiceSoap11Binding" maxBufferSize="214748 ...
- Luckily general gradient for spherical harmonics is defined
http://web4.cs.ucl.ac.uk/staff/j.kautz/publications/gradientSH_RS04.pdf
- oracle数据库中将clob字段内容利用java提取出至文本文档中
代码段: 1.执行clob转String public static String ClobToString(Clob sc) throws SQLException, IOException { S ...