mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Miscellaneous Math (https://www.mersenneforum.org/forumdisplay.php?f=56)
-   -   Thread for posting tiny primes (https://www.mersenneforum.org/showthread.php?t=13650)

science_man_88 2010-07-30 21:06

193 is the next one like 79.

science_man_88 2010-07-30 21:28

my original hypothesis was when I saw all the factors that (factor/p)%8=0 I assumed it was always true then I found 43 and it didn't work then I hypothesized it might be because it had 3 factors and noticed that if one of the lower 2 factors/p aren't divisible by 8 neither is the third but one of them can be divisible by 8 and not necessitate the third to be. that's about how far I've got into hypothesizing as i haven't figured out a pattern in the 4 and 5 factor ones.

3.14159 2010-07-30 22:40

[QUOTE=science_man_88]my original hypothesis was when I saw all the factors that (factor/p)%8=0 I assumed it was always true then I found 43 and it didn't work then I hypothesized it might be because it had 3 factors and noticed that if one of the lower 2 factors/p aren't divisible by 8 neither is the third but one of them can be divisible by 8 and not necessitate the third to be. that's about how far I've got into hypothesizing as i haven't figured out a pattern in the 4 and 5 factor ones.[/QUOTE]

Are you sure this isn't the law of small numbers at work again?

science_man_88 2010-07-30 22:46

[QUOTE=3.14159;223400]Are you sure this isn't the law of small numbers at work again?[/QUOTE]

[CODE]for(x=1,50,print(((factor(2^prime(x)-1)\prime(x))/8)%1==0","x))[/CODE]

is a code I tried to prove it it seems to work for every Mersenne number in a prime index any way to change it to check only factors less than 2^n-1 would help.

science_man_88 2010-07-30 23:58

[CODE]for(x=1,10,print((floor(factor(2^prime(x)-1)\prime(x)))%8))[/CODE]

yeah I figured I had an error there's one thing I have to ask. How do I check to see if it ever matches 0 0 for the response is it as simple as a if instead of the print ?

[CODE](21:06) gp > for(x=1,10,print((floor(factor(2^prime(x)-1)\prime(x)))%8))
Mat([1, 0])
Mat([2, 0])
Mat([6, 0])
Mat([2, 0])
[2, 0; 0, 0]
Mat([6, 0])
Mat([6, 0])
Mat([2, 0])
[2, 0; 0, 0]
[0, 0; 6, 0; 0, 0][/CODE]

I want to check for those 0, 0 parts

CRGreathouse 2010-07-31 03:48

First: We'll save the value like so:
[code]f=factor(2^prime(x)-1)\prime(x))%8;[/code]
so your code is now
[code]for(x=1,10,f=factor(2^prime(x)-1)\prime(x))%8;print(f))[/code]

[QUOTE=science_man_88;223404][CODE]for(x=1,10,print((floor(factor(2^prime(x)-1)\prime(x)))%8))[/CODE]

yeah I figured I had an error there's one thing I have to ask. How do I check to see if it ever matches 0 0 for the response is it as simple as a if instead of the print ?[/QUOTE]

I'm not quite sure what you want. If you want to check if it's *equal* to 0,0 then do
[code]if(f!=Mat([0,0]),print(f" is nonzero"))[/code]
rather than
[code]print(f)[/code]

If you want to check if any of the components are zero, it's best to make a special function for that:
[code]anyZero(M)={
for(i=1,#M[,1],
if(M[i,1] == 0 && M[i,2] == 0, return(1))
);
0
};
addhelp(anyZero, "anyZero(M): Returns 1 if any row of the n x 2 matrix M is all zeros, or 0 otherwise.")[/code]

Then you can do whatever you want with it:
[code]if(anyZero(f), zeroCount++, print(f" has no all-zero rows"))[/code]

But I'll admit I don't understand the purpose behind any of this.

science_man_88 2010-07-31 11:15

I'm trying to prove my hypothesis that most prime index Mersenne's that aren't prime have a factor that (factor\n)%8==0 and I have a phrase for the exceptions I'm finding to cover them.

CRGreathouse 2010-07-31 13:53

OK, so you need only check the first component.

science_man_88 2010-07-31 14:09

[QUOTE=CRGreathouse;223450]OK, so you need only check the first component.[/QUOTE]

yeah well I know printing takes a long time so I only want it to print prime(x) if (factor\n)%8==0 which I can't seem to pull off.

3.14159 2010-07-31 14:22

[QUOTE=science_man_88]I'm trying to prove my hypothesis that most prime index Mersenne's that aren't prime have a factor that (factor\n)%8==0 and I have a phrase for the exceptions I'm finding to cover them.
[/QUOTE]

Wait.. I seem to forget what the \ function is. Is it the floor function? By the way, are you still working on the trial factor code?

science_man_88 2010-07-31 14:23

[QUOTE=3.14159;223453]Wait.. I seem to forget what the \ function is. Is it the floor function? By the way" Are you still working on the trial factor code?[/QUOTE]

\ is floor and in a way but I'm trying to get more specific.


All times are UTC. The time now is 15:01.

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.