Ten digit numbers

无心人

不行的 haskell比他们都快 也简单 但都不如C/C++

无心人

2. #include <stdio.h>
3. #include <stdlib.h>
4. #include <gmp.h>
6. int test(mpz_t pow, int p)
7. {
8. unsigned tmp;
9. int d[10], i;
10. char t[250];
11. for (i = 0; i < 10; i ++) d[i] = 0;
13. mpz_get_str(t, 10, pow);
14. for (i = 0; i < strlen(t); i ++)
15. d[t[i] - '0']++;
16. for (i = 0; i < 10; i ++)
17. if (d[i] != p)
18. return 0;
20. return 1;
21. }
22. int main(void)
23. {
24. int p, b, i = 0;
25. mpz_t n, pow;
26. char start[16];
27. printf("输入方幂: ");
28. scanf("%d", &p);
29. printf("输入起始数字(最大16位): ");
30. scanf("%s", start);
31. mpz_init(pow);
32. mpz_init(n);
33. mpz_init(m);
34. mpz_init(t6);
35. mpz_set_ui(t6, 1000000);
36. mpz_set_str(n, start, 10);
37. b = 0;
38. printf("\n");
39. while (b == 0)
40. {
41. i ++;
42. mpz_pow_ui(pow, n, p);
43. if (test(pow, p) == 1)
44. {
45. gmp_printf("找到[%Zd, %Zd]\n", n, pow);
46. b = 1;
47. break;
48. }
49. mpz_add_ui(n, n, 1);
50. if (i >= 10000000)
51. {
52. i = 0;
53. gmp_printf("搜索到%Zd\n", n);
54. }
55. }
56. mpz_clear(n);
57. mpz_clear(pow);
58. mpz_clear(m);
59. mpz_clear(t6);
60. return 0;
61. }

复制代码

修改好了 不用预计算了 但这次时间达到了500秒

无心人

19的已经搜索到97 估计也是就1个

无心人

从8858667905到9873773268 只有一个 已经搜索完

mathe

19的确只有一个

无心人

2. #include <stdio.h>
3. #include <stdlib.h>
4. #include <gmp.h>
6. int test(mpz_t pow, int p)
7. {
8. unsigned tmp;
9. int d[10], i;
10. char t[250];
11. for (i = 0; i < 10; i ++) d[i] = 0;
13. mpz_get_str(t, 10, pow);
14. for (i = 0; i < strlen(t); i ++)
15. d[t[i] - '0']++;
16. for (i = 0; i < 10; i ++)
17. if (d[i] != p)
18. return 0;
20. return 1;
21. }
22. int main(void)
23. {
24. int p, b, i = 0;
25. mpz_t n, pow, t10, tp;
26. printf("输入方幂: ");
27. scanf("%d", &p);
28. mpz_init(pow);
29. mpz_init(n);
30. mpz_init(t10);
31. mpz_init(tp);
32. mpz_set_ui(t10, 1000000000);
33. mpz_mul_ui(t10, t10, 10);
34. mpz_set_ui(tp, 10);
35. mpz_pow_ui(tp, tp, 10*p - 1);
36. mpz_root(n, tp, p);
37. mpz_add_ui(n, n, 1);
38. b = 0;
39. gmp_printf("初始数字：%Zd\n", n);
40. while (mpz_cmp(n, t10) < 0)
41. {
42. mpz_pow_ui(pow, n, p);
43. if (test(pow, p) == 1)
44. gmp_printf("找到[%Zd, %Zd]\n", n, pow);
46. mpz_add_ui(n, n, 1);
47. }
48. mpz_clear(n);
49. mpz_clear(pow);
50. mpz_clear(t10);
51. mpz_clear(tp);
52. return 0;
53. }

复制代码

改成自动的了 只要输入20就能得到所有20次方的恰好有20个0123456789结果 附上程序，没有限制幂的大小 5300多秒扫描一个方幂 pn.rar (53.79 KB, 下载次数: 4)

2009-1-12 16:23 上传

点击文件名下载附件

无心人

输入方幂: 20 初始数字：8912509382 找到[8951993472, 109243381929665290144735387253484896485512153800629962749796403476740675298520463741968915306665585481026029884 06317917558817958274385371306360838853525309513410701102261427360428701722029347911499776] 20的找到一个最小的 经过Haskell核验， 是正确的

无心人

同时启动了20， 21 在考虑是否再启动几个

无心人

服务器上启动了22， 23

无心人

除非停电 明天给出20-23的结果 别人就不要测试20-23了 要测试的 下载我给出的可执行程序 算24以后的 明天写带日志文件的 一口气挂几个连续的幂出来