[原创]请教vfbpgyfk先生一个编程问题

白新岭

最后一组33生素数式不是33生素数，其中11个是素数，其中3个是3个素数的积，另外19个是2个素数的积，最后一组33生素数式为：510329,        510331,        510337,        510343,        510347,        510353,        510359,        510361,        510371,        510373,        510379,        510383,        510397,        510401,        510403,        510407,        510409,        510413,        510421,        510427,        510431,        510437,        510439,        510443,        510449,        510451,        510457,        510463,        510467,        510469,        510473,        510479,        510481,

白新岭

在素数式13中产生了40组跨度为152的最密33生素数式，到了17时就剩了18组14种最密33生素数式了，所以在小点的区间内得到基础数据是不稳定的，无法完成一直的连续膨胀，只能昙花一现就消失了。但是从33生素数式能在素数式17以后保持稳定增倍。对于33以下的k生素数式是可以找到其基础数据的。

白新岭

在用37素数式中最密33生素数式生成41素数式最密33生素数式时，因为其数值超出精确度（15位有效数字），或者超出MOD（）函数的计算范围，用程序获得的最密33生素数式的数量为163192组。为了验证其结果，我把周期304250263527210对41求余来代替其周期。（2*3*5*7*11*13*17*19*23*29*31*37*41=304250263527210），把步长弄错了，应是2*3*5*7*11*13*17*19*23*29*31*37=7420738134810

白新岭

在素数式17中产生了8种最密32生素数式（每种1组），到了素数式19就有4种夭折了，剩了4种，每种2组，最密32生素数式的间距为146

白新岭

在素数式17中产生了8种最密32生素数式（每种1组），到了素数式19就有4种夭折了，生了4种，每种2组，最密32生素数式的间距为146

白新岭

2010年8月8下午：18.47分  白新岭作如下分析                                                                       
合成方法：对于一组互质条件下，三元不定方程来说，每个条件有p-1种余数，余数的3维合成有                                                                       
(p-1)^3种方法，分别落到P种余数上，因为总方法为：P^3-3*P^2+3*P-1种，此数除P,少一种方法不能                                                                       
均分，当P≥3时，除不尽的占P-1种余数，大于或等于2，而且这些余数处于同一地位（在群中可以                                                                       
相互置换，都有与本身不同的逆元），所以除尽的（只有一类）与除不尽的扮演不同的角色，                                                                       
这样少的一种方法肯定是在整除类上，所以能整除的占P^2-3P+2种方法，不能整除的其余P-1类各占                                                                       
P^2-3P+3种方法，调节系数=P*对应的合成方法/总合成方法，即系数=P*(P^2-3P+3)/(P-1)^3,不能整                                                                       
除，能整除的系数=P*(P^2-3P+2)/(P-1)^3；整理后，系数=1+1/（P-1)^3,(不能整除）；                                                                       
系数=P(P-2)/(P-1)^2=(P^2-2P)/（P-1)^2=1-1/（P-1)^2,(能整除）。                                                                       
根据上述分析，20111不能整除2,3,5；20111能整除7，所以其综合系数=2.22167969                                                                       
渐进公式=调节系数*符合条件元素个数^m/n/(m-1)!=2.22167969*4596^3/20111/2=5362383组                                                                       
与实际组数（5366913）比较，差4530组，相对误差0.000844061，如果用20111*（1-1/2）*（1-1/3）*                                                                       
（1-1/5）*（1-1/7）=4596.8代替4596，则渐进公式的值=5365183，与实际相差1730组，相对                                                                       
误差0.000322345 ，所以用渐进公式获得的近似值与实际值相对误差可以趋近于0（随N的增大）。                                                                       
zhaolu先生给的27831，相差悬殊，不敢认同，即便改成排列数，27831*3！=27831*6=166986组，仅占实际的3.1%                                                                       
这应该是x+y+z=20111的正整解的问题（这里限定x，y，z不能是含有2,3
，5,7因子的自然数，其余的都可以取）

白新岭

2011年11月15日        中午从新建立
一般的系数形式        P^(k-1)*(P-k)/(P-1)^k
可以化简成        (1-k/P)/(1-1/P)^k
        只是分子部分的k要根据
        具体分析而定.
E----计算k生素数的系数

白新岭

温故而知新

白新岭

顶起此贴，有空时还学习一下编程问题

白新岭

1，求1到2000之间，能被3或7整除，并且含有数字5的个数，它们的和，平方根和（保留两位小数），以及他们的平方和。
这是网上的一道试题，也给了答案
clear
store 0 to N,S1,S2,S3
for i=1 to 2000
   if (mod(i,3)=0 or i%7=0)and '5'$str(i)
       N=N+1
       S1=S1+i
       S2=S2+sqrt(i)
       S3=S3+i^2
   endif
   next
   ?"个数：" ,N, "和：" ,S1, "平方根的和：" ,S2, "平方和：" ,S3
   return
上边是程序，可是刚开始的时候，不能正常运行，提示有不能识别的短语，在if行，我就把括号从新改了一下，原来写的时候可能没有注意到是否在英文的标点符号下写的括号，改后，那一行可以过去了；仅接着是找不到变量，还是同样的问题，也是出在标点符号上，因为store 0 to 后边的变量间的逗号是在中文情况下输进去的，所以找不到变量，改了后程序运行正常，得到了结果，所以在编写程序时一定是在英文的标点符号下编写，否则就会出现各种提示，一定注意，这种情况出错还不好找到原因（从程序上说）。