Submission #1154706

#	Time	Username	Problem	Language	Result	Execution time	Memory
1154706		tapilyoca	Cyberland (APIO23_cyberland)	C++20	8 / 100	31 ms	36676 KiB

cyberland

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

using ll = long long;
using vll = vector<ll>;
using ld = long double;


class comp{
    public:
        bool operator()(pair<pair<ll,ll>,ll> a, pair<pair<ll,ll>,ll> b){
            return a.first.first > b.first.first;
        }
};


double solve(int N, int M, int K, int H, std::vector<int> x, std::vector<int> y, std::vector<int> c, std::vector<int> arr) {
    
    vector<vector<pair<ll,ll>>> adj(N+1); // neigjbor, weight
    for(int i = 0; i < M; i++){
        ll a = x[i];
        ll b = y[i];
        ll w = c[i];
        adj[a].push_back({b,w});
        adj[b].push_back({a,w});
    }

    
    vector<bool> vis(N+1,0);
    vis[0] = 1;
    queue<ll> BFS_q;
    BFS_q.push(0);
    while(!BFS_q.empty()){
        ll at = BFS_q.front();
        BFS_q.pop();
        for(pair<ll,ll> p : adj[at]){
            ll node = p.first;
            if(vis[node]) continue;
            vis[node]  =1;
            if(node == H)
                continue;

            BFS_q.push(node);
        }
    }

    if(!vis[H]){
        cerr << "NOT EVEN CONNECTED " << endl;
        return -1;
    }

    // for(int i = 0; i < N; i++){
    //     for(int j = 0; j < M; j++){
    //         cout << adj[i][j].first << " ";
    //     } cout << endl;
    // }

    // dp-esque?
    // note that if K too big (over 30??) then we are within range lol
    K = min(30,K);

    vector<vector<ld>> dists(N+1, vector<ld>(K+1,-1));
    dists[H][0] = 0;

    // djk
    // store: dist, k, node
    priority_queue<pair<pair<ld,ll>,ll>, vector<pair<pair<ld,ll>,ll>>, comp> pq;

    pq.push({{ld(0),0},H});


    // suppose closest node is node i...
    // then let its neighbor be j
    // arr[j] = 1 -> normal, check dist[i] + weight(i,j)
    // arr[j] = 0 -> this is a starting point, so return 
    // arr[j] = 2 -> weird, check (dist[i] + weight(i,j))/2 if k is not maxed out

    // normal djk afterwards?

    while(!pq.empty()){
        pair<pair<ld,ll>,ll> front = pq.top();
        pq.pop();

        ll node = front.second;
        ld curDist = front.first.first;
        ll k    = front.first.second;
            
        for(pair<ll,ll> p : adj[node]){
            ll neighbor = p.first;
            if(!vis[neighbor]) continue;

            ld weight = ld(p.second);
            ld oldDist = dists[neighbor][k];

            if((arr[neighbor] != 2) or k == K or neighbor == 0){
                // normal case... so just do djk...
                if(curDist + weight < oldDist or oldDist == -1){
                    dists[neighbor][k] = curDist + weight;
                    pq.push({{dists[neighbor][k], k}, neighbor});
                }
            }

            if(arr[neighbor] == 0 or neighbor == 0){
                return dists[neighbor][k];
            }

            if(arr[neighbor] == 2 and k < K){
                if((curDist+weight)/ld(2) < oldDist or oldDist == -1){
                    dists[neighbor][++k] = (curDist + weight)/ld(2);
                    pq.push({{dists[neighbor][k], k}, neighbor});
                }
            }
        }
    }

    assert(false); // should never run?
}

• Subtask 1: 0 / 5
• Subtask 2: 8 / 8
• Subtask 3: 0 / 13
• Subtask 4: 0 / 19
• Subtask 5: 0 / 7
• Subtask 6: 0 / 16
• Subtask 7: 0 / 29
• Subtask 8: 0 / 3