Advertisement
MatsGranvik

Relation to square roots Linear programming

Jul 1st, 2019
290
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.83 KB | None | 0 0
  1. (*program 1*)
  2. a = 1;
  3. nn = 80;
  4. TableForm[
  5. L2 = Table[
  6. LinearProgramming[
  7. Table[1/n, {n, 1, k}], {Table[If[n == 1, k, 1], {n, 1, k}]}, {{1,
  8. 0}}, Table[
  9. If[n == 1, {-1, 1}, {-a*(n - 1), 0* (n - 1)}], {n, 1, k}]], {k,
  10. 1, nn}]];
  11. t1 = Table[Sum[If[L2[[n, k]] <= -1, 1, 0], {k, 2, n}], {n, 2, nn}]
  12. t2 = Table[Round[Sqrt[2/a*(n - 1)]], {n, 2, nn}]
  13. t1 - t2
  14. Count[%, 0] - Length[t2]
  15. (*end*)
  16.  
  17. (*program 2*)
  18. a = 2;
  19. nn = 80;
  20. TableForm[
  21. L2 = Table[
  22. LinearProgramming[
  23. Table[1/n, {n, 1, k}], {Table[If[n == 1, k, 1], {n, 1, k}]}, {{1,
  24. 0}}, Table[
  25. If[n == 1, {-1, 1}, {-a*(n - 1), 0* (n - 1)}], {n, 1, k}]], {k,
  26. 1, nn}]];
  27. t1 = Table[Sum[If[L2[[n, k]] <= -1, 1, 0], {k, 2, n}], {n, 2, nn}]
  28. t2 = Table[Round[Sqrt[2/a*(n - 1)]], {n, 2, nn}]
  29. t1 - t2
  30. Count[%, 0] - Length[t2]
  31. (*end*)
  32.  
  33. (*program 3*)
  34. a = 3;
  35. nn = 80;
  36. TableForm[
  37. L2 = Table[
  38. LinearProgramming[
  39. Table[1/n, {n, 1, k}], {Table[If[n == 1, k, 1], {n, 1, k}]}, {{1,
  40. 0}},
  41. Table[If[n == 1, {-1, 1}, {-a*(n - 1), 0* (n - 1)}], {n, 1,
  42. k}]], {k, 1, nn}]];
  43. t1 = Table[Sum[If[L2[[n, k]] <= -1, 1, 0], {k, 2, n}], {n, 2, nn}]
  44. t2 = Table[Round[Sqrt[2/a*(n - 1)]], {n, 2, nn}]
  45. t1 - t2
  46. Count[%, 0] - Length[t2]
  47. (*end*)
  48.  
  49.  
  50.  
  51. (*start*)
  52. Clear[nn];
  53. nn = 12
  54. (*Linear programming problem*)
  55. TableForm[
  56. L2 = Table[
  57. LinearProgramming[
  58. Table[1/n, {n, 1, k}], {Table[If[n == 1, k, 1], {n, 1, k}]}, {{1,
  59. 0}}, Table[
  60. If[n == 1, {-1, 1}, {-2 (n - 1), 0 (n - 1)}], {n, 1, k}]], {k,
  61. 1, nn}]];
  62. t1 = Table[Sum[L2[[n, k]], {k, 2, n}], {n, 1, nn}]
  63.  
  64. (*arithmetic sum*)
  65. TableForm[
  66. B = Table[
  67. Table[If[n <= b, 1, 0]*
  68. Sum[If[Mod[n, k] == 0,
  69. Sum[If[a >= n/k, -If[Mod[a, n/k] == 0, 1 - n/k, 1]*
  70. MoebiusMu[n/k], 0], {a, 1, b}], 0], {k, 1, nn}], {b, 1,
  71. nn}], {n, 1, nn}]]
  72. Total[B]
  73. (*end*)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement