logowanie

matematyka » forum » matematyka » temat

Liczby pierwsze na Z艂otej Spirali.

ostatnie wiadomo艣ci  |  regulamin  |  latex

AutorWiadomo艣膰

sylvi91
post贸w: 23
2019-07-02 04:24:21

Cze艣膰.
Zastanowi艂o mnie po艂膮czenie ci膮gu Fibonacciego ze zbiorem liczb pierwszych.
Uzna艂em, 偶e przeprowadz臋 eksperyment i wyrysuj臋 programowo spiral臋, a nast臋pnie nanios臋 liczby Fibonacciego oraz liczby Pierwsze.
Gdyby by艂a ju偶 taka spirala komu艣 znana to ca艂kiem dobrze si臋 sk艂ada, bo to co ja zauwa偶am mo偶e jest wiadome, ale nie by艂o dla mnie jeszcze przed eksperymentem.

Ot贸偶 liczby pierwsze na spirali wyst臋puj膮 w bardzo losowy spos贸b, to jest pewne.
Przybywa ich pomiedzy kolejnymi wyrazami ci膮gu Fibonacciego, to jest pewne.

Z obserwacji wynika te偶 co艣 takiego, 偶e:
Suma liczb pierwszych na spirali pomi臋dzy kolejnymi wyrazami ci膮gu Fibonacciego podzielona przez ilo艣膰 liczb pierwszych, a nast臋pnie podzielona przez wi臋ksz膮 liczb臋 Fibonacciego daje w przybli偶eniu po艂ow臋 Z艂otej liczby Phi.
Wyj膮tkiem jest pocz膮tek spirali, gdzie warto艣ci s膮 r贸wne 1.
Czyli zbi贸r jest w pewnym sensie zwi膮zany ze z艂ot膮 spiral膮 w ten spos贸b. Czy mo偶e to zbyt naci膮gany wniosek?

Prosz臋 zark膮膰 na obrazek okna aplikacji, w kt贸rej rysowa艂em spiral臋 i nanosi艂em liczby pierwsze.
Zrobi艂em kilka zrzut贸w.


1:

2:

3:



Na tych zrzutach ekranu powinno by膰 wida膰 spiral臋. Liczby pierwsze w zasadzie wida膰 tylko na pierwszym zrzucie. Dalej liczby zosta艂y wy艂膮czone, a jedynie ciut wieksze kropeczki od linii spirali s膮 troch臋 widoczne. Ale mamy dane o tych liczbach.
Istotn膮 rzecz膮 jest warto艣膰 przy zmiennej AVG/F, kt贸ra pokazuje 艣redni膮 liczb pierwszych podzielon膮 przez liczb臋 Fibonacciego, tak偶e ilo艣膰 liczb pierwszych QTY w przedziale pomi臋dzy liczbami Fibonacciego, czyli w polu danego kwadratu.

Gdyby kto艣 by艂 zainteresowany zabaw膮 z aplikacj膮, kt贸ra rysuje ta spiral臋, to mog臋 udost臋pni膰 w nast臋pnym wpisie. Jest to miniaturowy programik, ale troch臋 wymagaj膮cy pod wzgl臋dem mocy sprz臋towej. Mo偶e dzia艂a膰 zar贸wno na Windows jak i Linux.

To na razie tyle. Dzi臋ki za dotrwanie do ko艅ca i pozdrawiam.




sylvi91
post贸w: 23
2019-07-04 20:11:07


Za艂o偶enie z powy偶szej wiadomo艣ci zacz膮艂em sprawdza膰 dla wi臋kszych warto艣ci.
Mam wyprowdzon膮 formu艂臋 dla liczb pierwszych. Jest ona zupe艂nie nowa dla mnie... i pewnie dla wielu innych pasjonat贸w teorii liczb. Przynajmniej ja wcze艣niej o takiej nie s艂ysza艂em.

Oto jej posta膰:

$SP=1/2*Phi$

$SP= (((P1+P2+...+Pn)) / QP)/GF$

Gdzie:

$SF < Pn <= GF$

$P1, P2, Pn$ - s膮 to liczby pierwsze
$SF$ - mniejsza warto艣膰 ci膮gu Fibonacciego (dolna granica zbioru)
$GF$ - wieksza warto艣膰 ci膮gu Fibonacciego (g贸rna granica zbioru)
$SP$ - suma liczb pierwszych w podzbiorze
$QP$ - ilo艣膰 liczb pierwszych w podzbiorze
$Phi$ - Z艂ota Liczba Phi 1.618...

Przyk艂adowy listing oblicze艅 z programu do kalkulacji na liczbach pierwszych i Fibonacciego.

________________________________________________________________________
____________THAT APP IS CALCULATING THE PRIME NUMBERS SUBSETS___________
____________BETWEEN FIBONACCI AND FINDING APPROXIMATION_________________
____________TO THE GOLDEN NUMBER_PHI.___________________________________
_____________Author: Sylwester B aka Sylvi91____________________________
________________________________________________________________________
Fib(1) = 1 - Fib(2) = 1
Primes list:

Primes: Sum = 0 Qty = 0 Avg = 0.000000
Approximation to the The Golden Number Phi = 0.000000
________________________________________________________________________
Fib(2) = 1 - Fib(3) = 2
Primes list:
2
Primes: Sum = 2 Qty = 1 Avg = 2.000000
Approximation to the The Golden Number Phi = 2.000000
________________________________________________________________________
Fib(3) = 2 - Fib(4) = 3
Primes list:
3
Primes: Sum = 3 Qty = 1 Avg = 3.000000
Approximation to the The Golden Number Phi = 2.000000
________________________________________________________________________
Fib(4) = 3 - Fib(5) = 5
Primes list:
5
Primes: Sum = 5 Qty = 1 Avg = 5.000000
Approximation to the The Golden Number Phi = 2.000000
________________________________________________________________________
Fib(5) = 5 - Fib(6) = 8
Primes list:
7
Primes: Sum = 7 Qty = 1 Avg = 7.000000
Approximation to the The Golden Number Phi = 1.750000
________________________________________________________________________
Fib(6) = 8 - Fib(7) = 13
Primes list:
11 13
Primes: Sum = 24 Qty = 2 Avg = 12.000000
Approximation to the The Golden Number Phi = 1.846154
________________________________________________________________________
Fib(7) = 13 - Fib(8) = 21
Primes list:
17 19
Primes: Sum = 36 Qty = 2 Avg = 18.000000
Approximation to the The Golden Number Phi = 1.714286
________________________________________________________________________
Fib(8) = 21 - Fib(9) = 34
Primes list:
23 29 31
Primes: Sum = 83 Qty = 3 Avg = 27.000000
Approximation to the The Golden Number Phi = 1.588235
________________________________________________________________________
Fib(9) = 34 - Fib(10) = 55
Primes list:
37 41 43 47 53
Primes: Sum = 221 Qty = 5 Avg = 44.000000
Approximation to the The Golden Number Phi = 1.600000
________________________________________________________________________
Fib(10) = 55 - Fib(11) = 89
Primes list:
59 61 67 71 73 79 83 89
Primes: Sum = 582 Qty = 8 Avg = 72.000000
Approximation to the The Golden Number Phi = 1.617978
________________________________________________________________________
Fib(11) = 89 - Fib(12) = 144
Primes list:
97 101 103 107 109 113 127 131 137 139
Primes: Sum = 1164 Qty = 10 Avg = 116.000000
Approximation to the The Golden Number Phi = 1.611111
________________________________________________________________________
Fib(12) = 144 - Fib(13) = 233
Primes list:
149 151 157 163 167 173 179 181 191 193 197 199 211 223 227 229 233
Primes: Sum = 3223 Qty = 17 Avg = 189.000000
Approximation to the The Golden Number Phi = 1.622318
________________________________________________________________________
Fib(13) = 233 - Fib(14) = 377
Primes list:
239 241 251 257 263 269 271 277 281 283 293 307 311 313 317 331 337 347 349 353 359 367 373
Primes: Sum = 6989 Qty = 23 Avg = 303.000000
Approximation to the The Golden Number Phi = 1.607427
________________________________________________________________________
Fib(14) = 377 - Fib(15) = 610
Primes list:
379 383 389 397 401 409 419 421 431 433 439 443 449 457 461 463 467 479 487 491 499 503 509 521 523 541 547 557 563 569 571 577 587 593 599 601 607
Primes: Sum = 18165 Qty = 37 Avg = 490.000000
Approximation to the The Golden Number Phi = 1.606557
________________________________________________________________________

... i dalej bez listowania liczb pierwszych, a tylko z wyliczeniem sumy oraz ilo艣ci i 艣redniej i przyblizonej wartosci liczby Z艂otej.

Fib(15) = 610 - Fib(16) = 987
Primes: Sum = 43635 Qty = 55 Avg = 793.000000
Approximation to the The Golden Number Phi = 1.606890


Fib(16) = 987 - Fib(17) = 1597
Primes: Sum = 109567 Qty = 85 Avg = 1289.000000
Approximation to the The Golden Number Phi = 1.614277

Fib(19) = 4181 - Fib(20) = 6765
Primes: Sum = 1624217 Qty = 297 Avg = 5468.000000
Approximation to the The Golden Number Phi = 1.616556

Fib(27) = 196418 - Fib(28) = 31781
Primes: Sum = 2499980948 Qty = 9738 Avg = 256724.000000
Approximation to the The Golden Number Phi = 1.615577

....

Aktualnie ustawi艂em aplikacj臋 do wyci膮gania Z艂otej Liczby ze zbior贸w liczb pierwszych dla g贸rnego przedzia艂u do 44 liczby ci膮gu Fibonacciego, kt贸ra ma warto艣膰 701408733... komputer ju偶 oblicza kt贸r膮 godzin臋 i jeszcze pewnie sporo mu zosta艂o do ko艅ca.
Co do algorytmu to mo偶e nie jest najszybszy, bo stosuj臋 tablic臋 liczb pierwszych wyznaczon膮 metod膮 Sita Eratostenesa, a wi臋c staro偶ytny ju偶 algorytm... kt贸ry mimo to sprawdza si臋 w dobie komputer贸w.
Co do powi膮zania zbioru liczb pierwszych ze z艂ot膮 spiral膮 i z艂ot膮 liczb膮 to niewiele s艂ysza艂em do tej pory. A Wy?










masterkik
post贸w: 1
2020-05-16 20:08:40

Cze艣膰,

Szuka艂em troch臋 innych informacji, trafi艂em przypadkiem tu i postanowi艂em si臋 zarejestrowa膰 i odkopa膰 temat.

Zrobi艂em wizualizacj臋 (wykres excel)pierwszych ~65.5ty艣 liczb pierwszych na spirali r=e^fi (zapis biegunowy).
Po kr贸tkich poszukiwaniach uzna艂em 偶e mog艂em sobie r贸wnie dobrze te liczby zaznaczy膰 na osi liczbowej i porzuci艂em temat.

Ci膮g Fibonacciego to jeden (pewnie pierwszy historycznie) przyk艂ad z艂otych ci膮g贸w i ma swoje wady.
Do podobnych poszukiwa艅 raczej u偶ywa艂bym ci膮gu kolejnych pot臋g z艂otego Fi (kt贸ry spe艂nia swoj膮 drog膮 warunek Fn=Fn-1+Fn-2 a stosunek kolejnych wyraz贸w nie jest zbie偶ny do Fi tylko jest stale r贸wny Fi).

Dlaczego to mia艂by by膰 Fibonacci, a nie Lucca?
A mo偶e znajdziesz zasad臋 og贸ln膮. Powodzenia! :)


sylvi91
post贸w: 23
2020-10-13 18:09:36

呕adnego zwi膮zku liczb pierwszych z liczbami Lukasa i kolejnymi pot臋gami liczby Fi nie znalaz艂em jak dot膮d.

Natomiast podzbiory liczb pierwszych, kt贸rych granicami s膮 kolejne liczby Fibonacciego maj膮 taki zwi膮zek opisany wy偶ej.

Lub te偶 prawdopodobnie tak, je艣li nie pomyli艂em oznacze艅 Latexa.
$
\displaystyle{ SUMA = \sum_{i=1}^{n} a_{i} \approx \frac12 \Phi \approx \frac{ \frac{a_{1} + a_{2} + ... + a_{n}}{n}}{ fib(x) }}
$
$
\displaystyle{ fib(x-1) < a_{n} \le fib(x)}
$

Oto listing programu w C do obliczania tej zale偶no艣ci.
Aplikacja dzia艂a na 64 bitowych zmiennych.

#include <stdio.h>
#include <time.h>
#include <stdlib.h>
#include <string.h>


#define Phi 1.618033988749895

#define BITS 8
#define PIECES 7

unsigned long long fib(int n) {
if ((n == 1) || (n == 2))
return 1;
else
return fib(n - 1) + fib(n - 2);
}



// sprawdz czy x jest pierwsza czy nie
unsigned long long ifnotPrime(char prime[], unsigned long long x)
{

return (prime[x/BITS] & (1 << ((x >> 1) & PIECES) ) );

}




// zaznacz odpowiedni bit wedlug indexu
unsigned long long makeComposite(char prime[], unsigned long long x)
{
prime[x/BITS] |= (1 << ((x >> 1) & PIECES)); /// ok for test fast ok

return 0;

}




int main( int argc, const char* argv[] )
{

time_t start;

time_t end; // Get the system time


int x;
int y;
unsigned long long i,j,k;

unsigned long long sum_p=0; /// suma liczb pierwszych
unsigned long long qty_p=0; ///qty liczb pierwszych w podzbiorze
unsigned long long limit; /// limit of primes from fib

unsigned long long subset_start=0; /// konieczne do w艂a艣ciwej iteracji

unsigned long long subset_end=1;
unsigned long long tempNumber=0;


double approx_phi;
double avg_p=0;

qty_p=0;

if ( argc == 3 )
{
x = atoi( argv[1] );
y = atoi( argv[2] );


time(&start); // Get the system time





if ((y>x) && (x>1))


{




printf(\"________________________________________________________________________\n\");

printf(\"____________THAT APP IS CALCULATING THE PRIME NUMBERS SUBSETS___________\n\");

printf(\"____________BETWEEN FIBONACCI AND FINDING APPROXIMATION_________________\n\");

printf(\"____________TO THE GOLDEN NUMBER_PHI.___________________________________\n\");

printf(\"_____________Author: Sylwester B aka Sylvi91____________________________\n\");

printf(\"________________________________________________________________________\n\");

printf(\" Please wait for primes generator. Warning! For second argument f2 greater than 50 this may take couple minutes. Bigger f2 consume more and more memory.\n\");





limit = fib(y); /// tylko raz oblicz wartosc fibo


char * prime;
prime = (char*) malloc(sizeof(char)*limit/BITS); // alokacja pamieci

if (prime==NULL) exit (1);

/// memset(prime, 0, sizeof(prime)); // zainicjiuj zerowymi wartosciam



// 2 jest ignorowane
// petla od 3
// Eratostoenes .



for (i = 3; i * i <= limit; i += 2) {


if (!ifnotPrime(prime, i))
for ( j = i * i, k = i << 1; j < limit; j += k)
makeComposite(prime, j);
}


// drukuj na ekranie jesli chcesz
// printf(\"2 \");
/*
// drukuj pozostale
for (int i = 3; i <= limit; i += 2)
if (!ifnotPrime(prime, i))
{
printf(\"%d \", i);
}
printf(\"\n\");
*/




subset_start=fib(x-1); /// konieczne do rozpocz臋cia w艂a艣ciwej iteracji ale tylko raz

subset_end=fib(x); ///w p臋tli ju偶 tylko szybciej poprzez dodawanie liczb




for (i=x;i<y;i++) /// g艂贸wna p臋tla
{

/// fibo iteracja
tempNumber = subset_start;

subset_start = subset_end;

subset_end = tempNumber + subset_start;

/// zerowanie wartosci dla nowego zbioru
sum_p = 0;
qty_p = 0;
avg_p = 0;
approx_phi = 0;

/// wypisz
//printf(\"\n Fib(%llu) = %llu - Fib(%llu) = %llu \",i, subset_start,i+1,subset_end);

//printf(\"\n Primes list: \n\");



for (j = subset_start+1; j <= subset_end; j+=2) /// Liczy co drugi element aby omin膮膰 parzyste co daje 50% oszcz臋dno艣ci
{




if (j%2 == 0) /// je艣li j podzielne prze 2 to dodaj jeden

{
j++;
}


if (!ifnotPrime(prime, j)) /// prime czy nie

{



/// printf(\" %llu \",j); /// wypisz jak chcesz

sum_p += j;
qty_p += 1;



}





}

if (j>5)
{
avg_p = sum_p/qty_p;


///approx_phi = ((avg_p/subset_end) >> 1); /// przesuni臋cie bitowe zamiast mno偶enia odpoda dla double
approx_phi = ((avg_p/subset_end) *2);

printf(\" Fib(%llu) = %llu - Fib(%llu) = %llu SUM = %llu QTY = %llu AVG = %.2f Phi = %f\n\",i,subset_start,i+1, subset_end, sum_p, qty_p, avg_p, approx_phi);

}
else

{
printf(\" Fib(%llu) = %llu - Fib(%llu) = %llu SUM = %llu QTY = %u AVG = %.2f Phi = %f\n\",i,subset_start,i+1, subset_end, i,1, 1.0, 1.0);


}



/// wypisz jak chcesz
/*
printf(\"\n\");
printf(\" Primes: SUM = %llu QTY = %llu AVG = %.2f \n\", sum_p, qty_p, avg_p);



printf(\" Approximation to the The Golden Number Phi = %f \n\",approx_phi);
printf(\"________________________________________________________________________\n\");
*/

//




}


free(prime); // zwolnij pamiec

time(&end); // czas zakonczenia



double dif;




dif = difftime (end,start); // ustal roznice
printf (\"Your calculations took %.2lf seconds to run.\n\", dif ); // wypisz







}

else

printf(\" Simple usage: pn f1 f2 \n where f1 and f2 are Fibonacci words in sequence and f1 >= 1 and f2 > f1\n Warning! f2 should not exceed 50 otherwise you will reach memory limit.\n\");





} // if

else
printf(\" Simple usage: pn f1 f2 \n where f1 and f2 are Fibonacci words in sequence and f1 >= 1 and f2 > f1\n Warning! f2 should not exceed 50 otherwise you will reach memory limit.\n\");


return 0;

}

Wi臋cej wynik贸w.

Fib(30) = 832040 - Fib(31) = 1346269 SUM = 40227359343 QTY = 36981 AVG = 1087784 Approx. to Phi = 1.6159979915
Fib(31) = 1346269 - Fib(32) = 2178309 SUM = 101931466309 QTY = 57909 AVG = 1760200 Approx. to Phi = 1.6161159872
Fib(32) = 2178309 - Fib(33) = 3524578 SUM = 257870996074 QTY = 90550 AVG = 2847829 Approx. to Phi = 1.6159829631
Fib(33) = 3524578 - Fib(34) = 5702887 SUM = 654441400687 QTY = 142033 AVG = 4607671 Approx. to Phi = 1.6159082233
Fib(34) = 5702887 - Fib(35) = 9227465 SUM = 1661677489343 QTY = 222855 AVG = 7456316 Approx. to Phi = 1.6161136347
Fib(35) = 9227465 - Fib(36) = 14930352 SUM = 4221069024488 QTY = 349862 AVG = 12064954 Approx. to Phi = 1.6161647093
Fib(36) = 14930352 - Fib(37) = 24157817 SUM = 10735109882717 QTY = 549903 AVG = 19521824 Approx. to Phi = 1.6161910656
Fib(37) = 24157817 - Fib(38) = 39088169 SUM = 27324559743219 QTY = 865019 AVG = 31588392 Approx. to Phi = 1.6162635809
Fib(38) = 39088169 - Fib(39) = 63245986 SUM = 69592739215201 QTY = 1361581 AVG = 51111714 Approx. to Phi = 1.6162832531
Fib(39) = 63245986 - Fib(40) = 102334155 SUM = 177416424882449 QTY = 2145191 AVG = 82704255 Approx. to Phi = 1.6163568263
Fib(40) = 102334155 - Fib(41) = 165580141 SUM = 452491851513992 QTY = 3381318 AVG = 133821146 Approx. to Phi = 1.6163912555
Fib(41) = 165580141 - Fib(42) = 267914296 SUM = 1155101203883449 QTY = 5334509 AVG = 216533743 Approx. to Phi = 1.6164403784
Fib(42) = 267914296 - Fib(43) = 433494437 SUM = 2949939781199270 QTY = 8419528 AVG = 350368783 Approx. to Phi = 1.6164857174
Fib(43) = 433494437 - Fib(44) = 701408733 SUM = 7539299522749130 QTY = 13298630 AVG = 566923023 Approx. to Phi = 1.6165268447
Fib(44) = 701408733 - Fib(45) = 1134903170 SUM = 19277443065477372 QTY = 21014892 AVG = 917322966 Approx. to Phi = 1.6165660477
Fib(45) = 1134903170 - Fib(46) = 1836311903 SUM = 49319944945730044 QTY = 33227992 AVG = 1484289058 Approx. to Phi = 1.6165979816
Fib(46) = 1836311903 - Fib(47) = 2971215073 SUM = 126244882484469729 QTY = 52565409 AVG = 2401672219 Approx. to Phi = 1.6166263027
Fib(47) = 2971215073 - Fib(48) = 4807526976 SUM = 323315442948007805 QTY = 83198799 AVG = 3886059015 Approx. to Phi = 1.6166561454
Fib(48) = 4807526976 - Fib(49) = 7778742049 SUM = 828390738449096336 QTY = 131744274 AVG = 6287869015 Approx. to Phi = 1.6166801715
Fib(49) = 7778742049 - Fib(50) = 12586269025 SUM = 2123542722014291259 QTY = 208718785 AVG = 10174181121 Approx. to Phi = 1.6167112114
Fib(50) = 12586269025 - Fib(51) = 20365011074 SUM = 5445740631692717113 QTY = 330797447 AVG = 16462462697 Approx. to Phi = 1.6167398718
Fib(51) = 20365011074 - Fib(52) = 32951280099 SUM = 13971588279376661576 QTY = 524513152 AVG = 26637250612 Approx. to Phi = 1.6167657543
Fib(52) = 32951280099 - Fib(53) = 53316291173 SUM = 35859516892843327174 QTY = 831993816 AVG = 43100701235 Approx. to Phi = 1.6167929271
Fib(53) = 53316291173 - Fib(54) = 86267571272 SUM = 92072211900316816991 QTY = 1320232935 AVG = 69739369060 Approx. to Phi = 1.6168154043
Your calculations took 918.00 seconds to run.
64.5 % from 8GB memory usage.


strony: 1

Prawo do pisania przys艂uguje tylko zalogowanym u偶ytkownikom. Zaloguj si臋 lub zarejestruj

© 2019 Mariusz iwi駍ki      o serwisie | kontakt   drukuj