[字符串HASH && 阈值]Codeforces 741E.Arpa’s abnormal DNA and Mehrdad’s deep interest

传送门
这题写的心态都崩了…被卡常被卡HASH

首先可以把每个位置插入T后的字典序排序出来。
发现两个不同的位置i,j中插入T，两个字符串有一段前缀一段后缀是相同的，中间可以按T的端点分成3部分，只要比较这3部分的字典序就可以了。

我被卡常就是因为用了字符串hash来比较字典序…

那么写一个比较函数，然后直接用std::sort就好了

询问也很麻烦。
发现第二个限制可以写成$x+\lfloor{i\over k }\rfloor\times k\le i \le y+\lfloor{i\over k }\rfloor\times k$

这个形式就可以对k阈值，大于S的直接跳着做。
小于S的可以按k把询问分成 $\sqrt n$ 类，枚举k，再枚举模k的余数，用线段树记录区间最小值。枚举余数的部分处理的第二个限制，线段树处理的第一个限制。这样就可以了

好像不用双哈希，但是不能自然溢出

S取$\sqrt{{n\over\log n}}$复杂度最优

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <string>
#include <cstring>
#include <vector>
#include <cmath>
#include <ctime>
#define fi first
#define se second

using namespace std;

typedef unsigned long long ll;
typedef pair<ll,ll> pari;

const int N=100010;
const pari base={233,127};
const ll P=1e9+9;

int n,m,q,S,cnt,block,idx[N],lg2[N],rk[N],ans[N],mn[N<<2],tag[N<<2];
int st[N][20];
pari pw[N],hsa[N],hsb[N];
char a[N],b[N];

pari operator +(pari a,pari b){ return pari((a.fi+b.fi)%P,(a.se+b.se)%P); }
pari operator *(pari a,pari b){ return pari((a.fi*b.fi)%P,(a.se*b.se)%P); }
pari operator -(pari a,pari b){ return pari((a.fi+P-b.fi)%P,(a.se+P-b.se)%P); }
pari operator +(pari a,int b){ return pari((a.fi+b)%P,(a.se+b)%P); }
pari operator -(pari a,int b){ return pari((a.fi+P-b)%P,(a.se+P-b)%P); }

struct Qry{
  int l,r,x,y,g;
  Qry(int _l,int _r,int _x,int _y,int _g){l=_l;r=_r;x=_x;y=_y;g=_g;}
  Qry(){}
};
vector<Qry> Q[330];

inline int check(pari *a,pari *b,int a1,int b1,int a2,int b2){
  if(a1>b1) return 0;
  if(a[b1]-a[a1-1]*pw[b1-a1+1]==b[b2]-b[a2-1]*pw[b2-a2+1]) return 0;
  int l=1,r=b1-a1+1,mid,ret;
  while(l<=r){
    mid=l+r>>1;
    if(a[a1+mid-1]-a[a1-1]*pw[mid]==b[a2+mid-1]-b[a2-1]*pw[mid]) l=mid+1;
    else r=(ret=mid)-1;
  }
  return ret;
}

inline bool cmp(int x,int y){
  int ret=0,cur,nxt;
  if(x>y) swap(x,y),ret^=1;
  if(x+m>y){
    if(cur=check(hsa,hsb,x+1,y,1,y-x))
      return (b[cur]<a[x+cur])^ret;
    if(cur=check(hsb,hsb,y-x+1,m,1,m-y+x)) return (b[y-x+cur]<b[cur])^ret;
    if(cur=check(hsa,hsb,x+1,y,m-y+x+1,m)) return (a[x+cur]<b[m-y+x+cur])^ret;
    return ret^1;
  }
  else{
    if(cur=check(hsa,hsb,x+1,x+m,1,m)) return (b[cur]<a[x+cur])^ret;
    if(cur=check(hsa,hsa,x+1,y-m,x+m+1,y)) return (a[x+cur]<a[x+m+cur])^ret;
    if(cur=check(hsa,hsb,y-m+1,y,1,m)) return (a[y-m+cur]<b[cur])^ret;
    return ret^1;
  }
}

inline void Pre(){
  for(int i=2;i<=n;i++) lg2[i]=lg2[i-1]+((1<<lg2[i-1]+1)==i);
  int t=lg2[n];
  for(int i=0;i<=n;i++)
    st[i][0]=i,rk[idx[i]]=i;
  for(int k=1;k<=t;k++)
    for(int i=0;i+(1<<k)-1<=n;i++)
      if(rk[st[i][k-1]]<rk[st[i+(1<<k-1)][k-1]]) st[i][k]=st[i][k-1];
      else st[i][k]=st[i+(1<<k-1)][k-1];
}

inline int Query(int l,int r){
  int t=lg2[r-l+1];
  return rk[st[l][t]]<rk[st[r-(1<<t)+1][t]]?st[l][t]:st[r-(1<<t)+1][t];
}

inline void Push(int g){
  if(tag[g]==-1){
    mn[g<<1]=mn[g<<1|1]=-1; tag[g<<1]=tag[g<<1|1]=-1;
    tag[g]=0;
  }
}

inline void fix(int &x,int y){
  if(x==-1 || y==-1) x=x==-1?y:x;
  else if(rk[x]>rk[y]) x=y;
}

inline void Up(int g){
  mn[g]=-1;
  fix(mn[g],mn[g<<1]); fix(mn[g],mn[g<<1|1]);
}

void Modify(int g,int l,int r,int x){
  if(l==r){ mn[g]=x; return ; }
  int mid=l+r>>1; Push(g);
  if(x<=mid) Modify(g<<1,l,mid,x); else Modify(g<<1|1,mid+1,r,x);
  Up(g);
}

int Query(int g,int l,int r,int L,int R){
  if(l==L && r==R) return mn[g];
  int mid=L+R>>1; Push(g);
  if(r<=mid) return Query(g<<1,l,r,L,mid);
  else if(l>mid) return Query(g<<1|1,l,r,mid+1,R);
  int ret=Query(g<<1,l,mid,L,mid);
  fix(ret,Query(g<<1|1,mid+1,r,mid+1,R));
  return ret;
}

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}

inline void rea(char *a){
  char c=nc(); int len=0;
  for(;c>'z'||c<'a';c=nc());for(;c>='a'&&c<='z';a[++len]=c,c=nc());
  a[++len]=0;
}

inline void rea(int &x){
  char c=nc(); x=0;
  for(;c>'9'||c<'0';c=nc());for(;c>='0' && c<='9';x=x*10+c-'0',c=nc());
}

void print(int x){
  if(x>=10) print(x/10);
  putchar(x%10+'0');
}

int main(){
  rea(a); rea(b); rea(q);
  n=strlen(a+1); m=strlen(b+1);
  pw[0]={1,1}; for(int i=1;i<=max(n,m)+1;i++) pw[i]=pw[i-1]*base;
  for(int i=1;i<=n;i++)
    hsa[i]=hsa[i-1]*base+a[i]-'a';
  for(int i=1;i<=m;i++)
    hsb[i]=hsb[i-1]*base+b[i]-'a';
  for(int i=1;i<=n;i++) idx[i]=i;
  sort(idx,idx+1+n,cmp); Pre();
  S=max(0.8*sqrt(n/log2(n)),1.0);
  memset(ans,-1,sizeof(ans));
  for(int i=1;i<=q;i++){
    int l,r,k,x,y;
    rea(l); rea(r); rea(k); rea(x); rea(y);
    if(k<=S) Q[k].push_back(Qry(l,r,x,y,i));
    else{
      int cur=-1;
      for(;x<=r;x+=k,y+=k){
    if(y<l) continue; y=min(y,r);
    if(cur==-1 || rk[cur]>rk[Query(max(x,l),y)]) cur=Query(max(x,l),y);
      }
      ans[i]=cur;
    }
  }
  for(int i=1;i<=S;i++){
    if(Q[i].size()==0) continue;
    for(int j=0;j<i;j++){
      tag[1]=mn[1]=-1;
      for(int x=j;x<=n;x+=i) Modify(1,0,n,x);
      for(Qry x : Q[i]){
    if(j<x.x || j>x.y) continue;
    fix(ans[x.g],Query(1,x.l,x.r,0,n));
      }
    }
  }
  //printf("%d\n",clock());
  for(int i=1;i<=q;i++) printf("%d ",ans[i]);
  return 0;
}