题目61：找出唯一具有循环性质的6个四位数定形数集合的元素和

欧拉计划

Cyclical figurate numbers

Triangle, square, pentagonal, hexagonal, heptagonal, and octagonal numbers are all figurate (polygonal) numbers and are generated by the following formulae:


   
The ordered set of three 4-digit numbers: 8128, 2882, 8281, has three interesting properties.

1. The set is cyclic, in that the last two digits of each number is the first two digits of the next number (including the last number with the first).
2. Each polygonal type: triangle , square , and pentagonal , is represented by a different number in the set.
3. This is the only set of 4-digit numbers with this property.

Find the sum of the only ordered set of six cyclic 4-digit numbers for which each polygonal type: triangle, square, pentagonal, hexagonal, heptagonal, and octagonal, is represented by a different number in the set.

题目：

三角形数，四角形数，五角形数，六角形数，七角形数和八角形数都是定形数，他们分别由以下公式产生：



三个四位数形成的有序集合: 8128, 2882, 8281,有三个有趣的性质：

1. 这个集合是循环的：每个数的后两位数是下一个数的前两位数，包括第三个和第一个的关系。
2. 三种定形数中的每一种都能被这三个数中的一个不同的数代表：三角形数 , 四角形数 , 和五角形数 。
3. 这是唯一具有以上性质的四位数的集合。

找出唯一的一个六个四位数的循环集合，使得从三角形数到八角形数中的每一种都能由该集合中的一个不同的数代表。

求这个集合中元素之和。

迷雾少年

Mark

jerryxjr1220

1. ([8128, 2882, 8256, 5625, 2512, 1281], 28684)
   
2. Time:0.016s
   
3. 
   
4. import sys
   
5. import time
   
6. st = time.time()
   
7. def calc_numbers(start, end, proc):
   
8.     n = 1
   
9.     out = []
   
10.     while True:
    
11.         z = proc(n)
    
12.         if z >= start and z < end: out.append(z)
    
13.         if z >= end: return out
    
14.         n += 1
    
15. 
    
16. p3 = calc_numbers(1000, 10000, lambda n : n * (n + 1) / 2)
    
17. p4 = calc_numbers(1000, 10000, lambda n : n * n)
    
18. p5 = calc_numbers(1000, 10000, lambda n : n * (3 * n - 1) / 2)
    
19. p6 = calc_numbers(1000, 10000, lambda n : n * (2 * n - 1))
    
20. p7 = calc_numbers(1000, 10000, lambda n : n * (5 * n - 3) / 2)
    
21. p8 = calc_numbers(1000, 10000, lambda n : n * (3 * n - 2))
    
22. 
    
23. 
    
24. def find_loop(lists, start, end, found_numbers = []):
    
25.     if len(lists) == 1 and start * 100 + end in lists[0]:
    
26.         found_numbers.append(start * 100 + end)
    
27.         found_numbers.append(n)
    
28.         print (found_numbers, sum(found_numbers))
    
29.         ed = time.time()
    
30.         print ('Time:%.3fs'%(ed-st))
    
31.         sys.exit()
    
32.    
    
33.     for cur_list in lists:
    
34.         for c in cur_list:
    
35.             if c / 100 == start:
    
36.                 lists_copy = lists[0:]
    
37.                 lists_copy.remove(cur_list)
    
38.                 find_loop(lists_copy, c % 100, end, found_numbers + [c])
    
39. 
    
40. for n in p8:
    
41.     find_loop([p3, p4, p5, p6, p7], n % 100, n / 100)
    

复制代码

飘飞的白杨

1. import sys
   
2. from time import time
   
3. 
   
4. ttt = time()
   
5. 
   
6. def getlist(lam):
   
7.     '''获取x角形数的列表'''
   
8.     l = []
   
9.     for i in range(1, 1000):
   
10.         n = lam(i)
    
11.         if n < 1000: continue
    
12.         if n >= 10000: return l
    
13.         l.append(n)
    
14. 
    
15. def progrem61_2(numbers):
    
16.     if len(numbers) == 6 and numbers[-1] % 100 == numbers[0] // 100:
    
17.         print(numbers, sum(numbers), "\ntime:%.5fs"%(time()-ttt))
    
18.         sys.exit()
    
19.     for i in range(6):
    
20.         if flag[i]: continue
    
21.         flag[i] = True
    
22.         for j in p[i]:
    
23.             if j // 100 == numbers[-1] % 100:
    
24.                 progrem61_2(numbers+[j])
    
25.         flag[i] = False
    
26. 
    
27. def progrem61_1():
    
28.     flag[-1] = True
    
29.     # 八角形数最少，从八角形数中选择第一个数
    
30.     for i in p8:
    
31.         progrem61_2([i])
    
32. 
    
33. p3 = getlist(lambda n: n * (n + 1) // 2)
    
34. p4 = getlist(lambda n: n * n)
    
35. p5 = getlist(lambda n: n * (3 * n - 1) // 2)
    
36. p6 = getlist(lambda n: n * (2 * n - 1))
    
37. p7 = getlist(lambda n: n * (5 * n - 3) // 2)
    
38. p8 = getlist(lambda n: n * (3 * n - 2))
    
39. p = [p3, p4, p5, p6, p7, p8]
    
40. flag = [False] * 6 # 标志数组，用于判断是否遍历过
    
41. progrem61_1()

复制代码

guosl

1. /*
   
2. 答案：28684
   
3. 耗时：0.0000797秒
   
4. */
   
5. #include <iostream>
   
6. #include <vector>
   
7. #include <algorithm>
   
8. #include <omp.h>
   
9. using namespace std;
   
10. 
    
11. vector<int> v[6];//v[0],...,v[5]分别是三角数，...，八角数
    
12. int x[6];//记录查找的结果
    
13. int p[6] = { 0,1,2,3,4,5 };//改变表的查找次序用
    
14. 
    
15. void GenTriangle(void)//生成三角数表
    
16. {
    
17.   int k = 1;
    
18.   while (true)
    
19.   {
    
20.     int n = k * (k + 1) / 2;
    
21.     if (n > 9999)
    
22.       break;
    
23.     else if (n >= 1000)
    
24.       v[0].push_back(n);
    
25.     ++k;
    
26.   }
    
27. }
    
28. 
    
29. void GenSquare(void)//生成四角形数表
    
30. {
    
31.   int k = 1;
    
32.   while (true)
    
33.   {
    
34.     int n = k * k;
    
35.     if (n > 9999)
    
36.       break;
    
37.     else if (n >= 1000)
    
38.       v[1].push_back(n);
    
39.     ++k;
    
40.   }
    
41. }
    
42. 
    
43. void GenPentagonal(void)//生成五角形数表
    
44. {
    
45.   int k = 1;
    
46.   while (true)
    
47.   {
    
48.     int n = k * (3 * k - 1) / 2;
    
49.     if (n > 9999)
    
50.       break;
    
51.     else if (n >= 1000)
    
52.       v[2].push_back(n);
    
53.     ++k;
    
54.   }
    
55. }
    
56. 
    
57. void GenHexagonal(void)//生成六角形数表
    
58. {
    
59.   int k = 1;
    
60.   while (true)
    
61.   {
    
62.     int n = k * (2 * k - 1);
    
63.     if (n > 9999)
    
64.       break;
    
65.     else if (n >= 1000)
    
66.       v[3].push_back(n);
    
67.     ++k;
    
68.   }
    
69. }
    
70. 
    
71. void GenHeptagonal(void)//生成七角形数表
    
72. {
    
73.   int k = 1;
    
74.   while (true)
    
75.   {
    
76.     int n = k * (5 * k - 3) / 2;
    
77.     if (n > 9999)
    
78.       break;
    
79.     else if (n >= 1000)
    
80.       v[4].push_back(n);
    
81.     ++k;
    
82.   }
    
83. }
    
84. 
    
85. void GenOctagonal(void)//生成八角形数表
    
86. {
    
87.   int k = 1;
    
88.   while (true)
    
89.   {
    
90.     int n = k * (3 * k - 2);
    
91.     if (n > 9999)
    
92.       break;
    
93.     else if (n >= 1000)
    
94.       v[5].push_back(n);
    
95.     ++k;
    
96.   }
    
97. }
    
98. 
    
99. bool dfs(int k, int f, int b)//深度搜索找匹配链
    
100. {
     
101.   if (k == 6)
     
102.     return (f == b);//检查是否可以和第一个数匹配
     
103.   int a1 = b * 100, a2 = a1 + 99;
     
104.   //找出下一个表中可能匹配的范围
     
105.   auto itr1 = lower_bound(v[p[k]].begin(), v[p[k]].end(), a1);
     
106.   auto itr2 = upper_bound(v[p[k]].begin(), v[p[k]].end(), a2);
     
107.   for (auto itr = itr1; itr != itr2; ++itr)//具体查找
     
108.   {
     
109.     //找到可能的匹配
     
110.     x[k] = *itr;//记录查找的结果
     
111.     if (dfs(k + 1, f, x[k] % 100))//往下找
     
112.       return true;
     
113.   }
     
114.   return false;
     
115. }
     
116. 
     
117. int main(void)
     
118. {
     
119.   double t = omp_get_wtime();
     
120.   //生成各种数的表
     
121.   GenTriangle();
     
122.   GenSquare();
     
123.   GenPentagonal();
     
124.   GenHexagonal();
     
125.   GenHeptagonal();
     
126.   GenOctagonal();
     
127.   do
     
128.   {
     
129.     bool bFind = false;
     
130.     for (int i = 0; i < (int)v[p[0]].size(); ++i)
     
131.     {
     
132.       int a0 = v[p[0]][i];
     
133.       int b = a0 % 100;//后二位
     
134.       int f = a0 / 100;//前二位
     
135.       x[0] = v[p[0]][i];//记录查找的结果
     
136.       bFind = dfs(1, f, b);
     
137.       if (bFind)
     
138.         break;
     
139.     }
     
140.     if (bFind)
     
141.       break;
     
142.   } while (next_permutation(p, p + 6));//改变一下表的查找次序
     
143.   int nSum = 0;
     
144.   for (int i = 0; i < 6; ++i)
     
145.     nSum += x[i];
     
146.   t = omp_get_wtime() - t;
     
147.   cout << nSum << endl << t << endl;
     
148.   return 0;
     
149. }

复制代码

B1tetheDust

1. import time as t
   
2. 
   
3. start = t.perf_counter()
   
4. 
   
5. 
   
6. def divide_4_digit_num(a_num_list):
   
7.     new_list = []
   
8.     for each_num in a_num_list:
   
9.         new_list.append([str(each_num)[:2], str(each_num)[2:]])
   
10.     return new_list
    
11. 
    
12. 
    
13. def belong_to(a_num):
    
14.     range_list = [95, 163, 219, 267, 310, 350]
    
15.     for belong_i in range(6):
    
16.         if a_num <= range_list[belong_i]:
    
17.             return belong_i
    
18. 
    
19. 
    
20. def find_cyclic_num():
    
21.     Tri = [int(n * (n + 1) / 2) for n in range(45, 141)]  # length = 96
    
22.     Squ = [int(n ** 2) for n in range(32, 100)]  # length = 68
    
23.     Penta = [int(n * (3 * n - 1) / 2) for n in range(26, 82)]  # length = 56
    
24.     Hexa = [int(n * (2 * n - 1)) for n in range(23, 71)]  # length = 48
    
25.     Hept = [int(n * (5 * n - 3) / 2) for n in range(21, 64)]  # length = 43
    
26.     Octa = [int(n * (3 * n - 2)) for n in range(19, 59)]  # length = 40
    
27.     Tri_list, Squ_list = divide_4_digit_num(Tri), divide_4_digit_num(Squ)
    
28.     Penta_list, Hexa_list = divide_4_digit_num(Penta), divide_4_digit_num(Hexa)
    
29.     Hept_list, Octa_list = divide_4_digit_num(Hept), divide_4_digit_num(Octa)
    
30.     dict_list = {"Tri_list": Tri_list, "Squ_list": Squ_list, "Penta_list": Penta_list, "Hexa_list": Hexa_list,
    
31.                  "Hept_list": Hept_list, "Octa_list": Octa_list}
    
32.     total = []
    
33.     for lists in dict_list:
    
34.         total += dict_list[lists]
    
35.     len_total = len(total)
    
36.     for i00 in range(len_total):
    
37.         index_list = [belong_to(i00)]
    
38.         for i10 in range(len_total):
    
39.             if belong_to(i10) not in index_list and total[i00][1] == total[i10][0]:
    
40.                 index_list = [belong_to(i00), belong_to(i10)]
    
41.                 for i20 in range(len_total):
    
42.                     if belong_to(i20) not in index_list and total[i10][1] == total[i20][0]:
    
43.                         index_list = [belong_to(i00), belong_to(i10), belong_to(i20)]
    
44.                         for i30 in range(len_total):
    
45.                             if belong_to(i30) not in index_list and total[i20][1] == total[i30][0]:
    
46.                                 index_list = [belong_to(i00), belong_to(i10), belong_to(i20), belong_to(i30)]
    
47.                                 for i40 in range(len_total):
    
48.                                     if belong_to(i40) not in index_list and total[i30][1] == total[i40][0]:
    
49.                                         cur_num_list = [total[i00], total[i10], total[i20], total[i30], total[i40]]
    
50.                                         index_list = [belong_to(i00), belong_to(i10), belong_to(i20), belong_to(i30),
    
51.                                                       belong_to(i40)]
    
52.                                         for i50 in range(len_total):
    
53.                                             if belong_to(i50) not in index_list and \
    
54.                                                     total[i40][1] == total[i50][0] and \
    
55.                                                     total[i50][1] == cur_num_list[0][0]:
    
56.                                                 cur_num_list.append(total[i50])
    
57.                                                 return cur_num_list
    
58. 
    
59. 
    
60. res = find_cyclic_num()
    
61. print(res)
    
62. print(sum([int(i[0] + i[1]) for i in res]))
    
63. print("It costs %f s" % (t.perf_counter() - start))
    

复制代码

[['25', '12'], ['12', '81'], ['81', '28'], ['28', '82'], ['82', '56'], ['56', '25']]
28684
It costs 0.530015 s