PARI's commands

[science_man_88]

scratch f we can just use v=f to redefine it back to v after.

[CRGreathouse]

Originally Posted by science_man_88

CRG I got my code to work to find primes and fill a vector of 0's with them I'm thinking of how to make a function VecApp(v,x,f) that takes one vector(v) into another vector(f) with 1 extra element which is appended with value x maybe that will make it easier on me. I think this would need to say something like f[i] = v[i] f[i]+1 = x but I'm without reading them all I have no idea how to get i not to go over the amount of indexes of v.

concat([1,2,3],10)

[science_man_88]

I'm an idiot thanks.

[3.14159]

Originally Posted by CRGreathouse

Just look at my code for how to use the vector!

my(v=vectorsmall(sz, i, 1))?

What is supposed to be the substitution here?

[CRGreathouse]

Originally Posted by science_man_88

I'm an idiot thanks.

You were trying to write a function that already existed. That's pretty common. I wrote my own version of eint1, at the very least -- but probably others, too.

[3.14159]

A 1 in 10⁶⁵ chance:

I entered a(65), which is nextprime(random(10^65)). The result: 2.

[CRGreathouse]

Well, 3 in 10^65 anyway. Wow.

[3.14159]

Originally Posted by CRGreathouse

Well, 3 in 10^65 anyway. Wow.

It was in fact, too good to be true. It turned out PARI was messed for a while.

I entered a(7) and got 2, 6 times in a row. I tried a(8) through a(12), and got 2, about 20 times on end.

This would normally happen only about 1 in ≈10¹⁶³ times. I'm going with PARI being a bit buggy. Or did I have a supreme stroke of luck and just defy odds of 1 in 10¹⁶³?

Also: The lottery should be changed to guessing 10 random numbers between 1 and 10¹⁵.

(To ensure that nobody wins, )

[CRGreathouse]

You may have inadvertently overwritten a. Try ?a to see the current definition.

[3.14159]

Originally Posted by CRGreathouse

You may have inadvertently overwritten a. Try ?a to see the current definition.

Nope. I changed nothing. It was either incredible luck or the computer fucking up.

[science_man_88]

(15:27) gp > v=vector(1,n,0);a=1;for(x=1,10000,if(isprime(x),v[a]=x;concat(v[],0);a=a+1))
  ***   unused characters: v=vector(1,n,0);a=1;for(x=1,10000,if(is
                                              ^--------------------