Cod sursa(job #2194253)

Utilizator UnseenMarksmanDavid Catalin UnseenMarksman Data 12 aprilie 2018 17:54:57
Problema Arbore partial de cost minim Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 2.27 kb
#include <fstream>
#include <vector>
using namespace std;

ifstream fin("apm.in");
ofstream fout("apm.out");

const int N = 2e5+2;
const int inf = 1<<30;

int n, m, d[N], t[N], H[2*N+2], L, pos[N], result;
vector<pair<int,int>>G[N],R;

void introducere_in_apm(int x)
{
    for(unsigned int i=0; i<G[x].size(); i++)
    {
        int nod=G[x][i].first;
        int cost=G[x][i].second;
        d[nod]=min(d[nod],cost);
        if(d[nod]==cost)
        {
            t[nod]=x;
        }
    }
}

void upheap(int x)
{
    int poz=x;
    while(poz>1 && d[H[poz]]<d[H[poz/2]])
    {
        swap(H[poz],H[poz/2]);
        swap(pos[H[poz]],pos[H[poz/2]]);
        poz=poz/2;
    }
}

void downheap(int x)
{
    int poz=x;
    while((poz*2<=L && d[H[poz]]>d[H[poz*2]]) || (poz*2+1<=L && d[H[poz]]>d[H[poz*2+1]]))
    {
        if(d[H[poz*2]]<d[H[poz*2+1]] || poz*2+1>L)
        {
            swap(H[poz],H[poz*2]);
            swap(pos[H[poz]],pos[H[poz*2]]);
            poz=poz*2;
        }
        else
        {
            swap(H[poz],H[poz*2+1]);
            swap(pos[H[poz]],pos[H[poz*2+1]]);
            poz=poz*2+1;
        }
    }
}

void introducere_in_heap(int x)
{
    H[++L]=x;
    pos[x]=L;
    upheap(L);
}

int scoate_radacina_heap()
{
    int x=H[1];
    swap(H[1],H[L]);
    swap(pos[H[1]],pos[H[L]]);
    L--;
    pos[x]=0;
    downheap(1);
    return x;
}

int main()
{
    fin>>n>>m;
    for(int i=1; i<=m; i++)
    {
        int x,y,c;
        fin>>x>>y>>c;
        G[x].push_back(make_pair(y,c));
        G[y].push_back(make_pair(x,c));
    }
    for(int i=2; i<=n; i++) d[i]=inf;

    introducere_in_apm(1);

    for(int i=2; i<=n; i++)
        introducere_in_heap(i);

    for(int i=1; i<n; i++)
    {
        int x=scoate_radacina_heap();
        introducere_in_apm(x);
        result=result+d[x];
        R.push_back(make_pair(t[x],x));
        for(unsigned int j=0; j<G[x].size(); j++)
        {
            int nod=G[x][j].first;
            if(pos[nod])
            {
                upheap(pos[nod]);
            }
        }
    }
    fout<<result<<'\n'<<n-1<<'\n';
    for(unsigned int i=0; i<n-1; i++)
    {
        fout<<R[i].first<<' '<<R[i].second<<'\n';
    }
    return 0;
}