Запустимо пошук в глибину з довільної вершини дерева. Розглянемо ребро з вагою . Нехай кількість вершин у піддереві з коренем дорівнює . Тоді з одного боку ребра знаходиться вершин, а з іншого n - c вершин. Ребро входить в c \cdot (n - c) різних шляхів між парами вершин. Відповідно, внесок ребра в сумарну довжину всіх шляхів дорівнює c \cdot (n - c) \cdot w. Залишається перебрати всі ребра за допомогою пошуку в глибину та підсумувати їх внесок у загальну суму довжин шляхів.
Приклад
Графи, наведені у прикладах, мають вигляд:
Розглянемо внесок ребер у загальну суму довжин всіх шляхів у першому прикладі.
Ребро з вагою належить двом шляхам:
1 - 2;
2 - 1 - 3;
Його внесок у загальну суму дорівнює 1 \cdot 2 \cdot 1 = 2.
Ребро з вагою належить двом шляхам:
1 - 3;
2 - 1 - 3;
Його внесок у загальну суму дорівнює 2 \cdot 1 \cdot 3 = 6.
Сумарна довжина всіх шляхів дорівнює 2 + 6 = 8.
У другому прикладі розглянемо внесок ребра з вагою у загальну суму довжин усіх шляхів.
Кількість вершин у піддереві з коренем у вершині дорівнює c = 4 (включаючи саму вершину ). Тоді з іншого боку ребра знаходиться n - c = 6 - 4 = 2 вершини. Відповідно, кількість різних шляхів, що проходять через ребро , дорівнює c * (n - c) = 4 * 2 = 8. Дійсно, такими шляхами є всі пари вершин, в яких одна вершина лежить у піддереві вершини , а інша — поза ним:
Внесок ребра з вагою у суму довжин всіх шляхів становить
Реалізація алгоритму
Оголосимо значення модуля .
#define MOD 1000000000Вхідний зважений граф зберігаємо у списку суміжності .
vector<vector<pair<int,int> > > g;Функція dfs реалізує пошук в глибину з вершини і повертає кількість вершин у піддереві з коренем (включаючи саму вершину ). Підрахунок цих вершин ведеться у змінній . Ізначально вважаємо cnt = 1, оскільки саме піддерево вже містить вершину .
int dfs(int v, int p = -1)
{
int cnt = 1;
for (auto edge : g[v])
{
int to = edge.first;
int w = edge.second;Розглянемо ребро (v, to) з вагою . Визначаємо кількість вершин у піддереві з коренем у вершині . Тоді з одного боку ребра знаходяться вершин, а з іншого n - c вершин. Ребро (v, to) входить в c \cdot (n - c) різних шляхів між парами вершин.
Таким чином, внесок ребра (v, to) у сумарну довжину всіх шляхів дорівнює
if (to != p)
{
int c = dfs(to, v);
res = (res + 1LL * c * (n - c) * w) % MOD;
cnt += c;
}
}
return cnt;
}Основна частина програми. Читаємо вхідний зважений граф у список суміжності .
scanf("%d", &n);
g.resize(n + 1);
for (i = 1; i < n; i++)
{
scanf("%d %d %d", &u, &v, &d);
g[u].push_back({v, d});
g[v].push_back({u, d});
}Запускаємо пошук в глибину з вершини . Вершини графа нумеруються від до .
dfs(1);Виводимо відповідь.
printf("%lld\n",res);Java реалізація
import java.util.*;
import java.io.*;
class Edge
{
int to;
int weight;
public Edge(int to, int weight)
{
this.to = to;
this.weight = weight;
}
}
public class Main
{
static int MOD = 1000000000;
static ArrayList<Edge>[] g;
static int n, m;
static long res;
static int dfs(int v, int p)
{
int cnt = 1;
for(Edge edge : g[v])
{
int to = edge.to;
int w = edge.weight;
if (to != p)
{
int c = dfs(to, v);
res = (res + 1L * c * (n - c) * w) % MOD;
cnt += c;
}
}
return cnt;
}
public static void main(String[] args)
{
FastScanner con = new FastScanner(System.in);
n = con.nextInt();
g = new ArrayList[n+1]; // unchecked
for (int i = 0; i <= n; i++)
g[i] = new ArrayList<Edge>();
for (int i = 1; i < n; i++)
{
int u = con.nextInt();
int v = con.nextInt();
int d = con.nextInt();
g[u].add(new Edge(v,d));
g[v].add(new Edge(u,d));
}
dfs(1,-1);
System.out.println(res);
}
}
class FastScanner
{
BufferedReader br;
StringTokenizer st;
public FastScanner(InputStream inputStream)
{
br = new BufferedReader(new InputStreamReader(inputStream));
st = new StringTokenizer("");
}
public String next()
{
while (!st.hasMoreTokens())
{
try
{
st = new StringTokenizer(br.readLine());
} catch (Exception e)
{
return null;
}
}
return st.nextToken();
}
public int nextInt()
{
return Integer.parseInt(next());
}
}Python реалізація
Оголосимо значення модуля .
MOD = 1000000000Функція dfs реалізує пошук в глибину з вершини і повертає кількість вершин у піддереві з коренем (включаючи саму вершину ). Підрахунок цих вершин ведеться у змінній . Ізначально вважаємо cnt = 1, оскільки саме піддерево вже містить вершину .
def dfs(v, p = -1):
global res
cnt = 1
for to, w in g[v]:Розглянемо ребро (v, to) з вагою . Визначаємо кількість вершин у піддереві з коренем у вершині . Тоді з одного боку ребра знаходяться вершин, а з іншого n - c вершин. Ребро (v, to) входить в c \cdot (n - c) різних шляхів між парами вершин.
Таким чином, внесок ребра (v, to) у сумарну довжину всіх шляхів дорівнює
if to != p:
c = dfs(to, v)
res = (res + c * (n - c) * w) % MOD
cnt += c
return cntОсновна частина програми. Читаємо вхідні дані.
n = int(input())
g = [[] for _ in range(n + 1)]
res = 0
for _ in range(n - 1):
u, v, d = map(int, input().split())
g[u].append((v, d))
g[v].append((u, d))Запускаємо пошук в глибину з вершини . Вершини графа нумеруються від до .
dfs(1)Виводимо відповідь.
print(res)