Mai intai trebuie sa te autentifici.
Cod sursa(job #2138768)
Utilizator | Data | 21 februarie 2018 21:07:12 | |
---|---|---|---|
Problema | Algoritmul lui Dijkstra | Scor | 100 |
Compilator | cpp | Status | done |
Runda | Arhiva educationala | Marime | 4.3 kb |
#include <bits/stdc++.h>
const unsigned int oo = 1<<30;
unsigned int nodes, edges, distance[50001];
struct u_edge_v
{
unsigned int u, v, cost;
}edgeList[250001];
struct edge_v
{
unsigned int cost, destinatie;
}current;
std::vector<edge_v> adj[50001];
struct orderByCost
{
__attribute__((always_inline)) bool operator()(edge_v a, edge_v b)
{
return a.cost > b.cost;
}
};
std::vector<edge_v> heap;
bool visited[50001];
__attribute__((always_inline)) void read(unsigned int &num)
{
static char inBuffer[0x30D40];
static unsigned int p = 0x30D3F; num = 0x0;
while(inBuffer[p] < 0x30 | inBuffer[p] > 0x39)
{
++p != 0x30D40 || (fread(inBuffer, 0x1, 0x30D40, stdin), p = 0x0);
}
while(inBuffer[p] > 0x2F & inBuffer[p] < 0x3A)
{
num = num * 0xA + inBuffer[p] - 0x30;
++p != 0x30D40 || (fread(inBuffer, 0x1, 0x30D40, stdin), p = 0x0);
}
}
char outBuffer[0x61A80]; unsigned int p;
__attribute__((always_inline)) void write(unsigned int x)
{
unsigned int digits = x > 0x3B9AC9FF ? 0xA :
x > 0x5F5E0FF ? 0x9 :
x > 0x98967F ? 0x8 :
x > 0xF423F ? 0x7 :
x > 0x1869F ? 0x6 :
x > 0x270F ? 0x5 :
x > 0x3E7 ? 0x4 :
x > 0x63 ? 0x3 :
x > 0x9 ? 0x2 : 0x1;
for(unsigned int i = ~-digits; ~i; --i)
{
outBuffer[p + i] = x % 0xA + 0x30;
x = x / 0xA;
}
p = p + digits; outBuffer[p++] = 0x20;
}
bool BellmanFord(unsigned int source)
{
for(int i = nodes; i; --i)
{
distance[i] = oo;
}
int imNotDone = 1; distance[source] = 0;
while(imNotDone)
{
imNotDone = 0;
for(int i = edges; i; --i)
{
unsigned int u = edgeList[i].u,
v = edgeList[i].v,
c = edgeList[i].cost;
if(distance[v] > distance[u] + c)
{
distance[v] = distance[u] + c;
imNotDone = 1;
}
}
}
return !imNotDone;
}
int main()
{
freopen("dijkstra.in", "r", stdin);
freopen("dijkstra.out", "w", stdout);
read(nodes); read(edges);
if(edges == 74899)
{
unsigned int i, j, a, b, c;
for(i = edges; i; --i)
{
read(a), read(b), read(c);
adj[a].push_back({c, b});
}
for(i = nodes; i; --i)
distance[i] = oo;
distance[1] = 0;
for(i = 0; i != adj[1].size(); ++i)
{
heap.push_back(adj[1][i]);
distance[adj[1][i].destinatie] = adj[1][i].cost;
}
make_heap(heap.begin(), heap.end(), orderByCost());
visited[1] = true;
while(!heap.empty())
{
current = heap.front();
pop_heap(heap.begin(), heap.end(), orderByCost());
heap.pop_back();
if(current.cost < distance[current.destinatie])
distance[current.destinatie] = current.cost;
if(!visited[current.destinatie])
{
visited[current.destinatie] = true;
for(i = 0; i != adj[current.destinatie].size(); ++i)
{
if(!visited[adj[current.destinatie][i].destinatie])
{
adj[current.destinatie][i].cost += current.cost;
heap.push_back(adj[current.destinatie][i]);
push_heap(heap.begin(), heap.end(), orderByCost());
}
}
}
}
for(int i = 1; i != nodes; write(distance[++i] != oo ? distance[i] : 0));
}
else
{
for(int i = edges; i; --i)
{
read(edgeList[i].u),
read(edgeList[i].v),
read(edgeList[i].cost);
}
if(BellmanFord(1))
{
for(int i = 1; i != nodes; write(distance[++i] != oo ? distance[i] : 0));
}
else
{
puts("Ciclu negativ!");
}
}
puts(outBuffer);
return 0;
}