[Bzoj4641]基因改造

题意

给定字符串s和t,可以交换任意字符,求s的哪些连续子串与t匹配

---

题解

题目可以转化为

字符串匹配,匹配规则是模式串与文本串的子串不一样的情况相同 (e.g 3 1 3=2 1 2=1 3 1) ( e . g   3   1   3 = 2   1   2 = 1   3   1 )

既然是查看不一样的情况那么就记录一下相同元素上一次出现的位置(老套路?)

然后令 a[i]/b[i]=i a [ i ] / b [ i ] = i -上一次的位置,就可以直接 KMP K M P 了

注意几种特殊的合法情况

1: 1 : 当前字符 s s 是第一次出现

2: 2 : 当前字符位置 − − 上一次出现的位置大于等于当前匹配长度

#include
#include
#define fp(i,a,b) for(int i=a,I=b;i<=I;++i)
#define fd(i,a,b) for(int i=a,I=b;i>=I;--i)
#define file(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout)
template<class T>inline bool chkmax(T&a,const T&b){return a1:0;}
template<class T>inline bool chkmin(T&a,const T&b){return a>b?a=b,1:0;}
using namespace std;
char ss[1<<17],*A=ss,*B=ss;
inline char gc(){return A==B&&(B=(A=ss)+fread(ss,1,1<<17,stdin),A==B)?-1:*A++;}
template<class T>inline void sd(T&x){
    char c;T y=1;while(c=gc(),(c<48||571)if(c==45)y=-1;x=c^48;
    while(c=gc(),4758)x=(x<<1)+(x<<3)+(c^48);x*=y;
}
char sr[1<<21],z[20];int C=-1,Z;
template<class T>inline void we(T x){
    if(x<0)sr[++C]=45,x=-x;
    while(z[++Z]=x%10+48,x/=10);
    while(sr[++C]=z[Z],--Z);sr[++C]=32;
    if(C>1<<20)fwrite(sr,1,C+1,stdout),C=-1;
}
const int N=1e6+5;
typedef int arr[N];
int cas,n,m;arr a,b,fa,pos,nx;
inline bool ok(int x,int y){return b[x]==y||(b[x]==-1&&(x<=y||y==-1));}
inline void solve(){
    sd(n),sd(m);int x;pos[0]=0;
    memset(fa,-1,sizeof fa);
    fp(i,1,n)sd(x),a[i]=(fa[x]==-1)?-1:i-fa[x],fa[x]=i;
    memset(fa,-1,sizeof fa);
    fp(i,1,m)sd(x),b[i]=(fa[x]==-1)?-1:i-fa[x],fa[x]=i;
    for(int i=2,j=0;i<=m;++i){
        while(j&&!ok(j+1,b[i]))j=nx[j];
        if(ok(j+1,b[i]))++j;nx[i]=j;
    }
    for(int i=1,j=0;i<=n;++i){
        while(j&&!ok(j+1,a[i]))j=nx[j];
        if(ok(j+1,a[i]))++j;
        if(j==m)pos[++pos[0]]=i-m+1,j=nx[j];
    }
    we(pos[0]);sr[C]='\n';
    fp(i,1,pos[0])we(pos[i]);sr[++C]='\n';
}
int main(){
    #ifndef ONLINE_JUDGE
        file("reform");
    #endif
    sd(cas),sd(n);
    while(cas--)solve();
    fwrite(sr,1,C+1,stdout);
return 0;
}