#include <iostream>
#include <bits/stdc++.h>
#define inf 100000000
using namespace std;
class graph
{
// bitul 1 orientat sau neorientat, bitul 2 cost muchie
enum input {orientat = 1, neorientat = 2, cu_cost = 4, fillLat = 8};
struct muchie
{
int from, to, cost;
muchie(int _from, int _to, int _cost = 1):
from(_from), to(_to), cost(_cost)
{
}
};
int n,m;
int* extraArgs = nullptr;
vector<int>* la = nullptr;
vector<int>* lat = nullptr;
vector<pair<int,int>>* lac = nullptr;
vector<int> topo;
ifstream f;
ofstream g;
public:
void clearData();
void readData(string nume_fisier, input _tip_graf, int _extra_args_count = 0);
//////////////////////////////////////////////////////////////////////////////////
void bfs();
//////////////////////////////////////////////////////////////////////////////////
void dfs();
void dfsTopo(int from, int dist[]);
//////////////////////////////////////////////////////////////////////////////////
void ctc();
void dfsCtc(int nod, vector<vector<int>>& cc, int data[]);
//////////////////////////////////////////////////////////////////////////////////
void sortaret();
//////////////////////////////////////////////////////////////////////////////////
void critice();
void dfsCritice(int from, int parent, vector<pair<int,int>>& mcrit, int lowLink[], int ids[], int id);
//////////////////////////////////////////////////////////////////////////////////
void biconex();
void dfsBiconexe(int from, stack<pair<int,int>>& st, vector<vector<int>>& cb, int lowlink[], int ids[], int id);
//////////////////////////////////////////////////////////////////////////////////
void apm();
pair<int,int> getParentAndHeight(int x, int parents[]);
//////////////////////////////////////////////////////////////////////////////////
void disjoint();
//////////////////////////////////////////////////////////////////////////////////
void dijkstra();
//////////////////////////////////////////////////////////////////////////////////
void bellmanford();
//////////////////////////////////////////////////////////////////////////////////
void royfloyd();
};
void graph::clearData()
{
if (la != nullptr)
{
delete[] la;
la = nullptr;
}
if (lat != nullptr)
{
delete[] lat;
lat = nullptr;
}
if (lac != nullptr)
{
delete[] lac;
lac = nullptr;
}
if (extraArgs != nullptr)
{
delete[] extraArgs;
extraArgs = nullptr;
}
}
void graph::readData (string nume_fisier, input _tip_graf, int _extra_args_count)
{
f.open(nume_fisier);
f >> n >> m;
if (_extra_args_count)
{
extraArgs = new int[_extra_args_count];
for (int i = 0; i < _extra_args_count; i++)
f >> extraArgs[i];
}
int x,y,c;
if (_tip_graf & input::orientat)
{
if (_tip_graf & input::cu_cost)
{
lac = new vector<pair<int,int>>[n+1];
for (int i = 1; i <= m; i++)
{
f >> x >> y >> c;
lac[x].push_back({y,c});
}
}
else
{
if (_tip_graf & input::fillLat)
{
la = new vector<int>[n+1];
lat = new vector<int>[n+1];
for (int i = 1; i <= m; i++)
{
f >> x >> y;
// cout << x << ' ' << y << '\n';
la[x].push_back(y);
lat[y].push_back(x);
}
}
else
{
la = new vector<int>[n+1];
for (int i = 1; i <= m; i++)
{
f >> x >> y;
la[x].push_back(y);
}
}
}
}
else if (_tip_graf & input::neorientat)
{
if (_tip_graf & input::cu_cost)
{
lac = new vector<pair<int,int>>[n+1];
for (int i = 1; i <= m; i++)
{
f >> x >> y >> c;
lac[x].push_back({y,c});
lac[y].push_back({x,c});
}
}
else
{
la = new vector<int>[n+1];
for (int i = 1; i <= m; i++)
{
f >> x >> y;
la[x].push_back(y);
la[y].push_back(x);
}
}
}
f.close();
}
void graph::bfs()
{
ofstream g("bfs.out");
readData("bfs.in", input::orientat, 1);
int start = extraArgs[0];
int* dist = new int[n+1];
queue<int> q;
q.push(start);
dist[start] = 0;
while (q.size())
{
for (auto to: la[q.front()])
{
if (dist[to] == -1)
{
q.push(to);
dist[to] = dist[q.front()] + 1;
}
}
q.pop();
}
for (int i = 1; i <= n; i++)
g << dist[i] << ' ';
delete[] dist;
clearData();
}
void graph::dfsTopo(int from, int dist[])
{
dist[from] = 1;
for (auto to: la[from])
{
if (!dist[to])
dfsTopo(to, dist);
}
topo.push_back(from);
}
void graph::dfs()
{
ofstream g("dfs.out");
readData("dfs.in", input::neorientat);
int* dist = new int[n+1];
int cc = 0;
for (int i = 1; i <= n; i++)
{
if (!dist[i])
{
dfsTopo(i, dist);
cc++;
}
}
g << cc;
delete[] dist;
clearData();
}
void graph::dfsCtc(int nod, vector<vector<int>>& cc, int dist[])
{
dist[nod] = 1;
cc[cc.size() - 1].push_back(nod);
for (auto to: lat[nod])
{
if (!dist[to])
dfsCtc(to, cc, dist);
}
}
void graph::ctc()
{
ofstream g("ctc.out");
readData("ctc.in", (input)(input::orientat | input::fillLat));
int* dist = new int[n+1];
for (int i = 1; i <= n; i++)
{
if (!dist[i])
dfsTopo(i, dist);
}
reverse(topo.begin(), topo.end());
memset(dist, 0, sizeof(int) * (n+1));
vector<vector<int>> cc;
for (auto nod: topo)
{
if (!dist[nod])
{
cc.push_back(vector<int>());
dfsCtc(nod, cc, dist);
}
}
g << cc.size() << '\n';
for (auto& row: cc)
{
for (auto el: row)
g << el << ' ';
g << '\n';
}
delete[] dist;
clearData();
}
void graph::sortaret()
{
ofstream g("sortaret.out");
readData("sortaret.in", (input)(input::orientat));
int* dist = new int[n+1];
for (int i = 1; i <= n; i++)
{
if (!dist[i])
dfsTopo(i, dist);
}
reverse(topo.begin(), topo.end());
for (auto el: topo)
g << el << ' ';
delete[] dist;
clearData();
}
void graph::dfsCritice(int from, int parent, vector<pair<int,int>>& mcrit, int lowLink[], int ids[], int id)
{
lowLink[from] = id;
ids[from] = id;
id ++;
for (auto to: la[from])
{
if (ids[to] == -1)
{
dfsCritice(to, from, mcrit, lowLink, ids, id);
lowLink[from] = min(lowLink[from], lowLink[to]);
}
else if (to != parent)
lowLink[from] = min(lowLink[from], ids[to]);
}
if (lowLink[from] == ids[from] && parent != 0)
mcrit.push_back({from, parent});
}
void graph::critice()
{
ofstream g("critice.out");
readData("critice.in", (input)(input::neorientat));
int* ids = new int[n+1];
for (int i = 1; i <= n; i++)
ids[i] = -1;
int* lowLink = new int[n+1];
for (int i = 1; i <= n; i++)
lowLink[i] = 0;
vector<pair<int, int>> mcrit;
for (int i = 1; i <= n; i++)
{
if (ids[i] == -1)
{
dfsCritice(i, 0, mcrit, lowLink, ids, 0);
}
}
for (auto& per: mcrit)
g << per.first << ' ' << per.second << '\n';
delete[] ids;
delete[] lowLink;
}
void graph::dfsBiconexe(int from, stack<pair<int,int>>& st, vector<vector<int>>& cb, int lowLink[], int ids[], int id)
{
lowLink[from] = ids[from] = id++;
for (auto to: la[from])
{
if (ids[to] != -1)
lowLink[from] = min(lowLink[from], ids[to]);
else
{
st.push({from, to});
dfsBiconexe(to, st, cb, lowLink, ids, id);
lowLink[from] = min(lowLink[from], lowLink[to]);
if (lowLink[to] >= ids[from])
{
vector<int> comp;
while (st.top().first != from)
{
comp.push_back(st.top().second);
st.pop();
}
comp.push_back(st.top().first);
comp.push_back(st.top().second);
st.pop();
cb.push_back(comp);
}
}
}
}
void graph::biconex()
{
ofstream g("biconex.out");
readData("biconex.in", (input)(input::neorientat));
int* ids = new int[n+1];
for (int i = 1; i <= n; i++)
ids[i] = -1;
int* lowLink = new int[n+1];
for (int i = 1; i <= n; i++)
lowLink[i] = 0;
vector<vector<int>> cb;
stack<pair<int,int>> st;
for (int i = 1; i <= n; i++)
{
if (ids[i] == -1)
dfsBiconexe(i, st, cb, lowLink, ids, 0);
}
g << cb.size() << '\n';
for (auto& row: cb)
{
for (auto el: row)
g << el << ' ';
g << '\n';
}
delete[] ids;
delete[] lowLink;
}
pair<int,int> graph::getParentAndHeight(int x, int parents[])
{
vector<int> met;
int h = 0;
while (parents[x] != 0)
{
h++;
met.push_back(x);
x = parents[x];
}
for (auto el: met)
parents[el] = x;
return {x, h};
}
void graph::apm()
{
ofstream g("apm.out");
readData("apm.in", (input)(input::orientat | input::cu_cost));
int* parents = new int[n+1];
vector<muchie> muchii;
for (int i = 1; i <= n; i++)
{
for (auto per: lac[i])
muchii.push_back({i, per.first, per.second});
}
sort(muchii.begin(), muchii.end(), [](muchie& m1, muchie& m2){ return (m1.cost < m2.cost);});
vector<muchie> sol;
int s = 0;
for (auto& m: muchii)
{
auto [p1, h1] = getParentAndHeight(m.from, parents);
auto [p2, h2] = getParentAndHeight(m.to, parents);
if (p1 == p2)
continue;
if (h1 > h2)
{
parents[p1] = p2;
}
else
{
parents[p2] = p1;
}
sol.push_back(m);
s += m.cost;
}
g << s << '\n';
g << sol.size() << '\n';
for (auto& m: sol)
g << m.from << ' ' << m.to << '\n';
delete[] parents;
}
void graph::disjoint()
{
ifstream f("disjoint.in");
ofstream g("disjoint.out");
f >> n >> m;
int* parents = new int[n+1];
for (int i = 1; i <= m; i++)
{
int op, from, to;
f >> op >> from >> to;
auto [p1, h1] = getParentAndHeight(from, parents);
auto [p2, h2] = getParentAndHeight(to, parents);
if (op == 1)
{
if (h1 > h2)
parents[p2] = p1;
else
parents[p1] = p2;
}
else
{
if (p1 == p2)
g << "DA\n";
else
g << "NU\n";
}
}
delete[] parents;
}
void graph::dijkstra()
{
const int max_dist = 1111111111;
ofstream g("dijkstra.out");
readData("dijkstra.in", (input)(input::orientat | input::cu_cost));
int* dist = new int[n+1];
for (int i = 1; i <= n; i++)
dist[i] = max_dist;
int* viz = new int[n+1];
priority_queue<pair<int, int>> q; // cost negativ si nod
q.push({0, 1});
dist[1] = 0;
while (q.size())
{
int top = q.top().second;
if (viz[top])
continue;
viz[top] = 1;
for (auto& per: lac[top])
{
if (dist[per.first] > dist[top] + per.second)
{
dist[per.first] = dist[top] + per.second;
q.push({-per.second, per.first});
}
}
q.pop();
}
for (int i = 2; i <= n; i++)
{
if (dist[i] == max_dist)
{
dist[i] = 0;
}
g << dist[i] << ' ';
}
delete[] dist;
}
void graph::bellmanford()
{
const int max_dist = 1111111111;
ofstream g("bellmanford.out");
readData("bellmanford.in", (input)(input::orientat | input::cu_cost));
int* dist = new int[n+1];
for (int i = 1; i <= n; i++)
dist[i] = max_dist;
int* used = new int[n+1];
for (int i = 1; i <= n; i++)
used[i] = 0;
vector<int> to_check_now;
to_check_now.push_back(1);
dist[1] = 0;
for (int k = 1; k <= n && to_check_now.size(); k++)
{
vector<int> to_check_next;
for (int i = 1; i <= n; i++)
used[i] = 0;
for (auto el: to_check_now)
{
for (auto& m: lac[el])
{
if (dist[m.first] > dist[el] + m.second)
{
dist[m.first] = dist[el] + m.second;
if (!used[m.first])
{
used[m.first] = 1;
to_check_next.push_back(m.first);
}
}
}
}
to_check_next.swap(to_check_now);
}
for (auto el: to_check_now)
for (auto& m: lac[el])
if (dist[m.first] > dist[el] + m.second)
{
g << "Ciclu negativ!";
return;
}
for (int i = 2; i <= n; i++)
{
if (max_dist == dist[i])
dist[i] = 0;
g << dist[i] << ' ';
}
delete[] dist;
delete[] used;
}
void graph::royfloyd()
{
ifstream in("royfloyd.in");
ofstream out("royfloyd.out");
in >> n;
int minDist[100][100];
for (int i = 0; i < n; i++)
for (int j = 0; j < n; j++)
{
in >> minDist[i][j];
if (minDist[i][j] == 0 && i != j)
minDist[i][j] = inf;
}
for (int k = 0; k < n; k++)
for (int i = 0; i < n; i++)
for (int j = 0; j < n; j++)
//minDist[i][j] = min(min(minDist[i][j], minDist[i][k] + dist[k][j]), dist[i][k] + minDist[k][j]);
if (minDist[i][k] + minDist[k][j] < minDist[i][j])
minDist[i][j] = minDist[i][k] + minDist[k][j];
for (int i = 0; i < n; i++)
{
for (int j = 0; j < n; j++)
if (minDist[i][j] == inf)
{
out << 0 << ' ';
continue;
}
else
out << minDist[i][j] << ' ';
out << '\n';
}
}
int main()
{
graph g;
g.royfloyd();
return 0;
}