在VF语言环境下找出随机产生的m位的偶数所表示为的素数对

zhaolu48

“请求计算：
       1                  
————————，  
(k^n)ln(x/k^n)
这里 x=10^16，k=1.00008，n=1，2，……，400000，”
用如下程序计算
CLEAR
CLEAR ALL
SET TALK off
p=0
k=1.00008
a=1
b=16*log(10)
d=log(k)
c=0
s=0
for i=1 to 50
sa=0
FOR j=1 to 100
sb=0
FOR l=1 to 100
a=a*k
c=c+d
e=a*(b-c)
sb=sb+(1/e)
p=p+1
endf
sa=sa+sb
endf
s=s+sa
?p,s,1/e,sa,sb
endf
?p,s,1/e,sa,sb
SET TALK ON
用的是“n=1，2，……，500000”，这样结果能达到理想程度。
计算得ls=349.0483205148932

zhaolu48

请求计算：
      1
————————，  
(k^n)ln(x/k^n)
这里x = 1.00004×10^21，k = 1.00008，n = 1，2，……，500000，
用下面程序计算：
CLEAR
CLEAR ALL
SET TALK off
*这里x = 1.00004×10^21，k = 1.00008，n = 1，2，……，500000，
p=0
k=1.00008
a=1
b=21*log(10)+log(1.00004)
d=log(k)
c=0
s=0
for i=1 to 50
sa=0
FOR j=1 to 100
sb=0
FOR l=1 to 100
a=a*k
c=c+d
e=a*(b-c)
sb=sb+(1/e)
p=p+1
endf
sa=sa+sb
endf
s=s+sa
?p,s,1/e,sa,sb
endf
?p,s,1/e,sa,sb
SET TALK ON
计算得s=264.0908685808902

zhaolu48

“请求计算：
      1
————————，  
(k^n)ln(x/k^n)
这里 x = 1.000004×10^21，k = 1.000008，n = 1，2，……，5000000，”
用下面程序计算：
CLEAR
CLEAR ALL
SET TALK off
*这里 x = 1.000004×10^21，k = 1.000008，n = 1，2，……，5000000，
p=0
k=1.000008
a=1
b=21*log(10)+log(1.000004)
d=log(k)
c=0
s=0
for i=1 to 5
sa=0
FOR j=1 to 100
sb=0
FOR l=1 to 100
sc=0
FOR u=1 to 100
a=a*k
c=c+d
e=a*(b-c)
sc=sc+(1/e)
p=p+1
endf
sb=sb+sc
endf
sa=sa+sb
endf
s=s+sa
?p,s,1/e,sa,sb
endf
?p,s,1/e,sa,sb
SET TALK ON
得s=2640.908685824159

zhaolu48

              1
J(x) = ────────
      ln(x + n×10^8)  
设 x = 10^16，n = 1，2，…… ， 9,999,900,000,000
J(10^16) =
可能还是求和吧。用下面程序计算
x=10^16
a=10^8
s=0
b=0
p=0
for i=1 to 99999
sa=0
for j=1 to 10000
sb=0
for l=1 to 10000
b=b+a
c=log10(x+b)
sb=sb+1/c
p=p+1
endf
sa=sa+sb
endf
s=s+sa
?p,s
endf
?p,s
事实上是不可能的
因为当n=100000000时，就用时2分50秒，
那么当n=9999900000000时，用时要近197天。

zhaolu48

不过这个和相当于
a(n)=1/(36.84136149+n*0.000,000,01)
求s(n)
这个问题最好求教陆老师。
我也先做做看。如果做成，用DOC文档发上来。

zhaolu48

zhaolu48

能否计算：
               1
J(x) = ─────────
       ln(x + n×10^13)  
设 x = 10^16，n = 1，2，…… ，99,990,000  
a=10^16
b=10^13
s=0
m=0
p=0
for i=1 to 9999
sa=0
for j=1 to 100
sb=0
for l=1 to 100
p=p+1
m=m+b
c=log(a+m)
sb=sb+1/c
if mod(p,1000000)=0
?p,s
endi
endf
sa=sa+sb
endf
s=s+sa
endf
?p,s
s=2112512.895221372

白新岭

赵老师，能不能计算一下偶数104053950的素数对（无序数目）。

zhaolu48

白先生您好。
从找到的素数对看出，素数对的个数与这个大偶数的比约为0.02，
104053950*0.02约为2081000，每找1000个素数对用时约为66秒。
因此要找出全部素数对大约要38小时。
因此这样大的偶数要找出全部素数对是我很难完成的。
现在我只找出65711个素数对。
很对不起，完不成你交给我的任务。

zhaolu48

计算：
           1
J(x) = ─────
       ln(x + n)  
设 x = 100，n = 1，2，…… ，99,999,900
J(100) =
a=0
b=100
s=0
for i=1 to 100
sa=0
for j=1 to 1000
sb=0
for j=1 to 1000
a=a+1
c=log(a+b)
sb=sb+1/c
endf
sa=sa+sb
endf
s=s+sa
?a,s
endf
sc=0
for i=99999901 to 100000000
sc=sc+1/log(100+i)
endf
?sc
s=s-sc
?';s=';,s
s=9602.874582499207