趣味VisualFoxPro程序设计集锦.docx
《趣味VisualFoxPro程序设计集锦.docx》由会员分享，可在线阅读，更多相关《趣味VisualFoxPro程序设计集锦.docx（67页珍藏版）》请在冰豆网上搜索。
趣味VisualFoxPro程序设计集锦
-
趣味VisualFoxPro程序设计集锦
5．舍罕王的失算
settalkoff
input[请输入n:
]ton
|
t=1
s=1
fori=2ton
t=t*2
"
s=s+t
endfor
v=s/
h=（9/*v）^（1/3）
、
"总麦粒数约为:
"+str（s,10,3）+"粒."
"折合体积约为:
"+ltrim（str（v,15））+"立方米"
"堆成正圆锥高约为:
"+str（h,5）+"米"
return
]
程序运行结果与说明
请输入n:
64
总麦粒数约为:
+019
折合体积约为立方米
^
堆成正圆维高约为7193米
5-1失算的另一名题：
买马灯
input"请输入n:
"ton
t=1
\
s=1
fori=2ton
t=t*2
s=s+t
、
endfor
"共需花费"+ltrim（str（s,16））+"个包卢斯卡."
"合约"+str（s/400,10,4）+"卢卡."
return
~
程序运行结果与点评
请输入n:
24
共需花费1677215个外包卢斯卡
合约卢布
、
6分数不等式
"求n<1+1/2+1/3+1/4+...+1/msettalkoff
setdecito14
>
input"请输入n:
"ton
i=0
s=0
dowhiles（
i=i+1
s=s+1/i
enddo
c=i
%
dowhilesi=i+1
s=s+1/i
enddo
!
d=i-1
"满足不等式的m为:
"+ltrim（str（c,10））+"return
程序运行结果：
）
求n<1+1/2+1/3+1/4….+1/m请输入n:
10
满足不等式的m为：
12367请输入n:
15
】
满足不等式的m为：
18354217阶乘与阶乘和数
settalkoff
"三位阶乘和数有:
"
。
form=100to999
a=int（m/100）
b=mod（int（m/10）,10）
c=mod（m,10）
'
ifm=jc（a）+jc（b）+jc（c）
m
endif
endfor
~
return
funcjc
parax
p=1
*
fori=1tox
p=p*i
endfor
returnp
、
程序运行结果：
三位阶乘和数有：
145
可见145是唯一的三位阶乘和数
8综合高精度计算
《
clear
dimea（500）
a=0
"1:
计算阶乘n!
"
【
"2:
计算乘方m^n"
"3:
计算排列数A（n,m）"
input"选择（1-3）：
"toz
input"请输入正整数n:
"ton
^
ifz#1
input"请输入正整数m（m<=n）:
"tom
endif
d=500
;
t=1
ifz=3
t=n-m+1
endif
~
a
（1）=1
fori=tton
ifz=2
b=m
|
else
b=i
endif
f=0
*
forj=1tod
x=a（j）*b+f
f=int（x/10）
a（j）=x%10
—
endfor
endfor
j=d
dowhilea（j）=0
）
j=j-1
enddo
f=j
docase
）
casez=1
str（n,8）+"!
="
casez=2
str（m,6）+"^"+ltrim（str（n））+"="
$
casez=3
"A（"+ltrim（str（n））+","+ltrim（str（m））+"）="
endcase
d=0
-
forj=fto1step-1
d=d+1
str（a（j）,1）
ifd%50=0
；
space（10）
endif
endfor
"所得结果共"+ltrim（str（f））+"位"
|
return
运行程序，选择3，输入n=100，m=90，得
A（100,90）=0000000000000000000
所得结果共152位。
"
10最大公约数与最小公倍数
settalkoff
input"输入正整数a:
"toa
input"输入整数b:
"tob
、
st=ltrim（str（a,12））+","+ltrim（str（b,12））
m=a*b
r=mod（a,b）
dowhiler#0
-
a=b
b=r
r=mod（a,b）
enddo
"
"（"+st+"）="+ltrim（str（b,15））
"（"+st+"）="+ltrim（str（m/b））
return
运行程序，输入1104，1272得
？
（1104，1272）=24
{1104，1272}=58512
11水仙花数
settalkoff
.
input"输入正整数a:
"toa
input"输入整数b:
"tob
st=ltrim（str（a,12））+","+ltrim（str（b,12））
m=a*b
-
r=mod（a,b）
dowhiler#0
a=b
b=r
.
r=mod（a,b）
enddo
"（"+st+"）="+ltrim（str（b,15））
"（"+st+"）="+ltrim（str（m/b））
《
return
运行程序，得水仙花数：
153370371407
12勾股数
~
settalkoff
input"请输入区间下限值:
"tom1
inpu"请输入区间上限值:
"tom2
"区间["+ltrim（str（m1））+","+ltrim（str（m2））+"]中的勾股数组有:
"
[
n=0
fora=2tosqrt（m2）
forb=1toa-1
if（a+b）%2=0
b=b+1
endif
t=0
fork=2tob
（
ifb%k=0anda%k=0
t=1
endif
endfor
！
ift=0
d=a*a-b*b
e=2*a*b
f=a*a+b*b
·
fori=1tom2/f
x=i*d
y=i*e
z=i*f
~
ifxm2
loop
endif
str（x）+str（y）+str（z）
?
n=n+1
endfor
endif
endfor
《
endfor
"共"+ltrim（str（n））+"组勾股数."
return
运行程序，求得[10，30]中的勾股数：
]
区间（10，30）中的勾股数组有：
121620
152025
182430
~
102426
212029
13完全数
settalkoff
%
input[x=]tox
input[y=]toy
"区间["+ltrim（str（x））+","+ltrim（str（y））+"]中的完全数:
"
n=0
·
fora=xtoy
s=1
fork=2toa/2
ifa%k=0
&
s=s+k
endif
endfor
ifs=a
）
n=n+1
str（n,3）+":
"+ltrim（str（a））+"-1"
fork=2toa/2
ifa/k=int（a/k）
<
"+"+ltrim（str（k,5））
endif
endfor
endif
—
endfor
return
运行程序，输入X=2，Y=100，得
区间[2，1000]中的完全数：
'
1：
6=1+2+3
2：
28=1+2+4+7+14
3：
496=1+2=4+8+16+31+62+124+248
14相亲数
[
settalkoff
fori=11to9999
s=1
t=int（sqrt（i）+.01）
：
forj=2tot
ifi%j=0
s=s+j+i/j
endif
]
endfor
ifi=t*t
s=s-t
endif
]
ifis1=1
t=int（sqrt（s）+
forj=2tot
`
ifs%j=0
s1=s1+j+s/j
endif
endfor
）
ifs=t*t
s1=s1-t
endif
ifs1=i
·
"相亲数:
"+str（i）+str（s）
str（i,5）+"的真因数之和为:
1"
forj=2toi/2
ifi%j=0
—
"+"+ltrim（str（j））
endif
endfor
"="+ltrim（str（s））
>
str（s,5）+"的真因数之和为:
1"
forj=2tos/2
ifs%j=0
"+"+ltrim（str（j））
·
endif
endfor
"="+ltrim（str（i））
endif
endif
endfor
return
运行程序，得4位以内的相亲数：
（
相亲数：
220284
220的真因数之和为：
1+2+4+5+10+11+20+22+44+55+110=284
284的真因数之和为：
1+2+4+71+142=220
相亲数11841210
$
1184的真因数之和为：
1+2+4+8+16+32+37+74+148+296+592=1210
1210的真因数之和为：
1+2+5+10+11+22+55+110+121+242+605=1184
相亲数26202924
2620的真因数之和为：
1+2+4+5+10+20+131+262+524+655+1310=2924
￥
2924的真因数之和为：
1+2+4+17+34+43+68+86+172+731+1462=2620
相亲数50205564
5020的真因数之和为：
1+2+4+5+10+20+251+502+1004+1255+2510
=5564
、
5564的真因数之和为：
1+2+4+13+26+52+107+214+428+1391+2782
=5020
相亲数62326368
6232的真因数之和为：
1+2+4+8+19+38+41+76+82+152+164+328+779
+1558+3116=6368
6368的真因数之和为：
1+2+4+8+16+32+199+398+796+1592+3184
=6232
15守形数
>
settalkoff
input[x=]tox
input[y=]toy
fora=xtoy
'
s=a*a
b=1
k=a
dowhilek>0
。
b=b*10
k=int（k/10）
enddo
c=s%b
~
ifa=c
ltrim（str（a））+[^2=]+ltrim（str（s））
endif
endfor
~
return
运行程序，输入X=10，Y=10000得
25^2=625
76^2=5776
'
376^2=141376
625^2=390625
9376^2=
16素数
"
settalkoff
"本程序求区间[c,d]上的素数:
"
input"输入区间下限c:
"toc
input"输入区间上限d:
"tod
*
"区间["+ltrim（str（c））+","+ltrim（str（d））+"]上的素数有：
"
n=0
ifc%2=0
c=c+1
$
endif
""
fori=ctodstep2
t=0
…
forj=3tosqrt（i）step2
ifi%j=0
t=1
exit
[
endif
endfor
ift=0
str（i,6）
、
n=n+1
ifmod（n,10）=0
""
endif
"
endif
endfor
"共"+ltrim（str（n））+"个素数."
return
~
程序运行示例
本程序求区间[C，D]上的数数
输入区间下限C：
2000
输入区间上限D：
2100
"
2003201120172027202920392053206320692081
2083208720892099
共14个素数
17乌兰现象
<
settalkoff
clear
input"请输入数阵的阶数m:
"tom
dimex（（m+1）*（m+1））,y（（m+1）*（m+1））
{
x
（1）=0
y
（1）=0
n=1
t=int（m/2）
^
fori=1tot
n=n+1
x（n）=x（n-1）+1
y（n）=y（n-1）
（
dowhiley（n）n=n+1
x（n）=x（n-1）
y（n）=y（n-1）+1
】
enddo
dowhilex（n）>-i
n=n+1
x（n）=x（n-1）-1
`
y（n）=y（n-1）
enddo
dowhiley（n）>-i
n=n+1
'
x（n）=x（n-1）
y（n）=y（n-1）-1
enddo
dowhilex（n）|
n=n+1
x（n）=x（n-1）+1
y（n）=y（n-1）
enddo
（
endfor
str（m,2）+"阶数字方螺线为:
"
fork=1tom*m
r=0
;
forj=2tosqrt（k）
ifk%j=0
r=1
exit
；
endif
endfor
ifr=0andk#1
@25-y（k）*2,80+x（k）*5saystr（k,3）colorGR+/R
！
else
@25-y（k）*2,80+x（k）*5saystr（k,3）
endif
waitwindtime
;
endfor
return
运行程序，输入M=10，得10阶数字方螺线如下图：
100999897969594939291
-
65646362616059585790
66373635343332315689
67381716151413305588
68391854312295487
！
69401961211285386
70412078910275285
71422122232425265184
72434445464748495083
?
73747576777879808182
乌兰现象再现
18孪生素数
settalkoff
<
"求区间[c,d]上的孪生素数对"
input"请输入c（c>2）:
"toc
input"请输入d:
"tod
f=0
）
n=0
ifc%2=0
c=c+1
endif
|
fori=ctodstep2
t=0
forj=3tosqrt（i）step2
ifi%j=0
#
t=1
exit
endif
endfor
&
ift=0
ifi-f=2
"（"+ltrim（str（f））+","+ltrim（str（i））+"）"
n=n+1
$
endif
f=i
endif
endfor
￥
"共"+ltrim（str（n））+"对孪生素数."
return
运行程序，输入区间[101，200]，得
（101,103）,（107,109）,（137,139）,（149,151）,（179,181）
（
（191,193）,（197,199）
共7对孪生素数
19梅森尼数
settalkoff
/
clear
s=0
t=2
forn=2to50
`
t=t*2
m=t-1
x=0
forj=3tosqrt（m）+1step2
~
ifm%j=0
x=1
exit
endif
：
endfor
ifx=0
s=s+1
"2^"+ltrim（str（n））+"-1="+ltrim（str（m））
endif
endfor
"指数n于[2,50]中梅森尼数共有"+ltrim（str（s））+"个."
return
】
程序运行结果与讨论
2^2-1=3
2^3-1=7
2^5-1=31
]
2^7-1=127
2^13-1=8191
2^17-1=131071
2^19-1=524287
:
2^31-1=47
指数n于[2，50]中的梅森尼数共有8个
20金蝉素数
dimea（5）
；
"金蝉素数为:
"
fork=10001to99999step2
t=0
forj=3tosqrt（k）step2
|
ifk%j=0
t=1
exit
endif
（
endfor
ift=0
a
（1）=k%10
a（5）=int（k/10000）
。
d=int（k/10）%1000
forj=2tosqrt（d）
ifd%j=0
t=1
】
exit
endif
endfor
endif
-
ift=0
a
（2）=d%10
a（4）=int（d/100）
a（3）=int（d/10）%10
]
fori=1to4
forj=i+1to5
ifa（i）=a（j）
t=1
~
exit
endif
endfor
endfor
*
endif
ift=0
forj=1to5
ifa（j）%2=0ora（3）=1ora（3）=9
!
t=1
exit
endif
endfor
$
endif
ift=0
k
endif
!
endfor
return
程序运行，得5个金蝉素数：
1359753791795319157395713
;
在输出的这5个金蝉素数中，13597与79531是互逆的金蝉素数。
21素数多项式
settalkoff
m=0
-
input"请输入f:
"tof
"y=x^2-x+"+str（f,2）+","
"当x取值在[1,"+str（f-1,2）+"],y的素数分布:
"
""
！
forx=1tof-1
y=x*x-x+f
t=0
fork=3tosqrt（y）step2
:
ify%k=0
t=1
exit
endif
|
endfor
ift=0
"x="+str（x,2）+"时,"+str（y,4）+"素数:
"
m=m+1
~
ifm%4=0
""
endif
else
·
"x="+str（x,2）+"时,"+str（y,4）
"="+str（k,3）+"*"+str（y/
endif
endfor
）
return
请输入f:
41
y=x^2+41,当n取值在[1，40]，y的素数分布：
x=1时，41素数，x=2时，43素数，x=3时，47素数，
￥
x=4时，53素数，x=5时，61素数，x=6时，71素数，
x=7时，83素数，x=8时，97素数，x=9时，113素数，
x=10时，131素数，x=11时，151素数，x=12时，173素数，
x=13时，197素数，x=14时，223素数，x=15时，251素数，
`
x=16时，281素数，x=17时，313素数，x=18时，347素数，
x=19时，383素数，x=20时，421素数，x=21时，461素数，
x=22时，503素数，x=23时，547素数，x=24时，593素数，
x=25时，641素数，x=26时，691素数，x=27时，743素数，
》
x=28时，797素数，x=29时，853素数，x=30时，911素数，
x=31时，971素数，x=32时，1033素数，x=33时，1097素数，
x=34时，1163素数，x=35时，1231素数，x=36时，1301素数，
x=37时，1373素数，x=38时，1447素数，x=39时，1523素数，
，
x=40时，1601素数，
验证了当x=1,2,,,,40时，二次多项式x^2-x+41的值均为素数
22等差素数列
"求指定区间内等差素数列的最多项数:
"
;
input"请输入区间下限w0:
"tow0
input"请输入区间上限w:
"tow
dimea（（w-w0）/2）
u=0
（
ifw0%2=0
w0=w0+1
endif
fork=w0towstep2
>
t=0
forj=3tosqrt（k）step2
ifk%j=0
t=1
|
exit
endif
endfor
ift=0
（
u=u+1
a（u）=k
endif
endfor
】
max=0
forn=1tou-1
p=a（u）-a（n）
forj=n+1tou
?
d=a（j）-a（n）
ifd>p/3
exit
endif
}
h=j
m1=2
fori=j+1tou
ifa（i）-a（h）=d
…
h=i
m1=m1+1
endif
ifmax'
max=m1
m=a（n）
d1=d
endif
\
endfor
endfor
endfor
"区间["+ltrim（str（w0））+","+ltrim（str（w））+"]"
{
"内等差素数列最多有"+ltrim（str（max））+"项>"
""
fori=1tomax
ltrim（str（m+（i-1）*d1））+""
；
endfor
return
程序运行示例：
请输入区间下限w0:
3
>
请输入区限上限w:
2009
区间[3，2009]内等差素数列最多有9项：
19940961982910391249145916691879
请输入区间下限w0:
101
}
请输入区限上限w:
1000
区间[3，2009]内等差素数列最多有6项：
107137167197227257
23验证歌德巴赫想
@
settalkoff
"在区间[c,d]中验证哥德巴赫猜想:
"
input"请输入区间下限c="toc
input"请输入区间上限d="tod
（
ifc%2#0
c=
升级会员 