对于S较大的(非满分)字符串哈希解法,望诸君集思广益解决该算法正确性问题

#include<iostream> 
#include<algorithm>
#include<cmath>
using namespace std;
int ans[500000],cnt,d;
int ls[2000001],rs[2000001],v[2000001];
long long n,k,s,x,a[500000],b[500000],r[500000];
unsigned __int128 hs,h[2000001],br,pw[500001];
//base=31,usigned__int128自然溢出哈希
int nw(){
	d++;return d;
}//动态开点
//线段树单点修改
void upd(int ix,long long sl,long long sr,long long x,int c){
	if(sl!=sr){
		if(x<=(sl+sr)>>1){
			if(!ls[ix])ls[ix]=nw();
			upd(ls[ix],sl,(sl+sr)>>1,x,c);
		}else{
			if(!rs[ix])rs[ix]=nw();
			upd(rs[ix],((sl+sr)>>1)+1,sr,x,c);
		}v[ix]=v[ls[ix]]+v[rs[ix]];
		h[ix]=h[ls[ix]]*pw[v[rs[ix]]]+h[rs[ix]];
    //合并标记
	}else v[ix]++,h[ix]=c;
	return;
}//添加一个数
void del(int ix,long long sl,long long sr,long long x){
	if(sl!=sr){
		if(x<=(sl+sr)>>1){
			if(!ls[ix])ls[ix]=nw();
			del(ls[ix],sl,(sl+sr)>>1,x);
		}else{
			if(!rs[ix])rs[ix]=nw();
			del(rs[ix],((sl+sr)>>1)+1,sr,x);
		}v[ix]=v[ls[ix]]+v[rs[ix]];
		h[ix]=h[ls[ix]]*pw[v[rs[ix]]]+h[rs[ix]];
	}else v[ix]--,h[ix]=0;
	return;
}//删除一个数.
int main(){
	scanf("%lld%lld%lld",&n,&k,&s),d=1;
	for(int i=0;i<n;i++)scanf("%lld",&x),a[i]=n*(x-1)+i;
	for(int i=0;i<k;i++)scanf("%lld",&x),b[i]=k*(x-1)+i,r[i]=b[i];
	sort(r,r+k);
	for(int i=0;i<k;i++)b[i]=lower_bound(r,r+k,b[i])-r;
	for(int i=0;i<k;i++)r[b[i]]=i;
	for(int i=0;i<k;i++)hs=hs*31+r[i];
	pw[0]=1;for(int i=1;i<=k;i++)pw[i]=31*pw[i-1];
	for(int i=0;i<k;i++)br+=pw[i];//求S.
	for(int i=0;i<n;i++)r[i]=a[i];
	sort(r,r+n);
	for(int i=0;i<n;i++)a[i]=lower_bound(r,r+n,a[i])-r;
	for(int i=0;i<k-1;i++)upd(1,0,n-1,a[i],i);//先添加i<k-1的(a[i],i).
	for(int i=k-1;i<n;i++){
		upd(1,0,n-1,a[i],i);
		if(!((h[1]-(i-k+1)*br)^hs))ans[cnt]=i-k+2,cnt++;//!(a^b)就是a==b,但(仅限本人认为)常数小.
		del(1,0,n-1,a[i-k+1]);//一边添加一边删除.
	}printf("%d\n",cnt);
	for(int i=0;i<cnt;i++)printf("%d\n",ans[i]);
	return 0;
}