对于
t
t
的每个前缀,双指针扫一遍求出在 中出现的最左的位置,对每个后缀求出最右的位置。然后枚举
t
t
中分割的位置就好了。
注意判断 在一个子串中出现的情况。
UPD:
感谢 lifelikes 指出错误。
假如要求出最左的位置,由于要保证选的字符串长度至少为
k
k
,所以初始时指针在 。
但这样有一种情况,就是
t
t
串的一个前缀可以匹配 串以
k
k
为右端的一个子串,而此时指针已经移到了后面。一种解决办法是在这种情况出现时将指针移到 ,但这样无法保证复杂度。
我们先求出
t
t
中所有能匹配 串以
k
k
为右端的一个子串的位置,记最左边的位置为 ,最右边的位置为
second
s
e
c
o
n
d
,那么在
[first+1,second−1]
[
f
i
r
s
t
+
1
,
s
e
c
o
n
d
−
1
]
中的位置肯定不是最优的,所以分别做
[1,first],[second,m]
[
1
,
f
i
r
s
t
]
,
[
s
e
c
o
n
d
,
m
]
就好了。
#include<bits/stdc++.h>
using namespace std;
const int N=500010;
const int H=37;
const int P=1000000007;
int k,n,m;
int pos;
int f[N],g[N];
int hs[N],ht[N],p[N];
int x,y,lst,fst;
char s[N],t[N];
bool Check(int l,int r,int L,int R) {
int t=(hs[r]-1ll*hs[l-1]*p[r-l+1])%P;
t=(t+P)%P;
int T=(ht[R]-1ll*ht[L-1]*p[R-L+1])%P;
T=(T+P)%P;
return t==T;
}
bool Solve() {
pos=1;
for(int i=1;i<=m;i++) {
while(pos<i) ++pos;
while(pos<=n) {
if(Check(pos-i+1,pos,1,i)) break;
++pos;
}
}
if(max(pos,k)+k<=n) {
printf("Yes\n%d %d\n",max(1,pos-k+1),max(pos,k)+1);
return 1;
}
pos=n;
for(int i=m;i;i--) {
while(pos+m-i>n) --pos;
while(pos) {
if(Check(pos,pos+m-i,i,m)) break;
--pos;
}
}
if(n-max(n-pos+1,k)-k>0) {
printf("Yes\n%d %d\n",n-max(n-pos+1,k)-k,n-max(n-pos+1,k)+1);
return 1;
}
return 0;
}
void Solve1(int l,int r) {
int pos=k;
for(int i=l;i<=r;i++) {
while(pos<i) ++pos;
while(pos<=n) {
if(Check(pos-i+1,pos,1,i)) break;
++pos;
}
f[i]=pos<=n?pos:0;
}
}
void Solve2(int l,int r) {
int pos=n-k+1;
for(int i=r;i>=l;i--) {
while(pos+m-i>n) --pos;
while(pos) {
if(Check(pos,pos+m-i,i,m)) break;
--pos;
}
g[i]=pos;
}
}
int main() {
scanf("%d%d%d",&n,&m,&k);
scanf("%s%s",s+1,t+1);
p[0]=1;
for(int i=1;i<=n;i++) p[i]=1ll*p[i-1]*H%P;
for(int i=1;i<=n;i++) hs[i]=(1ll*hs[i-1]*H+s[i]-'a')%P;
for(int i=1;i<=m;i++) ht[i]=(1ll*ht[i-1]*H+t[i]-'a')%P;
if(m<=k&&Solve()) return 0;
lst=fst=1;
for(int i=1;i<=m&&i<=k;i++)
if(Check(k-i+1,k,1,i)) {
lst=i;
if(fst!=1) fst=i;
f[i]=k;
}
Solve1(1,fst-1);Solve1(lst+1,min(m,k));
lst=fst=m;
for(int i=m;i&&i>=m-k+1;i--)
if(Check(n-k+1,n-k+1+m-i,i,m)) {
lst=i;
if(fst!=m) fst=i;
g[i]=n-k+1;
}
Solve2(fst+1,m);Solve2(max(1,m-k+1),lst-1);
for(int i=1;i<m;i++)
if(f[i]&&g[i+1]&&f[i]<g[i+1]) {
x=f[i]-k+1;y=g[i+1];
break;
}
if(x) printf("Yes\n%d %d\n",x,y);else puts("No");
return 0;
}