Cod sursa(job #2789855)

Utilizator gruhtenZinnenberg Gruhten gruhten Data 28 octombrie 2021 01:33:18
Problema Algoritmul lui Dijkstra Scor 90
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 2.16 kb
#include<fstream>
#include<vector>
#include<queue>
#define oo 200000000
#define pb push_back
#define mp make_pair
#define dim 50004
using namespace std;
ifstream f("dijkstra.in");
ofstream g("dijkstra.out");
int d[dim],n,m,b[dim],H[dim],N=0;
bool w[dim];
int x,fiu,tata;

void urcam(int vecinul)
{
    tata=fiu/2;
    while(x<H[tata])
    {
        H[fiu]=H[tata];
        b[fiu]=b[tata];
        int v=b[tata];

        fiu=tata;
        tata=tata/2;
    }
    H[fiu]=x;
    b[fiu]=vecinul;
}

vector < pair < int, int > > a[dim];

inline void cit()
{
    f>>n>>m;
    int i,y,x,c;
    for(i = 1; i<= m; i ++)
    {
        f>>x>>y>>c;
        a[x].pb(mp(y,c));
    }
}

inline void dij(int i)
{
    int poz;

    for(int j = 1; j <= n; j ++)
        d[j] = oo;

    d[i] = 0;
    b[1]=i;

    w[i] = true;
    H[0]=-1;

    H[1]=i;
    N=1;

    while(N>0)
    {
        poz=b[1];//poz=nodul din varf
        x=H[N];
        int v=b[N];
        b[1]=v;
        --N; //stergem nodul din varful heap-ului
        w[poz] = false;

        tata=1;
        fiu=2;
        while(fiu<=N)
        {
            if(fiu<N)
                if(H[fiu]>H[fiu+1])
                    fiu=fiu+1;

            if(x>H[fiu])
            {
                H[tata]=H[fiu];
                b[tata]=b[fiu];
                tata=fiu;
                fiu*=2;
            }
            else
                fiu=N+1;
        }
        H[tata]=x;
        b[tata]=v;

        for(unsigned int j = 0; j < a[poz].size(); ++ j)
        {
            int vecin = a[poz][j].first;
            int cost = a[poz][j].second;
            if(d[vecin] > d[poz] + cost)
            {
                x=d[vecin]=d[poz]+cost;

                if(w[vecin] == false)
                {
                    ++N;
                    b[N]=vecin;
                    H[N]=x;
                    fiu=N;
                    urcam(vecin);
                }
            }
        }
    }
}

int main()
{
    int i;
    cit();
    dij(1);
    for(i = 2; i <= n; i ++)
       if(d[i] == oo)
        g<<"0 ";
       else
        g<<d[i]<<" ";
}