#define ll long long
ll prims(int n)
{
ll ans;
vector<bool> used (n);
#define INF 1000000000000LLvector<ll> min_e (n, INF), sel_e (n, -1);
min_e[0]=-1*INF;
ll dis=1;
for(int i=0;i<n;i++)
{
int v=-1;
for(int j=0;j<n;j++)
{
if (!used[j] && (v == -1 || min_e[j] < min_e[v]))
v = j;
}
used[v] = true;
if(sel_e[v]!=-1)
cout << v << " " << sel_e[v] << endl;
for (int to=0; to<n; ++to)
if (g[v][to] < min_e[to]) {
min_e[to] = g[v][to];
sel_e[to] = v;
}
}
for(int i=0;i<n;i++) cout<<i<<" "<<sel_e[i]<<" "<<g[i][sel_e[i]]<<endl;return dis;
}
Я пытаюсь применить алгоритм Прима для плотного неориентированного графа для отрицательных весов ребер, но я не могу понять, почему он дает неправильные результаты почти во всех случаях.
Я использую матрицу смежности g [N] [N] для хранения ребер.
На самом деле вывод для моего текущего кода — это минимальное связующее дерево с циклами. Почему не работает механизм проверки цикла?
На самом деле, проблема здесь:
for (int to=0; to<n; ++to)
if (g[v][to] < min_e[to]) {
min_e[to] = g[v][to];
sel_e[to] = v;
}
}
Вы должны только обновить sel_e
а также min_e
если to
еще не был
В противном случае рассмотрим этот случай:
0 -- 1 -- 2
где w({0, 1}) = 10
, а также w({1, 2} = 1)
, Вы бы установили sel_e[1] = 2
даже если вам нужно sel_e[1] = 0
,
Других решений пока нет …