传送门
这题写的心态都崩了…被卡常被卡HASH
首先可以把每个位置插入T后的字典序排序出来。
发现两个不同的位置i,j中插入T,两个字符串有一段前缀一段后缀是相同的,中间可以按T的端点分成3部分,只要比较这3部分的字典序就可以了。
我被卡常就是因为用了字符串hash来比较字典序…
那么写一个比较函数,然后直接用std::sort就好了
询问也很麻烦。
发现第二个限制可以写成
x+⌊ik⌋×k≤i≤y+⌊ik⌋×k
这个形式就可以对k阈值,大于S的直接跳着做。
小于S的可以按k把询问分成
n√
类,枚举k,再枚举模k的余数,用线段树记录区间最小值。枚举余数的部分处理的第二个限制,线段树处理的第一个限制。这样就可以了
好像不用双哈希,但是不能自然溢出
S取 nlogn−−−−√ 复杂度最优
#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;
}