Cod sursă (job #178739)

#include <cstdio>
#include <vector>
#include <set>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <fstream>

#define NMax 8005
#define INF 0x3f3f3f3f
using namespace std;
ifstream f("zapada.in");
ofstream g("zapada.out");

set<pair<int,int> > H;
vector<pair<int,int> > vecini[NMax];
int n,m,k,x,y,c;

void prim(int x0){
    int dist[NMax],viz[NMax];
    memset(dist,INF,sizeof(dist));
    memset(viz,0,sizeof(viz));

    dist[x0] = 0;
    H.insert(make_pair(0,x0));

    while(!H.empty()){
        int nod = (*H.begin()).second;
        H.erase(*H.begin());
        viz[nod] = 1;
        for(int i = 0; i < vecini[nod].size(); i++){
            if(vecini[nod][i].second < dist[vecini[nod][i].first] && !viz[vecini[nod][i].first]){
                dist[vecini[nod][i].first] = vecini[nod][i].second;
                H.insert(make_pair(dist[vecini[nod][i].first],vecini[nod][i].first));
            }
        }
    }
    int s = 0;
    for(int i = 1; i <= n; i++){
        s += dist[i];
    }
    g << s <<"\n";
}
int main()
{
    f >> n >> m >> k;
    for(int i = 1; i <= m; i++){
        f >> x >> y >> c;
        vecini[x].push_back(make_pair(y,c));
        vecini[y].push_back(make_pair(x,c));
    }
    for(int i = 0; i <= k; i++){
        prim(1);
        f >> x >> y >> c;
        vecini[x].push_back(make_pair(y,c));
        vecini[y].push_back(make_pair(x,c));
    }
    return 0;
}