技术宅的结界

 找回密码
 立即注册→加入我们

QQ登录

只需一步,快速开始

搜索
热搜: 下载 VB C 实现 编写
查看: 1302|回复: 11
收起左侧

经典算法求PI

[复制链接]

26

主题

114

帖子

1620

积分

用户组: 版主

UID
1821
精华
6
威望
67 点
宅币
1291 个
贡献
51 次
宅之契约
0 份
在线时间
258 小时
注册时间
2016-7-12
发表于 2017-4-18 18:03:16 | 显示全部楼层 |阅读模式

欢迎访问技术宅的结界,请注册或者登录吧。

您需要 登录 才可以下载或查看,没有帐号?立即注册→加入我们

x
蒲丰投针,其实这算法本质是个求定积分的算法
[C] 纯文本查看 复制代码
#include <windows.h>
#define N 10000
#define C 20
#define R 19

__int64 a=1000000, b=0, c=N*R, d, e, f[N*R+1], g;   
main()   
{    
for (;b-c;) f[b++]=C*N;
for(;d=0,g=c*2;c-=R,printf("%.6d",e+d/a),e=d%a)   
for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);   
system("pause");
return 0;   
}   

本帖被以下淘专辑推荐:

1044

主题

2345

帖子

5万

积分

用户组: 管理员

一只技术宅

UID
1
精华
218
威望
294 点
宅币
18326 个
贡献
37541 次
宅之契约
0 份
在线时间
1749 小时
注册时间
2014-1-26
发表于 2017-4-18 19:29:01 | 显示全部楼层
无力吐槽abcdefg,不过代码还真简短

0

主题

8

帖子

64

积分

用户组: 小·技术宅

UID
2357
精华
0
威望
1 点
宅币
54 个
贡献
0 次
宅之契约
0 份
在线时间
10 小时
注册时间
2017-3-28
发表于 2017-4-23 01:11:50 | 显示全部楼层
我今天在http://outofmemory.cn/code-snipp ... c-algorithm-PI-tell 上也看到一个求pi的算法了呢。不过好长,囧。

12

主题

50

帖子

386

积分

用户组: 中·技术宅

UID
2285
精华
0
威望
29 点
宅币
277 个
贡献
1 次
宅之契约
0 份
在线时间
30 小时
注册时间
2017-2-25
发表于 2017-5-4 12:45:50 | 显示全部楼层
N久前看到过,3行代码计算PI的东西,你这个有点长了。。。

26

主题

114

帖子

1620

积分

用户组: 版主

UID
1821
精华
6
威望
67 点
宅币
1291 个
贡献
51 次
宅之契约
0 份
在线时间
258 小时
注册时间
2016-7-12
 楼主| 发表于 2017-5-4 19:23:04 | 显示全部楼层
乘简 发表于 2017-5-4 12:45
N久前看到过,3行代码计算PI的东西,你这个有点长了。。。

当然有点长 因为对原代码进行了词义解释 方便计算精度的扩展

1044

主题

2345

帖子

5万

积分

用户组: 管理员

一只技术宅

UID
1
精华
218
威望
294 点
宅币
18326 个
贡献
37541 次
宅之契约
0 份
在线时间
1749 小时
注册时间
2014-1-26
发表于 2017-5-5 07:12:59 | 显示全部楼层
乘简 发表于 2017-5-4 12:45
N久前看到过,3行代码计算PI的东西,你这个有点长了。。。

代码长度不能决定其优秀程度

12

主题

50

帖子

386

积分

用户组: 中·技术宅

UID
2285
精华
0
威望
29 点
宅币
277 个
贡献
1 次
宅之契约
0 份
在线时间
30 小时
注册时间
2017-2-25
发表于 2017-5-26 16:58:34 | 显示全部楼层
本帖最后由 乘简 于 2017-5-26 17:03 编辑

原来代码是一样的。。。

12

主题

50

帖子

386

积分

用户组: 中·技术宅

UID
2285
精华
0
威望
29 点
宅币
277 个
贡献
1 次
宅之契约
0 份
在线时间
30 小时
注册时间
2017-2-25
发表于 2017-5-26 17:00:32 | 显示全部楼层
本帖最后由 乘简 于 2017-5-26 17:03 编辑
0xAA55 发表于 2017-5-5 07:12
代码长度不能决定其优秀程度


不好意思,发错了。。。

12

主题

50

帖子

386

积分

用户组: 中·技术宅

UID
2285
精华
0
威望
29 点
宅币
277 个
贡献
1 次
宅之契约
0 份
在线时间
30 小时
注册时间
2017-2-25
发表于 2017-5-26 17:02:15 | 显示全部楼层
0xAA55 发表于 2017-5-5 07:12
代码长度不能决定其优秀程度

不好意思,是我看错了。。。

1

主题

85

帖子

91

积分

用户组: 小·技术宅

UID
3026
精华
0
威望
1 点
宅币
3 个
贡献
0 次
宅之契约
0 份
在线时间
6 小时
注册时间
2017-10-31
发表于 2017-11-7 23:25:45 | 显示全部楼层
这就是精髓的

17

主题

184

帖子

1303

积分

用户组: 上·技术宅

UID
3808
精华
5
威望
31 点
宅币
972 个
贡献
60 次
宅之契约
0 份
在线时间
189 小时
注册时间
2018-5-6
发表于 2018-9-7 08:18:21 | 显示全部楼层
蒙特卡洛方法与数值模拟第一课就是讲的蒲丰式投针实验,
菜鸟一枚,直接指正,不必留情

0

主题

23

帖子

1198

积分

用户组: 上·技术宅

UID
3513
精华
0
威望
12 点
宅币
1151 个
贡献
0 次
宅之契约
0 份
在线时间
101 小时
注册时间
2018-3-1
发表于 2019-3-7 06:39:53 | 显示全部楼层
发帖辛苦,谢谢分享

本版积分规则

QQ|申请友链||Archiver|手机版|小黑屋|技术宅的结界 ( 滇ICP备16008837号|网站地图

GMT+8, 2019-11-14 14:48 , Processed in 0.111024 second(s), 36 queries , Gzip On.

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.

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