主题:求1000的阶乘的代码及思路?
gusoul1
[专家分:0] 发布于 2010-11-02 18:41:00
知道的说下
回复列表 (共10个回复)
沙发
yeg001 [专家分:14390] 发布于 2010-11-02 19:54:00
这个问题以前有人讨论过, 而且有代码. 可以搜索一下.
板凳
gusoul1 [专家分:0] 发布于 2010-11-02 22:48:00
我需要的是用Fortran语言写的
3 楼
jstzhurj [专家分:4680] 发布于 2010-11-03 15:10:00
program main
implicit none
integer :: i,k,t,max,x=1000
parameter (max=8192)
integer,dimension(max) :: z=0
z(1)=1
do k=2,x
do i=1,max
t=z(i)*k+t
z(i)=mod(t,10)
t=t/10
enddo
enddo
t=0
do i=max,1,-1
if(t/=0.or.z(i)/=0) then
write(*,100) z(i)
t=1
endif
enddo
100 format(i1,\)
end program main
402387260077093773543702433923003985719374864210714632543799910429938512398629020592044208486969404800479988610197196058631666872994808558901323829669944590997424504087073759918823627727188732519779505950995276120874975462497043601418278094646496291056393887437886487337119181045825783647849977012476632889835955735432513185323958463075557409114262417474349347553428646576611667797396668820291207379143853719588249808126867838374559731746136085379534524221586593201928090878297308431392844403281231558611036976801357304216168747609675871348312025478589320767169132448426236131412508780208000261683151027341827977704784635868170164365024153691398281264810213092761244896359928705114964975419909342221566832572080821333186116811553615836546984046708975602900950537616475847728421889679646244945160765353408198901385442487984959953319101723355556602139450399736280750137837615307127761926849034352625200015888535147331611702103968175921510907788019393178114194545257223865541461062892187960223838971476088506276862967146674697562911234082439208160153780889893964518263243671616762179168909779911903754031274622289988005195444414282012187361745992642956581746628302955570299024324153181617210465832036786906117260158783520751516284225540265170483304226143974286933061690897968482590125458327168226458066526769958652682272807075781391858178889652208164348344825993266043367660176999612831860788386150279465955131156552036093988180612138558600301435694527224206344631797460594682573103790084024432438465657245014402821885252470935190620929023136493273497565513958720559654228749774011413346962715422845862377387538230483865688976461927383814900140767310446640259899490222221765904339901886018566526485061799702356193897017860040811889729918311021171229845901641921068884387121855646124960798722908519296819372388642614839657382291123125024186649353143970137428531926649875337218940694281434118520158014123344828015051399694290153483077644569099073152433278288269864602789864321139083506217095002597389863554277196742822248757586765752344220207573630569498825087968928162753848863396909959826280956121450994871701244516461260379029309120889086942028510640182154399457156805941872748998094254742173582401063677404595741785160829230135358081840096996372524230560855903700624271243416909004153690105933983835777939410970027753472000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000Press any key to continue
4 楼
茶竹论道 [专家分:40] 发布于 2010-11-03 15:15:00
有人讨论过, 而且有代码. .[url=http://www.shenggang.com]磁力泵[/url]可以搜索一下
5 楼
jstzhurj [专家分:4680] 发布于 2010-11-03 17:25:00
来个效率高一点的
program main
implicit none
integer :: used,i,k,t,max,x=1000
parameter (max=1000)
integer,dimension(max) :: z=0
used=1
z(1)=1
do k=2,x
do i=1,used
t=z(i)*k+t
z(i)=mod(t,10000)
t=t/10000
if (used==i.and.z(i)/=0) used=used+1
enddo
enddo
do i=used,1,-1
write(*,100) z(i)
enddo
100 format(i4.4,\)
end program main
6 楼
dongyuanxun [专家分:7180] 发布于 2010-11-03 17:56:00
这种东西还是参考GMP比较好,虽然是C/C++的,自己可以混合编程或者把GMP封装成Fortran接口。
算1000!也就是几毫秒的事情
7 楼
gusoul1 [专家分:0] 发布于 2010-11-03 23:21:00
结果错误
8 楼
dongyuanxun [专家分:7180] 发布于 2010-11-03 23:25:00
[quote]结果错误[/quote]
你是说3L的结果?
我用GMPlib算1000!和他结果一样啊
9 楼
yrliu [专家分:750] 发布于 2010-11-06 19:37:00
计算机计算真是高深有趣哪:http://www.mathmu.cn/Doc.html
10 楼
yrliu [专家分:750] 发布于 2010-11-06 19:40:00
没发现3楼结果有错误,mathematica 计算也是这个
402387260077093773543702433923003985719374864210714632543799910429938512398629020592044208486969404800479988610197196058631666872994808558901323829669944590997424504087073759918823627727188732519779505950995276120874975462497043601418278094646496291056393887437886487337119181045825783647849977012476632889835955735432513185323958463075557409114262417474349347553428646576611667797396668820291207379143853719588249808126867838374559731746136085379534524221586593201928090878297308431392844403281231558611036976801357304216168747609675871348312025478589320767169132448426236131412508780208000261683151027341827977704784635868170164365024153691398281264810213092761244896359928705114964975419909342221566832572080821333186116811553615836546984046708975602900950537616475847728421889679646244945160765353408198901385442487984959953319101723355556602139450399736280750137837615307127761926849034352625200015888535147331611702103968175921510907788019393178114194545257223865541461062892187960223838971476088506276862967146674697562911234082439208160153780889893964518263243671616762179168909779911903754031274622289988005195444414282012187361745992642956581746628302955570299024324153181617210465832036786906117260158783520751516284225540265170483304226143974286933061690897968482590125458327168226458066526769958652682272807075781391858178889652208164348344825993266043367660176999612831860788386150279465955131156552036093988180612138558600301435694527224206344631797460594682573103790084024432438465657245014402821885252470935190620929023136493273497565513958720559654228749774011413346962715422845862377387538230483865688976461927383814900140767310446640259899490222221765904339901886018566526485061799702356193897017860040811889729918311021171229845901641921068884387121855646124960798722908519296819372388642614839657382291123125024186649353143970137428531926649875337218940694281434118520158014123344828015051399694290153483077644569099073152433278288269864602789864321139083506217095002597389863554277196742822248757586765752344220207573630569498825087968928162753848863396909959826280956121450994871701244516461260379029309120889086942028510640182154399457156805941872748998094254742173582401063677404595741785160829230135358081840096996372524230560855903700624271243416909004153690105933983835777939410970027753472000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000
我来回复