TravellingSalesAlg.java 14 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452
  1. package pathgame.algorithm;
  2. import pathgame.tilemap.TileMap;
  3. import java.util.ArrayList;
  4. import java.util.Collections;
  5. public class TravellingSalesAlg {
  6. public static int calcSalesPathLen(TileMap map) {
  7. DijkstraMagic dijkstra = new DijkstraMagic(map);
  8. ArrayList<ArrayList<SaleRoute>> salesPitch = dijkstra.getSalesPitch();
  9. //make minimum spanning tree
  10. ArrayList<TreeEdge> MSTree = makeMSTree(salesPitch);
  11. //find pairs-shortest path for verteces with odd degree and add those edges to MSTree
  12. ArrayList<TreeEdge> oddDegEdges = makeOddDegEdges(MSTree, salesPitch);
  13. MSTree.addAll(oddDegEdges);
  14. //finde euler tour
  15. ArrayList<Integer> eulerTour = getEulerTour(MSTree);
  16. //cut short
  17. cutShort(eulerTour);
  18. //TODO: write to logging file
  19. //calculate the total weight of the tour using the edge table (salesPitch)
  20. int tourWeight = calcTourWeight(eulerTour, salesPitch);
  21. System.out.println("min cost: " + tourWeight);
  22. return tourWeight;
  23. //brute force
  24. //return bruteForce(salesPitch);
  25. }
  26. private static int bruteForce(ArrayList<ArrayList<SaleRoute>> salesPitch) {
  27. Permutation permutation = new Permutation(salesPitch.size(), 1);
  28. OddDegreeList nodeList = new OddDegreeList();
  29. for(int i = 0; i < salesPitch.size()+1; i++) {
  30. nodeList.add(i);
  31. }
  32. while(!permutation.isOverFlowing()) {
  33. int newSum = bruteSumPermut(permutation, salesPitch, nodeList);
  34. if(newSum < permutation.getMinCost()) {
  35. permutation.setMinCost(newSum);
  36. }
  37. /*permutation.printPermut();
  38. System.out.println(newSum);
  39. System.out.println();*/
  40. permutation.increaseLowest();
  41. }
  42. permutation.makePermutMinimal();
  43. permutation.printPermut();
  44. System.out.println("min cost: " + permutation.getMinCost());
  45. return 0;
  46. }
  47. private static int bruteSumPermut(Permutation permutation, ArrayList<ArrayList<SaleRoute>> salesPitch, OddDegreeList nodeList) {
  48. int sum = 0;
  49. nodeList.resetUsed();
  50. int first = -1;
  51. int last = -1;
  52. for(int i = 0; i < permutation.size(); i++) {
  53. int pos = nodeList.getUnused(permutation.getValAtPos(i));
  54. //System.out.println("pos: " + pos);
  55. if(first == -1) {
  56. first = pos;
  57. }
  58. else {
  59. //System.out.println(last + " - " + pos);
  60. sum += bruteWeight(last, pos, salesPitch);
  61. //System.out.println(sum);
  62. }
  63. last = pos;
  64. nodeList.makePosUsed(pos);
  65. }
  66. //System.out.println(last + " - " + nodeList.getUnused(0));
  67. int pos = nodeList.getUnused(0);
  68. sum += bruteWeight(last, pos, salesPitch);
  69. //System.out.println(pos + " - " + first);
  70. sum += bruteWeight(pos, first, salesPitch);
  71. return sum;
  72. }
  73. private static int bruteWeight(int first, int sec, ArrayList<ArrayList<SaleRoute>> salesPitch) {
  74. if(first < sec) {
  75. return salesPitch.get(first).get(sec-first-1).getTotalCost();
  76. }
  77. else {
  78. return salesPitch.get(sec).get(first-sec-1).getTotalCost();
  79. }
  80. }
  81. private static int calcTourWeight(ArrayList<Integer> tour, ArrayList<ArrayList<SaleRoute>> salesPitch) {
  82. int totalWeight = 0;
  83. for(int i = 0; i < tour.size()-1; i++) {
  84. int startNode, endNode;
  85. if(tour.get(i) < tour.get(i+1)) {
  86. startNode = tour.get(i);
  87. endNode = tour.get(i+1) - startNode - 1;
  88. }
  89. else {
  90. startNode = tour.get(i+1);
  91. endNode = tour.get(i) - startNode - 1;
  92. }
  93. //System.out.println(startNode + "/" + endNode + ": " + salesPitch.get(startNode).get(endNode).getTotalCost());
  94. totalWeight += salesPitch.get(startNode).get(endNode).getTotalCost();
  95. }
  96. //System.out.println(totalWeight);
  97. return totalWeight;
  98. }
  99. private static void cutShort(ArrayList<Integer> eulerTour) {
  100. int counter = 2;
  101. while(counter < eulerTour.size()-1) {
  102. int current = eulerTour.get(counter);
  103. boolean found = false;
  104. for(int i = 0; i < counter; i++) {
  105. if(eulerTour.get(i) == current) {
  106. found = true;
  107. break;
  108. }
  109. }
  110. if(found) {
  111. eulerTour.remove(counter);
  112. }
  113. else {
  114. counter++;
  115. }
  116. }
  117. /*System.out.println("cut tour");
  118. for(int i = 0; i < eulerTour.size(); i++) {
  119. if(i != 0) {
  120. System.out.print(" - ");
  121. }
  122. System.out.print(eulerTour.get(i));
  123. }
  124. System.out.println();*/
  125. }
  126. private static ArrayList<Integer> getEulerTour(ArrayList<TreeEdge> graph) {
  127. ArrayList<Integer> tour = new ArrayList<>();
  128. while (graph.size() > 0) {
  129. /*System.out.println();
  130. for(int i = 0; i < graph.size(); i++) {
  131. System.out.println("'" + graph.get(i).getSrc() + " : " + graph.get(i).getDest());
  132. }
  133. System.out.println();*/
  134. if(tour.size() == 0) {
  135. tour = getSubtour(graph, graph.get(0).getSrc());
  136. }
  137. else {
  138. int start = -1;
  139. for(int e = 0; e < graph.size(); e++) {
  140. TreeEdge edge = graph.get(e);
  141. for(int tp = 0; tp < tour.size(); tp++) {
  142. if(edge.getSrc() == tour.get(tp)) {
  143. start = edge.getSrc();
  144. break;
  145. }
  146. else if(edge.getDest() == tour.get(tp)) {
  147. start = edge.getDest();
  148. break;
  149. }
  150. }
  151. if(start!= -1) {
  152. break;
  153. }
  154. }
  155. ArrayList<Integer> subTour = getSubtour(graph, start);
  156. mergeTours(tour, subTour);
  157. /*System.out.println("merged tour");
  158. for(int i = 0; i < tour.size(); i++) {
  159. if(i != 0) {
  160. System.out.print(" - ");
  161. }
  162. System.out.print(tour.get(i));
  163. }
  164. System.out.println();*/
  165. }
  166. }
  167. return tour;
  168. }
  169. private static ArrayList<Integer> getSubtour(ArrayList<TreeEdge> graph, int start) {
  170. ArrayList<Integer> tour = new ArrayList<>();
  171. tour.add(start);
  172. int pos = nextTourEdgePos(graph, start);
  173. int next = graph.get(pos).getOtherVertex(start);
  174. graph.remove(pos);
  175. tour.add(next);
  176. while (next != start) {
  177. pos = nextTourEdgePos(graph, next);
  178. next = graph.get(pos).getOtherVertex(next);
  179. graph.remove(pos);
  180. tour.add(next);
  181. }
  182. /*for(int i = 0; i < tour.size(); i++) {
  183. if(i != 0) {
  184. System.out.print(" - ");
  185. }
  186. System.out.print(tour.get(i));
  187. }
  188. System.out.println();*/
  189. return tour;
  190. }
  191. private static int nextTourEdgePos(ArrayList<TreeEdge> graph, int vertex) {
  192. for(int i = 0; i < graph.size(); i++) {
  193. if(graph.get(i).getSrc() == vertex || graph.get(i).getDest() == vertex) {
  194. return i;
  195. }
  196. }
  197. return -1;
  198. }
  199. private static void mergeTours(ArrayList<Integer> tour, ArrayList<Integer> subTour) {
  200. int mergeTo = subTour.get(0);
  201. int mergePos = -1;
  202. for(int i = 0; i < tour.size(); i++) {
  203. if (tour.get(i) == mergeTo) {
  204. mergePos = i;
  205. }
  206. }
  207. for(int i = subTour.size()-1; i > 0; i--) {
  208. tour.add(mergePos+1, subTour.get(i));
  209. }
  210. }
  211. private static ArrayList<TreeEdge> makeOddDegEdges(ArrayList<TreeEdge> msTree, ArrayList<ArrayList<SaleRoute>> salesPitch) {
  212. int numOfEdges[] = new int[salesPitch.size()];
  213. for (int i = 0; i < msTree.size(); i++) {
  214. numOfEdges[msTree.get(i).getSrc()]++;
  215. numOfEdges[msTree.get(i).getDest()]++;
  216. }
  217. OddDegreeList oddDegs = new OddDegreeList();
  218. for (int i = 0; i < numOfEdges.length; i++) {
  219. //System.out.println(numOfEdges[i]);
  220. if (numOfEdges[i] % 2 == 1) {
  221. oddDegs.add(i);
  222. }
  223. }
  224. Permutation permut = new Permutation(oddDegs.size(), 2);
  225. calcPairShortest(oddDegs, salesPitch, permut, 0, 0);
  226. permut.makePermutMinimal();
  227. //System.out.println(permut.tempCounter1);
  228. //System.out.println(permut.tempCounter2);
  229. //System.out.println("min cost: " + permut.getMinCost());
  230. //permut.printPermut();
  231. ArrayList<TreeEdge> oddEdges = new ArrayList<>();
  232. oddDegs.resetUsed();
  233. for(int i = 0; i < permut.size(); i++) {
  234. int offSet = permut.getValAtPos(i);
  235. addOddEdge(oddEdges, oddDegs, salesPitch, offSet);
  236. }
  237. addOddEdge(oddEdges, oddDegs, salesPitch, 0);
  238. return oddEdges;
  239. }
  240. private static void addOddEdge(ArrayList<TreeEdge> oddEdges, OddDegreeList oddDegs, ArrayList<ArrayList<SaleRoute>> salesPitch, int offSet) {
  241. int orig = oddDegs.getUnused(0);
  242. oddDegs.makeOffsetUsed(0);
  243. int dest = oddDegs.getUnused(offSet);
  244. oddDegs.makeOffsetUsed(offSet);
  245. oddEdges.add(new TreeEdge(orig, dest, salesPitch.get(orig).get(dest - orig - 1).getTotalCost()));
  246. }
  247. private static void calcPairShortest(OddDegreeList oddDegs, ArrayList<ArrayList<SaleRoute>> salesPitch, Permutation permut, int permutPos, int costSoFar) {
  248. while(true) {
  249. int offSet;
  250. if(permutPos == permut.size()) {
  251. offSet = 0;
  252. }
  253. else {
  254. offSet = permut.getValAtPos(permutPos);
  255. }
  256. int orig = oddDegs.getUnused(0);
  257. int dest = oddDegs.getUnused(1 + offSet);
  258. int edgeWeight = salesPitch.get(orig).get(dest - orig - 1).getTotalCost();
  259. int newCost = costSoFar + edgeWeight;
  260. //permut.tempCounter1++;
  261. /*System.out.println();
  262. System.out.println("permut pos: " + permutPos);
  263. System.out.println(orig + " : " + dest);
  264. System.out.println("newCost: " + newCost);
  265. permut.printPermut();*/
  266. if(newCost < permut.getMinCost()) {
  267. if(permutPos == permut.size()) {
  268. permut.setMinCost(newCost);
  269. /*System.out.println();
  270. System.out.println("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");
  271. System.out.println("New Min: " + newCost);
  272. System.out.println("!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!");*/
  273. }
  274. else {
  275. int used1 = oddDegs.makeOffsetUsed(0);
  276. int used2 = oddDegs.makeOffsetUsed(offSet);
  277. calcPairShortest(oddDegs, salesPitch, permut, permutPos+1, costSoFar+edgeWeight);
  278. oddDegs.makeUnused(used1);
  279. oddDegs.makeUnused(used2);
  280. }
  281. }
  282. if(permutPos == permut.size()) {
  283. //permut.tempCounter2++;
  284. break;
  285. }
  286. if(permut.isPosAtMax(permutPos)) {
  287. break;
  288. }
  289. permut.increaseAtPos(permutPos);
  290. /*permut.increaseAtPos(3);
  291. permut.printPermut();
  292. permut.increaseAtPos(2);
  293. permut.printPermut();
  294. permut.increaseAtPos(2);
  295. permut.printPermut();*/
  296. }
  297. }
  298. public static ArrayList<TreeEdge> makeMSTree(ArrayList<ArrayList<SaleRoute>> salesPitch) {
  299. ArrayList<TreeEdge> allEdges = new ArrayList<>();
  300. ArrayList<TreeEdge> msTree = new ArrayList<>();
  301. int vertNum = salesPitch.size();
  302. for (int orig = 0; orig < salesPitch.size(); orig++) {
  303. for (int dest = 0; dest < salesPitch.get(orig).size(); dest++) {
  304. allEdges.add(new TreeEdge(orig, dest + 1 + orig, salesPitch.get(orig).get(dest).getTotalCost()));
  305. //System.out.println(edge.getSrc() + " " + edge.getDest() + " " + edge.getWeight());
  306. //System.out.println(salesPitch.get(orig).get(dest).getTotalCost());
  307. }
  308. }
  309. //System.out.println(vertNum);
  310. Collections.sort(allEdges);
  311. while (msTree.size() < vertNum - 1) {
  312. if (notCycle(msTree, allEdges.get(0))) {
  313. msTree.add(allEdges.get(0));
  314. //System.out.println(allEdges.get(0).getSrc() + " " + allEdges.get(0).getDest() + " " + allEdges.get(0).getWeight());
  315. }
  316. allEdges.remove(0);
  317. }
  318. return msTree;
  319. }
  320. private static void qEdges(ArrayList<TreeEdge> tree, ArrayList<TreeEdge> edgeQ, int vertex) {
  321. for(int i = 0; i < tree.size(); i++) {
  322. TreeEdge edge = tree.get(i);
  323. if(edge.getSrc() == vertex || edge.getDest() == vertex) {
  324. if(!edge.isChecked()) {
  325. edgeQ.add(edge);
  326. //System.out.println("src: " + edge.getSrc() + "; dest: " + edge.getDest());
  327. edge.setChecked(true);
  328. }
  329. }
  330. }
  331. }
  332. private static boolean notCycle(ArrayList<TreeEdge> tree, TreeEdge additon) {
  333. ArrayList<TreeEdge> edgeQ = new ArrayList<>();
  334. int dest = additon.getDest();
  335. resetEdges(tree);
  336. qEdges(tree, edgeQ, additon.getSrc());
  337. while(edgeQ.size() > 0) {
  338. TreeEdge edge = edgeQ.get(0);
  339. edgeQ.remove(0);
  340. if(edge.getSrc() == dest || edge.getDest() == dest) {
  341. return false;
  342. }
  343. else {
  344. qEdges(tree, edgeQ, edge.getSrc());
  345. qEdges(tree, edgeQ, edge.getDest());
  346. }
  347. }
  348. return true;
  349. }
  350. private static void resetEdges(ArrayList<TreeEdge> tree) {
  351. for(int i = 0; i < tree.size(); i++) {
  352. tree.get(i).setChecked(false);
  353. }
  354. }
  355. }