#include using namespace std; typedef signed long long ll; #define _P(...) (void)printf(__VA_ARGS__) #define FOR(x,to) for(x=0;x<(to);x++) #define FORR(x,arr) for(auto& x:arr) #define FORR2(x,y,arr) for(auto& [x,y]:arr) #define ALL(a) (a.begin()),(a.end()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) template bool chmax(T &a, const T &b) { if(a bool chmin(T &a, const T &b) { if(a>b){a=b;return 1;}return 0;} //------------------------------------------------------- template struct FunctionGraph { int E[NV]; int N; int NG; int nex[61][NV]; // 2^i先の点 vector RE[NV]; // 逆向き辺 vector G[NV]; // 閉路 int GS[NV]; // 閉路長 int Gid[NV],Gindex[NV]; // 閉路IDとその中のindex。閉路外の場合、閉路IDは最寄りの閉路だがindexは-1 int GV[NV], D[NV]; //最寄りの閉路の点とそこまでの距離 void dfs(int cur,int id,int d) { GV[cur]=id; Gid[cur]=Gid[id]; D[cur]=d; FORR(e,RE[cur]) if(Gindex[e]==-1) dfs(e,id,d+1); } int move(int cur,ll step) { // step先に進む int i; FOR(i,60) if(step&(1LL< in(N); FOR(i,N) { RE[i].clear(); G[i].clear(); Gindex[i]=-2; nex[0][i]=E[i]; } NG=0; FOR(i,N) { RE[E[i]].push_back(i); in[E[i]]++; } queue Q; FOR(i,N) if(in[i]==0) Q.push(i); while(Q.size()) { int cur=Q.front(); Q.pop(); Gindex[cur]=-1; if(--in[E[cur]]==0) Q.push(E[cur]); } FOR(i,N) if(Gindex[i]==-2) { G[NG].push_back(i); while(1) { int x=G[NG].back(); Gid[x]=NG; Gindex[x]=G[NG].size()-1; x=E[x]; if(x==G[NG][0]) break; G[NG].push_back(x); } GS[NG]=G[NG].size(); NG++; } FOR(i,N) if(Gindex[i]>=0) dfs(i,i,0); FOR(j,60) FOR(i,N) nex[j+1][i]=nex[j][nex[j][i]]; } }; int N,L,R; FunctionGraph<101010> fg; map maxd; map did; int Q,S[101010],T[101010],B[101010]; int ret[101010]; vector cand[101010]; int dp[101010]; int tar[101010]; int cur[101010]; void solve() { int i,j,k,l,r,x,y; string s; cin>>N>>L>>R; FOR(i,N) { cin>>fg.E[i]; fg.E[i]--; } fg.build(N); cin>>Q; FOR(i,Q) { cin>>S[i]>>T[i]; S[i]--,T[i]--; ret[i]=-1; if(fg.Gid[S[i]]!=fg.Gid[T[i]] || fg.D[S[i]]fg.Gindex[T[i]]) B[i]+=fg.GS[fg.Gid[S[i]]]; cand[fg.GS[fg.Gid[S[i]]]].push_back(i); } } } FOR(i,N+1) if(cand[i].size()) { set alive,alivep; FOR(j,N+1) { dp[j]=-1; tar[j%i]=j; cur[j]=-1; if(j) { alive.insert(j); alivep.insert(j%i); alivep.insert(j%i+i); } } dp[0]=0; deque Q; Q.push_front(0); while(Q.size()) { int cur=Q.front(); Q.pop_front(); if(cur>=i&&(dp[cur-i]>dp[cur]||dp[cur-i]==-1)) { dp[cur-i]=dp[cur]; alive.erase(cur-i); Q.push_front(cur-i); } auto it=alive.lower_bound(cur+L); while(it!=alive.end()&&*it<=cur+R) { dp[*it]=dp[cur]+1; Q.push_back(*it); it=alive.erase(it); } int CL=cur+L; int CR=cur+R+1; if(CL<=N) CL=N+1; if(CR-CL>=i) { FORR(j,alivep) if(j=CR) break; y=*it%i; for(x=y;x<=N;x+=i) if(alive.count(x)) { dp[x]=dp[cur]+1; Q.push_back(x); alive.erase(x); } alivep.erase(y); alivep.erase(y+i); } } } FORR(q,cand[i]) { ret[q]=dp[B[q]]; } } FOR(i,Q) cout<