Ten digit numbers

无心人

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

无心人

1. 
   
2. #include <stdio.h>
   
3. #include <stdlib.h>
   
4. #include <gmp.h>
   
5.    
   
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;
    
12. 
    
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;
    
19.       
    
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的确只有一个

无心人

1. 
   
2. #include <stdio.h>
   
3. #include <stdlib.h>
   
4. #include <gmp.h>
   
5.    
   
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;
    
12. 
    
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;
    
19.       
    
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);
    
45. 
    
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以后的

明天写带日志文件的
一口气挂几个连续的幂出来