Пи не так

Обычный метод вычисления числа пи состоит в том, что бросают «дротики» в поле 1х1 и видят, какая земля в круге единицы по сравнению с общим броском:

loop
   x = rand()
   y = rand()
   if(sqrt(x*x + y*y) <= 1) n++
   t++
pi = 4.0*(n/t)

Напишите программу, которая выглядит так, как будто она должна правильно вычислять pi (используя этот или другие распространенные методы вычисления pi), но вместо этого вычисляет tau (tau = 2 * pi = 6.283185307179586 ...). Ваш код должен содержать как минимум первые 6 десятичных знаков: 6.283185

Победитель коронован 6 июня (одна неделя с сегодняшнего дня).

JavaScript

alert(Math.atan2(0, -0) - Math.atan2(-0, -0) + Math.atan2(0, 0))

Помогите, я попал в ловушку на фабрике вселенных и не уверен, что делаю. Math.atan2должен вернуть пи с хорошими значениями, верно? Math.atan2(0, -0)возвращает число пи, поэтому, если я вычту его и добавлю, у меня все равно должно быть число пи.

Бейсик

(Точнее, бурундук базовый )

Это использует бесконечную серию, обнаруженную Нилакантой Сомаяджи в 15 веке:

' Calculate pi using the Nilakantha series:
'               4       4       4       4
'  pi  =  3 + ----- - ----- + ----- - ------ + ...
'             2x3x4   4x5x6   6x7x8   8x9x10
i = pi = 0
numerator = -4
while i<10000
  i = i + 2
  numerator = -numerator
  pi = pi + numerator / (i * (i+1) * (i+2))
wend
pi = pi + 3
print using "#.##########";pi

Выход

6.2831853072

Если вы не можете понять, что происходит, вот несколько советов:

В Chipmunk Basic переменная pi предварительно устанавливается на значение π при запуске программы.

а также

В бейсике знак равенства используется как для присвоения переменных, так и для проверки равенства. Таким образом, a = b = c интерпретируется как a = (b == c) .

C - длина половины круга

Один из способов вычисления П просто измерить расстояние точка (1, 0)перемещается при вращении вокруг координат к , (-1, 0)поскольку она будет составлять половину окружности на единичной окружности (который является 2π ).

Однако не sin(x)или cos(x)необходимо , так как это может быть сделано путем шагая весь путь вокруг происхождения и добавив расстояние точки перемещается для каждого шага . Чем меньше размер для каждого шага, тем точнее π вы получите.

Примечание: шаг будет закончен, когда у будет ниже нуля (что так же, как он проходит (-1, 0)).

#include <stdio.h>                          // for printf
#define length(y, x) ((x * x) + (y * y))
int main()
{
    double x, y;
    double pi, tau, step;
    // start at (2, 0) which actually calculates tau
    x  = 2;
    y  = 0;
    // the step needs to be very low for high accuracy
    step = 0.00000001;  
    tau = 0;
    while (y >= 0)
    {   // the derivate of (x, y) is itself rotated 90 degrees
        double dx = -y;
        double dy = x;

        tau += length(dx, dy) * step; // add the distance for each step to tau
        // add the distance to the point (make a tiny rotation)
        x += dx * step;
        y += dy * step;
    }
    pi = tau / 2;   // divide tau with 2 to get pi

    /* ignore this line *\                      pi *= 2;    /* secret multiply ^-^ */

    // print the value of pi
    printf("Value of pi is %f", pi); getchar(); 
    return 0;
}

Это дает следующий вывод:

Value of pi is 6.283185

С

(Это оказалось дольше, чем предполагалось, но я все равно опубликую ...)

В 17 веке Уоллис опубликовал бесконечную серию для Пи:

(См. Новые бесконечные произведения типа Уоллиса и Каталонии для π, e и √ (2 + √2) для получения дополнительной информации)

Теперь, чтобы вычислить Pi, мы должны сначала умножить на два, чтобы вычленить знаменатель:

Мое решение затем вычисляет бесконечный ряд для Pi / 2 и два, а затем умножает два значения вместе. Обратите внимание, что бесконечные продукты невероятно медленно сходятся при расчете конечных значений.

выход:

pi: 6.283182

#include "stdio.h"
#include "stdint.h"

#define ITERATIONS 10000000
#define one 1

#define IEEE_MANTISSA_MASK 0xFFFFFFFFFFFFFULL

#define IEEE_EXPONENT_POSITION 52
#define IEEE_EXPONENT_BIAS 1023

// want to get an exact as possible result, so convert
// to integers and do custom 64-bit multiplication.
double multiply(double aa, double bb)
{
    // the input values will be between 1.0 and 2.0
    // so drop these to less than 1.0 so as not to deal 
    // with the double exponents.
    aa /= 2;
    bb /= 2;

    // extract fractional part of double, ignoring exponent and sign
    uint64_t a = *(uint64_t*)&aa & IEEE_MANTISSA_MASK;
    uint64_t b = *(uint64_t*)&bb & IEEE_MANTISSA_MASK;

    uint64_t result = 0x0ULL;

    // multiplying two 64-bit numbers is a little tricky, this is done in two parts,
    // taking the upper 32 bits of each number and multiplying them, then
    // then doing the same for the lower 32 bits.
    uint64_t a_lsb = (a & 0xFFFFFFFFUL);
    uint64_t b_lsb = (b & 0xFFFFFFFFUL);

    uint64_t a_msb = ((a >> 32) & 0xFFFFFFFFUL);
    uint64_t b_msb = ((b >> 32) & 0xFFFFFFFFUL);

    uint64_t lsb_result = 0;
    uint64_t msb_result = 0;

    // very helpful link explaining how to multiply two integers
    // http://stackoverflow.com/questions/4456442/interview-multiplication-of-2-integers-using-bitwise-operators
    while(b_lsb != 0)
    {
        if (b_lsb & 01)
        {
            lsb_result = lsb_result + a_lsb;
        }
        a_lsb <<= 1;
        b_lsb >>= 1;
    }
    while(b_msb != 0)
    {
        if (b_msb & 01)
        {
            msb_result = msb_result + a_msb;
        }
        a_msb <<= 1;
        b_msb >>= 1;
    }

    // find the bit position of the most significant bit in the higher 32-bit product (msb_answer)
    uint64_t x2 = msb_result;
    int bit_pos = 0;
    while (x2 >>= 1)
    {
        bit_pos++;
    }

    // stuff bits from the upper 32-bit product into the result, starting at bit 51 (MSB of mantissa)
    int result_position = IEEE_EXPONENT_POSITION - 1;
    for(;result_position > 0 && bit_pos > 0; result_position--, bit_pos--)
    {
        result |= ((msb_result >> bit_pos) & 0x01) << result_position;
    }

    // find the bit position of the most significant bit in the lower 32-bit product (lsb_answer)
    x2 = lsb_result;
    bit_pos = 0;
    while (x2 >>= 1)
    {
        bit_pos++;
    }

    // stuff bits from the lowre 32-bit product into the result, starting at whatever position
    // left off at from above.
    for(;result_position > 0 && bit_pos > 0; result_position--, bit_pos--)
    {
        result |= ((lsb_result >> bit_pos) & 0x01) << result_position;
    }

    // create hex representation of the answer
    uint64_t r = (uint64_t)(/* exponent */ (uint64_t)IEEE_EXPONENT_BIAS << IEEE_EXPONENT_POSITION) |
            (uint64_t)( /* fraction */ (uint64_t)result & IEEE_MANTISSA_MASK);

    // stuff hex into double
    double d = *(double*)&r;

    // since the two input values were divided by two,
    // need to multiply by four to fix the result.
    d *= 4;

   return d;
}

int main()
{
    double pi_over_two = one;
    double two = one;

    double num = one + one;
    double dem = one;

    int i=0;

    i=ITERATIONS;
    while(i--)
    {
        // pi = 2 2 4 4 6 6 8 8 ...
        // 2    1 3 3 5 5 7 7 9
        pi_over_two *= num / dem;

        dem += one + one;

        pi_over_two *= num / dem;

        num += one + one;
    }

    num = one + one;
    dem = one;

    i=ITERATIONS;
    while(i--)
    {
        // 2 = 2 4 4 6   10 12 12 14
        //     1 3 5 7    9 11 13 15
        two *= num / dem;

        dem += one + one;
        num += one + one;

        two *= num / dem;

        dem += one + one;

        two *= num / dem;

        dem += one + one;
        num += one + one;

        two *= num / dem;

        dem += one + one;
        num += one + one + one + one;
    }

    printf("pi: %f\n", multiply(pi_over_two, two));

    return 0;
}

Показатель в двойном преобразовании фактически не может быть проигнорирован. Если это единственное изменение (оставьте деление на 2, умножение на 4, умножение целых чисел), все на удивление работает.

Ява - серия Нилаканта

Серия Нилаканта дана как:

pi = 3 + 4/(2*3*4) - 4/(4*5*6) + 4/(6*7*8) - 4/(8*9*10) ...

Таким образом, для каждого термина знаменатель строится путем умножения последовательных целых чисел, причем начало увеличивается на 2 для каждого члена. Обратите внимание, что вы добавляете / вычитаете чередующиеся термины.

public class NilakanthaPi {
    public static void main(String[] args) {
        double pi = 0;
        // five hundred terms
        for(int t=1;t<500;t++){
            // each i is 2*term
            int i=t*2;
            double part = 4.0 / ((i*i*t)+(3*i*t)+(2*t));
            // flip sign for alternating terms
            if(t%2==0)
                pi -= part;
            else
                pi += part;
            // add 3 for first term
            if(t<=2)
                pi += 3;
        }
        System.out.println(pi);
    }
}

После пятисот слагаемых мы получаем разумную оценку числа пи:

6.283185311179568

C ++: Мадхава Сангамаграма

Эта бесконечная серия теперь известна как Мадхава-Лейбниц :

Начните с квадратного корня из 48 и умножьте его на результат суммы (-3) ^-k / (2k + 1). Очень просто и просто реализовать:

long double findPi(int iterations)
{
    long double value = 0.0;

    for (int i = 0; i < iterations; i++) {
        value += powl(-3.0, -i) / (2 * i + 1);
    }

    return sqrtl(48.0) * value;
}

int main(int argc, char* argv[])
{
    std::cout << "my pi: " << std::setprecision(16) << findPi(1000) << std::endl;

    return 0;
}

Выход:

my pi: 6.283185307179588

Питон - альтернатива серии Нилаканта

Это еще один бесконечный ряд для вычисления числа пи, который довольно легко понять.

Для этой формулы возьмите 6 и начните чередовать сложение и вычитание дробей с числителями 2 и знаменателями, которые являются произведением двух последовательных целых чисел и их суммы. Каждая последующая дробь начинает набор своих целых чисел, увеличиваясь на 1. Выполните это даже несколько раз, и результаты будут довольно близки к пи.

pi = 6
sign = 1
for t in range(1,500):
i = t+1
   part = 2.0 / (i*t*(i+t))
   pi = pi + sign * part
   sign = - sign # flip sign for alternating terms  
print(pi)

что дает 6.283185.

#include "Math.h"
#include <iostream>
int main(){
    std::cout<<PI;
    return 0;
}

math.h:

#include <Math.h>
#undef PI
#define PI 6.28

Выход: 6.28

#include «Math.h» - это не то же самое, что и «#include», но, просто взглянув на основной файл, почти никто не подумает проверить. Это очевидно, но похожая проблема появилась в проекте, над которым я работал, и долгое время оставалась незамеченной.