数学研发网设为首页收藏本站
切换到窄版

 找回密码
 欢迎注册
数学研发论坛»论坛 › 【算法研发】 › 编程擂台 › 立方数数字和问题
返回列表 发新帖
查看: 10441|回复: 5

[讨论] 立方数数字和问题

[复制链接]
无心人
发表于 2011-9-12 19:41:50 | 显示全部楼层 |阅读模式

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要 登录 才可以下载或查看，没有账号？欢迎注册

×
沿袭平方数数字和问题的描述 http://bbs.emath.ac.cn/thread-642-1-1.html 考虑平方数 先给出递归公式 T(0) = 0 N(0) = 0 F(0) = 1 T(i+1)=T(i)+F(i) N(i+1)=N(i)+6 F(i+1)=F(i)+N(i) 则，T(0)...T(i)形成一系列的平方数
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复

使用道具 举报

无心人
 楼主| 发表于 2011-9-12 20:58:05 | 显示全部楼层
代码
  2. #include <stdio.h>
  3. #include <stdlib.h>
  5. #define M 6
  6. #define D6 1000000
  8. int D[D6];
  9. int T[M];
  10. int N[M];
  11. int F[M];
  12. int dT, dN, dF;
  13. int mask[480];
  15. void init(void)
  16. {
  17. int i, j, k, l, m, n;
  18. for (i = 0; i <= 9; i ++)
  19. for (j = 0; j <= 9; j ++)
  20. for (k = 0; k <= 9; k ++)
  21. for (l = 0; l <= 9; l ++)
  22. for (m = 0; m <= 9; m ++)
  23. for (n = 0; n <= 9; n ++)
  24. D[i*100000+j*10000+k*1000+l*100+m*10+n] = i + j + k + l + m + n;
  26. for (i = 0; i < M; i ++)
  27. T[i] = N[i] = F[i] = 0;
  28. F[0] = 1;
  29. dT = dN = dF = 1;
  30. for (i = 0; i < 480; i ++)
  31. mask[i] = 0;
  32. }
  34. void newT( void )
  35. {
  36. int i, t = 0;
  37. for (i = 0; i < dF; i ++)
  38. {
  39. t += T[i] + F[i];
  40. t >= D6 ? (T[i] = t - D6, t = 1) : (T[i] = t, t = 0);
  41. }
  42. i = dF;
  43. while ((t > 0) && (i < dT))
  44. {
  45. t += T[i];
  46. t >= D6 ? (T[i] = t - D6, t = 1) : (T[i] = t, t = 0);
  47. i ++;
  48. }
  49. if (t > 0)
  50. T[dT++] = 1;
  51. }
  54. void newN( void )
  55. {
  56. int i = 0, t = 6;
  57. while ((t > 0) && (i < dN))
  58. {
  59. t += N[i];
  60. t >= D6 ? (N[i] = t - D6, t = 1) : (N[i] = t, t = 0);
  61. i ++;
  62. }
  63. if (t > 0)
  64. N[dN ++] = 1;
  65. }
  67. void newF( void )
  68. {
  69. int i, t = 0;
  70. for (i = 0; i < dN; i ++)
  71. {
  72. t += F[i] + N[i];
  73. t >= D6 ? (F[i] = t - D6, t = 1) : (F[i] = t, t = 0);
  74. }
  75. i = dN;
  76. while ((t > 0) && (i < dF))
  77. {
  78. t += F[i];
  79. t >= D6 ? (F[i] = t - D6, t = 1) : (F[i] = t, t = 0);
  80. i ++;
  81. }
  82. if (t > 0)
  83. F[dF++] = 1;
  84. }
  86. void checkT( void )
  87. {
  88. int i, t = 0;
  89. for (i = 0; i < dT; i ++)
  90. t += D[T[i]];
  91. if (mask[t] == 0)
  92. {
  93. printf("Find %d: ", t);
  94. i = dT - 1;
  95. printf("%d", T[i--]);
  96. while (i >= 0)
  97. printf("%06d", T[i--]);
  98. printf("\n");
  99. mask[t] = 1;
  100. }
  101. }
  103. void checkVar( int * v, int len)
  104. {
  105. if (len <= 0)
  106. return ;
  107. printf("%d", v[-- len]);
  108. while ( len > 0)
  109. printf("%06d", v[-- len]);
  110. printf("\n");
  111. }
  113. int main( void )
  114. {
  115. long long i = 0;
  116. init( );
  117. while ( i <= 100000000 )
  118. {
  119. //printf("%d: %d %d %d T, N, F====\n", i, dT, dN, dF);
  120. //checkVar(T, dT);
  121. //checkVar(N, dN);
  122. //checkVar(F, dF);
  123. newT( );
  124. newN( );
  125. newF( );
  126. checkT( );
  127. i ++;
  128. }
  129. }
复制代码
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

无心人
 楼主| 发表于 2011-9-12 20:59:13 | 显示全部楼层
10^24内结果
  2. Find 1: 1
  3. Find 8: 8
  4. Find 9: 27
  5. Find 10: 64
  6. Find 18: 729
  7. Find 19: 2197
  8. Find 17: 2744
  9. Find 28: 6859
  10. Find 26: 17576
  11. Find 27: 19683
  12. Find 35: 148877
  13. Find 36: 287496
  14. Find 37: 438976
  15. Find 44: 778688
  16. Find 45: 2299968
  17. Find 46: 3869893
  18. Find 53: 43986977
  19. Find 54: 75686967
  20. Find 55: 174676879
  21. Find 62: 596947688
  22. Find 63: 796597983
  23. Find 64: 1693669888
  24. Find 73: 7598896696
  25. Find 71: 9649992689
  26. Find 72: 56888939736
  27. Find 80: 78898389569
  28. Find 81: 197747699976
  29. Find 82: 677298787768
  30. Find 89: 1778597976896
  31. Find 98: 1999899757799
  32. Find 90: 4329949978899
  33. Find 91: 6999657683689
  34. Find 100: 27969886988875
  35. Find 99: 68986579789197
  36. Find 107: 254996898985997
  37. Find 108: 299982589968888
  38. Find 109: 994878767299789
  39. Find 116: 2679768989786699
  40. Find 118: 9949889468976787
  41. Find 117: 18696769447799997
  42. Find 125: 39677989979796875
  43. Find 126: 129937879999769499
  44. Find 127: 188979889884278779
  45. Find 135: 849894379779889989
  46. Find 136: 968397868897889977
  47. Find 134: 1737876998969969768
  48. Find 144: 9688689596689799688
  49. Find 143: 19786889787799686488
  50. Find 145: 48879314999968898899
  51. Find 154: 67988999959868987875
  52. Find 153: 85897988298989489979
  53. Find 152: 179999887787788988375
  54. Find 162: 2979909779796589978989
  55. Find 161: 4479959888985687898688
  56. Find 163: 8783893698974994898999
  57. Find 170: 67899899689199886577697
  58. Find 171: 68885699985378975978999
  59. Find 172: 84997849968878977898776
  60. Find 179: 529967889999499896479969
  61. Find 181: 688859969398989879749896
复制代码
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

无心人
 楼主| 发表于 2011-9-12 21:59:59 | 显示全部楼层
四次方递归公式 Q(0)=0 F(0)=1 N(0)=14 M(0)=36 Q(i+1)=Q(i)+F(i) F(i+1)=F(i)+N(i) N(i+1)=N(i)+M(i) M(i+1)=M(i)+24 Q(i)组成四次方序列
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

sprime
发表于 2012-2-10 13:52:18 | 显示全部楼层
多谢分享！
毋因群疑而阻独见　　毋任己意而废人言
毋私小惠而伤大体　　毋借公论以快私情
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-11-21 21:00 , Processed in 0.026780 second(s), 17 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表