Cod sursa(job #391510)

Utilizator alexandru92alexandru alexandru92 Data 5 februarie 2010 19:55:17
Problema Arbore partial de cost minim Scor 90
Compilator cpp Status done
Runda Arhiva educationala Marime 2.32 kb
/* 
 * File:   main.cpp
 * Author: virtualdemon
 *
 * Created on February 5, 2010, 7:18 PM
 */
#include <fstream>
#define Nmax 200000
#define inf 0x3f3f3f3f

/*
 *
 */
using namespace std;
struct list
{
    int vertex, cost;
    list *next;
} *L[Nmax];
int N;
int Heap[Nmax], Position[Nmax], cmin[Nmax], apm[Nmax];
inline void push( int x, int y, int c )
{
    list *q;
    q=new list;
    q->vertex=y; q->cost=c;
    q->next=L[x];
    L[x]=q;
}
void DownHeap( int k )
{
    int son;
    while( true )
    {
        son=2*k;
        if( son > N )
            return;
        if( son < N && cmin[Heap[son]] > cmin[Heap[son+1]])
            ++son;
        if( cmin[Heap[son]] >= cmin[Heap[k]] )
            return;
        swap( Heap[k], Heap[son] );
        swap( Position[Heap[k]], Position[Heap[son]] );
        k=son;
    }
}
void UpHeap( int k )
{
    int key=cmin[Heap[k]], f=k/2;
    while( k > 1 && key < cmin[Heap[f]] )
    {
        swap( Heap[k], Heap[f] );
        swap( Position[Heap[k]], Position[Heap[f]] );
        k=f;
        f/=2;
    }
}
void cut_root()
{
    swap( Heap[1], Heap[N] );
    swap( Position[Heap[1]], Position[Heap[N]] );
    Position[Heap[N]]=-1;
    --N;
    DownHeap(1);
}
void insert( int vertex, int value )
{
    Heap[++N]=vertex;
    cmin[vertex]=value;
    Position[vertex]=N;
    UpHeap(N);
}
int main( void )
{
    list *p;
    int n, m, x, y, c, v, s=0;
    ifstream in("apm.in");
    in>>n>>m;
    for( ; m; --m )
    {
        in>>x>>y>>c;
        x-=1; y-=1;
        push( x, y, c );
        push( y, x, c );
        if( !x )
            insert( y, c );
        else if( !y )
                insert( x, c );
    }
    Position[0]=-1;
    for( m=1; m < n; ++m )
    {
        v=Heap[1];
        s+=cmin[v];
        cut_root();
        for( p=L[v]; p; p=p->next )
        {
            if( -1 == Position[p->vertex] )
                continue;
            if( 0 == Position[p->vertex] )
              apm[p->vertex]=v,  insert( p->vertex, p->cost );
            else if( cmin[p->vertex] > p->cost )
                 {
                    apm[p->vertex]=v;
                    cmin[p->vertex]=p->cost;
                    UpHeap( Position[p->vertex] );
                }
        }
    }
    ofstream out("apm.out");
    out<<s<<' '<<(n-1)<<'\n';
    for( x=1; x < n; ++x )
        out<<(x+1)<<' '<<(apm[x]+1)<<'\n';
    return 0;
}