PARI's commands

[CRGreathouse]

Originally Posted by 3.14159

If there is a k < b^b restriction, b ≥ 1.

That was my third suggestion.

Originally Posted by 3.14159

If there's no k < b^b restriction, b > 1.

That would be a fourth way to do it. It has relative density ≈ 0.528 in the primes.

[science_man_88]

Originally Posted by 3.14159

If there is no k < b^b restriction, and b = 1 is allowed, it is precisely the same as A000040.

I will begin work on the sequence with the k < b^b restriction.

sorry CRG I see what he means today lol k*2^1+1 can give all primes.

[3.14159]

Also: Is there a command or script that allows printing in numerical order?

[CRGreathouse]

Originally Posted by 3.14159

If there is no k < b^b restriction, and b = 1 is allowed, it is precisely the same as A000040.

Neither of my first two sequences restrict k to be less than b^b and both allow b = 1, but neither is A000040. I invite you to re-read #1209.

[science_man_88]

Originally Posted by 3.14159

Also: Is there a command or script that allows printing in numerical order?

maybe vecsort could help you in pari but I'm not sure how to use it.

[CRGreathouse]

Originally Posted by 3.14159

Also: Is there a command or script that allows printing in numerical order?

If you put the numbers in a vector, you can use vecsort to order them. vecsort(v) returns a sorted version of v, while vecsort(v,,8) returns a sorted version of v with the duplicates removed.

[CRGreathouse]

Originally Posted by science_man_88

maybe vecsort could help you in pari but I'm not sure how to use it.

Bingo! sm88 beat me to it.

[science_man_88]

Originally Posted by CRGreathouse

Bingo! sm88 beat me to it.

calls for a celebration I think that's a first lol.

[CRGreathouse]

Originally Posted by science_man_88

sorry CRG I see what he means today lol k*2^1+1 can give all primes.

You can see that I understand that, since I mention A000040 explicitly in post #1209. But none of my three sequences are equal to A000040, as even a cursory examination will show.

[CRGreathouse]

Originally Posted by science_man_88

calls for a celebration I think that's a first lol.

I'll lift a glass in your honor. (Literally, why not!)

[science_man_88]

Originally Posted by CRGreathouse

You can see that I understand that, since I mention A000040 explicitly in post #1209. But none of my three sequences are equal to A000040, as even a cursory examination will show.

oh wait k*b^b+1 would give k*1+1 for b=1 still all the primes I think.

doh I fell to his thinking lol.