圆周率（@陶远航）

METALARISING

import time
from math import *
from decimal import *
getcontext().prec=10

def ft(value):
    value=int(value)
    if value==0 or value==1:
        return 1
    total=1
    for i in range(1,value+1):
        total*=i
    return total

def rj(frequency):
    getcontext().prec=frequency*8+2
    k=0
    k=Decimal(k)
    result=0
    num1=2*(Decimal(2)**Decimal(0.5))/9801

    while k!=Decimal(frequency):
        num2=Decimal(ft(4*k))*(Decimal(1103)+Decimal(26390)*k)
        num3=Decimal((ft(k)))**Decimal(4)*(Decimal(396)**(Decimal(4)*k))
        num4=num2/num3
        result+=num4
        k+=Decimal(1)

    ending=Decimal(1)/(Decimal(result)*num1)
    return ending

def improvement(frequency):
    t1=time.time()
    getcontext().prec=frequency*15+2
    q=str(time.time()-t1)
    print(q)
    time.sleep(0.5)
    k=0
    k=Decimal(k)
    result=0
    t2=time.time()
    kf=Decimal(10005)**Decimal(0.5)
    w=str(time.time()-t2)
    print(w)
    time.sleep(0.5)
    t3=time.time()
    num1=Decimal(426880)*kf
    e=str(time.time()-t3)
    print(e)
    time.sleep(0.2)
    t=time.time()
    while k!=Decimal(frequency):
        lt=time.time()
        num2=Decimal(ft(6*k))*(Decimal(13591409)+Decimal(545140134)*k)
        num3=Decimal(ft(3*k))*Decimal(ft(k)**3)*(Decimal(-640320)**(3*k))
        result+=num2/num3
        print('次数：'+str(k+1)+'位数：'+str((k+1)*15))
        k+=1
    ending=num1/Decimal(result)
    return ending

x=input('1：限定位数，2：自动增加')
if x=='1':
    c=input('1, 2')
    if c !='1' and c !='2':
        c='3'

    while True:
        print('-'*80)
        if c=='1':
            print('1')
        elif c=='2':
            print('2')
        else:
            print('3')

        d=int(input('小数点后几位'))

        while True:
            if d>5:
                break
            elif d<=0:
                print('??')
            elif d<=5:
                print('这都不会？')

            d=int(input('重写！'))

        result=0
        old_time=time.time()
        if c=='1' or c=='3':
            result=rj(floor(d/8))
            print('1ended:'+str(result))
            print('1time:'+str(time.time()-old_time)+'\n')
            old_time=time.time()
        if c=='2' or c=='3':
            result=improvement(floor(d/15))
            print('2end:'+str(result))
            print('2time:'+str(time.time()-old_time)+'\n')

FishC

这段代码实现了计算圆周率的功能，主要使用了泰勒级数和阶乘的概念。可以分为两个部分：rj函数和improvement函数。

rj函数用来计算圆周率的近似值，输入参数为要计算的位数frequency。在函数内部，首先通过getcontext()函数设置计算精度，然后使用循环计算出近似值的部分，并累加到result变量中。最后通过一些数学运算得到最终的结果。

improvement函数也是用来计算圆周率的近似值，与rj函数不同的是，improvement函数不需要输入位数，它会自动增加计算位数直到达到指定的精度。同样，该函数也是通过循环计算出近似值的部分，并累加到result变量中。最后通过一些数学运算得到最终的结果。

可以根据需要选择调用rj函数或improvement函数来计算圆周率的近似值。

以下是一个示例代码：

1. import time
   
2. from math import *
   
3. from decimal import *
   
4. getcontext().prec=10
   
5. 
   
6. def ft(value):
   
7.     value=int(value)
   
8.     if value==0 or value==1:
   
9.         return 1
   
10.     total=1
    
11.     for i in range(1,value+1):
    
12.         total*=i
    
13.     return total
    
14. 
    
15. def rj(frequency):
    
16.     getcontext().prec=frequency*8+2
    
17.     k=0
    
18.     k=Decimal(k)
    
19.     result=0
    
20.     num1=2*(Decimal(2)**Decimal(0.5))/9801
    
21. 
    
22.     while k!=Decimal(frequency):
    
23.         num2=Decimal(ft(4*k))*(Decimal(1103)+Decimal(26390)*k)
    
24.         num3=Decimal((ft(k)))**Decimal(4)*(Decimal(396)**(Decimal(4)*k))
    
25.         num4=num2/num3
    
26.         result+=num4
    
27.         k+=Decimal(1)
    
28. 
    
29.     ending=Decimal(1)/(Decimal(result)*num1)
    
30.     return ending
    
31. 
    
32. def improvement(frequency):
    
33.     t1=time.time()
    
34.     getcontext().prec=frequency*15+2
    
35.     q=str(time.time()-t1)
    
36.     print(q)
    
37.     time.sleep(0.5)
    
38.     k=0
    
39.     k=Decimal(k)
    
40.     result=0
    
41.     t2=time.time()
    
42.     kf=Decimal(10005)**Decimal(0.5)
    
43.     w=str(time.time()-t2)
    
44.     print(w)
    
45.     time.sleep(0.5)
    
46.     t3=time.time()
    
47.     num1=Decimal(426880)*kf
    
48.     e=str(time.time()-t3)
    
49.     print(e)
    
50.     time.sleep(0.2)
    
51.     t=time.time()
    
52.     while k!=Decimal(frequency):
    
53.         lt=time.time()
    
54.         num2=Decimal(ft(6*k))*(Decimal(13591409)+Decimal(545140134)*k)
    
55.         num3=Decimal(ft(3*k))*Decimal(ft(k)**3)*(Decimal(-640320)**(3*k))
    
56.         result+=num2/num3
    
57.         print('次数：'+str(k+1)+'位数：'+str((k+1)*15))
    
58.         k+=1
    
59.     ending=num1/Decimal(result)
    
60.     return ending
    
61. 
    
62. x=input('1：限定位数，2：自动增加')
    
63. if x=='1':
    
64.     c=input('1, 2')
    
65.     if c !='1' and c !='2':
    
66.         c='3'
    
67. 
    
68.     while True:
    
69.         print('-'*80)
    
70.         if c=='1':
    
71.             print('1')
    
72.         elif c=='2':
    
73.             print('2')
    
74.         else:
    
75.             print('3')
    
76. 
    
77.         d=int(input('小数点后几位'))
    
78. 
    
79.         while True:
    
80.             if d>5:
    
81.                 break
    
82.             elif d<=0:
    
83.                 break
    
84.             elif d>0 and d<=5:
    
85.                 if c=='1':
    
86.                     print(rj(int(d)))
    
87.                     break
    
88.                 elif c=='2':
    
89.                     print(improvement(int(d)))
    
90.                     break
    
91.                 else:
    
92.                     print("输入错误，请重新输入")
    
93.                     break

复制代码

这样，用户可以选择是限定位数计算还是自动增加计算，并可以输入要计算的位数。最后，根据用户的选择，调用不同的函数来计算圆周率的近似值。

以上回复来自 FishC AI 小助理，如未能正确解答您的问题，请继续追问。