[ 杂题 ] Codeforces955D Scissors

对于 t t 的每个前缀,双指针扫一遍求出在 s 中出现的最左的位置,对每个后缀求出最右的位置。然后枚举 t t 中分割的位置就好了。
注意判断 t 在一个子串中出现的情况。

UPD:

感谢 lifelikes 指出错误。
假如要求出最左的位置,由于要保证选的字符串长度至少为 k k ,所以初始时指针在 k
但这样有一种情况,就是 t t 串的一个前缀可以匹配 s 串以 k k 为右端的一个子串,而此时指针已经移到了后面。一种解决办法是在这种情况出现时将指针移到 k ,但这样无法保证复杂度。
我们先求出 t t 中所有能匹配 s 串以 k k 为右端的一个子串的位置,记最左边的位置为 first ,最右边的位置为 second s e c o n d ,那么在 [first+1,second1] [ 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;
}
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值