mersenneforum.org a new Deterministic primality testing
 Register FAQ Search Today's Posts Mark Forums Read 2012-12-07, 07:43   #23
LaurV
Romulan Interpreter

"name field"
Jun 2011
Thailand

100110010100002 Posts Quote:
 Originally Posted by LaurV It turned out there is no big deal to implement this. No idea if pari has one already.
Well, I just figured out that the a t_COMPLEX type in pari accepts coefficients of the t_MOD type. This voids the requirement for a function, as one can define "base=Mod(3,n)+Mod(2,n)*I" for a given n, then pari takes care by itself to do modular things with the coefficients. Although not so fast like my last implementation, it is however much faster then the initial version, and it would allow somebody to test by himself, using nothing except pari default functions and types (no "unknown" complex_pow_mod function).

So, one can modify the forstep parameters to go as high as he likes, eventually going 4 by 4 to be faster, and/or cancelling 5+12i test (I use the step 2 by 2 and testing both cases to show that the thing is working, you can also remove the primality test and it will print primes, or use ispseudoprime(n,1) to make it much faster, but here be careful that you can lose solutions if the random base for RM test is unlucky selected).

Code:
gp > nr=0; cnt=0; forstep(n=11,10^7,2,if(((z=(Mod(3,n)+Mod(2,n)*I)^(n+1))==13||z==5+12*I)&&!isprime(n),nr++; print(n" : "n%4"     ")); if(n
>cnt,cnt+=10^6; printf("...%d...%c",n,13))); nr
1105 : 1
2465 : 1
10585 : 1
29341 : 1
...<snip many lines deleted>...
9131401 : 1
9638785 : 1
9890881 : 1
%6 = 69
gp >
There should be no counterexample 3 (mod 4) to 10^9. 2012-12-07, 11:45 #24 fivemack (loop (#_fork))   Feb 2006 Cambridge, England 22×32×179 Posts random thoughts comprised entirely of inutility My suspicion is that the '==3 mod 4' part means that counterexamples are large because they have to be a product of three primes ==3 mod 4. Looking at the way that pseudoprimes work in base 2: Code: for k in [1000000..2000000] do b:=IntegerRing(k)!2; t:=b^k; if (t eq 2 and not IsPrime(k)) then print k, Factorisation(k); for j in PrimeDivisors(k) do print " ",j,Order(FiniteField(j)!2); end for; end if; end for; eg 251*4051 is a base-2 pseudoprime because the order of 2 mod 251 and mod 4051 are both 50, and 50 divides 251*4051-1. what happens if you look at products of three primes-congruent-to-3-mod-4 where (3+2i) has small order mod each p_i? Code: TT:=[]; S:=-1; for t in [1..200000] do S:=4+S; while (not IsPrime(S)) do S:=4+S; end while; Q:=FiniteField(S^2); QQ:=S^2-1; i:=Sqrt(Q!-1); if (Order(3+2*i) lt 30*QQ/S) then print S,QQ/Order(3+2*i),Factorisation(Order(3+2*i)); Append(~TT,S); end if; end for; print TT; P:=PolynomialRing(Integers()); print "begin"; for a in TT do for b in TT do for c in TT do if (a gt b and b gt c) then So:=a*b*c; S:=IntegerRing(So); Sx:=PolynomialRing(S); Sx1:=quo; k:=Coefficient((3+2*q)^So,0); if (k lt 10000 or k gt So-10000) then print a,b,c,So, k; end if; end if; end for; end for; end for; print "end"; (http://magma.maths.usyd.edu.au/calc/) 2012-12-07, 12:00 #25 fivemack (loop (#_fork))   Feb 2006 Cambridge, England 22×32×179 Posts hmm, this is possibly quite an interesting additional question: How big a set of 4k+3-primes can you find where 3+2i has the same order mod each prime? How about the same odd-part-of-order? The birthday paradox isn't really on our side because the orders are around-p^2 rather than around-p. think of something like 1011863 990288 23 [ <47, 1>, <21529, 1> ] 47 23 43994 43056 21529 23 43994 43056 1012441 1008000 120 [ <241, 1>, <4201, 1> ] 241 120 8437 8400 4201 120 8437 8400 base-53 pseudoprimes the top lines are N, phi(N), order of 53 mod N, factorisation of N the indented lines are p, O=order of 53 mod p, (N-1)/O, phi(N)/O so having two primes where 53 has the same order was at least a good start 2012-12-07, 16:01 #26 fivemack (loop (#_fork))   Feb 2006 Cambridge, England 11001001011002 Posts ... ok, so the numbers where the order of 2 mod p is 47 are 2351 and 4513, which are among the factors of 2^47-1. So maybe we want to be looking at the factorisation into primes of Z[i] of (3+2*i)^k-1. But, hmm, it looks as if the norms of such numbers with k prime only have prime factors which are congruent to 1 mod 4 ... so we're getting 'accidental' factors for composite k, and they're a lot rarer and a lot smaller. Last fiddled with by fivemack on 2012-12-12 at 14:48 2012-12-08, 11:22 #27 LaurV Romulan Interpreter   "name field" Jun 2011 Thailand 980810 Posts There are no counterexample 3 (mod 4) for b=3+2i, up to n=10^10 (and a bit upper). I stoped it. 2012-12-10, 02:45 #28 wsc812   Dec 2012 China 11 Posts how to prove has no Positive integer root, is variable , is constant and are parameter ,and all character represents Positive integer.it's a question about the conjecture. my proof on mathoverflow : http://mathoverflow.net/questions/11...ese-conditions Last fiddled with by wsc812 on 2012-12-10 at 02:54 2012-12-10, 18:45 #29 JFranke   Nov 2012 Bonn University 108 Posts Pomerance's reasoning for BPSW works here as well I think Alex Kruppa is correct about this being essentially a combination of a (p-1) and a (p+1)-test. Moreover, the examples (if any) constructed by C. Pomerance's heuristic are all expamples of pseudoprimes for the test considered in this thread. 2012-12-11, 02:51   #30
ixfd64
Bemusing Prompter

"Danny"
Dec 2002
California

32×269 Posts Quote:
 Originally Posted by JFranke I think Alex Kruppa is correct about this being essentially a combination of a (p-1) and a (p+1)-test. Moreover, the examples (if any) constructed by C. Pomerance's heuristic are all expamples of pseudoprimes for the test considered in this thread.
I see you've finally registered here. Welcome! 2012-12-11, 09:16 #31 akruppa   "Nancy" Aug 2002 Alexandria 2,467 Posts Welcome on mersenneforum, Jens! 2012-12-13, 10:18 #32 wsc812   Dec 2012 China 11 Posts I have found the Complete proof for N=4k+3!!!!! 2012-12-13, 13:32   #33
Uncwilly
6809 > 6502

"""""""""""""""""""
Aug 2003
101×103 Posts

1010310 Posts Quote:
 Originally Posted by wsc812 I have found the Complete proof for N=4k+3!!!!!
Big blue letters means it must be true. Thread Tools Show Printable Version Email this Page Similar Threads Thread Thread Starter Forum Replies Last Post a1call Miscellaneous Math 194 2018-03-19 05:54 lukerichards Software 8 2018-01-24 22:30 1260 Software 17 2015-08-28 01:35 CRGreathouse Computer Science & Computational Number Theory 18 2013-06-08 19:12 jasong Math 1 2007-11-06 21:46

All times are UTC. The time now is 05:09.

Mon Nov 29 05:09:50 UTC 2021 up 128 days, 23:38, 0 users, load averages: 1.86, 1.78, 1.49
{"ft_lang_label":"__label__en","ft_lang_prob":0.7439555,"math_prob":0.8229692,"size":5164,"snap":"2021-43-2021-49","text_gpt3_token_len":1784,"char_repetition_ratio":0.094767444,"word_repetition_ratio":0.021765416,"special_character_ratio":0.40782338,"punctuation_ratio":0.16527545,"nsfw_num_words":0,"has_unicode_error":false,"math_prob_llama3":0.9676691,"pos_list":[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98,99,100,101,102,103,104,105,106,107,108],"im_url_duplicate_count":[null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null,null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2021-11-29T05:09:50Z\",\"WARC-Record-ID\":\"<urn:uuid:dc78e8d3-281c-4728-a39b-3cf351163cdc>\",\"Content-Length\":\"78985\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:b02a0d01-47d3-44f2-9dda-724fbcb3aa68>\",\"WARC-Concurrent-To\":\"<urn:uuid:f737f320-aee4-477a-91d7-23fd6c7d041c>\",\"WARC-IP-Address\":\"216.92.123.112\",\"WARC-Target-URI\":\"https://www.mersenneforum.org/showthread.php?s=926ad1d287d6e9db8f43bc7ecb84e76f&p=321552\",\"WARC-Payload-Digest\":\"sha1:KKUNAOFJGK7Q65PHO72W4MVEQPDX6SLU\",\"WARC-Block-Digest\":\"sha1:GJ55XPRS7E2P63LJCB6FQ6BJBSJ736IQ\",\"WARC-Identified-Payload-Type\":\"application/xhtml+xml\",\"warc_filename\":\"/cc_download/warc_2021/CC-MAIN-2021-49/CC-MAIN-2021-49_segments_1637964358688.35_warc_CC-MAIN-20211129044311-20211129074311-00600.warc.gz\"}"}