klx's profile - activity

https://doc.sagemath.org/html/en/reference/arithmetic_curves/sage/schemes/elliptic_curves/ell_point.html

I want to implement Extended GCD with SageMath so that the internal can be printed. The below is a modification of this. This should be straightforward, however, I might miss something;

def extended_euclides(a,b):

    s = 0; old_s = 1
    t = 1; old_t = 0
    r = b; old_r = a

    while r != 0:
        quotient,rem = old_r.quo_rem(r)

        old_r, r = r, old_r - quotient*r
        old_s, s = s, old_s - quotient*s
        old_t, t = t, old_t - quotient*t
    return [old_r, old_s, old_t]

K.<a> =  GF(2^8, modulus=x^8+x^4+x^3+x+1)
print(K)
print("m = ", K.modulus())

g = a^7 + a^5 + a^4 + a
h = a^4 + a^2 + 1

r,s,t = extended_euclides(g,h)

print("The GCD of \n \t [ {} ] and  \n \t [ {} ]  is \n\t [ {} ]".format(g,h,r))
print("Its Bézout coefficients are {} and {}".format(s,t))
assert r == g.gcd(h), "The gcd should be {}!".format(g.gcd(h))
assert r == s*g + t*h, "The cofactors are wrong!"

In the end, I've got the assertion error AssertionError: The gcd should be 1!

Any idea to solve this?

I see. I rather consider calculating it over the finite field instead of polynomials. Since the polynomials are always l

ExtGCD in Finite Fields I want to implement Extended GCD with SageMath so that the internal can be printed. The below is

ExtGCD in Finite Fields I want to implement Extended GCD with SageMath so that the internal can be printed. The below is

This should be the code. You have forgotten the ** instead of ^ from sage.all import * F = GF(2) R = F['x']; (x,) = R.

This should be the code. You have forgotten the ** instead of ^ from sage.all import * F = GF(2) R = F['x']; (x,) = R.

Yes, I want to the reverse, use the SageMath in Python.

I've seen this on ask.sagemath and stackoverflow

I want to use this in Python

k = GF(2)
R.<x> = k[]
k.extension(x^1000 + x^5 + x^4 + x^3 + 1, 'a')

The python code

from sage.all import *

F = GF(2)
R.<x> = k[]
K = F.extension(x^4 + x + 1, 'a')

print(K)

the R.<x> = k[] fails...

Is there a way to do this in python?

My final aim is finding the multiplicative inverse of an element using python with the sagemath import.

This should be the code. You have forgotten the ** instead of ^ from sage.all import * F = GF(2) R = F['x']; (x,) = R.

Using SageMath Finite Field Extension on Python. Yes, I want to the reverse, use the SageMath in Python. I've seen this

Using SageMath Finite Field Extension on Python. Yes, I want to the reverse, use the SageMath in Python. I've seen this

Also, I've seen this answer of yours and the function is not working, at least for me. Something has changed with Python

Also, I've seen this answer of yours and the function is not working, at least for me.

Also, I've seen this answer of yours and the function is not working, at least for me.

I guess substitute_expression is for something else, however, I couldn't find it in the docs.

I need to change Y into 1/2 *( y - a1 * x - a3) so that I can change variables in Elliptic curves.

var("X,Y,Z,x,y,a1,a2,a3,a4,a5,a6")
F = Y^2  *Z + a1 * X*Y *Z + a3 *Y *Z^2 - X^3 - a2 *X^2 * Z + a4*X*Z^ 2 + a6*Z^3
print(F)
F.substitute_expression(Y == 1/2 *( y - a1 * x - a3))  
print(F)

This produces an error on the 4th line;

AttributeError: 'sage.symbolic.expression.Expression' object has no attribute 'substitute_expression'

What is the proper way to change the variables?

Change of variables in symbolic computation I need to change Y into 1/2 *( y - a1 * x - a3) so that I can change variabl

def probit(p):
    return (sqrt(2)*erfinv(2*p-1)).n()

print(probit(0.025)) #Wikipedia Verification
print(probit(0.975))

P1 = [] 
for n in xsrange(1,100):
    P1.append( ( n , probit(n/100.0)) )

G1 = point2d(P1,color='blue',pointsize=10)
G1.show()

https://doc.sagemath.org/html/en/refe...

I've used the erfinv

def probit(p):
    return e^(sqrt(2)*erfinv(2*p-1)).n()

print(probit(0.025))
print(probit(0.975))

and outputs

0.140863494093217
7.09907138423134

Wikipedia states that

$\operatorname{probit}(0.025) = -1.96 = -\operatorname{probit}(0.975)$

Interestingly

print(e^(sqrt(2)*erfinv(1/2)).n())

outputs

1.96303108415826

Am I missing a point here?

[k_1 + 11x2, k_2 + 11x2, k_3, k_4 + 1, x1 + 23*x2, 2]. I've got this as an experiment, what is the value of k_3 and the last 2?

No automatic way form basis to eq.sub?

Sorry, $i$ was a typo that takes your time. Is this the only way? I was expecting that the symbolic package can try solve by Gaussian Elimination.

To solve the Hill Cipher a have setup symbolic equations;

R = IntegerModRing(26)
a,b,c,d,x1,x2 = var('k_1,k_2,k_3,k_4,x1,x2')
K = Matrix(SR,[[a,b],[c,d]]); K
print(K)

pl = vector(SR,[1,29])
ci = vector(SR,[8,21])
eq1,eq2 = (K*pl-ci)

pl = vector(SR,[25,22])
ci = vector(SR,[x1,i])
eq3,eq4 = (K*pl-ci)


pl = vector(SR,[7,16])
ci = vector(SR,[x2,19])
eq5,eq6 = (K*pl-ci)

solve((eq1,eq2,eq3,eq4,eq5,eq6), [a,b,c,d,x1,x2])

There is a problem that The symbolic equations are not working modulo 26. For example, the first equation is

k_1 + 29*k_2 - 8

So how can I properly turn them into modulo 26 and solve the equations in modulo 26?

I want to draw a line passing through a point P and Q on an elliptic curve.

E = EllipticCurve([0,0,0,4,20])
print(E)
E.plot(xmin=-20, xmax=20, ymin=-20, ymax=20)
plotE = E.plot()
P=E(1,5)
R = - (P + P)
plotE += line([P.xy(),R.xy()],color='red')

Using this only draw a line segment from P to Q. How can I draw a line? I've looked on the doc but did not see an option for this. Any solution?

Sorry, cells mixed! I have put it into a curve defined over GF. Yes, it is working. Thanks.

PR = (lift(vR - vP)).normalized() is needed and after that TypeError: unsupported operand parent(s) for -: 'Vector space of dimension 2 over Finite Field of size 29' and 'Vector space of dimension 2 over Symbolic Ring' on the last line