Skip to content

Instantly share code, notes, and snippets.

@evinism
Last active September 9, 2022 20:04
Show Gist options
  • Save evinism/23d5c199779118d6013b20f88ed2853f to your computer and use it in GitHub Desktop.
Save evinism/23d5c199779118d6013b20f88ed2853f to your computer and use it in GitHub Desktop.
Jones-Sato-Wada-Wiens prime polynomial encoded in Python
# Below is the Jones-Sato-Wada-Wiens prime polynomial encoded in Python
# Ultimately, expressing it in this format is somewhat useless because of the
# size of numbers involved.
# https://thomas1111.wordpress.com/2015/02/01/what-do-solutions-of-the-jones-sato-wada-wiens-polynomial-look-like/
def jones_sato_wada_wiens(a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u, v, w, x, y, z):
return (k * 1) * (
1
-(w*z+h+j-q)**2
-((g*k+2*g+k+1)*(h+j)+h-z)**2
-(2*n+p+q+z-e)*2
-(16*((k+1)**3)*(k+2)*((n+1)**2)+1-f**2)**2
-((e**3)*(e+2)*((a+1)**2)+1-o**2)**2
-((a**2-1)*y**2+1-x**2)**2
-(16*r**2*y**4*(a**2-1)+1-u**2)**2
-(((a+u**2*(u**2-a))**2-1)*(n+4*d*y)**2+1-(x+c*u)**2)**2
-(n+l+v-y)**2
-((a**2-1)*l**2+1-m**2)**2
-(a*i+k+1-l-i)*2
-(p+l*(a-n-1)+b*(2*a*n+2*a-n**2-2*n-2)-m)**2
-(q+y*(a-p-1)+s*(2*a*p+2*a-p**2-2*p-2)-x)**2
-(z+p*l*(a-p)+t*(2*a*p-p**2-1)-p*m)**2
)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment