New Prime Unless Someone Else found it:)

Stangerzv · 2013-06-07 07:23:18

Consider this equation.

Where, n is an integer, Pt is a prime number and Ps is the resulting Prime.

For n=1

Ps=5, 7, 11, 23, 47, 59, 83, 107, 167, 179, 227, 263, 347, 359, 383 for Pt<200 A well known prime (http://oeis.org/A005385)

For n=2

A well known prime (http://oeis.org/A155841)

For n=3

Ps=17, 263 for Pt<200

For n=4

Ps=79, 1049, 17179869269, 30354201441027016733116592294117482916287606860189680019559568902170379456331383469 for Pt<200

For n=5

Ps=37, 78167, 11920928955078263, 186264514923095703299, 4656612873077392578311 for Pt<200

For n=6

Ps=789730223053602977, 293242067884135544935936513642647623193965101483 for Pt<200

For n=7

Ps=367 for Pt<200

For n=8

Ps=144115188075856043, 154742504910672534362390789, 12259964326927110866866776217202473468949912977468817957, 363419362147803445274 661903944002267176820680343659030140745099590319644056698961663095525356881782780381260803133088966767300814308669 for Pt<200

For n=9

Ps=101, 4783039, 2541865828459 for Pt<200

For n=10

Ps=1033, 100000000000000000000253, 100000000000000000000000000319, 100000000000000000000000000000000000000000000000000000000000000000
000000000000000000000000000000000000000000001243 for Pt<200

Last edited by Stangerzv (2013-06-07 12:48:17)

Stangerzv · 2013-06-07 12:59:34

Some values for larger n:

For n=100000000

Ps={1000000000000000300000003}, {10000000000000000000000000000000000000000000000000000000000000000
000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000002300000023} for Pt<200

phrontister · 2013-06-08 07:02:41

Hi Stangerzv,

For n=156000000000000000000000

At Pt=197, Ps=a 4570-digit prime:

"PrimeQ[ps]" in Mathematica says the number is prime.

Do you have any means of testing this number for primality? I tried the online Elliptic Curve Method at http://www.alpertron.com.ar/ECM.HTM, but gave up because it was taking forever.

Stangerzv · 2013-06-08 08:07:04

Hi phrontister

I would try to run it on alpertron and maybe after 1-2 days I would get the answer.

phrontister · 2013-06-08 12:48:08

Hi Stangerzv,

Overnight I found a larger prime than in my last post, this time for Pt=199:

For n=2340240773870584944725754

At Pt=199, Ps=a 4919-digit prime:

"PrimeQ[ps]" in Mathematica says the number is prime.

I'll try the alpertron test on it. Pity it's so slow, though...their Java app time doesn't compare at all with Mathematica's PrimeQ function, which only takes 20 seconds to determine that Ps is prime.

Last edited by phrontister (2013-06-08 12:51:29)

Stangerzv · 2013-06-08 13:42:13

Hi phrontister

There is a primality software (open source) used by gimps to find the largest mersenne's prime. There is USD150,000 award for those who could find the 100,000,000 millions digit prime by electronic frontier. Since you are good at programming why not you try this software. I am looking for to find non-mersenne's prime. So far all biggest primes are mersenne numbers.

You can read about the software here http://en.wikipedia.org/wiki/Prime95

Last edited by Stangerzv (2013-06-08 13:42:55)

Stangerzv · 2013-06-09 03:58:28

After 20h 14 mins 37s

{1110 569841 961736 521093 527020 321061 997637 323487 171183 916164 743035
890337 478305 626254 080227 459404 322209 157812 216425 509264 920965 225954
216452 814660 943212 732459 541390 782750 026232 133733 246934 155146 590841
811410 546270 679392 624264 117413 197437 484536 121492 118007 586229 922433
542359 267360 140967 814731 499057 585417 184901 893328 122999 044651 625494
638152 554456 608373 144186 988936 455348 214706 248421 696806 985296 814269
741497 769100 215242 886224 200235 468455 936000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000000 000000
000000 000000 000000 000000 000000 000000 000000 000000 000000 000030 732000
000000 000000 000197 is prime }

Last edited by Stangerzv (2013-06-09 04:32:05)

Stangerzv · 2013-06-09 04:14:51

Hi phrontister

I think mathematica's PrimeQ function is the same with alpertron.

Stangerzv · 2013-06-09 15:03:01

This one also a prime

phrontister · 2013-06-10 00:06:16

Yes, I saw that from M's solution but I only posted the largest Pt<200 for n=156000000000000000000000.

Alpertron is still busy thinking (and has been for 35 hours) about the 4919-digit prime from my post #5. If that number is prime it will be quite a while yet for it to find that out, as it's only up to 22% completion.

phrontister · 2013-06-10 00:54:07

Stangerzv wrote:

There is a primality software (open source) used by gimps to find the largest mersenne's prime...Since you are good at programming why not you try this software. I am looking for to find non-mersenne's prime. So far all biggest primes are mersenne numbers.
You can read about the software here http://en.wikipedia.org/wiki/Prime95

Thanks for that, Stangerzv, but I only have a passing interest in primes. In fact, I don't even know what a Mersenne prime is, and although I've heard the term I've never been curious enough about it to find out what it is.

I've just enjoyed tackling the challenges you've posted because of the puzzle-solving aspect and the practical opportunity they give me to learn some more about Mathematica.

Stangerzv · 2013-06-10 02:01:47

I would run it on my computer as mine is core i5. I am learning mathematica/matlab/mathcad while looking for a solution too. Regarding Mersenne, he was a bishop not even a mathematician but his work in mathematics was/is known by many people and Mersenne's prime is the most famous and there are group of people working everyday on grid computing to get the latest largest prime (Mersenne's prime 2^p-1). They have more than 700,000 cpus in their grid computing group linking by online network. Finding the largest prime having more than 20 millions digits on my computer maybe would take more than 50-100 years:)

You can read about Mersenne Prime Search group here http://www.mersenne.org/

phrontister · 2013-06-10 02:15:24

Thanks for the Mersenne info and link...I'll check it out.

Alpertron is now up to 23% after 37.4 hours. I'll keep it running.

phrontister · 2013-06-10 03:38:03

Hi Stangerzv,

I'll keep it running.

I've changed my mind and have put it out of its misery (it was still showing 23%!!), as I need all my computer's resources for rendering a video.

I'll just accept Mathematica's word for it that the 4919-digit prime from my post #5 is prime, as its assessment of the 4570-digit prime you tested in alpertron is accurate.

Last edited by phrontister (2013-06-10 03:49:27)

Stangerzv · 2013-06-10 04:44:34

Hi phrontister

Don't worry I would get the result by tomorrow, now it is 11% after 2hours 41mins.

Stangerzv · 2013-06-12 00:20:27

After 22h 1m 53s

{30 398195 986259 526179 302633 268366 341776 895486 915223 261429 574719
664291 540660 885243 535785 201913 587451 484005 348788 275448 583272 837731
792151 379645 657141 289600 115327 481119 099892 899538 205041 658220 871131
340323 267671 045321 059211 448641 665276 423953 359467 071362 517048 972796
783208 151823 033442 569886 693949 797990 780274 538453 391806 301941 108378
980830 100565 371775 006873 093696 868688 016144 255230 898989 809567 050225
323275 607703 687762 522940 558683 847549 109358 444056 315584 569759 894895
853117 791274 592614 104409 709799 600529 731948 419730 544427 679776 748268
520895 415201 403259 373393 153404 826778 407181 864405 394786 619383 497885
310639 267898 281896 648296 344821 942234 374165 992885 817128 361179 945381
633928 341963 008409 672119 238504 813798 748701 731076 524233 995979 037135
112013 125089 076959 564708 774147 064248 510993 307813 545556 125029 063630
979454 419686 668368 980592 341743 713160 789823 477957 223366 959589 127709
726671 285337 561116 082429 436247 887131 802447 128809 437369 699808 236416
869486 011279 516397 879073 978825 828138 275898 121489 545055 153789 495298
287782 609756 668631 885428 887207 234165 953658 525093 608212 872829 817616
897872 825591 133767 686913 805120 000233 349286 160789 671527 704607 300416
169569 144323 105042 722426 422443 812131 001906 946132 244641 402701 017848
279768 066740 446802 951019 403917 681739 034847 620165 718209 525443 343139
990870 233188 838990 636652 079277 747702 635000 175449 494652 945197 943767
800497 721393 243171 942432 510181 204627 959054 426470 046416 189604 755487
709391 175160 202859 354352 214522 663126 515831 156580 635237 816371 163629
650321 830465 779640 038989 294097 778448 983549 727413 438363 623539 370582
701096 997246 992442 856422 252514 631447 199461 167972 699692 981088 802617
305892 999886 424840 830016 850395 159985 210600 576785 246050 555630 226881
712346 253724 229575 965637 179055 416760 125602 635786 564289 272398 311639
028493 706691 689470 093804 171563 621093 318073 508780 679537 427403 712935
093488 480149 502882 082750 299299 725620 143238 273604 513362 005393 974880
929832 950837 286295 307803 731414 081797 117488 488706 987012 363659 584698
739201 491322 776827 131638 691915 759098 574221 662323 775660 641820 337557
686123 062467 230719 267108 461294 830515 256708 547936 649101 436256 916287
080022 020968 131833 980438 505859 535946 297895 549573 937851 837803 655367
697435 471373 171128 921516 034700 453070 373119 269987 600813 348147 747593
283219 410435 013102 213184 321847 764294 578982 939482 044762 240513 284360
623095 443017 092684 298495 473181 118909 854583 257084 018997 763005 405168
013211 908698 960659 621483 771757 189548 366850 247144 152532 922212 552972
387618 341615 708110 771317 189713 480846 699806 630387 874550 725610 040679
303451 761912 175720 986272 701549 952653 700735 209388 322018 064817 675310
571550 429833 111120 270042 713505 180791 801165 545179 105960 209966 187825
384551 764250 202712 137210 086553 788264 072549 979957 134099 387576 751785
656961 618444 978687 644963 282700 236760 268530 813604 637432 889435 141910
107761 217097 745618 737802 251574 941169 702762 672146 191375 554781 718798
719500 549245 167103 818228 263089 484305 073753 513235 067911 978072 959448
440133 927168 244613 389571 779876 216366 883576 761521 484438 644774 245631
549366 898853 586859 668221 910698 051419 941601 662416 303895 250089 832984
955670 326840 927157 666612 650040 422830 592515 448006 983014 518098 753757
104947 828737 750364 746282 527567 869064 402576 740405 757797 696596 963260
802508 214432 038301 422340 567553 576111 543089 820519 951849 274054 626187
460049 311505 246714 739190 280672 270949 744017 044845 289032 834946 719634
390567 982358 791849 976707 376730 687178 297852 371555 749822 644519 085644
277283 055547 762841 669186 454173 175865 276058 277318 449344 985025 615384
291174 386504 196879 071112 929462 190389 321942 856491 261908 453971 618585
435754 337361 642244 843513 252322 509462 561211 765564 700485 447676 627376
896098 654723 181064 453043 806826 307951 985131 419444 853531 689850 236956
825982 762512 269374 066056 601454 401163 993905 448595 539906 335674 564507
937127 693051 800933 695191 662372 348054 433853 883236 436485 176331 476912
446036 921791 363131 021340 816378 316515 046696 660063 019154 345408 452124
187344 947776 153698 261613 730692 792790 661817 052457 928552 089672 111169
147298 314157 842582 625250 408112 480450 717761 715699 140084 323461 830034
496611 046620 493693 910232 931337 516751 995329 131684 082187 241259 074941
022802 224708 954193 173782 411179 962034 379634 641463 350412 717873 680926
597003 889443 360565 775542 201391 517773 334467 564744 023543 161321 857218
864656 613333 518105 194103 940791 094326 276032 341799 194121 875486 989244
082751 966811 821503 325974 021431 513664 467835 214978 413081 402109 595899
245138 514624 660938 781678 082931 249557 706993 240863 820581 029575 720513
198614 561137 615471 027973 021007 753543 721433 248088 967933 511039 898937
868991 798669 539796 986371 473504 160946 531543 950916 111373 827191 505785
252257 254867 433843 702102 261556 841244 568009 610774 554275 584364 251770
447599 303535 487158 565497 429864 491429 806394 934395 523682 076989 422612
896413 148657 746269 349931 009912 648755 694640 950868 726110 354183 805200
411322 132619 458254 964666 590219 477107 536034 919657 928757 213804 672352
904141 085034 826411 761525 903750 141990 018812 772382 964912 889304 998920
175919 700123 231270 291957 937589 534003 059920 375834 576770 532782 791000
934283 410807 996100 667964 338757 346589 is prime }

bobbym · 2013-06-12 00:28:48

Hi;

Yes, that checks out as prime. Congratulations!

phrontister · 2013-06-12 01:07:22

Hi Stangerzv,

Thanks for that...and I got my video done. Nice to have that large prime confirmed by a program that seems to be well regarded in the prime community.

After 22h 1m 53s...

Mathematica took just 21.344 seconds.

bobbym · 2013-06-12 01:09:33

You could probably beat that too with an AKS algorithm. Using a webpage has got to be slower.

phrontister · 2013-06-12 01:54:35

Hi Bobby,

I did a quick research on AKS but now know less about it than before I started...far too many multi-syllable words for me.

I found this question and two responses, though (for whatever they're worth):-

Question:

A few years ago, it was proven that PRIMES is in P. Are there any algorithms implementing their primality test in Python? I wanted to run some benchmarks with a naive generator and see for myself how fast it is. I'd implement it myself, but I don't understand the paper enough yet to do that.

Response #1:

AKS is of great theoretical importance but its performance is terrible (Miller-Rabin is much better).

Response #2

Quick answer: no, the AKS test is not the fastest way to test primality. There are much much faster primality tests that either assume the (generalized) Riemann hypothesis and/or are randomized. (E.g. Miller-Rabin is fast and simple to implement.) The real breakthrough of the paper was theoretical, proving that a deterministic polynomial-time algorithm exists for testing primality, without assuming the GRH or other unproved conjectures.

bobbym wrote:

Using a webpage has got to be slower.

The alpertron website says, "Since all calculations are performed in your computer, you can disconnect it from the Internet while the factorization is in progress."

I know you know I haven't got much idea at all about any of this, so please forgive me if I've taken a wrong turn.

bobbym · 2013-06-12 02:01:18

Hi;

The Miller Rabin is a probabilistic test and as such is not certain. But if I tell you that this number is prime with only a chance of .000000000000000001% of not being a prime, then you can feel very safe about it. They are right about AKS, it takes about an hour to prove the primality of a 1000 digit number on an old desktop.

I know about Dario's site, I use the other features and stick to Miller Rabin for primality. About the GRH? That is very difficult to understand.

Math Is Fun Forum

#1 2013-06-07 07:23:18

New Prime Unless Someone Else found it:)

#2 2013-06-07 12:59:34

Re: New Prime Unless Someone Else found it:)

#3 2013-06-08 07:02:41

Re: New Prime Unless Someone Else found it:)

#4 2013-06-08 08:07:04

Re: New Prime Unless Someone Else found it:)

#5 2013-06-08 12:48:08

Re: New Prime Unless Someone Else found it:)

#6 2013-06-08 13:42:13

Re: New Prime Unless Someone Else found it:)

#7 2013-06-09 03:58:28

Re: New Prime Unless Someone Else found it:)

#8 2013-06-09 04:14:51

Re: New Prime Unless Someone Else found it:)

#9 2013-06-09 15:03:01

Re: New Prime Unless Someone Else found it:)

#10 2013-06-10 00:06:16

Re: New Prime Unless Someone Else found it:)

#11 2013-06-10 00:54:07

Re: New Prime Unless Someone Else found it:)

#12 2013-06-10 02:01:47

Re: New Prime Unless Someone Else found it:)

#13 2013-06-10 02:15:24

Re: New Prime Unless Someone Else found it:)

#14 2013-06-10 03:38:03

Re: New Prime Unless Someone Else found it:)

#15 2013-06-10 04:44:34

Re: New Prime Unless Someone Else found it:)

#16 2013-06-12 00:20:27

Re: New Prime Unless Someone Else found it:)

#17 2013-06-12 00:28:48

Re: New Prime Unless Someone Else found it:)

#18 2013-06-12 01:07:22

Re: New Prime Unless Someone Else found it:)

#19 2013-06-12 01:09:33

Re: New Prime Unless Someone Else found it:)

#20 2013-06-12 01:54:35

Re: New Prime Unless Someone Else found it:)

#21 2013-06-12 02:01:18

Re: New Prime Unless Someone Else found it:)

Board footer