已屏蔽 原因:{{ notice.reason }}已屏蔽
{{notice.noticeContent}}
~~空空如也

贴一下我写的程序...

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
总主题 15 帖总回复 815 楼拥有证书:进士 机友 笔友
注册于 2019-12-22 09:55最后登录 2024-04-28 22:00
主体类型:个人
所属领域:无
认证方式:身份证号
IP归属地:浙江

个人简介


UNITED INTERGALACTIC


Shanghai Jiao Tong University
Artificial Intelligence

AI

Physics

Informatics

Aerospace

Inertial Navigation

Energetic Materials


Visit Website

文件下载
加载中...
{{errorInfo}}
{{downloadWarning}}
你在 {{downloadTime}} 下载过当前文件。
文件名称:{{resource.defaultFile.name}}
下载次数:{{resource.hits}}
上传用户:{{uploader.username}}
所需积分:{{costScores}},{{holdScores}}下载当前附件免费{{description}}
积分不足,去充值
文件已丢失

当前账号的附件下载数量限制如下:
时段 个数
{{f.startingTime}}点 - {{f.endTime}}点 {{f.fileCount}}
视频暂不能访问,请登录试试
仅供内部学术交流或培训使用,请先保存到本地。本内容不代表科创观点,未经原作者同意,请勿转载。
音频暂不能访问,请登录试试
投诉或举报
加载中...
{{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}}
加载中...
详情
详情
推送到专栏从专栏移除
设为匿名取消匿名
查看作者
回复
只看作者
加入收藏取消收藏
收藏
取消收藏
折叠回复
置顶取消置顶
评学术分
鼓励
设为精选取消精选
管理提醒
编辑
通过审核
评论控制
退修或删除
历史版本
违规记录
投诉或举报
加入黑名单移除黑名单
查看IP
{{format('YYYY/MM/DD HH:mm:ss', toc)}}
ID: {{user.uid}}