Advertisement
MatsGranvik

Combinatorial approach to linear programming problem

Feb 2nd, 2020
220
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.64 KB | None | 0 0
  1. (*start*)
  2. Clear[A, B, a];
  3. nn = 7;
  4. B = Table[
  5. Tuples[Table[
  6. Table[If[k == 1, m, n], {n, -(k - 1), k - 1}], {k, 1, m}]], {m,
  7. 1, nn}];
  8. B1 = Table[
  9. DeleteCases[
  10. Table[If[Total[B[[n, k]]] == 1, B[[n, k]], a], {k, 1,
  11. Length[B[[n]]]}], a], {n, 1, nn}];
  12.  
  13. Clear[T, n, k, a];
  14. a[n_] := If[n < 1, 0, Sum[d MoebiusMu@d, {d, Divisors[n]}]]
  15. TableForm[
  16. M = Table[
  17. Table[Sum[If[n >= k, a[GCD[n, k]], 0], {n, 1, m}], {k, 1,
  18. nn}], {m, 1, nn}]];
  19. M1 = Table[Table[M[[n, k]], {k, 1, n}], {n, 1, nn}] ;
  20. Table[Flatten[Position[B, M1[[n]]]][[2]], {n, 1, nn}]
  21. Table[Flatten[Position[B1, M1[[n]]]][[2]], {n, 1, nn}]
  22. (*end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement