53302023-04-25 22:29:09CattVázsony vonatjegyet vásárolcpp17Time limit exceeded 55/1001.567s156084 KiB
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const ll maxN = 2e5 + 5;


vector<vector<pair<ll, ll> > > g;
vector<ll> dist;
vector<ll> p(maxN), s(maxN, 1);
vector<set<ll> > to(maxN);
vector<bool> volt(maxN, 0);

ll root(ll x) {
    if(p[x] == x) return x;

    return p[x] = root(p[x]);
}

void connect(ll x, ll y) {
    x = root(x);
    y = root(y);

    if(x == y) return;
    if(to[x].size() < to[y].size()) swap(x, y);

    p[y] = x;
    for(ll sz : to[y]) {
        to[x].insert(root(sz));
    }
    s[x] += s[y];
}

void dijkstra(ll start) {
    priority_queue<pair<ll, ll> > pq;
    pq.push({0, start});

    while (!pq.empty()) {
        auto curr = pq.top();
        pq.pop();

        ll w = -curr.first, x = curr.second;
        if(dist[x] != -1) continue;
        dist[x] = w;
        for(pair<ll, ll> sz : g[x]) {
            if(dist[sz.first] != -1) continue;
            pq.push({-w-sz.second, sz.first});
        }
    }
}

void dfs(int x) {
    volt[x] = 1;
    for(int sz : to[x]) {
        if(volt[root(sz)]) continue;
        dfs(root(sz));
    }
}

int main() {

    ios::sync_with_stdio(0);cin.tie(0); cout.tie(0);
    
    ll n,m,a,b;
    cin >> n >> m >> a >> b;
    g.resize(n+1);
    dist.resize(n+1, -1);

    for(ll i = 0; i < m; i++) {
        ll x,y,z;
        cin >> x >> y >> z;

        g[x].push_back({y, z});
        g[y].push_back({x, z});
    }

    dijkstra(b);
    for(int i = 1; i <= n; i++) p[i] = i;
    for(ll i = 1; i <= n; i++) {
        for(auto sz : g[i]) {
            if(dist[sz.first] == dist[i]) {
                connect(sz.first, i);
            }
        }
    }
    for(ll i = 1; i <= n; i++) {
        for(auto sz : g[i]) {
            if(dist[sz.first] < dist[i]) {
                to[root(i)].insert(root(sz.first));
            }
            else if(dist[sz.first] > dist[i]) {
                to[root(sz.first)].insert(root(i));
            }
        }
    }

    dfs(root(a));
    vector<ll> be(n+1, 0);
    for(ll i = 1; i <= n; i++) {

        if(!volt[i] || root(i) != i) continue;
        for(ll sz : to[i]) be[root(sz)]++;
    }

    vector<ll> tav(n+1, -1), parent(n+1, -1);
    
    tav[root(a)] = 0;
    queue<ll> sor;
    sor.push(root(a));
    while(!sor.empty()) {
        ll cur = sor.front();
        sor.pop();

        for(ll sz : to[cur]) {
            be[sz]--;
            if(tav[sz] < tav[cur] + s[cur]) {
                tav[sz] = tav[cur] + s[cur];
                parent[sz] = cur;
            }
            if(be[sz] == 0) {
                sor.push(sz);
            }
        }
    }

    set<ll> mo = {root(b)};
    ll cur = root(b);
    while(cur != root(a)) {
        cur = parent[cur];
        mo.insert(root(cur));
    }
    
    vector<ll> ki;
    for(ll i = 1; i <= n; i++) {
        if(mo.find(root(i)) != mo.end()) ki.push_back(i);
    }

    cout << ki.size() << "\n";
    for(ll sz : ki) cout << sz << ' ';
    
    return 0;
}
SubtaskSumTestVerdictTimeMemory
subtask10/0
1Accepted14ms26868 KiB
2Accepted201ms59660 KiB
subtask220/20
3Accepted87ms43160 KiB
4Accepted759ms119800 KiB
5Accepted551ms99492 KiB
6Accepted379ms82336 KiB
7Accepted556ms97800 KiB
8Accepted368ms80784 KiB
9Accepted1.049s145388 KiB
subtask315/15
10Accepted280ms77448 KiB
11Accepted29ms31676 KiB
12Accepted25ms31580 KiB
13Accepted64ms40452 KiB
14Accepted123ms55724 KiB
subtask420/20
15Accepted92ms49216 KiB
16Accepted250ms63476 KiB
17Accepted63ms45284 KiB
18Accepted97ms55640 KiB
19Accepted50ms40116 KiB
20Accepted94ms48172 KiB
subtask50/45
21Accepted12ms29152 KiB
22Time limit exceeded1.567s96692 KiB
23Accepted90ms43384 KiB
24Accepted397ms83504 KiB
25Accepted312ms73740 KiB
26Accepted587ms98408 KiB
27Accepted305ms75448 KiB
28Accepted1.2s156084 KiB
29Accepted197ms61144 KiB
30Accepted1.095s132140 KiB
31Accepted303ms65840 KiB
32Accepted85ms42608 KiB
33Accepted382ms74516 KiB
34Accepted321ms67788 KiB
35Accepted321ms67776 KiB
36Accepted301ms67768 KiB
37Accepted326ms67772 KiB