猜猜我的手机号码？

winxos

直接试除还是挺快的: #include using namespace std; using namespace LiDIA; #define HPRIME 1000 int primelist[HPRIME*HPRIME]; int pc; void initp() { int i,j; primelist[0]=primelist[1]=1; ... mathe 发表于 2009-8-26 11:27

mathe的using namespace LiDIA; 程序中使用了 LiDIA的什么？怎么没看出来

gxqcn

可以把那行去掉。 因为 mathe 曾用 LiDIA 进行因数分解，改变算法后忘了去掉它了。 另外，我用 VS2008 编译 47# 的代码，需要将第 36 行代码由

1. printf("%lld=%lld*%lld\n",r,primelist[i],r/primelist[i]);

复制代码

改成

1. printf("%lld=%d*%d\n",r,primelist[i],r/primelist[i]);

复制代码

才能正确输出。

winxos

52# gxqcn 原来如此。

〇〇

47# mathe 我用你的算法改写的vc版本，怎么出来252个结果？

2. #include <stdio.h>
3. #include <conio.h>
4. //#include <windows.h>
5. #include "2M.txt"
6. void output(int *p); //输出函数
7. void permutation(int *p);//全排列函数
8. void test(__int64 r);
10. int a[]={0,1,2,5,6,7};//待排列数组
11. int N=sizeof(a)/sizeof(int); //数组长度
12. int *p=a; //数组指针
13. int n=0;
14. int m[720];
15. void test(__int64 r)
17. {
19. int i;
21. for(i=0;i<141800;i++){///skip prime 2 and 3
23. if(prime[i]>=1900000)return;
25. if(r%prime[i]==0)break;
27. }
29. if(prime[i]<10000)return;
31. printf("%I64d=%d*%d\n",r,prime[i],r/prime[i]);
33. }
36. void main()
37. {
38. output(p);
39. permutation(p); //调用全排列函数
40. __int64 r=189e8;
41. for(int i=0;i<720;i++)
42. test(r+m[i]*100+43);
44. //output(p); //输出排列结果
45. //getch();
46. }
48. void permutation(int *p)
49. {
50. bool flag=false; //标志排列是否结束
51. int temp; //临时变量,用来交换数组元素
53. while(1)
54. {
55. for(int i=0;i<N-1;i++)
56. {
57. flag=true; //如已排列完,falg为真
59. if(p[i]<p[i+1]) //如果发现第一个倒序的数组元素,p[i+1]被称为倒序元素
60. {
61. flag=false; //说明未排列完
63. int j=0;
64. while(p[j++]>=p[i+1]); //找到第一个大于或等于倒序元素的数组元素
65. j-=1;
67. temp=p[i+1]; //交换这两个元素
68. p[i+1]=p[j];
69. p[j]=temp;
71. j=0;
72. while(j<(i+1)/2) //倒置倒充元素原位置以前的所有元素
73. {
74. temp=p[j];
75. p[j]=p[i-j];
76. p[i-j]=temp;
78. j++;
79. }
80. break; //跳出for循环,
81. }
82. }
83. if(flag) //如果排列完,结束循环
84. break;
85. output(p); //输出排列结果
86. }
87. }
89. void output(int *p)
90. {
91. int s=0;
92. for(int i=0;i<N;i++)
93. s=s*10+p[i];
94. m[n++]=s;
95. //printf("%06d",s);
96. //printf("\n");
98. }

复制代码

mathe

你的素数列表有问题.我的要求的是2开始的素数列表

〇〇

55# mathe 为什么不是10000开始的质数呢？不是说其中一个为5位数？

gxqcn

因为若用试除法， 得先排除有小于10000的小因子的情形， 如果你的素数表缩水了，将无法进行该过程。

〇〇

57# gxqcn 你的意思是说 我的结果中有的除数是质数，商是合数？

gxqcn

right!

〇〇

我明白了，如果一个电话号码没有比10000小的质因数，而又有一个比10000大的质因数，那么除以其得到的商一定也是质数 这里用到了电话号码<19e10的条件