제출 #1013795

#	제출 시각	아이디	문제	언어	결과	실행 시간	메모리
1013795		bobbilyking	가장 긴 여행 (IOI23_longesttrip)	C++17	0 / 100	212 ms	344 KiB

longesttrip

#include "longesttrip.h"
using namespace std;
#include <bits/stdc++.h>

#define F(i, l, r) for (int i = (l); i < (r); ++i)
#define A(a) (a).begin(), (a).end()

int edge[300][300];
int parent[300];
int find(int i) {
    return i == parent[i] ? i : parent[i] = find(parent[i]);
}

int p[300];

bool query(vector<int> A, vector<int> B) {
    // for (auto x: A) cout << x << endl;
    // for (auto x: B) cout << x << endl;
    
    if (A.size() == 1 and B.size() == 1) {
        if (edge[A[0]][B[0]]) return 1;
    }
    vector<int> OA = A, OB = B;
    for (auto &x: A) x = p[x];
    for (auto &x: B) x = p[x];

    bool res = are_connected(A, B);

    if (res) {
        if (A.size() == 1 and B.size() == 1) {
            auto x = OA[0], y = OB[0];
            edge[x][y] = edge[y][x] = 1;
            parent[find(x)] = find(y);
        }
    }

    // cout << "FUCK " << endl;

    return res; 
}

vector<int> longest_trip(int N, int D)
{
    // random_device rd;
    mt19937_64 rng(39);
    iota(p, p + N, 0ll);
    iota(parent, parent + N, 0ll);
    shuffle(p, p + N, rng);
    D = 1;

    vector<int> hamil;

    auto gen = [&]() {
        set<int> v;
        vector<bool> seen(N+1);
        for (auto x: hamil) seen[x] = 1;
        F(i, 0, N) if (!seen[i]) v.insert(i);
        return v;
    };

    if (query({0}, {1})) {
        hamil = {0, 1};
    } else if (query({1}, {2})) {
        hamil = {1, 2};
    } else {
        hamil = {0, 2};
    }
    {
        int x = hamil[0], y=  hamil[1];
        edge[x][y] = edge[y][x] = 1;
        parent[x] = parent[y] = y;
    }

    vector<int> prev = hamil;
    // for (auto x: hamil) cout << x << " "; cout << endl;

    // cout << "FUCKY" << endl;

    while (hamil.size() != N) { // this will reach N, with the exception of one case 
        bool can_do_any = hamil.size() == 2;

        // cout << "ITERS " << endl;

        if (!can_do_any) {

            // cout << "HOW " << endl;
            // cout << hamil[0] << " " << hamil.back() << endl;
            // for (auto x: prev) cout << x << " "; cout << endl;

            if (query({hamil[0]}, {hamil.back()})) {
                can_do_any = true;
            } else {
                // Pick any node not in hamil, this must be connected to either the front or the back.
                auto i = *gen().begin();

                // cout << "HAHA " << endl;
                if (query({hamil[0]}, {i})) hamil.insert(hamil.begin(), i);
                else hamil.push_back(i);
                continue;
            }
        }

        // Okay, we want to find *any* connection between hamil and others, if exists already. 
        auto bad = gen();

        bool seen_edge = false;
        F(i, 0, hamil.size()) for (auto y: bad) if (edge[hamil[i]][y]) {
            hamil.insert(hamil.begin() + i, y);
            // for (auto x: hamil) cout << x << " "; cout << endl;

            // we dont do this if it came from a 2 array tho 
            if (hamil.size() != 3) while (hamil.back() != y) { // really dumb rotate 
                auto v = hamil.back(); hamil.pop_back(); hamil.insert(hamil.begin(), v);
            }
            seen_edge = true; 
            
            // cout << "GOTO END " << endl;
            goto end;
        }
        end: if (seen_edge) continue;

        // Okay, now we want to do the "pick one from end, and two disjoint components" strategy.
        while (true) {
            int state = 0;
            // cout << "WHAT " << endl;
            for (auto x: bad) for (auto y: bad) if (find(x) != find(y)) {
                if (query({x}, {y})) {
                    state = 1; goto end2;
                } else { // one MUST be connected to one of the endpoints
                    state = 2; 
                    if (query({x}, {hamil[0]})) {
                        hamil.insert(hamil.begin(), x);
                    } else {
                        hamil.push_back(x);
                    }
                    goto end2;
                }
            }
            end2:
            if (state == 1) continue;
            if (state == 2) break; 
            // OTHERWISE, RARE ONCE IN AN ITERATION CASE: WE ASK IF THERE'S ANY CONNECTION ACROSS BOTH COMPS.

            // IF NOT, THEN INSTANTLY KNOW TWO K_N COMPS.

            if (!query(hamil, vector<int>(A(bad)))) {
                if (hamil.size() < bad.size()) {
                    hamil = vector<int>(A(bad));
                }
                return hamil;
            }
            
            // ELSE, FIND ANY SPANNING EDGE WITH BSEARCH, TAKE IT. BREAK, WE WILL FIND THE SPANNING EDGE AGAIN NEXT LOOP. 

            vector<int> cur(A(bad));
            while (cur.size() > 1) {
                vector<int> shit(cur.begin(), cur.begin() + cur.size()/2);
                if (query(hamil, shit)) cur = shit;
                else cur = vector<int>(cur.begin() + cur.size() / 2, cur.end());
            }

            int x = cur[0];

            cur = hamil;
            while (cur.size() > 1) {
                vector<int> shit(cur.begin(), cur.begin() + cur.size()/2);
                if (query({x}, shit)) cur = shit;
                else cur = vector<int>(cur.begin() + cur.size() / 2, cur.end());
            }

            int y = cur[0];

            edge[x][y] = edge[y][x] = 1;
        }
    }


    return hamil;
}

컴파일 시 표준 에러 (stderr) 메시지

longesttrip.cpp: In function 'std::vector<int> longest_trip(int, int)':
longesttrip.cpp:79:25: warning: comparison of integer expressions of different signedness: 'std::vector<int>::size_type' {aka 'long unsigned int'} and 'int' [-Wsign-compare]
   79 |     while (hamil.size() != N) { // this will reach N, with the exception of one case
      |            ~~~~~~~~~~~~~^~~~
longesttrip.cpp:5:40: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
    5 | #define F(i, l, r) for (int i = (l); i < (r); ++i)
      |                                        ^
longesttrip.cpp:107:9: note: in expansion of macro 'F'
  107 |         F(i, 0, hamil.size()) for (auto y: bad) if (edge[hamil[i]][y]) {
      |         ^

• Subtask 0: 0 / 0
• Subtask 1: 0 / 5
• Subtask 2: 0 / 10
• Subtask 3: 0 / 25
• Subtask 4: 0 / 60