Submission #876881

# Submission time Handle Problem Language Result Execution time Memory
876881 2023-11-22T12:55:08 Z Ice_man Jail (JOI22_jail) C++14
49 / 100
80 ms 62036 KB
#include <iostream>
#include <vector>
#include <cstring>
#include <stack>

#define maxn 120005
#define maxlog 20
#define pb(x) push_back(x)
#define control cout<<"passed"<<endl;

#pragma GCC optimize("O3" , "Ofast" , "unroll-loops")
#pragma GCC target(avx2)

using namespace std;

int n, m,  q;
vector <int> v[maxn];
int start[maxn], _final[maxn];
vector <int> new_tree[maxn];
int bin_lift[maxn][maxlog];
int pom1[maxn * 12];
int pom2[maxn * 12];
int chain[maxn];

void read()
{
    cin >> n;

    for(int i = 1; i <= n; i++)
    {
        v[i].clear();
        for(int j = 0; j < maxlog; j++) bin_lift[i][j] = 0;
    }

    for(int i = 1; i <= n * 12 - 1; i++)
    {
        pom1[i] = 0;
        pom2[i] = 0;
        ///v[i].clear();
        new_tree[i].clear();
    }

    int from, to;
    for(int i = 1; i <= (n - 1); i++)
    {
        cin >> from >> to;

        v[from].pb(to);
        v[to].pb(from);
    }

    cin >> m;

    for(int i = 1; i <= m; i++)
    {
        pom1[i] = 0;
        pom2[i] = 0;
    }

    for(int i = 1; i <= m; i++) cin >> start[i] >> _final[i];

}

int tree[12 * maxn];

void build_seg(int node , int l , int r)
{
    if(l == r)
    {
        tree[l] = node;
        return;
    }

    int mid = (r - l) / 2 + l;

    build_seg(node * 2 , l , mid);
    build_seg(node* 2 + 1 , mid + 1 , r);

    new_tree[node].pb(node * 2);
    new_tree[node].pb(node * 2 + 1);

    new_tree[node * 2 + 4 * n].pb(node + 4 * n);
    new_tree[node * 2 + 1 + 4 * n].pb(node + 4 * n);
}

void _add(int node , int l , int r , int ql , int qr , int val)
{

    if(ql <= l && r <= qr)
    {
        new_tree[node + n * 4].pb(val + n * 8);

        new_tree[val + n * 8].pb(node);

        return;
    }

    if(qr < l || r < ql || ql > qr) return;

    int mid = (r - l) / 2 + l;

    _add(node * 2 , l , mid , ql , qr , val);
    _add(node * 2 + 1 , mid + 1 , r , ql , qr , val);

}


///int bin_lift[maxn][maxlog];
int depth[maxn];
int subtree[maxn];
int help[maxn * 12];

void calc_bin(int node, int parent)
{
    
    subtree[node] = 1;
    
    for(int &nb : v[node])
    {
        if(nb != parent)
        {

            bin_lift[nb][0] = node;
            for(int i = 1; i < maxlog; i++) bin_lift[nb][i] = bin_lift[bin_lift[nb][i - 1]][i - 1];

            depth[nb] = depth[node] + 1;
            calc_bin(nb, node);

            subtree[node] += subtree[nb];

            if(subtree[v[node][0]] < subtree[nb] || v[node][0] == parent) swap(v[node][0] , nb);

        }
    }

}


int get_lca(int a, int b)
{
    if(depth[a] < depth[b]) swap(a, b);

    ///cout << a << "-" << b << endl;

    int levels = depth[a] - depth[b];
    ///cout << "** " << levels << endl;;
    for(int i = maxlog - 1; i >= 0; i--)
    {
        if(levels >> i & 1)
        {
            ///cout  << ">< "<< (1 << i) << endl;
            a = bin_lift[a][i];
            ///levels -= (1 << i);
            ///cout << (1 << i) << endl;
        }
    }

    ///cout << a << "-" << b << endl;
    if(a == b) return a;

    for(int i = maxlog - 1; i >= 0; i--)
    {
        if(!bin_lift[a][i]) continue;
        if(!bin_lift[b][i]) continue;

        if(bin_lift[a][i] != bin_lift[b][i])
        {
            a = bin_lift[a][i];
            b = bin_lift[b][i];
        }

    }

    return bin_lift[a][0];

}

int heads[maxn];

void add_edges(int node , int from , int to)
{
    bool lamp1 = true;
    bool lamp2 = true;

    while(heads[from] != heads[to])
    {
        ///cout << from << " " << to << endl;
        if(heads[from] == heads[to]) break;

        if(depth[heads[to]] > depth[heads[from]])
        {
            swap(from , to);
            swap(lamp1 , lamp2);
        }

        _add(1 , 1 , n , chain[heads[from]] , chain[from] - lamp1 , node);

        lamp1 = false;
        from = bin_lift[heads[from]][0];
    }


    if(chain[from] > chain[to])
    {
        swap(from , to);
        swap(lamp1 , lamp2);
    }

    _add(1 , 1 , n , chain[from] + lamp1 , chain[to] - lamp2 , node);

}


bool check(int a, int b, int c)
{
    int lca = get_lca(a, b);
    ///int pom = get_lca(lca , c);

    if(get_lca(lca, c) == lca && get_lca(c, a) == c) return true;
    if(get_lca(lca, c) == lca && get_lca(c, b) == c) return true;

    return false;

}

int cur_chain = 0;
///int heads[maxn];

void dfs_hld(int node , int parent , int leader)
{

    cur_chain++;

    chain[node] = cur_chain;
    heads[node] = leader;

    if(v[node][0] != parent) dfs_hld(v[node][0] , node , leader);

    for(int nb : v[node]) if(nb != parent && nb != v[node][0])
    {
        dfs_hld(nb , node , nb);
    }

}



//vector <int> new_tree[maxn];

void solve()
{

    calc_bin(1, 0);

    ///control

    /**for(int a = 1; a <= m; a++)
    {
        for(int b = 1; b <= m; b++)
        {
            if(a == b) continue;
            ///int lca1 = get_lca(start[a], _final[a]);

            if(check(start[a], _final[a], start[b]) == true) new_tree[b].pb(a);
            if(check(start[a], _final[a], _final[b]) == true) new_tree[a].pb(b);

            ///cout << check(start[a] , _final[a] , start[b]) << " " <<  check(start[a] , _final[a] , _final[b]) << endl;

        }
    }*/

    cur_chain = 0;
    ///_time = 0;
    dfs_hld(1 , 0 , 1);

    /**cout << "--------------" << endl;
    for(int i = 1; i <= n; i++) cout << heads[i] << " ";
    cout << endl << "--------------" << endl;*/


    ///control

    build_seg(1 , 1 , n);

    ///cout << "sd;inkmsad" << endl;
    ///control

    for(int i = 1; i <= n * 12; i++) help[i] = 0;

    for(int i = 1; i <= m; i++)
    {
        new_tree[tree[chain[_final[i]]]].pb(8 * n + i);
        new_tree[i + n * 8].pb(n * 4 + tree[chain[start[i]]]);

        help[_final[i]] = i;
        add_edges(i , start[i] , _final[i]);

        ///cout << "passed " << i << endl;

    }


    ///cout << "-----------------" << endl;
    ///control

    for(int i = 1; i <= m; i++) if(help[start[i]])
    {
        new_tree[i + 8 * n].pb(help[start[i]] + 8 * n);
    }

    /**for(int i = 1; i <= n; i++)
    {
       cout << i << ": ";
       for(int nb : new_tree[i]) cout << nb << " ";
       cout << endl;
    }*/

}

stack <int> s;
///int pom1[maxn] , pom2[maxn];
int _time = 0;
int br = 0;
bool ans = true;

void count_ssc(int node, int parent)
{
    s.push(node);
    _time++;

    pom1[node] = _time;
    pom2[node] = _time;

    for(int nb : new_tree[node])
    {
        if(pom2[nb] == -1) continue;
        if(pom2[nb])
        {
            pom1[node] = min(pom1[node], pom1[nb]);
        }
        else
        {
            count_ssc(nb, node);
            pom1[node] = min(pom1[node], pom1[nb]);
        }
    }

    int h;

    if(pom1[node] == pom2[node])
    {
        br++;

        h = 0;

        while(!s.empty() && s.top() != node)
        {
            pom2[s.top()] = -1;

            h += s.top() > 8 * n? 1 : 0;

            s.pop();
        }

        h += node > 8 * n? 1 : 0;

        pom2[node] = -1;
        if(h > 1) ans = false;

        s.pop();
    }


}

int _pom[maxn];

void combine()
{
    cin >> q;
    while(q--)
    {
        br = 0;
        _time = 0;
        ans = true;

        read();

        ///control

        solve();

        ///control

        /**for(int i = 1; i <= n; i++)
        {
            cout << i << ": ";
            for(int power = 0; power < maxlog; power++)
            {
                cout << bin_lift[i][power] << " ";
            }
            cout << endl;
        }*/

        ///cout << "-> " << get_lca(1 , 5) << endl;
        while(s.size()) s.pop();
        for(int i = 1; i <= 8 * n + m; i++) if(!pom2[i]) count_ssc(i, -1);

        if(ans == true) cout << "Yes" << endl;
        else cout << "No" << endl;

    }
}

int main()
{
    /**ios_base::sync_with_stdio(false);
    cin.tie(nullptr);*/

    combine();

    return 0;

}
/**

3
3
1 2
2 3
2
2 1
3 2
7
1 2
2 3
3 4
4 5
5 6
6 7
3
1 3
4 2
2 5
8
1 2
2 3
3 4
4 5
5 6
6 7
7 8
4
1 5
2 6
3 7
4 8


*/

Compilation message

jail.cpp:12:20: warning: '#pragma GCC option' is not a string [-Wpragmas]
   12 | #pragma GCC target(avx2)
      |                    ^~~~
# Verdict Execution time Memory Grader output
1 Correct 3 ms 18524 KB Output is correct
2 Correct 4 ms 18548 KB Output is correct
3 Correct 3 ms 18520 KB Output is correct
4 Correct 28 ms 18524 KB Output is correct
5 Correct 58 ms 18524 KB Output is correct
6 Correct 6 ms 18520 KB Output is correct
7 Correct 6 ms 18524 KB Output is correct
8 Correct 8 ms 18612 KB Output is correct
9 Correct 78 ms 21364 KB Output is correct
10 Runtime error 30 ms 62032 KB Execution killed with signal 11
11 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 3 ms 18524 KB Output is correct
2 Correct 3 ms 18524 KB Output is correct
3 Correct 6 ms 18524 KB Output is correct
4 Correct 6 ms 18524 KB Output is correct
5 Correct 5 ms 18524 KB Output is correct
6 Correct 5 ms 18524 KB Output is correct
7 Correct 5 ms 18524 KB Output is correct
8 Correct 6 ms 18524 KB Output is correct
9 Correct 6 ms 18524 KB Output is correct
10 Correct 5 ms 18524 KB Output is correct
11 Correct 5 ms 18524 KB Output is correct
12 Correct 4 ms 18524 KB Output is correct
13 Correct 4 ms 18568 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 18524 KB Output is correct
2 Correct 3 ms 18524 KB Output is correct
3 Correct 6 ms 18524 KB Output is correct
4 Correct 6 ms 18524 KB Output is correct
5 Correct 5 ms 18524 KB Output is correct
6 Correct 5 ms 18524 KB Output is correct
7 Correct 5 ms 18524 KB Output is correct
8 Correct 6 ms 18524 KB Output is correct
9 Correct 6 ms 18524 KB Output is correct
10 Correct 5 ms 18524 KB Output is correct
11 Correct 5 ms 18524 KB Output is correct
12 Correct 4 ms 18524 KB Output is correct
13 Correct 4 ms 18568 KB Output is correct
14 Correct 3 ms 18520 KB Output is correct
15 Correct 3 ms 18496 KB Output is correct
16 Correct 6 ms 18672 KB Output is correct
17 Correct 6 ms 18772 KB Output is correct
18 Correct 8 ms 18520 KB Output is correct
19 Correct 4 ms 18776 KB Output is correct
20 Correct 6 ms 18524 KB Output is correct
21 Correct 6 ms 18524 KB Output is correct
22 Correct 6 ms 18528 KB Output is correct
23 Correct 3 ms 18520 KB Output is correct
24 Correct 3 ms 18524 KB Output is correct
25 Correct 5 ms 18524 KB Output is correct
26 Correct 4 ms 18524 KB Output is correct
27 Correct 6 ms 18524 KB Output is correct
28 Correct 3 ms 18524 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 18524 KB Output is correct
2 Correct 3 ms 18524 KB Output is correct
3 Correct 6 ms 18524 KB Output is correct
4 Correct 6 ms 18524 KB Output is correct
5 Correct 5 ms 18524 KB Output is correct
6 Correct 5 ms 18524 KB Output is correct
7 Correct 5 ms 18524 KB Output is correct
8 Correct 6 ms 18524 KB Output is correct
9 Correct 6 ms 18524 KB Output is correct
10 Correct 5 ms 18524 KB Output is correct
11 Correct 5 ms 18524 KB Output is correct
12 Correct 4 ms 18524 KB Output is correct
13 Correct 4 ms 18568 KB Output is correct
14 Correct 3 ms 18520 KB Output is correct
15 Correct 3 ms 18496 KB Output is correct
16 Correct 6 ms 18672 KB Output is correct
17 Correct 6 ms 18772 KB Output is correct
18 Correct 8 ms 18520 KB Output is correct
19 Correct 4 ms 18776 KB Output is correct
20 Correct 6 ms 18524 KB Output is correct
21 Correct 6 ms 18524 KB Output is correct
22 Correct 6 ms 18528 KB Output is correct
23 Correct 3 ms 18520 KB Output is correct
24 Correct 3 ms 18524 KB Output is correct
25 Correct 5 ms 18524 KB Output is correct
26 Correct 4 ms 18524 KB Output is correct
27 Correct 6 ms 18524 KB Output is correct
28 Correct 3 ms 18524 KB Output is correct
29 Correct 7 ms 18688 KB Output is correct
30 Correct 8 ms 18532 KB Output is correct
31 Correct 6 ms 18524 KB Output is correct
32 Correct 6 ms 18524 KB Output is correct
33 Correct 6 ms 18608 KB Output is correct
34 Correct 7 ms 18524 KB Output is correct
35 Correct 6 ms 18524 KB Output is correct
36 Correct 6 ms 18524 KB Output is correct
37 Correct 5 ms 18524 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 3 ms 18524 KB Output is correct
2 Correct 3 ms 18524 KB Output is correct
3 Correct 6 ms 18524 KB Output is correct
4 Correct 6 ms 18524 KB Output is correct
5 Correct 5 ms 18524 KB Output is correct
6 Correct 5 ms 18524 KB Output is correct
7 Correct 5 ms 18524 KB Output is correct
8 Correct 6 ms 18524 KB Output is correct
9 Correct 6 ms 18524 KB Output is correct
10 Correct 5 ms 18524 KB Output is correct
11 Correct 5 ms 18524 KB Output is correct
12 Correct 4 ms 18524 KB Output is correct
13 Correct 4 ms 18568 KB Output is correct
14 Correct 3 ms 18520 KB Output is correct
15 Correct 3 ms 18496 KB Output is correct
16 Correct 6 ms 18672 KB Output is correct
17 Correct 6 ms 18772 KB Output is correct
18 Correct 8 ms 18520 KB Output is correct
19 Correct 4 ms 18776 KB Output is correct
20 Correct 6 ms 18524 KB Output is correct
21 Correct 6 ms 18524 KB Output is correct
22 Correct 6 ms 18528 KB Output is correct
23 Correct 3 ms 18520 KB Output is correct
24 Correct 3 ms 18524 KB Output is correct
25 Correct 5 ms 18524 KB Output is correct
26 Correct 4 ms 18524 KB Output is correct
27 Correct 6 ms 18524 KB Output is correct
28 Correct 3 ms 18524 KB Output is correct
29 Correct 7 ms 18688 KB Output is correct
30 Correct 8 ms 18532 KB Output is correct
31 Correct 6 ms 18524 KB Output is correct
32 Correct 6 ms 18524 KB Output is correct
33 Correct 6 ms 18608 KB Output is correct
34 Correct 7 ms 18524 KB Output is correct
35 Correct 6 ms 18524 KB Output is correct
36 Correct 6 ms 18524 KB Output is correct
37 Correct 5 ms 18524 KB Output is correct
38 Correct 78 ms 21264 KB Output is correct
39 Runtime error 33 ms 62036 KB Execution killed with signal 11
40 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 3 ms 18524 KB Output is correct
2 Correct 3 ms 18524 KB Output is correct
3 Correct 3 ms 18524 KB Output is correct
4 Correct 4 ms 18336 KB Output is correct
5 Correct 16 ms 18524 KB Output is correct
6 Correct 4 ms 18520 KB Output is correct
7 Correct 4 ms 18524 KB Output is correct
8 Correct 3 ms 18524 KB Output is correct
9 Correct 4 ms 18524 KB Output is correct
10 Correct 3 ms 18524 KB Output is correct
11 Correct 3 ms 18524 KB Output is correct
12 Correct 6 ms 18632 KB Output is correct
13 Correct 48 ms 18524 KB Output is correct
14 Correct 80 ms 18580 KB Output is correct
15 Correct 61 ms 18524 KB Output is correct
16 Runtime error 30 ms 62036 KB Execution killed with signal 11
17 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 3 ms 18524 KB Output is correct
2 Correct 4 ms 18548 KB Output is correct
3 Correct 3 ms 18520 KB Output is correct
4 Correct 28 ms 18524 KB Output is correct
5 Correct 58 ms 18524 KB Output is correct
6 Correct 6 ms 18520 KB Output is correct
7 Correct 6 ms 18524 KB Output is correct
8 Correct 8 ms 18612 KB Output is correct
9 Correct 78 ms 21364 KB Output is correct
10 Runtime error 30 ms 62032 KB Execution killed with signal 11
11 Halted 0 ms 0 KB -