找回密码
 欢迎注册
楼主: northwolves

[求助] 六次等幂和

[复制链接]
发表于 2026-4-29 15:50:06 | 显示全部楼层
n=44, m=88
76 73 72 71 70 69 68 65 57 56 53 52 50 49 48 47 46 43 42 41 39 37 35 32 31 30 29 27 25 21 20 18 17 16 15 14 13 10 9 7 6 4 2 1

评分

参与人数 1威望 +8 金币 +8 贡献 +8 经验 +8 鲜花 +8 收起 理由
northwolves + 8 + 8 + 8 + 8 + 8 赞一个!

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2026-4-29 20:36:50 | 显示全部楼层
n=53, m=93
88 77 70 68 66 65 64 63 62 58 57 56 55 54 53 52 51 50 49 48 46 45 44 42 40 38 36 35 34 33 31 30 29 28 26 23 22 21 19 17 15 14 13 12 11 9 8 7 6 5 4 3 1

n=54,m=105
99 88 86 77 70 69 68 67 66 62 60 59 58 57 56 55 54 53 52 51 50 49 48 46 45 44 42 41 40 39 38 36 35 34 33 32 31 29 28 27 24 22 21 18 17 15 13 12 11 8 6 4 3 1

n=55,m=99
90 82 76 75 74 73 72 71 70 69 68 67 65 62 58 57 56 54 53 52 51 50 49 48 47 45 43 41 40 39 38 37 36 35 34 32 31 30 29 28 26 25 24 19 17 16 15 14 10 8 6 4 3 2 1

n=56,m=97
81 80 79 78 76 75 72 71 70 68 67 65 64 63 62 61 60 58 56 54 53 51 49 48 47 46 45 42 41 40 39 38 37 36 34 32 30 28 27 25 21 20 18 17 15 14 12 10 9 8 6 5 4 3 2 1

n=57,m=94
84 75 74 73 72 71 70 69 68 67 64 62 61 60 59 58 57 53 52 50 48 47 46 45 44 42 40 39 38 37 36 35 34 32 31 30 29 28 27 26 25 24 23 20 19 17 16 14 12 10 9 8 6 5 4 2 1

n=58, m=93
80 78 76 74 72 70 68 67 65 64 63 62 61 60 58 57 56 54 53 51 48 46 45 44 42 40 39 38 37 36 35 32 30 28 27 26 25 24 23 21 20 19 18 17 16 15 14 13 12 11 10 9 7 6 5 4 3 1


毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2026-4-30 07:39:48 | 显示全部楼层
n=59 FOUND 95 81 79 78 75 74 73 72 70 67 65 63 62 61 60 59 57 56 55 54 53 50 48 47 46 45 44 43 39 38 35 34 32 31 30 29 28 27 26 25 24 23 21 19 18 17 16 15 14 13 11 10 9 8 7 5 4 3 2 1
n=60 FOUND 94 84 77 74 73 71 70 69 67 66 65 64 63 62 61 60 59 57 56 55 53 52 51 50 49 48 47 45 43 42 41 38 35 34 32 31 30 29 28 27 26 25 24 23 20 19 18 16 15 14 13 12 11 10 9 8 7 5 4 2 1
n=61 FOUND 96 86 81 78 75 72 71 69 68 66 64 63 60 59 57 56 55 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 34 32 31 30 29 28 27 26 24 23 22 21 20 19 18 17 15 13 12 11 9 8 7 6 5 3 1
n=62 FOUND 95 85 81 78 74 71 69 66 65 64 63 62 61 60 58 57 55 53 52 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 30 29 27 26 25 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 7 6 5 4 3 2
n=63 FOUND 96 85 79 77 76 75 72 71 68 66 65 64 63 62 61 59 57 56 55 53 52 50 49 48 47 46 45 44 43 42 40 39 38 36 35 33 32 31 29 28 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 4 2 1
n=64 FOUND 100 88 82 81 79 77 76 74 72 70 69 68 66 65 62 61 59 57 56 55 54 53 52 51 50 49 48 46 45 44 43 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 25 24 23 22 20 18 15 14 13 12 11 9 8 7 6 5 3 2 1
n=65 FOUND 107 99 88 84 79 78 77 75 74 72 71 70 69 68 67 66 65 64 63 62 61 59 58 56 55 54 53 52 51 48 47 46 45 44 43 42 40 39 38 34 33 32 29 28 27 26 25 24 23 22 21 19 18 17 16 15 14 13 11 10 9 7 5 4 3 2
n=66 FOUND 110 91 89 87 86 85 84 83 82 81 80 79 78 76 74 73 72 71 69 68 65 64 63 62 61 60 59 58 57 54 53 52 51 50 48 47 46 43 42 40 39 38 37 36 35 34 31 30 28 26 23 20 19 18 17 16 15 14 13 10 9 8 7 6 5 2 1

评分

参与人数 1威望 +8 金币 +8 贡献 +8 经验 +8 鲜花 +8 收起 理由
northwolves + 8 + 8 + 8 + 8 + 8 赞一个!

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2026-4-30 10:56:01 | 显示全部楼层
我的计算进度:
{67} >= 101
{68} >= 101
{69} >= 103
{70} >= 104
{71} >= 105
{72} >= 101
{73} >= 107
{74} >= 108
{75} >= 108
{76} >= 109
{77} >= 110
{78} >= 112
{79} >= 113
{80} >= 115
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2026-4-30 15:12:04 | 显示全部楼层
{80} >= 116,卡很久了
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2026-4-30 16:43:19 来自手机 | 显示全部楼层
12和11比较困难,感觉要为它们写定制代码才行。而可以利用的信息是左边模8和11都是0或1,方程右边每项同样都是0或1,区分0和1的数目情况,可以分成简单几类,每类都可以先搜索数目少的情况。看上去模11应该更加适合分类。
另外6次幂情况6个数,我写了个针对性版本,穷举到100000左右都没有找到解,而计算速度显示后面没增加10000左右的范围,计算时间需要翻倍,所以几乎很难再穷举下去了。当然还有一个问题时由于找不到任何候选解来建议,代码是否存在bug也比较难以确认。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2026-4-30 19:06:28 | 显示全部楼层
k=6,n<=6的定制程序可以轻松搜索完m<=100000的范围(前面代码搜索到5000左右就卡住了),可以并没有找到新的解,主要方案AI汇总如下
# k=6, n≤6 定制搜索代码总结

## 核心思想

对于本原解,m 不是 2, 3, 7 的倍数。利用 42=2×3×7 的因子分解,将 m^6 的分解问题转化为:先确定 1-3 个"特殊数"(满足特定模条件),再搜索剩余的 42 倍数。

## 五种 Case 分类

根据 x2(非2倍数), x3(非3倍数), x7(非7倍数) 三者是否合并:

| Case | 确定数 | 约束 | 余下搜索 | 候选数 |
|------|--------|------|---------|--------|
| 1: x2=x3=x7 | 1个x | x^6≡m^6 (mod 42^6) | ≤5个42倍数 | 144 |
| 2: x2=x3≠x7 | x23, x7 | x23:7倍数, x7:6倍数 | ≤4个42倍数 | 24×6 |
| 3: x2=x7≠x3 | x27, x3 | x27:3倍数, x3:14倍数 | ≤4个42倍数 | 24×6 |
| 4: x3=x7≠x2 | x37, x2 | x37:2倍数, x2:21倍数 | ≤4个42倍数 | 36×4 |
| 5: 全不同 | x2,x3,x7 | x2:21倍, x3:14倍, x7:6倍 | ≤3个42倍数 | 4×6×6 |

## 候选模数

| 变量 | 模数 | 值 | 候选残基数 |
|------|------|-----|-----------|
| Case1 x | 64×729×117649 | 5,489,031,744 | 144 |
| Case2 x23 | 64×729×7 | 326,592 | 24 |
| Case2 x7 | 6×117649 | 705,894 | 6 |
| Case3 x27 | 64×117649×3 | 22,588,608 | 24 |
| Case3 x3 | 14×729 | 10,206 | 6 |
| Case4 x37 | 729×117649×2 | 171,532,242 | 36 |
| Case4 x2 | 21×64 | 1,344 | 4 |
| Case5 | 同Case4 x2 + Case3 x3 + Case2 x7 | | 4×6×6 |

## 预计算表格

**mod 64 表**:8 个特解 `min_sol64[8]`,通过遍历 x=1..15(奇数)构建。4 个解为 `x0, 32-x0, x0+32, 64-x0`。

**mod 729 表**:81 个特解 `min_sol729[81]`,通过遍历 x=1..364(非3倍数)构建。6 个解为 `x0, x0+243, x0+486, 729-x0, 729-x0-243, 729-x0-486`。

**mod 117649 表**:16807 个特解 `min_sol7[16807]`,通过遍历 x=1..58824(非7倍数)构建。6 个解通过乘以基解 `{1, 34967, 34968, 82681, 82682, 117648}` 再取模得到。

## 搜索裁剪

`searchUpTo_count` 搜索不超过 h_max 个 42 倍数,裁剪条件:
- **mod 7/8/9 裁剪**:rem % d > h_max 时剪枝(每个数的 6 次方 mod d 只能贡献 0 或 1)
- **下界裁剪**:rem < T[h_max] 时无解
- **上下界优化**:maxX = rem^(1/6), minX = (rem/h_max)^(1/6)

## 两个版本

| | l128dec_k6_lm.cpp | gmp2dec_k6_lm.cpp |
|---|---|---|
| 数据类型 | `__int128_t` | `mpz_class` (GMP) |
| m 范围 | ≤ ~38万(T数组溢出限制) | 理论无限 |
| 速度 | 快 | 较慢(GMP 开销) |
| 多线程 | OpenMP `schedule(dynamic)` | 同 |

## 已修复的 Bug(从 l128 同步到 GMP)

1. `initMod64Tables`:遍历范围 `x <= 7` → `x <= 15`
2. `getX64Candidates`:解公式 `x0+16, x0+32, x0+48` → `x0, 32-x0, x0+32, 64-x0`
3. `searchUpTo_count`:删除 `multiple` 参数,改用下界 `lb = (rem/h_max)^(1/6)` 优化
4. `searchCase` 参数类型:`ll` → `long`(模数值在 long 范围内)
5. 删除多线程中的 `ensureT` 调用(避免竞争)

## 性能参考(l128 版本)

| m 范围 | 时间 |
|--------|------|
| 10,000 | 0.13s |
| 20,000 | 0.79s |
| 40,000 | 3.8s |
| 60,000 | 12.3s |
| 100,000 | ~5min |

对应代码
  1. #include <iostream>
  2. #include <vector>
  3. #include <cmath>
  4. #include <algorithm>
  5. #include <cstdio>
  6. #include <gmpxx.h>
  7. using namespace std;

  8. typedef __int128_t ll;

  9. ll power6(ll x) {
  10.   ll x2 = x * x;
  11.   ll x4 = x2 * x2;
  12.   return x4 * x2;
  13. }

  14. ostream& operator<<(ostream& os, ll n) {
  15.     if (n == 0) return os << '0';
  16.     if (n < 0) { os << '-'; n = -n; }
  17.     char buf[64]; int i = 0;
  18.     while (n > 0) { buf[i++] = char('0' + (n % 10)); n /= 10; }
  19.     while (i--) os << buf[i];
  20.     return os;
  21. }

  22. long power6_mod(ll x, long mod) {
  23.   x %= mod;
  24.   long x2 = (x * x) % mod;
  25.   long x4 = (x2 * x2) % mod;
  26.   return (x4 * x2) % mod;
  27. }

  28. ll ext_gcd(ll a, ll b, ll& x, ll& y) {
  29.   if (b == 0) { x = 1; y = 0; return a; }
  30.   ll g = ext_gcd(b, a % b, y, x);
  31.   y -= a / b * x;
  32.   return g;
  33. }

  34. ll crt_merge(ll r1, ll m1, ll r2, ll m2) {
  35.   ll x, y;
  36.   ll g = ext_gcd(m1, m2, x, y);
  37.   if ((r2 - r1) % g != 0) return -1;
  38.   ll M = m1 / g * m2;
  39.   ll t1 = ((r2 - r1) / g) % (m2 / g);
  40.   if (t1 < 0) t1 += (m2 / g);
  41.   ll t2 = (x % (m2 / g));
  42.   if (t2 < 0) t2 += (m2 / g);
  43.   ll ans = (r1 + t1 * t2 % (m2 / g) * m1) % M;
  44.   if (ans < 0) ans += M;
  45.   return ans;
  46. }

  47. vector<ll> T;
  48. unsigned long T_max = 0;

  49. void initT(unsigned long maxVal) {
  50.   T.resize(maxVal + 1);
  51.   T[0] = 0;
  52.   ll sum = 0;
  53.   for (unsigned long i = 1; i <= maxVal; i++) {
  54.     sum += power6((ll)i);
  55.     T[i] = sum;
  56.   }
  57.   T_max = maxVal;
  58. }

  59. void ensureT(unsigned long maxVal) {
  60.   if (maxVal <= T_max) return;
  61.   T.resize(maxVal + 1);
  62.   ll sum = (T_max > 0) ? T[T_max] : 0;
  63.   for (unsigned long i = T_max + 1; i <= maxVal; i++) {
  64.     sum += power6((ll)i);
  65.     T[i] = sum;
  66.   }
  67.   T_max = maxVal;
  68. }

  69. int n;
  70. long M64 = 64, M729 = 729, M7 = 117649, M42_6;

  71. int min_sol64[8];
  72. int min_sol729[81];
  73. int min_sol7[16807];
  74. long sol7_base[6] = {1, 34967, 34968, 82681, 82682, 117648};

  75. void initMod64Tables() {
  76.   for (int x = 1; x <= 15; x += 2) {
  77.     long m6 = power6_mod(x, 64);
  78.     int idx = (m6 - 1) / 8;
  79.     min_sol64[idx] = x;
  80.   }
  81. }

  82. void initMod729Tables() {
  83.   for (int x = 1; x < 365; x++) {
  84.     if (x % 3 == 0) continue;
  85.     long m6 = power6_mod(x, 729);
  86.     if (m6 % 9 == 1) {
  87.       int idx = (m6 - 1) / 9;
  88.       if (min_sol729[idx] == 0) min_sol729[idx] = x;
  89.     }
  90.   }
  91. }

  92. void initMod7Tables() {
  93.   for (int x = 1; x < 58825; x++) {
  94.     if (x % 7 == 0) continue;
  95.     long m6 = power6_mod(x, M7);
  96.     if (m6 % 7 == 1) {
  97.       int idx = (m6 - 1) / 7;
  98.       if (idx >= 0 && idx < 16807 && min_sol7[idx] == 0) min_sol7[idx] = x;
  99.     }
  100.   }
  101. }

  102. vector<long> getX64Candidates(long m6) {
  103.   vector<long> xs;
  104.   if (m6 == 0) {
  105.     for (long x = 2; x < 64; x += 2) xs.push_back(x);
  106.   } else if (m6 % 8 == 1) {
  107.     int idx = (m6 - 1) / 8;
  108.     long x0 = min_sol64[idx];
  109.     xs.push_back(x0);
  110.     xs.push_back(32 - x0);
  111.     xs.push_back(x0 + 32);
  112.     xs.push_back(64 - x0);
  113.   }
  114.   return xs;
  115. }

  116. vector<long> getX729Candidates(long m6) {
  117.   vector<long> xs;
  118.   if (m6 == 0) {
  119.     for (long x = 3; x < 729; x += 3) xs.push_back(x);
  120.   } else if (m6 % 9 == 1) {
  121.     int idx = (m6 - 1) / 9;
  122.     long x0 = min_sol729[idx];
  123.     if (x0 == 0) return xs;
  124.     xs.push_back(x0);
  125.     xs.push_back(x0 + 243);
  126.     xs.push_back(x0 + 486);
  127.     xs.push_back(729 - x0);
  128.     xs.push_back(729 - x0 - 243);
  129.     xs.push_back(729 - x0 - 486);
  130.   }
  131.   return xs;
  132. }

  133. vector<long> findX7Candidates(long m6) {
  134.   vector<long> xs;
  135.   if (m6 == 0) {
  136.     for (long x = 7; x < M7; x += 7) xs.push_back(x);
  137.   } else if (m6 % 7 == 1) {
  138.     int idx = (m6 - 1) / 7;
  139.     long x0 = min_sol7[idx];
  140.     if (x0 == 0) return xs;
  141.     for (int i = 0; i < 6; i++) {
  142.       long x = (x0 * sol7_base[i]) % M7;
  143.       xs.push_back(x);
  144.     }
  145.   }
  146.   return xs;
  147. }

  148. // Case 1: x2=x3=x7=x
  149. vector<long> generateCase1Candidates(long m6_mod64, long m6_mod729, long m6_mod7) {
  150.   vector<long> cands64 = getX64Candidates(m6_mod64);
  151.   vector<long> cands729 = getX729Candidates(m6_mod729);
  152.   vector<long> cands7 = findX7Candidates(m6_mod7);

  153.   vector<long> residues;

  154.   for (long r64 : cands64) {
  155.     for (long r729 : cands729) {
  156.       ll r1 = crt_merge(r64, M64, r729, M729);
  157.       if (r1 < 0) continue;
  158.       for (long r7 : cands7) {
  159.         ll r2 = crt_merge(r1, M64 * M729, r7, M7);
  160.         if (r2 < 0) continue;
  161.         residues.push_back((long)r2);
  162.       }
  163.     }
  164.   }
  165.   sort(residues.begin(), residues.end(), greater<ll>());
  166.   return residues;
  167. }

  168. // Case 2: x23 (7 multiple) + x7 (6 multiple)
  169. vector<long> generateCase2x23Candidates(long m6_mod64, long m6_mod729) {
  170.   vector<long> cands64 = getX64Candidates(m6_mod64);
  171.   vector<long> cands729 = getX729Candidates(m6_mod729);

  172.   ll M_x23 = M64 * M729;
  173.   vector<long> residues;

  174.   for (long r64 : cands64) {
  175.     for (long r729 : cands729) {
  176.       ll r = crt_merge(r64, M64, r729, M729);
  177.       if (r < 0) continue;
  178.       ll r7 = crt_merge(r, M_x23, 0, 7);
  179.       if (r7 < 0) continue;
  180.       residues.push_back((long)r7);
  181.     }
  182.   }
  183.   sort(residues.begin(), residues.end(), greater<ll>());
  184.   return residues;

  185. }

  186. vector<long> generateCase2x7Candidates(long m6_mod7) {
  187.   vector<long> cands7 = findX7Candidates(m6_mod7);

  188.   vector<long> residues;

  189.   for (long r7 : cands7) {
  190.     ll r = crt_merge(0, 6, r7, M7);
  191.     if (r < 0) continue;
  192.     residues.push_back((long)r);
  193.   }
  194.   sort(residues.begin(), residues.end(), greater<ll>());
  195.   return residues;

  196. }

  197. // Case 3: x27 (3 multiple) + x3 (14 multiple)
  198. vector<long> generateCase3x27Candidates(long m6_mod64, long m6_mod7) {
  199.   vector<long> cands64 = getX64Candidates(m6_mod64);
  200.   vector<long> cands7 = findX7Candidates(m6_mod7);

  201.   ll M_x27 = M64 * M7;
  202.   vector<long> residues;

  203.   for (long r64 : cands64) {
  204.     for (long r7 : cands7) {
  205.       ll r = crt_merge(r64, M64, r7, M7);
  206.       if (r < 0) continue;
  207.       ll r3 = crt_merge(r, M_x27, 0, 3);
  208.       if (r3 < 0) continue;
  209.       residues.push_back((long)r3);
  210.     }
  211.   }
  212.   sort(residues.begin(), residues.end(), greater<ll>());
  213.   return residues;
  214. }

  215. vector<long> generateCase3x3Candidates(long m6_mod729) {
  216.   vector<long> cands729 = getX729Candidates(m6_mod729);

  217.   vector<long> residues;

  218.   for (long r729 : cands729) {
  219.     ll r = crt_merge(0, 14, r729, M729);
  220.     if (r < 0) continue;
  221.     residues.push_back((long)r);
  222.   }

  223.   sort(residues.begin(), residues.end(), greater<ll>());
  224.   return residues;

  225. }

  226. // Case 4: x37 (2 multiple) + x2 (21 multiple)
  227. vector<long> generateCase4x37Candidates(long m6_mod729, long m6_mod7) {
  228.   vector<long> cands729 = getX729Candidates(m6_mod729);
  229.   vector<long> cands7 = findX7Candidates(m6_mod7);

  230.   ll M_x37 = M729 * M7;
  231.   vector<long> residues;

  232.   for (long r729 : cands729) {
  233.     for (long r7 : cands7) {
  234.       ll r = crt_merge(r729, M729, r7, M7);
  235.       if (r < 0) continue;
  236.       ll r2 = crt_merge(r, M_x37, 0, 2);
  237.       if (r2 < 0) continue;
  238.       residues.push_back((long)r2);
  239.     }
  240.   }
  241.   sort(residues.begin(), residues.end(), greater<ll>());
  242.   return residues;

  243. }

  244. vector<long> generateCase4x2Candidates(long m6_mod64) {
  245.   vector<long> cands64 = getX64Candidates(m6_mod64);

  246.   vector<long> residues;

  247.   for (long r64 : cands64) {
  248.     ll r = crt_merge(0, 21, r64, M64);
  249.     if (r < 0) continue;
  250.     residues.push_back((long)r);
  251.   }

  252.   sort(residues.begin(), residues.end(), greater<ll>());
  253.   return residues;
  254. }

  255. // Case 5: x2 (21), x3 (14), x7 (6) all different
  256. vector<long> generateCase5x2Candidates( long m6_mod64) {
  257.   return generateCase4x2Candidates( m6_mod64);  // Same logic
  258. }

  259. vector<long> generateCase5x3Candidates( long m6_mod729) {
  260.   return generateCase3x3Candidates(m6_mod729);  // Same logic
  261. }

  262. vector<long> generateCase5x7Candidates( long m6_mod7) {
  263.   return generateCase2x7Candidates( m6_mod7);  // Same logic
  264. }

  265. bool searchUpTo_count(ll rem, int h_max, vector<ll>& xs, bool first_round=true) {
  266.   if (rem == 0&&xs.size()>1) return true;
  267.   if (h_max == 0) return false;
  268.   if ((int)xs.size() >= n) return false;

  269.   ll r8 = rem % 8;
  270.   ll r9 = rem % 9;
  271.   ll r7v = rem % 7;

  272.   if (r8 > h_max) return false;
  273.   if (r9 > h_max) return false;
  274.   if (r7v > h_max) return false;

  275.   if (h_max > 0 && (unsigned long)h_max <= T_max && rem < T[(unsigned long)h_max]) return false;

  276.   ll prevX = first_round ? rem : xs.back();
  277.   double ub = pow((double)rem, 1.0/6);
  278.   long maxX = ceil(ub) + 1;
  279.   while (maxX > 0 && power6(maxX) > rem) maxX--;
  280.   if (maxX >= prevX) maxX = prevX - 1;
  281.   double lb = pow((double)rem/h_max, 1.0/6);

  282.   long minX = (long)ceil(lb);
  283.   if (minX > maxX) return false;

  284.   for (long x = maxX; x >= minX; x --) {
  285.     ll x6 = power6(x);
  286.     if (rem < x6) continue;

  287.     xs.push_back(x);

  288.     if (searchUpTo_count(rem - x6, h_max - 1, xs, false)) return true;

  289.     if (rem == x6) return true;

  290.     xs.pop_back();
  291.   }

  292.   return false;
  293. }

  294. bool searchCase(ll m, int case_num, long m6_mod64, long m6_mod729, long m6_mod7) {
  295.   ll mk = power6(m);
  296.   ll diff;
  297.   ll rem;
  298.   vector<ll> xs;

  299.   if (case_num == 1) {
  300.     vector<long> x_cands = generateCase1Candidates((long)m6_mod64, m6_mod729, m6_mod7);
  301.     if (x_cands.empty()) return false;

  302.     for (ll x_s : x_cands) {
  303.       ll M6=5489031744L;
  304.       for(ll x=x_s;x<m;x+=M6){
  305.         ll x6 = power6(x);
  306.         diff = mk - x6;
  307.         if(diff<=0)break;
  308.         if (diff % M42_6 != 0) {
  309.           cerr<<"Error1\n";
  310.           continue;
  311.         }
  312.         rem = diff / M42_6;

  313.         xs.clear();
  314.         xs.push_back(x);
  315.         if (searchUpTo_count(rem, n - 1, xs)) {
  316. #pragma omp critical
  317.           {
  318.             cout << "FOUND Case1 m=" << m << ": ";
  319.             for (ll xi : xs) cout << xi << " ";
  320.             cout << endl;
  321.             cout.flush();
  322.           }
  323.           return true;
  324.         }
  325.       }
  326.     }
  327.   }
  328.   else if (case_num == 2) {
  329.     vector<long> x23_cands = generateCase2x23Candidates(m6_mod64, m6_mod729);
  330.     vector<long> x7_cands = generateCase2x7Candidates(m6_mod7);
  331.     if (x23_cands.empty() || x7_cands.empty()) return false;
  332.     ll M23=326592;
  333.     for (ll x23_s : x23_cands) {
  334.       for(ll x23 = x23_s;x23<m;x23+=M23){
  335.         ll x23_6 = power6(x23);
  336.         for (ll x7_s : x7_cands) {
  337.           ll M7=705894L;
  338.           for(ll x7=x7_s;x7<m;x7+=M7){
  339.             if (x7 == x23) continue;
  340.             ll x7_6 = power6(x7);
  341.             diff = mk - x23_6 - x7_6;
  342.             if(diff<0)break;
  343.             if (diff % M42_6 != 0){
  344.               cerr<<"Error2\n";
  345.               continue;
  346.             }
  347.             rem = diff / M42_6;

  348.             xs.clear();
  349.             xs.push_back(x23);
  350.             xs.push_back(x7);
  351.             if (searchUpTo_count(rem, n - 2, xs)) {
  352. #pragma omp critical
  353.               {
  354.                 cout << "FOUND Case2 m=" << m << ": ";
  355.                 for (ll xi : xs) cout << xi << " ";
  356.                 cout << endl;
  357.                 cout.flush();
  358.               }
  359.               return true;
  360.             }
  361.           }
  362.         }
  363.       }
  364.     }
  365.   }
  366.   else if (case_num == 3) {
  367.     vector<long> x27_cands = generateCase3x27Candidates(m6_mod64, m6_mod7);
  368.     vector<long> x3_cands = generateCase3x3Candidates(m6_mod729);
  369.     if (x27_cands.empty() || x3_cands.empty()) return false;

  370.     for (ll x27_s : x27_cands) {
  371.       ll M27=22588608L;
  372.       for(ll x27=x27_s;x27<m;x27+=M27){
  373.         ll x27_6 = power6(x27);
  374.         for (ll x3_s : x3_cands) {
  375.           ll M3=10206;
  376.           for(ll x3=x3_s;x3<m;x3+=M3){
  377.             if (x3 == x27) continue;
  378.             ll x3_6 = power6(x3);
  379.             diff = mk - x27_6 - x3_6;
  380.             if(diff<0)break;
  381.             if (diff % M42_6 != 0) {
  382.               cerr<<"error3\n";
  383.               continue;
  384.             }
  385.             rem = diff / M42_6;

  386.             xs.clear();
  387.             xs.push_back(x27);
  388.             xs.push_back(x3);
  389.             if (searchUpTo_count(rem, n - 2, xs)) {
  390. #pragma omp critical
  391.               {
  392.                 cout << "FOUND Case3 m=" << m << ": ";
  393.                 for (ll xi : xs) cout << xi << " ";
  394.                 cout << endl;
  395.                 cout.flush();
  396.               }
  397.               return true;
  398.             }
  399.           }
  400.         }
  401.       }
  402.     }
  403.   }
  404.   else if (case_num == 4) {
  405.     vector<long> x37_cands = generateCase4x37Candidates(m6_mod729, m6_mod7);
  406.     vector<long> x2_cands = generateCase4x2Candidates(m6_mod64);
  407.     if (x37_cands.empty() || x2_cands.empty()) return false;

  408.     for (ll x37_s : x37_cands) {
  409.       ll M37=171532242L;
  410.       for(ll x37=x37_s; x37<m; x37+=M37){
  411.         ll x37_6 = power6(x37);
  412.         for (ll x2_s : x2_cands) {
  413.           ll M2=1344;
  414.           for(ll x2=x2_s; x2<m;x2+=M2){
  415.             if (x2 == x37) continue;
  416.             ll x2_6 = power6(x2);
  417.             diff = mk - x37_6 - x2_6;
  418.             if(diff<0)break;
  419.             if (diff % M42_6 != 0) {
  420.               cerr<<"error4\n";
  421.               continue;
  422.             }
  423.             rem = diff / M42_6;

  424.             xs.clear();
  425.             xs.push_back(x37);
  426.             xs.push_back(x2);
  427.             if (searchUpTo_count(rem, n - 2, xs)) {
  428. #pragma omp critical
  429.               {
  430.                 cout << "FOUND Case4 m=" << m << ": ";
  431.                 for (ll xi : xs) cout << xi << " ";
  432.                 cout << endl;
  433.                 cout.flush();
  434.               }
  435.               return true;
  436.             }
  437.           }
  438.         }
  439.       }
  440.     }
  441.   }
  442.   else if (case_num == 5) {
  443.     vector<long> x2_cands = generateCase5x2Candidates(m6_mod64);
  444.     vector<long> x3_cands = generateCase5x3Candidates(m6_mod729);
  445.     vector<long> x7_cands = generateCase5x7Candidates(m6_mod7);
  446.     if (x2_cands.empty() || x3_cands.empty() || x7_cands.empty()) return false;

  447.     for (ll x2_s : x2_cands) {
  448.       ll M2=1344;
  449.       for(ll x2=x2_s;x2<m;x2+=M2){
  450.         ll x2_6 = power6(x2);
  451.         for (ll x3_s : x3_cands) {
  452.           ll M3=10206;
  453.           for(ll x3=x3_s;x3<m;x3+=M3){
  454.             if (x3 == x2) continue;
  455.             ll x3_6 = power6(x3);
  456.             if(mk<x2_6+x3_6)break;
  457.             for (ll x7_s : x7_cands) {
  458.               ll M7=705894;
  459.               for(ll x7=x7_s;x7<m;x7+=M7){
  460.                 if (x7 == x3 || x7==x2) continue;
  461.                 ll x7_6 = power6(x7);
  462.                 diff = mk - x2_6 - x3_6 - x7_6;
  463.                 if(diff<0)break;
  464.                 if (diff % M42_6 != 0) {
  465.                   cerr<<"error5\n";
  466.                   continue;
  467.                 }
  468.                 rem = diff / M42_6;

  469.                 xs.clear();
  470.                 xs.push_back(x2);
  471.                 xs.push_back(x3);
  472.                 xs.push_back(x7);
  473.                 if (searchUpTo_count(rem, n - 3, xs)) {
  474. #pragma omp critical
  475.                   {
  476.                     cout << "FOUND Case5 m=" << m << ": ";
  477.                     for (ll xi : xs) cout << xi << " ";
  478.                     cout << endl;
  479.                     cout.flush();
  480.                   }
  481.                   return true;
  482.                 }
  483.               }
  484.             }
  485.           }
  486.         }
  487.       }
  488.     }
  489.   }
  490.   return false;
  491. }

  492. bool is_primitive_candidate(ll m) {
  493.   return (m % 2 != 0 && m % 3 != 0 && m % 7 != 0);
  494. }

  495. int main(int argc, char* argv[]) {
  496.   if (argc < 4) {
  497.     cerr << "Usage: " << argv[0] << " n m_start m_end" << endl;
  498.     return 1;
  499.   }

  500.   time_t start=time(NULL);

  501.   n = atoi(argv[1]);
  502.   long m_start = atol(argv[2]);
  503.   long m_end = atol(argv[3]);

  504.   M42_6 = power6(42);

  505.   printf("k=6 n=%d m_start=%ld m_end=%ld\n", n, (long)m_start, (long)m_end);
  506.   printf("42^6 = %ld\n", (long)M42_6);

  507.   initMod64Tables();
  508.   initMod729Tables();
  509.   initMod7Tables();

  510.   long T_needed = m_end + 100;
  511.   initT((unsigned long)T_needed);

  512.   printf("initT done, T_max=%lu\n", T_max);

  513.   bool found=false;
  514.   long minm=-1;

  515. #pragma omp parallel for schedule(dynamic)
  516.   for (long m = m_start; m <= m_end; m++) {
  517.     if (!is_primitive_candidate(m)) continue;

  518.     long m6_mod64 = power6_mod(m, 64);
  519.     long m6_mod729 = power6_mod(m, 729);
  520.     long m6_mod7 = power6_mod(m, M7);

  521.     for (int case_num = 1; case_num <= 5; case_num++) {
  522.       if(searchCase(m, case_num, m6_mod64, m6_mod729, m6_mod7)){
  523. #pragma omp critical
  524.         {
  525.           {
  526.             found=true;
  527.             if(minm<0||m<minm){
  528.               minm=m;
  529.             }
  530.           }
  531.         }
  532.       }
  533.     }
  534.     if(m%100==1)
  535. #pragma omp critical
  536.     {
  537.       long diff = time(NULL) -start;
  538.       cerr<<"m="<<m<<" processed cost "<<diff<<"s\n";
  539.       cerr.flush();
  540.     }
  541.   }

  542.   if(!found){
  543.     cout << "NO_SOLUTION in range" << endl;
  544.   }else{
  545.     cout << "Found best m = "<<minm<<endl;
  546.   }
  547.   return 1;
  548. }
复制代码
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2026-4-30 19:14:09 | 显示全部楼层
前面的总结不好,找了个中间版本的(可能包含部分错误)
k=6 大 m 搜索方案总结报告
一、问题背景
对于 k=6, n≤6 的本原解搜索,当 m 较大时(如 m > 5000),传统方法(模 7/8/9 裁剪)效率不足。需要利用 42=2×3×7 的因子分解特性进行更高效的剪枝。
二、数学基础
本原解必要条件
- m 不是 2, 3, 7 的倍数
- 对于每个解 x_i,要求 x_i^6 ≡ m^6 (mod 64/729/117649)
关键恒等式
对于 x^6 mod d:
- mod 64: 奇数 x → x^6 ≡ 1, 9, 17, 25, 33, 41, 49, 57(共8种,每种4解)
- mod 729: 非3倍数 x → x^6 ≡ 1 (mod 9),共81种,每种6解
- mod 117649: 非7倍数 x → x^6 ≡ 1 (mod 7),共16807种,每种6解
三、五种 Case 分类
根据 x2, x3, x7 是否相等,分为 5 种情况:
Case        描述        确定变量        倍数约束        余下搜索
1        x2=x3=x7        1个x        无        ≤5个42倍数
2        x2=x3≠x7        x23, x7        x23:7倍数, x7:6倍数        ≤4个42倍数
3        x2=x7≠x3        x27, x3        x27:3倍数, x3:14倍数        ≤4个42倍数
4        x3=x7≠x2        x37, x2        x37:2倍数, x2:21倍数        ≤4个42倍数
5        全不同        x2,x3,x7        x2:21, x3:14, x7:6        ≤3个42倍数
四、各 Case 的候选模数
Case 1: x2=x3=x7=x
模        值        候选数
M_total        64×729×117649 = 5489031744        144
最小候选        ~5.5×10^9        需 m > 5.5e9
Case 2: x23 (7倍数) + x7 (6倍数)
变量        模数        值        候选数        m范围影响
x23        M64×M729×7        326592        24        模<38万,需多重候选
x7        6×M7        705894        6        模>38万,单候选
Case 3: x27 (3倍数) + x3 (14倍数)
变量        模数        值        候选数        m范围影响
x27        M64×M7×3        22588608        24        模>38万,单候选
x3        14×M729        10206        6        模<38万,需多重候选
Case 4: x37 (2倍数) + x2 (21倍数)
变量        模数        值        候选数        m范围影响
x37        M729×M7×2        171532242        36        模>38万,单候选
x2        21×M64        1344        4        模<38万,需多重候选
Case 5: x2, x3, x7 全不同
变量        模数        值        候选数        m范围影响
x2        21×M64        1344        4        模<38万,需多重候选
x3        14×M729        10206        6        模<38万,需多重候选
x7        6×M7        705894        6        模>38万,单候选
五、模数总结表
模数        值        适用Case        38万范围内
42^6        5489031744        Case1全部候选        > m,单候选
326592        3.3×10^5        Case2 x23        < m,需循环生成
705894        7.1×10^5        Case2/5 x7        > m,单候选
22588608        2.3×10^7        Case3 x27        > m,单候选
10206        1.0×10^4        Case3/5 x3        < m,需循环生成
171532242        1.7×10^7        Case4 x37        > m,单候选
1344        1.3×10^3        Case4/5 x2        < m,需循环生成
六、候选生成策略
当模数 < m 时
需要生成所有候选:对每个 residue r,生成 x = r + k×mod for x < m
例如 Case 4 x2(模1344):
- residue 可能有 4 种
- 对于 m = 100000,每个 residue 生成约 100000/1344 ≈ 74 个候选
- 总候选数约 4×74 = 296 个
当模数 > m 时
每个 residue 只有一个候选(如果 residue < m)
七、搜索复杂度对比
m范围        单线程时间        适合的Case
60000        27秒        Case2-5有效
380000        预估~3分钟        Case2-5有效
>5.5e9        Case1才有效        Case1开始有候选
八、优化建议
1. 多进程策略:按 m 范围分块,每块一个进程
2. 候选预计算:mod64/mod729 表格查询 O(1),mod7 表格大小16807
3. 剪枝优化:
   - rem % 8/9/7 > h_max 时剪枝
   - rem < Th_max 时剪枝(需T数组不溢出)
4. 输出处理:使用 cout.setf(ios::unitbuf) 或 cout << flush 确保即时输出
九、当前测试结果
- m < 60000: 27秒,无解
- m < 380000: 多进程版本已编译,待测试
- Case 1-5 候选验证: 符合预期(144, 24×6, 24×6, 36×4, 4×6×6)
十、下一步
1. 测试 m < 380000 完整范围
2. 如果无解,考虑更大的 m 范围(Case 1 开始有效时)
3. 考虑 n=5 或更小 n 的情况
==================================
再添加一个初始描述内容:
gmp2dec_fork2.cpp对于k=6的场景使用了模8,模7和模9三种操作,但是对于n<=6加超大m如m>5000后搜索速度还是太慢了。我们现在目标是定制修改gmp2dec_fork2,处理k=6,n<=6的场景,希望能够处理超大的m. 对于每个m,如果是本源解,我们直到m必须不是2,3,7的倍数。然后我们可以首先找到一个小于m不是2的倍数的数x2使得m^6-x2^6是2^6=64的倍数。同样我们需要找一个小于m不是3的倍数的数x3使得m^6-x3^6是3^6的倍数,同样还需要找一个x7使得m^6-x7^6是7^6的倍数。这里x2,x3,x7之间可以有重叠,也可以完全不同。而除了x2,x3,x7以外,余下最多还有5个数,要求它们必须同时是2*3*7=42的倍数。而同样,如果x2和x3不同,那么我们要求x2是3的倍数,x3是2的倍数。同样如果x2和x7不同,x2必须是7的倍数,x7必须是2的倍数等等。当然还可以有x2,x3,x7都相同的情况,对应m^6-x2^6是42^6的倍数。所以根据三者是否两两相等共有8大类。其中每一类对应有1-3个已经确定的数。比如x2=x3!=x7, 那么我们已经确定(x2,x7)两个数,于是余下就变成我们求不超过4个数,使得它们的6次方之和等于(m^6-x2^6-x7^6)/42^6, 这时可以继续使用原先gmp2dec_fork2中递归搜索部分的逻辑去搜索,区别是原先搜索固定数目的数,现在改为不超过4个数,任何时候达到目的就报告结果。请帮忙设计方案

评分

参与人数 1威望 +8 金币 +8 贡献 +8 经验 +8 鲜花 +8 收起 理由
northwolves + 8 + 8 + 8 + 8 + 8 膜拜大神!赞一个!

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2026-5-2 16:26:50 | 显示全部楼层
先记录一下现在找到的7次方结果(A395180)

  1. n=6, m>=2168

  2. n=7,m=568
  3. 525 439 430 413 266 258 127

  4. n=8,m=102
  5. 90 85 83 64 58 53 35 12

  6. n=9, m=62
  7. 59 50 41 33 27 22 20 14 6

  8. n=10, m=89
  9. 77 75 70 65 52 48 35 29 12 4

  10. n=11, m=54
  11. 49 47 39 27 24 21 5 4 3 2 1

  12. n=12, m=57
  13. 55 43 37 33 27 26 23 18 17 15 8 7

  14. n=13, m=49
  15. 46 36 34 33 32 30 26 24 12 10 9 7 2

  16. n=14, m=47
  17. 43 40 35 22 21 20 18 17 11 9 8 7 4 2

  18. n=15, m=40
  19. 39 28 25 22 21 20 18 17 16 14 12 9 5 3 1

  20. n=16, m=45
  21. 37 36 35 34 33 31 25 24 20 19 17 14 5 4 3 2

  22. n=17, m=43
  23. 38 34 33 30 29 27 26 22 19 18 17 15 13 9 4 2 1

  24. n=18, m=43
  25. 37 36 32 31 29 27 24 22 19 18 16 13 11 8 6 4 3 1

  26. n=19, m=44
  27. 38 36 33 31 30 29 27 22 21 20 18 17 14 13 11 7 6 4 3

  28. n=20, m=46
  29. 42 36 34 30 29 28 27 25 22 20 19 16 15 14 7 6 5 4 2 1

  30. n=21, m=45
  31. 40 36 33 32 31 27 26 24 22 19 16 15 12 11 10 8 7 5 4 2 1

  32. n=22, m=44
  33. 37 34 33 32 31 30 29 26 25 24 23 22 21 20 16 15 13 9 8 7 5 4

  34. n=23, m=47
  35. 40 36 35 34 33 32 31 30 26 25 23 21 20 17 14 12 11 8 7 6 3 2 1

  36. n=24, m=46
  37. 39 36 35 34 32 30 29 28 25 24 22 20 19 18 16 14 12 9 8 6 4 3 2 1

  38. n=25, m=48
  39. 42 38 37 35 31 29 28 27 25 23 22 20 18 14 13 12 11 9 8 7 6 5 4 3 1

  40. n=26, m=49
  41. 41 40 39 36 31 30 29 28 27 25 23 21 19 18 17 14 13 12 10 9 8 7 6 4 3 1

  42. n=27, m=47
  43. 43 35 33 32 31 29 28 27 26 25 24 22 20 19 18 17 16 13 11 9 8 7 6 4 3 2 1

  44. n=28, m=47
  45. 41 39 33 32 31 30 29 28 25 24 22 21 20 18 16 15 14 13 12 10 8 7 6 5 4 3 2 1

  46. n=29, m=52
  47. 43 42 40 38 36 33 32 31 30 29 26 25 24 22 21 20 19 15 13 12 11 8 7 6 5 4 3 2 1

  48. n=30, m=55
  49. 46 43 42 41 37 36 34 33 32 30 29 28 27 26 25 24 23 18 17 16 14 13 12 10 8 7 5 4 3 2

  50. n=31, m=52
  51. 42 41 39 38 37 35 33 32 31 30 29 28 27 26 25 23 22 21 20 18 17 14 12 10 9 8 5 4 3 2 1

  52. n=32, m=57
  53. 46 44 43 42 40 38 37 36 34 33 32 31 30 29 28 27 26 25 24 23 20 19 18 17 16 15 12 9 8 5 4 2

  54. n=33, m=57
  55. 45 44 43 42 40 39 38 37 35 34 33 31 30 27 26 25 24 23 22 21 20 19 18 17 14 13 12 11 9 8 7 5 1

  56. n=34, m=57
  57. 47 44 43 41 40 38 37 35 34 33 32 31 30 29 27 26 25 21 20 17 16 15 14 13 12 11 10 8 7 5 4 3 2 1

  58. n=35, m=60
  59. 48 46 45 44 42 40 39 38 37 36 35 34 32 31 30 29 28 27 26 25 23 21 20 19 17 16 14 13 11 10 9 6 4 3 2

  60. n=36, m=61
  61. 47 46 45 44 43 42 41 40 39 38 37 36 34 33 32 31 30 28 27 26 24 22 21 19 17 16 13 12 11 10 9 8 7 6 5 4

  62. n=37, m=61
  63. 49 47 46 44 43 42 39 38 37 36 35 34 33 32 31 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 11 10 7 6 5 2 1

  64. n=38, m=64
  65. 55 51 49 46 43 40 37 36 35 33 32 31 30 29 28 27 26 24 23 22 21 20 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 1

  66. n=39, m=65
  67. 57 50 48 46 43 41 39 38 37 36 35 34 33 31 30 28 27 25 24 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 6 5 4 3 2

  68. n=40, m=65
  69. 54 49 48 47 46 45 44 41 38 37 35 33 32 31 30 29 28 27 26 25 24 22 21 20 19 18 17 16 13 12 11 10 9 8 7 6 5 3 2 1

  70. n=41, m=69
  71. 55 52 51 50 49 48 45 44 42 41 40 39 38 37 36 35 34 33 31 30 28 26 24 23 22 21 20 17 16 14 13 12 11 10 8 7 6 5 3 2 1

  72. n=42, m=69
  73. 56 54 53 51 46 45 44 43 41 40 39 38 36 34 33 32 31 29 28 27 26 24 23 21 20 18 17 16 15 14 13 12 11 10 9 7 6 5 4 3 2 1

  74. n=43, m=70
  75. 58 54 52 50 48 47 44 43 42 41 40 39 38 37 35 33 32 31 30 29 28 27 26 25 24 23 21 20 18 17 16 15 13 12 11 9 8 7 6 5 4 3 1

  76. n=44, m=69
  77. 53 52 51 50 49 48 47 45 44 43 42 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 22 21 19 18 17 16 15 14 12 11 10 9 8 7 6 5 3 1

  78. n=45, m=72
  79. 61 56 52 51 50 49 48 43 41 40 39 38 35 34 33 32 31 30 29 28 27 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 7 6 5 4 3 2 1

  80. n=46, m=75
  81. 62 56 54 52 51 50 49 48 47 46 45 43 42 41 40 39 38 36 35 34 33 32 30 29 28 27 26 25 24 23 20 18 17 15 14 11 10 9 8 7 6 5 4 3 2 1

  82. n=47, m=77
  83. 59 58 57 56 55 54 53 49 48 46 45 44 43 42 41 39 37 36 35 33 32 30 29 28 27 26 25 24 23 21 20 19 18 17 16 15 14 12 11 10 8 7 6 5 3 2 1

  84. n=48, m=77
  85. 61 58 57 56 55 52 51 49 48 47 44 43 42 41 40 39 37 36 35 34 33 32 31 30 29 28 27 26 24 22 21 20 19 18 17 16 15 14 13 12 11 10 9 6 5 4 3 1

  86. n=49, m=78
  87. 61 60 57 56 55 54 51 50 49 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 24 23 22 21 20 19 18 17 16 15 14 13 11 7 6 5 3 2 1

  88. n=50, m=78
  89. 62 59 57 55 54 53 52 51 50 49 46 45 43 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 14 13 12 10 9 8 7 6 5 4 1

  90. n=51, m=82
  91. 67 65 59 58 57 56 54 49 48 47 46 45 44 43 42 41 40 39 36 35 34 32 31 30 29 28 27 26 25 24 23 22 21 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 2 1

  92. n=52, m=84
  93. 69 63 60 59 58 57 55 53 51 50 49 48 47 46 45 44 43 42 40 39 38 37 36 35 34 33 32 31 30 29 28 27 25 24 23 22 21 19 18 17 16 14 12 11 10 9 8 7 6 5 4 1

  94. n=53, m=86
  95. 69 63 62 61 60 57 56 55 54 53 52 51 50 49 48 46 45 44 43 40 39 38 37 36 35 34 32 31 30 29 28 27 26 24 23 22 21 19 18 16 15 14 13 11 10 9 8 6 5 4 3 2 1

  96. n=54, m=83
  97. 65 61 60 59 58 56 55 54 53 52 50 49 48 47 46 45 44 42 41 40 39 37 35 34 33 32 31 30 29 27 26 25 24 23 22 21 20 19 17 16 15 14 13 12 11 10 9 8 7 5 4 3 2 1

  98. n=55, m=87
  99. 67 65 62 61 60 59 58 56 55 54 53 52 51 50 49 48 47 46 45 44 42 41 40 39 38 37 35 34 33 32 31 30 28 27 26 25 24 23 22 21 20 18 15 14 13 12 11 10 7 6 5 4 3 2 1

  100. n=56, m=89
  101. 69 68 65 64 63 62 59 56 54 53 52 51 50 49 48 47 46 44 43 42 41 40 37 36 35 34 33 32 31 30 29 28 26 25 24 23 22 21 20 19 18 16 15 13 12 11 10 9 8 7 6 5 4 3 2 1

  102. n=57, m=90
  103. 69 66 65 64 63 62 61 59 58 56 54 53 52 51 50 49 47 45 44 43 42 41 40 39 38 37 36 34 33 32 31 29 28 26 25 24 23 22 21 19 18 17 16 15 14 13 12 11 10 9 8 7 5 4 3 2 1

  104. n=58, m=93
  105. 80 75 65 62 60 58 57 55 54 53 52 51 50 49 48 46 45 44 43 42 40 39 38 37 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  106. n=59, m=93
  107. 76 66 65 64 63 62 61 60 59 58 57 55 54 53 52 51 50 49 48 47 46 45 44 41 40 39 38 37 36 35 34 33 32 31 29 28 27 26 25 23 22 21 20 19 18 17 16 15 14 12 11 10 9 8 6 4 3 2 1

  108. n=60, m=95
  109. 75 71 67 66 65 64 62 61 60 59 58 57 56 55 54 53 51 50 49 46 45 44 43 42 40 38 37 36 35 33 32 31 30 29 28 27 26 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 5 4 3 2 1

  110. n=61, m=94
  111. 73 71 67 66 64 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 38 37 35 34 33 30 28 27 26 25 24 23 22 21 20 19 18 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  112. n=62, m=97
  113. 76 74 73 68 67 66 62 61 60 59 57 55 54 53 52 51 50 49 48 47 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 6 4 3 2 1

  114. n=63, m=98
  115. 74 73 69 68 67 66 65 64 63 62 61 59 58 57 56 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 34 33 32 31 30 29 28 27 25 24 23 21 20 19 18 17 15 13 12 11 9 8 7 6 5 4 3 2 1

  116. n=64, m=101
  117. 80 79 73 70 69 68 67 66 62 61 58 57 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3

  118. n=65, m=100
  119. 77 76 73 71 70 69 64 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 21 20 19 18 17 16 15 14 13 12 10 9 7 6 5 4 3 1

  120. n=66, m=104
  121. 82 74 73 72 71 70 69 67 66 65 64 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 46 44 42 41 40 39 38 37 36 34 33 32 31 30 28 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 1

  122. n=67, m=104
  123. 81 78 75 73 72 69 67 66 65 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 34 33 32 31 29 28 27 26 25 24 23 22 21 20 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 1

  124. n=68, m=106
  125. 82 78 77 76 75 72 70 69 66 65 64 62 60 59 57 56 55 54 53 52 51 50 49 48 47 46 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  126. n=69, m=108
  127. 82 81 79 78 75 72 71 70 68 67 64 63 61 60 59 58 57 56 55 54 53 52 51 50 49 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 15 12 11 10 9 8 7 6 5 4 3 2 1

  128. n=70, m=110
  129. 86 82 80 78 73 72 71 70 68 67 66 65 64 63 62 61 60 59 58 57 56 55 53 52 51 50 49 48 45 44 43 42 41 40 39 38 37 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 13 12 11 10 9 8 7 6 4 3 2 1

  130. n=71, m=109
  131. 88 82 76 74 73 72 71 70 67 66 65 64 63 62 61 60 59 58 57 56 54 53 52 51 50 49 48 47 46 44 43 42 41 40 39 38 37 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  132. n=72, m=113
  133. 82 81 80 79 77 76 75 74 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 54 53 52 51 50 48 47 46 45 43 42 41 40 39 38 37 36 35 34 33 31 30 29 28 27 26 25 24 22 21 20 18 17 16 15 14 13 12 11 10 9 8 7 6 4 3 2 1

  134. n=73, m=113
  135. 88 81 80 78 77 76 74 73 71 70 69 68 67 66 65 63 62 61 60 59 58 57 56 55 54 52 51 50 49 47 46 45 44 43 42 41 40 39 38 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 6 4 3 2 1

  136. n=74, m=116
  137. 86 85 82 81 79 77 76 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 43 42 41 40 39 38 37 34 33 31 30 29 28 27 26 25 24 21 20 19 18 17 16 15 14 13 12 11 10 8 7 6 5 4 3 2 1

  138. n=75, m=117
  139. 91 87 82 81 79 78 75 74 73 72 71 70 69 68 67 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 40 39 38 37 35 34 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  140. n=76, m=118
  141. 94 90 83 81 80 76 75 74 73 72 71 70 69 68 67 66 64 63 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 36 35 34 33 32 31 30 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  142. n=77, m=120
  143. 90 87 84 83 82 80 79 78 77 76 74 73 72 71 70 68 67 66 64 63 62 61 60 59 57 56 55 54 53 51 50 49 48 47 46 45 44 43 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  144. n=78, m=121
  145. 90 87 86 85 81 80 79 78 77 76 75 74 72 71 70 69 68 67 66 65 64 62 61 60 59 58 57 56 55 54 53 52 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 20 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3

  146. n=79, m=123
  147. 89 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 66 65 64 63 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 21 20 19 18 17 16 15 14 13 12 9 8 7 6 5 4 2

  148. n=80, m=125
  149. 94 88 87 86 85 84 83 81 80 79 78 77 76 75 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 36 35 34 33 32 31 30 29 28 27 26 25 23 22 21 20 18 17 16 14 13 12 11 10 9 8 7 6 5 4 2 1

  150. n=81, m=127
  151. 96 92 88 87 86 85 84 83 80 79 78 77 76 75 73 72 71 70 69 68 67 66 65 63 62 61 59 58 57 56 55 54 53 51 50 49 48 47 46 45 44 43 42 41 40 39 37 36 35 34 33 32 31 30 29 28 27 26 25 23 22 21 20 19 18 17 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  152. n=82, m=127
  153. 97 90 88 87 86 84 83 82 81 80 79 78 77 76 75 73 72 71 70 68 67 66 63 61 60 59 58 57 56 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  154. n=83, m=130
  155. 103 97 93 92 85 84 82 81 80 79 78 77 76 75 74 72 71 70 69 67 66 65 64 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 9 8 7 6 5 4 3 2 1

  156. n=84, m=131
  157. 96 94 92 90 88 87 86 85 84 83 82 81 80 78 76 75 74 72 71 70 69 68 67 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 31 30 29 28 27 26 25 24 23 22 20 18 17 16 15 14 13 12 11 9 8 7 6 5 4 3 2 1

  158. n=85, m=133
  159. 98 93 92 91 90 89 88 87 86 84 82 81 80 79 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 56 54 53 52 51 50 49 48 47 46 45 43 42 41 40 39 38 37 36 35 33 32 31 30 29 28 27 26 25 24 23 22 20 19 18 17 16 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  160. n=86, m=134
  161. 101 96 95 91 89 88 87 86 85 84 83 82 81 79 78 77 76 74 73 72 71 70 69 68 67 66 65 64 63 62 59 58 57 55 54 53 52 51 50 49 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 9 8 7 6 5 4 3 2 1

  162. n=87, m=135
  163. 102 97 94 93 92 89 88 86 85 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 44 43 42 41 40 39 38 37 36 35 34 32 31 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 11 10 9 8 7 6 5 4 3 2 1

  164. n=88, m=136
  165. 98 97 96 94 93 92 88 87 86 85 84 83 82 81 80 79 77 76 75 74 73 72 71 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 27 26 25 24 23 21 20 19 18 17 16 15 14 13 12 11 10 8 7 6 5 4 3 2 1

  166. n=89, m=140
  167. 105 98 96 95 94 93 92 90 89 87 86 85 84 83 82 81 80 79 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 57 56 55 54 53 52 51 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 33 32 31 30 29 28 27 26 25 24 23 21 20 19 18 17 16 15 14 13 12 11 10 8 7 6 4 3 2 1

  168. n=90, m=142
  169. 107 101 98 97 96 95 93 91 90 89 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 69 68 67 66 65 64 63 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 26 25 24 23 22 21 20 19 18 17 16 15 13 12 11 10 9 6 5 4 3 2 1

  170. n=91, m=142
  171. 104 100 99 98 96 94 93 92 90 89 88 86 85 84 83 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 57 56 55 54 53 52 51 50 49 48 46 45 44 42 41 40 39 38 37 36 35 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 9 8 7 6 5 4 3 2 1

  172. n=92, m=144
  173. 111 100 98 97 96 95 94 93 92 91 90 88 87 86 84 83 82 79 78 77 76 75 74 73 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 47 46 45 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 8 7 6 5 4 3 2 1

  174. n=93, m=146
  175. 108 103 102 100 97 96 95 94 93 92 90 89 88 87 86 85 84 83 81 80 79 76 75 74 73 72 71 70 69 68 67 66 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 11 10 9 8 7 6 5 4 3 2 1

  176. n=94, m=146
  177. 106 104 103 102 101 96 95 93 92 91 90 89 88 86 85 83 82 81 80 79 78 77 76 74 73 72 71 70 69 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 46 45 44 43 42 41 40 39 38 37 36 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  178. n=95, m=148
  179. 113 104 102 101 100 97 96 95 93 92 91 89 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 12 11 10 9 8 6 5 4 3 2 1

  180. n=96, m=150
  181. 114 106 102 100 99 98 97 96 95 94 92 91 90 89 88 87 86 85 84 83 82 81 80 79 77 76 75 73 72 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 5 4 3 2 1

  182. n=97, m=151
  183. 116 110 105 104 100 98 96 95 94 93 92 91 89 88 87 86 85 83 82 81 80 79 78 77 76 75 74 73 72 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 1

  184. n=98, m=153
  185. 107 106 105 104 103 102 100 98 97 96 95 94 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 60 59 58 57 56 55 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 8 7 6 5 4 2 1

  186. n=99, m=154
  187. 112 109 104 103 102 101 100 99 98 97 96 95 93 92 91 90 88 87 86 85 84 83 82 81 80 79 78 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 53 52 51 50 49 48 47 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 22 21 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  188. n=100, m=157
  189. 111 108 107 106 105 104 103 102 101 100 99 97 96 95 93 92 91 90 88 87 86 85 84 82 81 80 78 77 76 75 74 72 71 70 69 68 67 66 65 64 63 62 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 12 11 10 9 8 7 6 5 4 3 2 1

  190. n=101, m=159
  191. 114 109 108 106 105 104 103 102 101 100 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 65 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 43 42 40 39 38 37 36 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 17 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  192. n=102, m=160
  193. 114 112 111 108 105 104 103 102 101 100 99 98 97 96 95 94 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 73 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 50 49 48 47 46 44 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 11 10 9 8 7 6 5 4 3 2 1

  194. n=103, m=162
  195. 118 113 110 109 107 106 104 103 102 101 100 99 98 97 96 95 94 93 92 91 89 88 87 86 85 84 83 82 80 79 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 29 28 27 26 25 23 22 21 20 19 18 17 16 15 14 13 12 10 9 8 7 6 5 4 3 2 1

  196. n=104, m=164
  197. 120 114 112 110 109 108 107 105 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 43 41 40 39 38 37 35 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  198. n=105, m=164
  199. 121 113 112 111 109 107 106 105 104 103 102 101 100 99 97 95 93 92 91 90 89 88 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 26 25 24 23 22 21 20 19 17 16 15 14 13 12 11 10 9 8 7 5 4 3 2 1

  200. n=106, m=166
  201. 119 114 113 112 111 110 109 106 105 104 103 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 54 53 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 32 31 30 29 28 27 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 1

  202. n=107, m=168
  203. 119 116 115 114 111 110 109 107 106 105 104 103 102 101 99 98 97 96 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 31 30 29 28 27 26 25 24 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  204. n=108, m=171
  205. 122 117 116 115 114 113 112 111 110 108 106 105 104 101 100 99 98 97 96 95 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 77 76 75 74 73 72 71 70 69 68 67 66 65 64 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 14 13 12 11 9 8 7 6 5 4 3 2 1

  206. n=109, m=170
  207. 123 122 117 116 114 112 108 107 106 105 104 103 102 101 100 99 98 97 95 94 93 91 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 6 5 4 3 2 1

  208. n=110, m=174
  209. 124 119 118 116 115 114 113 111 110 109 108 107 105 104 103 102 101 100 99 98 97 95 94 93 92 90 89 88 87 86 85 84 83 82 81 80 79 78 77 75 74 73 72 71 70 69 68 67 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 19 18 17 16 15 13 12 11 10 9 8 7 6 5 4 3 2 1

  210. n=111, m=175
  211. 124 118 117 116 115 114 113 112 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 66 65 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 47 46 45 44 43 42 41 39 38 37 36 35 34 33 32 31 30 29 28 27 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 9 8 7 6 5 4 3 2 1

  212. n=112, m=176
  213. 129 119 118 117 116 115 114 112 111 110 109 108 107 106 105 103 102 101 99 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 66 65 64 63 62 61 60 59 58 57 56 55 54 53 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 28 27 26 25 24 23 22 21 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  214. n=113, m=177
  215. 126 121 119 118 117 116 115 113 112 110 109 108 107 106 105 104 103 102 101 100 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 46 45 44 43 42 41 39 38 37 36 35 33 32 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  216. n=114, m=180
  217. 125 124 122 121 120 119 116 115 113 112 111 110 109 108 107 105 104 103 102 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 70 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 19 18 17 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  218. n=115, m=182
  219. 133 130 123 121 120 119 118 117 115 113 112 111 110 108 107 105 104 102 101 100 99 98 97 96 95 93 92 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  220. n=116, m=183
  221. 133 127 124 123 120 119 118 116 115 114 113 111 110 109 108 107 105 104 103 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 82 81 80 79 78 77 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 14 13 12 11 10 9 8 7 6 5 4 3 2

  222. n=117, m=183
  223. 129 128 126 125 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 96 95 94 93 92 91 90 89 88 87 85 84 83 82 81 80 79 78 77 76 75 74 73 72 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 19 18 17 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  224. n=118, m=184
  225. 133 128 126 125 123 120 119 118 115 114 113 112 111 110 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  226. n=119, m=186
  227. 131 128 125 124 123 121 120 119 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 100 99 98 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 32 31 30 29 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  228. n=120, m=190
  229. 138 131 127 126 125 124 122 121 120 119 118 116 115 114 113 112 111 110 107 106 104 103 102 100 99 98 97 96 95 94 93 92 90 89 88 87 86 85 84 83 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  230. n=121, m=192
  231. 143 137 130 128 126 123 121 120 119 118 116 115 114 113 112 111 110 109 108 107 106 105 103 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 12 11 10 9 8 7 6 5 4 3 2 1

  232. n=122, m=194
  233. 140 136 132 131 130 127 123 122 121 120 119 116 115 114 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 54 53 51 50 49 48 47 46 45 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2

  234. n=123, m=194
  235. 141 138 131 130 126 124 123 121 120 119 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 63 62 61 60 59 58 57 56 55 53 52 51 50 49 48 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 4 3 2 1

  236. n=124, m=196
  237. 145 144 131 128 127 126 125 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 101 100 99 98 97 96 95 94 93 92 91 90 88 87 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  238. n=125, m=199
  239. 139 134 132 131 130 129 128 127 126 125 124 121 120 119 118 117 116 115 114 113 112 111 109 108 107 105 104 103 102 101 100 99 98 96 95 94 93 92 91 90 89 88 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  240. n=126, m=199
  241. 141 136 134 133 131 130 128 127 126 125 124 123 121 120 118 117 115 114 112 111 108 107 106 105 104 103 102 101 100 99 98 97 96 95 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  242. n=127, m=202
  243. 151 141 138 132 130 129 128 127 126 124 123 122 121 120 119 118 116 114 113 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 21 20 19 18 17 16 15 13 12 11 10 9 8 7 6 5 4 3 2 1

  244. n=128, m=203
  245. 146 139 136 133 131 130 129 128 127 126 124 123 122 121 119 118 117 116 115 114 113 112 111 110 109 108 107 106 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 63 62 61 60 59 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 5 4 3 2 1

  246. n=129, m=205
  247. 146 143 137 135 134 133 132 129 128 127 126 124 123 122 121 119 118 117 116 115 114 113 112 111 110 109 108 106 105 104 103 102 101 100 97 96 95 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 12 11 10 9 8 7 6 5 4 3 2 1

  248. n=130, m=207
  249. 141 140 138 137 136 135 134 133 132 130 128 126 125 124 123 121 120 119 118 117 116 114 113 112 111 110 108 107 106 105 104 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 79 78 77 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  250. n=131, m=208
  251. 141 140 138 137 136 135 134 133 131 128 127 126 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 105 104 103 102 101 100 99 98 97 96 95 92 91 90 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  252. n=132, m=210
  253. 148 145 141 139 137 136 135 132 131 130 129 126 125 124 123 122 121 120 119 118 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 88 87 86 85 84 83 82 81 80 79 78 76 75 74 73 72 71 70 69 68 67 66 65 64 63 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 2 1

  254. n=133, m=212
  255. 150 148 143 142 141 140 138 134 133 129 128 127 126 125 123 122 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 17 16 15 14 13 12 11 10 8 7 6 5 4 3 2 1

  256. n=134, m=214
  257. 155 142 141 139 138 137 136 135 134 133 131 130 128 127 126 125 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 86 85 84 83 82 81 80 79 78 77 76 75 74 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 16 14 13 12 10 9 8 7 6 5 4 3 2 1

  258. n=135, m=215
  259. 148 147 146 145 143 139 138 135 134 131 130 129 128 127 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 14 13 12 11 10 9 8 7 6 5 4 2 1

  260. n=136, m=217
  261. 156 152 148 144 143 137 136 135 133 131 130 129 128 127 126 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 71 70 69 67 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  262. n=137, m=218
  263. 150 147 144 143 142 140 139 138 137 136 135 134 133 131 130 126 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 2 1

  264. n=138, m=220
  265. 154 152 150 146 145 142 140 139 136 135 134 133 131 130 128 127 126 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 106 105 104 103 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  266. n=139, m=222
  267. 155 153 147 143 142 141 140 139 138 137 136 135 134 133 132 131 130 129 128 126 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  268. n=140, m=224
  269. 154 151 148 147 146 145 144 142 139 137 136 135 134 133 132 131 130 129 128 127 126 125 124 123 121 120 119 118 117 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 38 37 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  270. n=141, m=226
  271. 153 152 150 148 147 146 143 142 140 139 138 137 135 134 133 132 131 130 129 128 127 126 125 124 123 122 121 120 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 7 5 4 3 2 1

  272. n=142, m=227
  273. 159 155 153 150 148 146 144 142 141 139 138 137 136 135 134 133 131 129 128 127 126 125 124 123 122 121 120 119 118 116 115 114 113 112 111 110 109 108 107 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  274. n=143, m=228
  275. 154 152 149 148 147 146 145 144 142 141 140 138 137 136 135 133 132 131 130 129 128 127 126 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 48 47 46 45 44 43 42 41 39 38 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  276. n=144, m=231
  277. 153 152 151 150 149 148 147 146 145 144 143 141 139 138 137 136 135 134 133 131 130 129 128 127 126 125 124 123 122 121 120 118 117 116 115 114 113 112 111 110 109 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 18 17 16 15 14 13 12 11 10 9 7 6 5 3 2 1

  278. n=145, m=231
  279. 155 153 152 151 150 147 146 145 144 143 141 140 138 137 136 135 134 133 132 131 130 129 128 127 126 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 79 78 77 76 75 74 73 72 71 70 69 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  280. n=146, m=233
  281. 168 155 153 151 150 148 147 146 145 142 141 140 139 138 137 136 135 134 133 132 131 130 129 128 126 125 124 123 122 120 119 118 117 116 115 114 113 112 111 110 109 108 107 105 104 103 102 101 100 99 98 97 96 94 93 92 91 90 89 88 87 86 85 84 83 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  282. n=147, m=236
  283. 164 159 156 155 153 152 150 148 146 145 144 143 142 139 138 137 136 135 133 132 131 130 129 128 127 126 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 7 6 5 4 3 2 1

  284. n=148, m=238
  285. 168 166 158 155 153 152 151 149 146 145 144 143 142 141 140 139 138 136 135 133 131 130 129 128 127 126 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 24 23 22 21 20 19 18 17 16 14 13 12 11 10 9 8 7 6 5 3 2 1

  286. n=149, m=240
  287. 163 158 157 156 155 152 151 150 149 147 146 145 144 142 141 140 139 138 137 136 135 134 133 132 131 130 129 128 127 126 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 66 65 63 62 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  288. n=150, m=241
  289. 173 168 167 157 153 152 151 149 147 146 145 144 143 142 141 140 139 137 136 135 134 133 131 130 129 128 127 126 124 123 122 121 120 119 118 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 78 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1

  290. n=151, m=243
  291. 169 161 160 159 158 155 154 153 152 150 149 148 147 145 143 142 140 138 137 136 135 134 133 132 131 130 129 128 126 125 124 123 122 121 120 119 118 117 116 115 114 113 112 111 110 109 108 107 106 105 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 87 86 85 84 83 82 81 80 79 77 76 75 74 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 44 43 42 41 40 39 38 37 36 35 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 6 5 4 3 2 1
复制代码

点评

mathe版主7次幂使用的哪个模加速的?  发表于 2026-5-3 07:37
我试了一下,7次幂模29也不是很好使  发表于 2026-5-2 22:24

评分

参与人数 1威望 +8 金币 +8 贡献 +8 经验 +8 鲜花 +8 收起 理由
northwolves + 8 + 8 + 8 + 8 + 8 神马都是浮云

查看全部评分

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2026-5-3 08:13:27 | 显示全部楼层
这个计算难度更大一些:

A395519

a(n) is the smallest number k such that k^n is the sum of 2n distinct positive n-th powers.

DATA
3, 9, 13, 31, 32, 46, 47, 113, 47, 140

OFFSET
1,1

EXAMPLE
3^1 = 1^1 + 2^1;

9^2 = 2^2 + 4^2 + 5^2 + 6^2;

13^3 = 1^3 + 5^3 + 6^3 + 7^3 + 8^3 + 10^3;

31^4 = 8^4 + 10^4 + 14^4 + 16^4 + 17^4 + 18^4 + 20^4 + 26^4;

32^5 = 2^5 + 3^5 + 7^5 + 8^5 + 9^5 + 11^5 + 15^5 + 17^5 + 19^5 + 31^5;

46^6 = 1^6 + 3^6 + 5^6 + 7^6 + 10^6 + 13^6 + 14^6 + 19^6 + 21^6 + 28^6 + 40^6 + 41^6;

47^7 = 2^7 + 4^7 + 7^7 + 8^7 + 9^7 + 11^7 + 17^7 + 18^7 + 20^7 + 21^7 + 22^7 + 35^7 + 40^7 + 43^7;

113^8 = 108^8 + 86^8 + 84^8 + 79^8 + 74^8 + 56^8 + 52^8 + 48^8 + 38^8 + 36^8 + 34^8 + 28^8 + 26^8 + 12^8 + 8^8 + 6^8;

47^9 = 1^9 + 2^9 + 4^9 + 7^9 + 11^9 + 14^9 + 15^9 + 18^9 + 26^9 + 27^9 + 30^9 + 31^9 + 32^9 + 33^9 + 36^9 + 38^9 + 39^9 + 43^9;

140^10 = 132^10 + 121^10 + 110^10 + 103^10 + 102^10 + 96^10 + 86^10 + 77^10 + 66^10 + 55^10 + 54^10 + 52^10 + 44^10 + 38^10 + 34^10 + 33^10 + 24^10 + 19^10 + 13^10 + 7^10;

目前其它幂次计算进度:

{11} >= 73
{12} >= 134
{13} >= 73
{14} >= 75
{15} >= 75
{16} >= 244
{17} >= 79
{18} >= 312
{19} >= 81
{20} >= 81
{21} >= 83
{22} >= 86
{23} >= 87
{24} >= 88
{25} >= 89
{26} >= 90
{27} >= 92
{28} >= 93
{29} >= 94
{30} >= 95
{31} >= 95
{32} >= 96
{33} >= 96
{34} >= 96
{35} >= 96
{36} >=201
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2026-5-14 01:04 , Processed in 0.035046 second(s), 17 queries .

Powered by Discuz! X3.5

© 2001-2026 Discuz! Team.

快速回复 返回顶部 返回列表