IntegerExponent Equivalent in Pari-GP

[a1call]

Hi,

I can't find an equivalent to IntegerExponent function in Pari-GP.

Is there one?, or should One write his own routine?

Thanks in advance


https://reference.wolfram.com/langua...rExponent.html

http://math.mit.edu/~brubaker/PARI/PARIrefcard.pdf

[paulunderwood]

?logint

logint(x,b,&z): return the largest integer e so that b^e <= x, where the
parameters b > 1 and x > 0 are both integers. If the parameter z is present,
set it to b^e.

Example:

? logint(139,2)
7

Of course, you will have take the remainder when divided by the base to see if it divides exactly.

[a1call]

Thank you paulunderwood,
That was fast.
For the record that function is missing from the reference card that I linked.
That saves me quite a bit of time.
Thank you again.

[science_man_88]

Originally Posted by paulunderwood

?logint



Example:

? logint(139,2)
7

Of course, you will have take the remainder when divided by the base to see if it divides exactly.

only in later versions of PARI-GP than that refcard URL is good for. 2.75 is the earliest I know of it from ( only because I looked at my old installs). edit: of course before this all that's needed is floor(log(x)/log(y)) where y is the number you want to take the log base as. but built in's are usually faster. see the PARI commands thread for some other things that may help. of course depending on what you plan on doing with the number digits() may also come in handy or even fromdigits() in newer versions. however on reading the other link I think you technically want valuation() not logint as logint doesn't care about divisibility valuation does.

(18:29) gp > ??valuation
valuation(x,p):

Computes the highest exponent of p dividing x. If p is of type integer, x must be an integer, an intmod whose modulus is divisible by p, a fraction, a q-adic number with q = p, or a polynomial or power series in which case the valuation is the minimum of the valuation of the coefficients.

If p is of type polynomial, x must be of type polynomial or rational function, and also a power series if x is a monomial. Finally, the valuation of a vector, complex or quadratic number is the minimum of the component valuations.

If x = 0, the result is LONG_MAX (2^{31}-1 for 32-bit machines or 2^{63}-1 for 64-bit machines) if x is an exact object. If x is a p-adic numbers or power series, the result is the exponent of the zero. Any other type combinations gives an error.

The library syntax is long gvaluation(GEN x, GEN p).

[a1call]

Originally Posted by science_man_88

only in later versions of PARI-GP than that refcard URL is good for. 2.75 is the earliest I know of it from ( only because I looked at my old installs). edit: of course before this all that's needed is floor(log(x)/log(y)) where y is the number you want to take the log base as. but built in's are usually faster. see the PARI commands thread for some other things that may help. of course depending on what you plan on doing with the number digits() may also come in handy or even fromdigits() in newer versions. however on reading the other link I think you technically want valuation() not logint as logint doesn't care about divisibility valuation does.

Thank you science_man_88,
That would-have-been/is as good in saving me time.

[science_man_88]

Originally Posted by a1call

Thank you science_man_88,
That would-have-been/is as good in saving me time.

of course if it's used for factoring factor() can do it for all primes you might be able to save yourself time with that also then of course there's factorint() (though that one doesn't have an upper limit to my knowledge and only factors using integers. if you want a list of divisors that are powers there's divsors() and ispower() ( or even fordiv in some versions) if you want to check if something is square free there's an issquare free function etc. edit: and of course if something is squarefree it is also power-free since all powers greater than squares that could divide it are not squarefree.

[paulunderwood]

Originally Posted by science_man_88

??valuation

[a1call]

Originally Posted by science_man_88

of course if it's used for factoring factor() can do it for all primes you might be able to save yourself time with that also then of course there's factorint() (though that one doesn't have an upper limit to my knowledge and only factors using integers. if you want a list of divisors that are powers there's divsors() and ispower() ( or even fordiv in some versions) if you want to check if something is square free there's an issquare free function etc.

Thank you,
My current use is not for factoring. Just need to formulate the highest integer power to run my WDP code on Pari-gp.
The only thing I need to figure out now is how to use variables.

[science_man_88]

Originally Posted by a1call

Thank you,
My current use is not for factoring. Just need to formulate the highest integer power to run my WDP code on Pari-gp.
The only thing I need to figure out now is how to use variables.

probably would depend on scope and if they are predefined values at first. not to mention if you want to do things in parallel global variables would be NOT a great choice. I'm not even sure what WDP is. one thing that comes to mind is in theory primes are the simplest case if they can't divide by prime power's they don't divide by any power other than 1^n which is always 1.

[a1call]

Well variables seem straight forward.
Is there a way to save/load scripts/routines/programs?

Thank you.


ETA: WDP = Wolfram Development Platform:

https://develop.wolframcloud.com/app/

[science_man_88]

Originally Posted by a1call

Well variables seem straight forward.
Is there a way to save/load scripts/routines/programs?

Thank you.


ETA: WDP = Wolfram Development Platform:

https://develop.wolframcloud.com/app/

you can have a codes.gp load on start. you'll need to setup your GPRC file to do that though or copy a text file in each time. anything done earlier in a session will still be accessible and in some versions you may even be able to get it from older log files as I can hit back ( or up in my case I think) and see me looking up about fromdigits etc. ( though I don't think that's worth the risk)