結果

問題 No.2704 L to R Graph
ユーザー maksimmaksim
提出日時 2024-03-29 22:43:49
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 202 ms / 3,000 ms
コード長 4,460 bytes
コンパイル時間 2,450 ms
コンパイル使用メモリ 191,724 KB
実行使用メモリ 41,216 KB
最終ジャッジ日時 2024-09-30 16:35:10
合計ジャッジ時間 7,694 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 7 ms
9,856 KB
testcase_01 AC 6 ms
9,728 KB
testcase_02 AC 6 ms
9,856 KB
testcase_03 AC 6 ms
9,856 KB
testcase_04 AC 7 ms
9,856 KB
testcase_05 AC 7 ms
9,856 KB
testcase_06 AC 7 ms
9,856 KB
testcase_07 AC 7 ms
9,728 KB
testcase_08 AC 118 ms
34,048 KB
testcase_09 AC 115 ms
34,048 KB
testcase_10 AC 114 ms
34,048 KB
testcase_11 AC 196 ms
39,808 KB
testcase_12 AC 177 ms
36,352 KB
testcase_13 AC 192 ms
40,448 KB
testcase_14 AC 183 ms
37,120 KB
testcase_15 AC 176 ms
37,120 KB
testcase_16 AC 181 ms
37,120 KB
testcase_17 AC 176 ms
36,992 KB
testcase_18 AC 178 ms
37,632 KB
testcase_19 AC 202 ms
40,064 KB
testcase_20 AC 181 ms
37,120 KB
testcase_21 AC 183 ms
41,216 KB
testcase_22 AC 163 ms
41,216 KB
testcase_23 AC 171 ms
40,832 KB
testcase_24 AC 173 ms
40,832 KB
testcase_25 AC 179 ms
36,096 KB
testcase_26 AC 180 ms
36,096 KB
testcase_27 AC 93 ms
34,944 KB
testcase_28 AC 96 ms
34,944 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
#define int long long
#define app push_back
#define all(x) (x).begin(),(x).end()
#ifdef LOCAL
#define debug(...) [](auto...a){ ((cout << a << ' '), ...) << endl;}(#__VA_ARGS__, ":", __VA_ARGS__)
#else
#define debug(...)
#endif
#ifdef LOCAL
#define __int128 long long
#endif // LOCAL
const int inf=1e18;
const int maxn=2e5+5;
int a[maxn];int col[maxn];bool used[maxn];int h[maxn];int de[maxn];int u[maxn][20];bool cy[maxn];int di[maxn];
int me[maxn];
int mo[maxn];
int f[maxn];
vector<int> ps;
vector<int> mem[maxn];
int32_t main()
{
    ios_base::sync_with_stdio(false);cin.tie(0);cout.tie(0);
    int n,l,r;cin>>n>>l>>r;
    set<int> unu;
    f[0]=0;fill(f+1,f+maxn,inf);
    for(int i=1;i<=n;++i) unu.insert(i);
    queue<int> q;q.push(0);
    while(!q.empty())
    {
        int x=q.front();q.pop();
        set<int>::iterator it=unu.lower_bound(x+l);
        while(it!=unu.end() && (*it)<=x+r)
        {
            f[*it]=f[x]+1;q.push(*it);
            unu.erase(it);
            it=unu.lower_bound(x+l);
        }
    }
    for(int i=1;i<=n;++i) {cin>>a[i];u[i][0]=a[i];}
    for(int j=1;j<20;++j) {for(int i=1;i<=n;++i) {u[i][j]=u[u[i][j-1]][j-1];}}
    for(int i=1;i<=n;++i) {cy[u[i][19]]=true;}
    for(int i=1;i<=n;++i)
    {
        if(cy[i] && !used[i])
        {
            int x=i;int c=0;
            while(true)
            {
                col[x]=i;
                used[x]=true;
                h[x]=c;
                ++c;
                x=a[x];
                if(x==i) break;
            }
            de[i]=c;
            ps.app(de[i]);
        }
    }
    for(int i=1;i<=n;++i)
    {
        if(cy[i]) {me[i]=i;continue;}
        di[i]=0;
        int x=i;
        for(int j=19;j>=0;--j)
        {
            if(!cy[u[x][j]])
            {
                di[i]+=(1<<j);x=u[x][j];
            }
        }
        me[i]=a[x];
        di[i]++;
        mo[de[col[me[i]]]]=max(mo[de[col[me[i]]]],di[i]);
    }
    sort(all(ps));ps.erase(unique(all(ps)),ps.end());
    auto ok=[&](int p,int x)->bool
    {
        return ((l+2*p-1-x)/p!=(r+2*p-x)/p);
    };
    for(int p:ps)
    {
        vector<int> ans(mo[p]+p+1);fill(all(ans),inf);ans[0]=0;
        set<int> unuse;for(int i=1;i<ans.size();++i)  {unuse.insert(i);unuse.insert(i+p);}
        queue<int> q;q.push(0);
        while(!q.empty())
        {
            int x=q.front();q.pop();
            int wis=x+l;for(int i=20;i>=0;--i) {if(wis-p*(1<<i)>=((int) (ans.size()))) wis-=p*(1<<i);}
            if(wis>=ans.size()) wis-=p;
            debug(x,wis);
            set<int>::iterator it=unuse.lower_bound(wis);
            while(it!=unuse.end() && (*it)<=x+r && ok(p,(*it)-x))
            {
                int v=(*it);
                unuse.erase(it);
                if(v>=ans.size()) {v-=p;}
                if(ans[v]==inf) {ans[v]=ans[x]+1;q.push(v);}
                if(unuse.count(v)) unuse.erase(v);
                it=unuse.lower_bound(wis);
            }
        }
        for(int i=ans.size()-1-p;i>=0;--i) {ans[i]=min(ans[i],ans[i+p]);}
        debug(ans.size());
        for(int x:ans) debug(x);
        mem[p]=ans;
    }
    auto up=[&](int x,int h)
    {
        for(int i=19;i>=0;--i)
        {
            if(h & (1<<i))
            {
                x=u[x][i];
            }
        }
        return x;
    };
    int q1;cin>>q1;
    while(q1--)
    {
        int x,y;cin>>x>>y;
        if(col[me[x]]!=col[me[y]])
        {
            cout<<(-1)<<'\n';
            continue;
        }
        if(!cy[y])
        {
            if(di[x]<di[y])
            {
                cout<<(-1)<<'\n';
                continue;
            }
            if(x==y)
            {
                cout<<0<<'\n';continue;
            }
            int h=di[x]-di[y];
            if(up(x,h)!=y)
            {
                cout<<(-1)<<'\n';continue;
            }
            cout<<(f[h]==inf ? -1 : f[h])<<'\n';
            continue;
        }
        int p=de[col[me[x]]];
        int o=((h[y]-h[me[x]])%p+p)%p;
        o+=di[x];
        if(o<mem[p].size())
        {
            cout<<(mem[p][o]==inf ? -1 : mem[p][o])<<'\n';
        }
        else
        {
            for(int i=20;i>=0;--i)
            {
                if(o-p*(1<<i)>=((int) (mem[p].size()))) {o-=p*(1<<i);}
            }
            o-=p;
            cout<<(mem[p][o]==inf ? -1 : mem[p][o])<<'\n';
        }
    }
    return 0;
}
/*
2 101 101
2 1
1
1 2
*/
0