Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- SeedRandom[1]
- arrayR = Append[Thread[RandomSample[Tuples[Range[1000], {2}], 5000] ->
- RandomReal[1, 5000]], {_, _} -> 0];
- weightedGraph=WeightedAdjacencyGraph[SparseArray[Most@arrayR, {1000, 1000}, ∞]];
- weG = weightedGraph;
- resMatrix = AdjacencyMatrix[weG];
- Table[minST = FindSpanningTree[weG];
- adjM1 = AdjacencyMatrix[minST];
- adjWG1 = AdjacencyMatrix[weG];
- weG = WeightedAdjacencyGraph[
- SparseArray[
- Most@ArrayRules[
- Chop[(adjWG1 - adjM1)]*WeightedAdjacencyMatrix[weG]], {1000,
- 1000}, [Infinity]]];, {3}];
- Chop[resMatrix - AdjacencyMatrix[weG]]
Add Comment
Please, Sign In to add comment