Submission #949724

# Submission time Handle Problem Language Result Execution time Memory
949724 2024-03-19T15:26:26 Z andrei_boaca Tourism (JOI23_tourism) C++17
5 / 100
4196 ms 56616 KB
#include <bits/stdc++.h>

using namespace std;
typedef pair<int,int> pii;
int bucksize=320;
int n,m,q,c[100005],niv[100005],in[100005],out[100005],lin[500005],rmq[21][500005],loga[500005],timp,par[100005];
int sol[100005],nivmin,cnt;
set<pii> s;
vector<int> muchii[100005];
void dfs(int nod)
{
    timp++;
    lin[timp]=nod;
    in[nod]=timp;
    bool have=0;
    for(int i:muchii[nod])
        if(i!=par[nod])
        {
            par[i]=nod;
            niv[i]=niv[nod]+1;
            dfs(i);
            timp++;
            lin[timp]=nod;
            have=1;
        }
    if(!have)
    {
        timp++;
        lin[timp]=nod;
    }
    out[nod]=timp;
}
struct date
{
    int l,r,ind;
};
int getbuck(int x)
{
    return x/bucksize+(x%bucksize!=0);
}
bool byr(date a, date b)
{
    return a.r<b.r;
}
vector<date> qr[100005];
bool isancestor(int a,int b)
{
    return in[a]<=in[b]&&out[a]>=out[b];
}
int getmin(int l,int r)
{
    int lg=loga[r-l+1];
    int a=rmq[lg][r-(1<<lg)+1];
    int b=rmq[lg][l];
    if(niv[a]<niv[b])
        return a;
    return b;
}
int LCA(int a,int b)
{
    if(in[a]>in[b])
        swap(a,b);
    if(isancestor(a,b))
        return a;
    int l=out[a];
    int r=in[b];
    return getmin(l,r);
}
void addnode(int nod)
{
    if(s.find({in[nod],nod})!=s.end())
        return;
    s.insert({in[nod],nod});
    if(s.size()==1)
    {
        nivmin=niv[nod];
        cnt=nivmin+1;
    }
    int nmin=1e9;
    int nmax=0;
    auto it=s.find({in[nod],nod});
    bool have=0;
    if(it!=s.begin())
    {
        int st=(*prev(it)).second;
        int x=LCA(st,nod);
        nmin=min(nmin,niv[x]);
        nmax=max(nmax,niv[x]);
    }
    it++;
    if(it!=s.end())
    {
        int dr=(*it).second;
        if(isancestor(nod,dr))
            have=1;
        int x=LCA(nod,dr);
        nmin=min(nmin,niv[x]);
        nmax=max(nmax,niv[x]);
    }
    if(nivmin>nmin)
    {
        if(!have)
            cnt+=niv[nod]-nmin;
        nivmin=min(nivmin,nmin);
    }
    else
    {
        if(!have)
            cnt+=niv[nod]-nmax;
    }
}
int main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(0);
    cin>>n>>m>>q;
    for(int i=1;i<n;i++)
    {
        int a,b;
        cin>>a>>b;
        muchii[a].push_back(b);
        muchii[b].push_back(a);
    }
    for(int i=1;i<=m;i++)
        cin>>c[i];
    dfs(1);
    for(int i=1;i<=timp;i++)
    {
        for(int bit=0;bit<=20;bit++)
            if((1<<bit)<=i)
                loga[i]=bit;
    }
    for(int i=1;i<=timp;i++)
        rmq[0][i]=lin[i];
    for(int j=1;j<=18;j++)
        for(int i=1;i<=timp;i++)
        {
            rmq[j][i]=rmq[j-1][i];
            int poz=i+(1<<(j-1));
            if(poz<=timp&&niv[rmq[j-1][poz]]<niv[rmq[j][i]])
                rmq[j][i]=rmq[j-1][poz];
        }
    int bmax=0;
    for(int i=1;i<=q;i++)
    {
        int l,r;
        cin>>l>>r;
        int b=getbuck(l);
        bmax=max(bmax,b);
        qr[b].push_back({l,r,i});
    }
    for(int b=1;b<=bmax;b++)
        if(!qr[b].empty())
        {
            sort(qr[b].begin(),qr[b].end(),byr);
            reverse(qr[b].begin(),qr[b].end());
            while(!qr[b].empty()&&getbuck(qr[b].back().r)==b)
            {
                int l=qr[b].back().l;
                int r=qr[b].back().r;
                int ind=qr[b].back().ind;
                qr[b].pop_back();
                s.clear();
                nivmin=niv[c[l]];
                cnt=nivmin+1;
                s.insert({in[c[l]],c[l]});
                for(int i=l+1;i<=r;i++)
                {
                    int nod=c[i];
                    addnode(nod);
                }
                sol[ind]=cnt-nivmin;
            }
            s.clear();
            int poz=(b-1)*bucksize;
            cnt=0;
            nivmin=0;
            while(!qr[b].empty())
            {
                int l=qr[b].back().l;
                int r=qr[b].back().r;
                int ind=qr[b].back().ind;
                qr[b].pop_back();
                while(poz<r)
                {
                    poz++;
                    int nod=c[poz];
                    addnode(nod);
                }
                int last=(b-1)*bucksize;
                int auxniv=nivmin;
                int auxcnt=cnt;
                vector<pii> toerase;
                for(int i=l;i<=last;i++)
                {
                    int nod=c[i];
                    if(s.find({in[nod],nod})!=s.end())
                        continue;
                    toerase.push_back({in[nod],nod});
                    addnode(nod);
                }
                for(pii p:toerase)
                    s.erase(p);
                sol[ind]=cnt-nivmin;
            }
        }
    for(int i=1;i<=q;i++)
        cout<<sol[i]<<'\n';
    return 0;
}

Compilation message

tourism.cpp: In function 'int main()':
tourism.cpp:191:21: warning: unused variable 'auxniv' [-Wunused-variable]
  191 |                 int auxniv=nivmin;
      |                     ^~~~~~
tourism.cpp:192:21: warning: unused variable 'auxcnt' [-Wunused-variable]
  192 |                 int auxcnt=cnt;
      |                     ^~~~~~
# Verdict Execution time Memory Grader output
1 Correct 6 ms 47704 KB Output is correct
2 Correct 6 ms 47708 KB Output is correct
3 Correct 10 ms 43864 KB Output is correct
4 Correct 10 ms 45660 KB Output is correct
5 Correct 9 ms 43612 KB Output is correct
6 Correct 9 ms 45660 KB Output is correct
7 Correct 8 ms 43612 KB Output is correct
8 Correct 8 ms 45660 KB Output is correct
9 Correct 10 ms 43728 KB Output is correct
10 Correct 10 ms 43612 KB Output is correct
11 Correct 12 ms 47708 KB Output is correct
12 Correct 16 ms 45804 KB Output is correct
13 Correct 17 ms 45660 KB Output is correct
14 Correct 17 ms 47708 KB Output is correct
15 Correct 11 ms 45660 KB Output is correct
16 Correct 9 ms 43612 KB Output is correct
17 Correct 9 ms 45836 KB Output is correct
18 Correct 12 ms 43612 KB Output is correct
19 Correct 10 ms 45660 KB Output is correct
20 Correct 10 ms 45660 KB Output is correct
21 Correct 10 ms 43608 KB Output is correct
22 Correct 10 ms 43608 KB Output is correct
23 Correct 10 ms 45820 KB Output is correct
24 Correct 10 ms 45660 KB Output is correct
25 Correct 13 ms 47704 KB Output is correct
26 Correct 12 ms 47704 KB Output is correct
27 Correct 6 ms 43608 KB Output is correct
28 Correct 6 ms 45660 KB Output is correct
29 Correct 7 ms 45660 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 6 ms 47704 KB Output is correct
2 Correct 6 ms 47708 KB Output is correct
3 Correct 10 ms 43864 KB Output is correct
4 Correct 10 ms 45660 KB Output is correct
5 Correct 9 ms 43612 KB Output is correct
6 Correct 9 ms 45660 KB Output is correct
7 Correct 8 ms 43612 KB Output is correct
8 Correct 8 ms 45660 KB Output is correct
9 Correct 10 ms 43728 KB Output is correct
10 Correct 10 ms 43612 KB Output is correct
11 Correct 12 ms 47708 KB Output is correct
12 Correct 16 ms 45804 KB Output is correct
13 Correct 17 ms 45660 KB Output is correct
14 Correct 17 ms 47708 KB Output is correct
15 Correct 11 ms 45660 KB Output is correct
16 Correct 9 ms 43612 KB Output is correct
17 Correct 9 ms 45836 KB Output is correct
18 Correct 12 ms 43612 KB Output is correct
19 Correct 10 ms 45660 KB Output is correct
20 Correct 10 ms 45660 KB Output is correct
21 Correct 10 ms 43608 KB Output is correct
22 Correct 10 ms 43608 KB Output is correct
23 Correct 10 ms 45820 KB Output is correct
24 Correct 10 ms 45660 KB Output is correct
25 Correct 13 ms 47704 KB Output is correct
26 Correct 12 ms 47704 KB Output is correct
27 Correct 6 ms 43608 KB Output is correct
28 Correct 6 ms 45660 KB Output is correct
29 Correct 7 ms 45660 KB Output is correct
30 Incorrect 12 ms 43892 KB Output isn't correct
31 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 6 ms 45660 KB Output is correct
2 Correct 7 ms 45732 KB Output is correct
3 Correct 7 ms 45660 KB Output is correct
4 Incorrect 3826 ms 56616 KB Output isn't correct
5 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 6 ms 47704 KB Output is correct
2 Incorrect 58 ms 50620 KB Output isn't correct
3 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 7 ms 47704 KB Output is correct
2 Correct 7 ms 47704 KB Output is correct
3 Correct 7 ms 47708 KB Output is correct
4 Incorrect 4196 ms 55488 KB Output isn't correct
5 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 6 ms 47704 KB Output is correct
2 Correct 6 ms 47708 KB Output is correct
3 Correct 10 ms 43864 KB Output is correct
4 Correct 10 ms 45660 KB Output is correct
5 Correct 9 ms 43612 KB Output is correct
6 Correct 9 ms 45660 KB Output is correct
7 Correct 8 ms 43612 KB Output is correct
8 Correct 8 ms 45660 KB Output is correct
9 Correct 10 ms 43728 KB Output is correct
10 Correct 10 ms 43612 KB Output is correct
11 Correct 12 ms 47708 KB Output is correct
12 Correct 16 ms 45804 KB Output is correct
13 Correct 17 ms 45660 KB Output is correct
14 Correct 17 ms 47708 KB Output is correct
15 Correct 11 ms 45660 KB Output is correct
16 Correct 9 ms 43612 KB Output is correct
17 Correct 9 ms 45836 KB Output is correct
18 Correct 12 ms 43612 KB Output is correct
19 Correct 10 ms 45660 KB Output is correct
20 Correct 10 ms 45660 KB Output is correct
21 Correct 10 ms 43608 KB Output is correct
22 Correct 10 ms 43608 KB Output is correct
23 Correct 10 ms 45820 KB Output is correct
24 Correct 10 ms 45660 KB Output is correct
25 Correct 13 ms 47704 KB Output is correct
26 Correct 12 ms 47704 KB Output is correct
27 Correct 6 ms 43608 KB Output is correct
28 Correct 6 ms 45660 KB Output is correct
29 Correct 7 ms 45660 KB Output is correct
30 Incorrect 12 ms 43892 KB Output isn't correct
31 Halted 0 ms 0 KB -