新版本公告
~~空空如也
加载中
加载中
表情图片
评为精选
鼓励
加载中...
分享
加载中...
文件下载
加载中...

贴一下我写的程序...

import multiprocessing
import time

def calc(cpu_number, cpu_total, result, total):
    sumlimit = 1000
    for n in range(cpu_number, total, cpu_total):
        s1 = 0
        s2 = 0
        sum = 0
        for k in range(0, n):
            s1 += (pow(16, n-k) % (k * 8 + 1)) / (k * 8 + 1)
        for k in range(n+1, sumlimit):
            s2 += (pow(16, n-k)) / (k * 8 + 1)
        sum += (s1 + s2) * 4
        s1 = 0
        s2 = 0
        for k in range(0, n):
            s1 += (pow(16, n-k) % (k * 8 + 4)) / (k * 8 + 4)
        for k in range(n+1, sumlimit):
            s2 += (pow(16, n-k)) / (k * 8 + 4)
        sum -= (s1 + s2) * 2
        s1 = 0
        s2 = 0
        for k in range(0, n):
            s1 += (pow(16, n-k) % (k * 8 + 5)) / (k * 8 + 5)
        for k in range(n+1, sumlimit):
            s2 += (pow(16, n-k)) / (k * 8 + 5)
        sum -= (s1 + s2)
        s1 = 0
        s2 = 0
        for k in range(0, n):
            s1 += (pow(16, n-k) % (k * 8 + 6)) / (k * 8 + 6)
        for k in range(n+1, sumlimit):
            s2 += (pow(16, n-k)) / (k * 8 + 6)
        sum -= (s1 + s2)
        ret = (sum % 1) * 16
        if ret < 0:
            ret += 16
        result[n] = int(ret)

def main():
    print('Using CPU')
    startt = time.time()
    grid = 4
    block = 1024
    total = grid * block
    cpu_n = multiprocessing.cpu_count()
    cpu_result = multiprocessing.Array('d', total)
    process = []
    for i in range(0, cpu_n):
        print('Starting CPU ' + str(i))
        p = multiprocessing.Process(target=calc, args=(i, cpu_n, cpu_result, total))
        process.append(p)
        p.start()
    for i in process:
        i.join()
    f = open('pi_bin_CPU.txt''w')
    for i in range(0, total):
        f.write((bin(int(cpu_result[i]))).replace('0b','').zfill(4) + ' ')
    f.close()
    final = 3.140625
    for i in range(2, total):
        try:
            final = final + cpu_result[i] / pow(16, i+1)
        except:
            break
    print(final)
    print('Time elapsed:' + str(time.time() - startt))

if __name__ == "__main__":
    main()

自学了python一两天,程序写不好见谅

文号 / 885743

百炼成钢
名片专栏发私信
学术分 0
总主题 17 帖总回复 848 楼拥有证书:进士 机友 笔友
注册于 2019-12-22 09:55最后登录 2025-02-20 00:34
主体类型:个人
所属领域:无
认证方式:身份证号
IP归属地:浙江

个人简介


UNITED INTERGALACTIC


Shanghai Jiao Tong University
Artificial Intelligence

AI

Physics

Informatics

Aerospace

Inertial Navigation

Energetic Materials


Visit Website

文件下载
加载中...
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
投诉或举报
加载中...
{{tip}}
请选择违规类型:
{{reason.type}}

空空如也

插入资源
全部
图片
视频
音频
附件
全部
未使用
已使用
正在上传
空空如也~
上传中..{{f.progress}}%
处理中..
上传失败,点击重试
等待中...
{{f.name}}
空空如也~
(视频){{r.oname}}
{{selectedResourcesId.indexOf(r.rid) + 1}}
处理中..
处理失败
插入表情
我的表情
共享表情
Emoji
上传
注意事项
最大尺寸100px,超过会被压缩。为保证效果,建议上传前自行处理。
建议上传自己DIY的表情,严禁上传侵权内容。
点击重试等待上传{{s.progress}}%处理中...已上传,正在处理中
空空如也~
处理中...
处理失败
加载中...
草稿箱
加载中...
此处只插入正文,如果要使用草稿中的其余内容,请点击继续创作。
{{fromNow(d.toc)}}
{{getDraftInfo(d)}}
标题:{{d.t}}
内容:{{d.c}}
继续创作
删除插入插入
插入公式
评论控制
加载中...
文号:{{pid}}
笔记
{{note.content}}
{{n.user.username}}
{{fromNow(n.toc)}} {{n.status === noteStatus.disabled ? "已屏蔽" : ""}} {{n.status === noteStatus.unknown ? "正在审核" : ""}} {{n.status === noteStatus.deleted ? '已删除' : ''}}
  • 编辑
  • 删除
  • {{n.status === 'disabled' ? "解除屏蔽" : "屏蔽" }}
我也是有底线的