1000 cifre di pi in python

Question

Ho pensato a questo problema e non riesco a capirlo. Forse puoi aiutarmi. Il problema è che il mio codice non funziona per produrre 1000 cifre di pi nel linguaggio di codifica python.

Ecco il mio codice:

def make_pi(): 
    q, r, t, k, m, x = 1, 0, 1, 1, 3, 3 
    while True: 
     if 4 * q + r - t < m * t: 
      yield m 
      q, r, t, k, m, x = (10*q, 10*(r-m*t), t, k, (10*(3*q+r))//t - 10*m, x) 
     else: 
      q, r, t, k, m, x = (q*k, (2*q+r)*x, t*x, k+1, (q*(7*k+2)+r*x)//(t*x), x+2) 

digits = make_pi() 
pi_list = [] 
my_array = [] 
for i in range(1000): 
    my_array.append(str("hello, I'm an element in an array \n")) 
big_string = "".join(my_array) 

print "here is a big string:\n %s" % big_string 

So che questo codice può essere fissato al lavoro, ma non sono sicuro cosa sistemare ... L'istruzione print dicendo qui è una grande corda e la my_array.append(str("hello, im an element in an array \n)) è solo un riempitivo per ora. So come tutto il codice è usato per funzionare, ma come ho detto prima, non riesco a farlo sparare quel codice.

Answer 1

Esegui questo

def make_pi(): 
    q, r, t, k, m, x = 1, 0, 1, 1, 3, 3 
    for j in range(1000): 
     if 4 * q + r - t < m * t: 
      yield m 
      q, r, t, k, m, x = 10*q, 10*(r-m*t), t, k, (10*(3*q+r))//t - 10*m, x 
     else: 
      q, r, t, k, m, x = q*k, (2*q+r)*x, t*x, k+1, (q*(7*k+2)+r*x)//(t*x), x+2 


digits = make_pi() 
pi_list = [] 
my_array = [] 

for i in make_pi(): 
    my_array.append(str(i)) 

my_array = my_array[:1] + ['.'] + my_array[1:] 
big_string = "".join(my_array) 
print "here is a big string:\n %s" % big_string 

e leggere su yield operatore da qui: What does the "yield" keyword do in Python?

ecco la risposta:

3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679821480865132823066470938446095505822317253594081284811174502841027019385211055596446229489549303819644288109756659334461284756482337 

Answer 2

Non ho familiarità con il tuo algoritmo. È un'implementazione di BBP?

In ogni caso, il tuo make_pi è un generatore. Provare a utilizzare in un ciclo for:

for digit in make_pi(): 
    print digit 

Si noti che questo ciclo è infinito: make_pi() mai tiri StopIteration

Answer 3

Questo fare quello che vuoi?

i = 0; 
pi_str = "" 
for x in make_pi(): 
    pi_str += str(x) 
    i += 1 
    if i == 1001: 
     break 

print "pi= %s.%s" % (pi_str[0],pi_str[1:]) 

Answer 4

Se non si desidera implementare il tuo algoritmo, puoi usare mpmath.

try: 
    # import version included with old SymPy 
    from sympy.mpmath import mp 
except ImportError: 
    # import newer version 
    from mpmath import mp 
mp.dps = 1000 # set number of digits 
print(mp.pi) # print pi to a thousand places 

Reference

Aggiornamento:. Codice supporta installazioni più vecchie e più recenti del SymPy (vedi comment) *

Answer 5

Qui è possibile controllare se i tuoi risultati del programma corrette 1000 digit: http://spoj.com/CONSTANT

Naturalmente puoi anche usare diff o tc, ma dovresti copiare questi 1000 cifre da qualche parte e lì basta inviare il programma e verificare se il punteggio è più grande di 999.

Si può provare a stampare ancora più cifre lì e quindi ottenere più punti. Forse ti piacerebbe.

Answer 6

Ecco un modo diverso che ho trovato qui ->Python pi calculation? per approssimare pitone sulla base della formula fratelli Chudnovsky per la generazione di Pi che ho sightly modificato per il mio programma.

def pifunction(): 
    numberofdigits = int(input("please enter the number of digits of pi that you want to generate")) 
    getcontext().prec = numberofdigits 

def calc(n): 
    t = Decimal(0) 
    pi = Decimal(0) 
    deno = Decimal(0) 
    k = 0 
    for k in range(n): 
     t = (Decimal(-1)**k)*(math.factorial(Decimal(6)*k))*(13591409+545140134*k) 
     deno = math.factorial(3*k)*(math.factorial(k)**Decimal(3))*(640320**(3*k)) 
     pi += Decimal(t)/Decimal(deno) 
    pi = pi * Decimal(12)/Decimal(640320**Decimal(1.5)) 
    pi = 1/pi 
    return str(pi) 
print(calc(1)) 

Spero che questo aiuti in quanto è possibile generare un numero qualsiasi di cifre di pi che si desidera generare.

Answer 7

la formula di wallis può arrivare a 3.141592661439964 ma è necessario un modo più efficiente per risolvere questo problema.

https://www.youtube.com/watch?v=EZSiQv_G9HM

and now my code 
x = 2 
y = 3 
summing = 4 
for count in range (0,100000000): 
    summing *= (x/y) 
    x += 2 
    summing *= (x/y) 
    y += 2 
print (summing)