Judith Gatz on Tue, 15 Jan 2002 10:05:16 +0100


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

easy precision problem ( I hope ;)


Hi!

I am a real pari-gp-beginner, so I hope it´s easy for you:
I want to make a number of exact factorisations like:

? for(x=65,75,print(factor(truncate(precision(exp(x),1000)),0)))

the results cant be exact, there are too many 2^x:
[2, 1; 3, 2; 941605135783518730078768673, 1]
[2, 2; 19, 2; 43, 1; 1217, 1; 119723, 1; 5092511140946699, 1]
[2, 3; 401, 1; 1213, 1; 20929, 1; 1537753234015292017, 1]
[2, 3; 42553450624146756517211339837, 1]
[2, 4; 3, 1; 7, 1; 13, 1; 211853977234152646525695923, 1]
[2, 5; 5, 1; 7, 1; 19, 1; 23, 1; 433, 1; 677, 1; 17532252168972131581, 1]
[2, 9; 17, 1; 43, 1; 18269256662969027516768267, 1]
[2, 8; 5, 1; 43, 1; 239, 1; 7927, 1; 39181, 1; 4549272969690959, 1]
[2, 11; 3, 1; 4637, 1; 1773408868724932787900621, 1]
[2, 11; 67059715797860165419149830557, 1]
[2, 14; 7, 1; 181, 1; 17984136432019657849799891, 1]
 ^^^^

before that, I tried to increase the precision with:  (I don´t know if this
makes sence...)
? primelimit=1000000000000000000000000000000000000000000000000000000000000
%1 = 1000000000000000000000000000000000000000000000000000000000000
? parisize=100000000000000
%2 = 100000000000000
? realprecision= 500
%3 = 500

Can I do this in a better way?
I am mainly interested in the set of resulted primenumbers of the series
[exp(n)], n= 1,2,3,...


Thank You very much,
regards, Judith