连续自然数构成的素数

无心人

定时以一个小时保存一次状态 和状态的MD5编码

mathe

原帖由 无心人 于 2008-4-23 17:42 发表 你加上定时存盘和冗余纠错 我给你挂服务器上三个月

逃跑

无心人

这么大的工作量 无法保证在运行时机器断电 不做我说的那两点 最大运行时间只能是24小时 做不到7天

无心人

如果降序排列 在82产生素数 82818079...987654321是素数

无心人

http://mathworld.wolfram.com/ConsecutiveNumberSequences.html 似乎目前没找到素数

数论爱好者

无心人上面给出的外国网址，我去看了一下，他们给出的通项式 1/81*(10^(n+1)-9n-10) 当n在1至9时，符合我的连续自然数形式，当n=10以后，不满足连续自然数。 如n=10时，计算值1234567900，不是我要的12345678910 n=11时，计算值12345679011，不是我要的1234567891011 n=17，12345679012345677 请无心人验证我的计算是否正确 说明这个通项式不是连续自然数的通项式

gxqcn

文中说过的，该公式仅适用于“up to n=9”

mathe

而对于n在10和99之间，100和999之间等等，都有类似的公式，只是一些系数不同。不过这个公式好像用处不大

gxqcn

1. #include <stdlib.h>
2. #include <stdio.h>
3. #include <assert.h>
5. #include "../../../HugeCalc_API/CppAPI/Include/HugeCalc.h" // 公共接口
6. #include "../../../HugeCalc_API/CppAPI/Include/HugeInt.h" // 10进制系统
8. #pragma message( "automatic link to ../../../HugeCalc_API/CppAPI/Lib/HugeCalc.lib" )
9. #pragma comment( lib, "../../../HugeCalc_API/CppAPI/Lib/HugeCalc.lib" )
12. #define NUMB 82
13. #define DIGITS_OF_NUMB 2 // digits of NUMB
14. #define POW10_DIGITS 100 // 10^DIGITS_OF_NUMB
16. int main( void )
17. {
18. const UINT32 len = ( NUMB + 1 )*DIGITS_OF_NUMB - ( POW10_DIGITS - 1 ) / 9;
20. char *pBuf, *p;
21. UINT32 i, d=0;
23. p = pBuf = new char[len + 1];
24. assert( NULL != pBuf );
26. memset( pBuf, 0, ( len + 1 )*sizeof( char ));
27. #if 0
28. for ( i = 1; i <= NUMB; ++i )
29. #else
30. for ( i = NUMB; i > 0; --i )
31. #endif
32. {
33. d = sprintf( p+=d, "%u", i );
34. }
36. assert( strlen(pBuf) == len );
38. HugeCalc::EnableTimer();
40. printf( CHugeInt( pBuf ).IsPrime() \
41. ? "Is a prime number!\r\n" \
42. : "Is not a prime number.\r\n" );
44. HugeCalc::EnableTimer( FALSE );
45. printf( "computation took %s\r\n", HugeCalc::GetTimerStr( FT_HHMMSS_ms ));
47. delete []pBuf;
49. return 0;
50. }

复制代码

无心人

玩素性检测看来不行了 可以考虑玩分解 不过分解代价很高