Advertisement
Guest User

Untitled

a guest
May 22nd, 2019
160
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 5.59 KB | None | 0 0
  1. n = 480 (*n==8,48,480,5760,...*)
  2.  
  3. If[n == 8,
  4. Print["n=8"];
  5. v = {11, 2, 1, 8, 7, 14, 13, 4};
  6. penultimateValue = v[[n - 1]]; (*13*)
  7. constant1 = Total[v]; (*(60) sum of all the values*)
  8. constant2 =
  9. constant1/(n/2); (*(15) sum of pairwise differences (OEIS A060753)*)
  10.  
  11.  
  12. constant3 =
  13. 7; (*absolute value of min value of pairwise differences, and also \
  14. the first value minus last value, and also the sum of pairwise \
  15. differences*)
  16. constant4 = 35; (*sum of absolute values of pairwise differences*)
  17. constant5 = 16;(*16=2 x 2 x 2 x 2*)
  18. constant6 = 9; (*max value of pairwise differences*)
  19. constant7 =
  20. 1; (*value from the sequence 1,3,15,45,45,405,4455,... given by \
  21. A005867(n-1)/2^A057773(n-1). used for constant8=constant1/constant7*)
  22.  
  23.  
  24. constant8 = constant1/constant7 (*(60)*)
  25. ]
  26.  
  27. If[n == 48,
  28. Print["n=48"];
  29. v = {27, -18, 1, 4, 23, 26, 13, 32, 19, 22, 41, 44, 31, 18, 37, 24,
  30. 27, 46, 33, 36, 23, -6, -3, 16, 19, 38, 41, 12, -1, 2, -11, 8,
  31. 11, -2, 17, 4, -9, -6, 13, 16, 3, 22, 9, 12, 31, 34, 53, 8};
  32. penultimateValue = v[[n - 1]]; (*53*)
  33. constant1 = Total[v]; (*(840)*)
  34. constant2 = constant1/(n/2); (*(35)*)
  35. constant3 = 19;
  36. constant4 = 641;
  37. constant5 = 192;(*192=2 x 2 x 2 x 2 x 2 x 2 x 3*)
  38. constant6 = 45;
  39. constant7 = 3;
  40. constant8 = constant1/constant7 (*(280)*)
  41. ]
  42.  
  43. If[n == 480,
  44. Print["n=480"];
  45. v = {61, -54, -41, 4, 17, -2, 43, 24, 37, 82, 95, 76, 57, 102, 83,
  46. 96, 141, 122, 135, 116, 65, 78, 123, 136, 181, 194, 47, 60, 41, 86,
  47. 3, 48, 29, 10, 23, 68, 81, 62, 107, 24, 69, 82, 127,
  48. 12, -71, -26, -13, 32, 45, 26, 71, 52, 65, 46, 27, 8, 53, 34, 47,
  49. 92, 73, 86, 67, 16, 29, 74, 87, 68,
  50. 17, -2, -85, -40, -27, -46, -1, -20, -39, -26, 19, 32, 13, 58, 39,
  51. 52, 97, 78, -5, 40, -43, 2, 15, 60, 73, 54, 3, 16, 61, 74, -41,
  52. 4, -15, -2, 43, 24, 37, 18, -97, -52, -39, 6,
  53. 19, -32, -51, -38, -57, -12, 1, -18, 27, 8, -75, -30, -17, -36,
  54. 9, -10, 3, 48, 61, 106, 23, 68, -15, 30, 43, 24, 5, 50, 31, 12, 25,
  55. 6, -13, 32, 13, 26, 71, 52, 65, 46, -5, 8, 53, 34,
  56. 47, -4, -23, -10, -29, 16, 29,
  57. 10, -41, -60, -47, -2, -85, -40, -59, -46, -1, 12, 57, -26,
  58. 19, -64, -19, -6, 39, 52, 1, -18, -5, 40, 53, 34, 15, 60, 41, 54,
  59. 3, 16, -3, -54, -41, 4, 17, 62, 75, 24, 5, 18, -1, -20, -39,
  60. 6, -13, -32, -19, 26, 39, 20, 65, 46, 59, 104, 117, 162, 79, 124,
  61. 41, 86, 67, 80, 61, 106, 87, 100, 145, 158, 139, 120, 69, 82, 127,
  62. 108, 25, -26, -13, 32, 45, 90, 103, 52, -31, -50, -5,
  63. 8, -43, -62, -81, -68, -23, -10, -29, 16, -3, 10, -9,
  64. 36, -47, -2, -85, -40, -27, 18, 31, 12, 57, 38, 51, 96, 109, 90,
  65. 71, 116, 97, 78, 59, 72, 53, 2, 15, 60, 73, 118, 131, 80, 61, 74,
  66. 23, 36, 17, 62, 43, 24, 37, 82, 95, 76, 25, 38, 83, 96, 141, 58,
  67. 103, 20, 65, 78, 123, 136, 117, 162, 79, 124, 137, 118, 67, 48, 61,
  68. 106, 87, 100, 81, 30, 43, 24, 69, 82, 31, 12, 25, 6, 51, 64, 45,
  69. 90, 71, 52, 65, 46, 27, 72, 53, 34, 47, 92, 9, 54, -29, 16, 29, 74,
  70. 87, 68, 113, 94, 107, 152, 69, 50, 95, 76, 89, 134, 115, 128, 109,
  71. 58, 71, 116, 129, 174, 59, 40, 53, 34, 79, 92, 73, 118, 3, 16, 61,
  72. 74, 23, 4, 17, 62, 75, 120, 37, 82, -1, -20, 25, 38, 19, 64, 45,
  73. 58, 103, 116, 97, 78, 123, 104, 117, 162, 79, 60, 9, -10, 3, 48,
  74. 61, 10, -9, 4, -15, 30, 43, 24, 69, 50, 31, 12, 25, 6, 51, 32, 45,
  75. 90, 103, 148, 65, -50, -5, 8, 53, -30, 15, -4, 9, 54, 67, 48, 29,
  76. 74, -9, 36, 17, 30, -117, -104, -59, -46, -1,
  77. 12, -39, -58, -45, -64, -19, -6, -25, 20, 1, -18, -5, 40, 53, 34,
  78. 79, 60, 73, 118, 131, 16};
  79. penultimateValue = v[[n - 1]]; (*131*)
  80. constant1 = Total[v]; (*(18480)*)
  81. constant2 = constant1/(n/2);(*(77)*)
  82. constant3 = 45;
  83. constant4 = 15615;
  84. constant5 =
  85. 3840;(*3840=2 x 2 x 2 x 2 x 2 x 2 x 2 x 2 x 3 x 5, (3x5=15)*)
  86. constant6 = 147;
  87. constant7 = 15;
  88. constant8 = constant1/constant7(*(1232)*)
  89. ]
  90.  
  91. If[n == 5760,
  92. Print["n=5760"];
  93. v = {1, 2, 3, 4, 5}; (*add the correct 5760 values here later*)
  94. penultimateValue = v[[n - 1]]; (*2263*)
  95. constant1 = Total[v]; (*(2882880)*)
  96. constant2 = constant1/(n/2); (*(1001)*)
  97. constant3 = 617;
  98. constant4 = 2523883;
  99. constant5 =
  100. 552960;(*552960=2 x 2 x 2 x 2 x 2 x 2 x 2 x 2 x 2 x 2 x 2 x 2 x 3 x \
  101. 3 x 3 x 5, (3x3x3x5=135)*)
  102. constant6 = 3223;
  103. constant7 = 45;
  104. constant8 = constant1/constant7(*(64064)*)
  105. ]
  106.  
  107.  
  108. Total[Total[
  109. Table[v[[i]] - Sort[Table[v[[i]], {i, 1, n}]], {i, 1, n}]]] == 0
  110. (*possibly always true*)
  111.  
  112. Abs[Total[Table[v[[i]], {i, 1, n/2}]] -
  113. Total[Table[v[[i]], {i, n/2 + 1, n}]]] == constant5
  114.  
  115. v[[n - 1]] == penultimateValue
  116. v[[n - 1]] > 0
  117.  
  118. Table[v[[i]] + v[[n - i + 1]] == constant2, {i, 1, n/2}]
  119.  
  120. Total[Table[v[[i]], {i, 1, n}]] == constant1
  121.  
  122. Total[Table[v[[i]], {i, 1, n}]] ==
  123. Total[Table[v[[i]], {i, n/2 - n/4 + 1, n/2 + n/4}]]*2
  124. Total[Table[v[[i]], {i, n/2 - n/4 + 1, n/2 + n/4}]] ==
  125. Total[Table[v[[i]], {i, n/2 - n/8 + 1, n/2 + n/8}]]*2
  126.  
  127. (*for n=480 etc: (can add more of these)*)
  128. If[n >= 480,
  129. Total[Table[v[[i]], {i, n/2 - n/8 + 1, n/2 + n/8}]] ==
  130. Total[Table[v[[i]], {i, n/2 - n/16 + 1, n/2 + n/16}]]*2
  131. ]
  132.  
  133. Table[v[[i]] - v[[i + 1]] == v[[n - i]] - v[[n - i + 1]], {i, 1,
  134. n/2 - 1}]
  135.  
  136. v[[n/2]] - v[[n/2 + 1]] == v[[n/2 - 2]] - v[[n/2 - 1]]
  137. v[[n/2]] - v[[n/2 + 1]] == v[[n/2 + 2]] - v[[n/2 + 3]]
  138. v[[1]] - v[[n]] == constant3
  139. v[[n]] - v[[1]] == -constant3
  140.  
  141. Total[Abs[Table[v[[i]] - v[[i + 1]], {i, 1, n - 1}]]] == constant4
  142.  
  143. Total[Table[v[[i]] - v[[i + 1]], {i, 1, n - 1}]] == constant3
  144.  
  145. Abs[Min[Table[v[[i]] - v[[i + 1]], {i, 1, n - 1}]]] == constant3
  146.  
  147. Min[Table[v[[i]] - v[[i + 1]], {i, 1, n - 1}]] == -constant3
  148.  
  149. Max[Table[v[[i]] - v[[i + 1]], {i, 1, n - 1}]] == constant6
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement