Physics
Maple provides a state-of-the-art environment for algebraic computations in Physics, with emphasis on ensuring that the computational experience is as natural as possible. The theme of the Physics project for Maple 2016 has been the consolidation of the new functionality introduced during the last three releases, through more than 300 enhancements across the board, together with signification enhancements and new functionality in General Relativity. Most notably, for the 100th anniversary of the presentation of the theory of relativity, in Maple 2016 we are bringing to completion the digitalization of the solutions to Einstein's equations collected, from more than 4000 papers, in the classic book "Exact Solutions to Einstein's Equations" by H. Stephani - et al.
In addition, Maple 2016 implements new general functionality that is relevant within and beyond Physics, including a new Factor command with the ability to perform factorization in expressions involving products of noncommutative operators and the ability to compute with differential operators algebraically, that is, operating with them using multiplication to express application. Combining these two developments, it is now possible to factorize algebraic expressions involving differential operators and so, for example, solve partial differential equations through factorization.
Taking all together, the more than 350 enhancements throughout the entire package increased robustness, versatility and functionality, extending again the range of Physics-related algebraic computations that can be done using computer algebra software and in a natural way.
As part of its commitment to providing the best possible environment for algebraic computations in Physics, Maplesoft launched a Maple Physics: Research and Development web site with Maple 18, which enabled users to download research versions, ask questions, and provide feedback. The results from this accelerated exchange with people around the world have been incorporated into the Physics package in Maple 2016.
Completion of the Database of Solutions to Einstein's Equations
Operatorial Algebraic Expressions Involving the Differential Operators ∂μ, 𝒟μ and ∇ (nabla)
Factorization of Expressions Involving Noncommutative Operators
Tensors in Special and General Relativity
Vectors Package
The Physics Library
Redesigned Functionality and Miscellaneous
A database of solutions to Einstein's equations was added to the Maple library in Maple 15 with a selection of metrics from "Stephani, H.; Kramer, D.; MacCallum, M.; Hoenselaers, C.; and Herlt, E., Exact Solutions to Einstein's Field Equations". More metrics from this book were added for Maple 16, Maple 17, and Maple 18, up to 225 metrics. In Maple 2016, for the 100th anniversary of the presentation of the theory of Relativity by A. Einstein, we brought to completion the digitalization of all the 971 metric solutions collected from more than 4000 papers presented in the General Relativity literature. All these metrics can be loaded or searched using g_ (the Physics command representing the spacetime metric that also sets the metric to your choice in one go) or using the command DifferentialGeometry:-Library:-MetricSearch. When a metric is loaded, all the General Relativity tensors (Christoffel, Ricci, Einstein, and Riemann) are automatically computed on background and available for further computations. These metrics can also be changed in various ways and, automatically, all the General Relativity tensors are recomputed on background. You can work with these metrics, the ones in the database or any other one you set, using the Tetrads formalism too, with the commands of the Physics:-Tetrads package.
In the Maple PDEtools package, you have the mathematical tools - including a complete symmetry approach - to work with the underlying partial differential equations. By combining the functionality of the Physics:-Tetrads package, the Physics:-TransformCoordinates command, and the ability to compute Riemann and Weyl invariants, you can also formulate and, depending on the the metrics also resolve, the equivalence problem; that is: to answer whether or not, given two metrics, they can be obtained from each other by a transformation of coordinates, as well as compute the transformation.
Examples
Load Physics, set the metric to Schwarzschild (and everything else automatically) in one go
restart; withPhysics:
g_sc
Systems of spacetime Coordinates are: X=r,θ,φ,t
Default differentiation variables for d_, D_ and dAlembertian are: X=r,θ,φ,t
The Schwarzschild metric in coordinates r,θ,φ,t
Parameters: m
That is all you need to do: although the strength in Physics is to compute with tensors using indicial notation, by setting the metric as in all of the tensor components of the General Relativity tensors are also derived on the fly and ready for use. For instance these are the definition in terms of Christoffel symbols, and the covariant components of the Ricci tensor for Schwarzschild solution
Riccidefinition
Rμ,ν=∂α⁡Γ⁢αμ,ν⁢αμ,ν−∂ν⁡Γ⁢αμ,α⁢αμ,α+Γ⁢βμ,ν⁢βμ,ν⁢Γ⁢αβ,α⁢αβ,α−Γ⁢βμ,α⁢βμ,α⁢Γ⁢αν,β⁢αν,β
Ricci
Rμ,ν=0000000000000000
These are the 16 Riemann invariants for Schwarzschild solution, using the formulas by Carminati and McLenaghan
Riemanninvariants
r0=0,r1=0,r2=0,r3=0,w1=6⁢m2r6,w2=6⁢m3r9,m1=0,m2=0,m3=0,m4=0,m5=0
The related Weyl scalars in the context of the Newman-Penrose formalism; the definition is in terms of the Weyl tensor and the tetrad of tensors lμ,nμ,mμ,m&conjugate0;μ of the Newman-Penrose formalism
Weylscalarsdefinition
ψ__0=−C⁢μ,ν,α,β⁢μ,ν,α,β⁢lμ⁢mν⁢lα⁢mβ,ψ__1=−C⁢μ,ν,α,β⁢μ,ν,α,β⁢lμ⁢nν⁢lα⁢mβ,ψ__2=−C⁢μ,ν,α,β⁢μ,ν,α,β⁢lμ⁢mν⁢m&conjugate0;α⁢nβ,ψ__3=−C⁢μ,ν,α,β⁢μ,ν,α,β⁢lμ⁢nν⁢m&conjugate0;α⁢nβ,ψ__4=−C⁢μ,ν,α,β⁢μ,ν,α,β⁢nμ⁢m&conjugate0;ν⁢nα⁢m&conjugate0;β
Weylscalars
ψ__0=0,ψ__1=0,ψ__2=−mr3,ψ__3=0,ψ__4=0
These are the matrix components of the Christoffel symbols of the second kind (that describe, in coordinates, the effects of parallel transport in curved surfaces), when the first of its three indices is equal to 1; contravariant indices are prefixed by ~
Christoffel~1,alpha,beta,matrix
Γ⁢1α,β⁢1α,β=mr⁢−r+2⁢m0000−r+2⁢m0000−r+2⁢m⁢sin⁡θ20000−2⁢m2+m⁢rr3
In Physics, the Christoffel symbols of the first kind are represented by the same object (one command, Christoffel, not two) just by taking the first index covariant, as we do when computing with paper and pencil
Christoffel1,alpha,beta,matrix
Γ1,α,β=m−r+2⁢m20000r0000r⁢sin⁡θ20000−mr2
One could query the database, directly from the spacetime metrics, for example about the solutions (metrics) to Einstein's equations related to Levi-Civita, the Italian mathematician
g_civi
____________________________________________________________
12,16,1=Authors=Bertotti (1959),Kramer (1978),Levi-Civita (1917),Robinson (1959),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=Homogeneous
12,18,1=Authors=Bertotti (1959),Kramer (1978),Levi-Civita (1917),Robinson (1959),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=Homogeneous
12,19,1=Authors=Bertotti (1959),Kramer (1978),Levi-Civita (1917),Robinson (1959),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=Homogeneous,Comments=_lambda=_zeta
22,7,1=Authors=Levi-Civita (1917), Frehland (1971),PrimaryDescription=Vacuum,SecondaryDescription=Cylindrically-Symmetric,Comments=Locally static, Weyl class_m=0,1 - flat, _m=1/2, 2, -1 - PetrovType D
Each triad of numbers indicates the chapter and equation number with which each of these metrics appear in "Exact Solutions to Einstein's Field Equations". For example, [12, 16, 1] is the metric found in Chapter 12 with equation number 16.1. These solutions can be set in one go from the metrics command, just by indicating the triad of numbers as follows
g_12,16,1
Systems of spacetime Coordinates are: X=t,x,θ,φ
Default differentiation variables for d_, D_ and dAlembertian are: X=t,x,θ,φ
The Bertotti (1959), Kramer (1978), Levi-Civita (1917), Robinson (1959) metric in coordinates t,x,θ,φ
Parameters: k,κ0,β
Resetting the signature of spacetime from "- - - +" to `- + + +` in order to match the signature in the database of metrics:
Automatically, everything gets set accordingly; these are the contravariant components of the related Ricci tensor
Ricci~
R⁢μ,ν⁢μ,ν=1k4⁢sinh⁡x20000−1k400001k400001k4⁢sin⁡θ2
One works with the Newman-Penrose formalism frequently using tetrads (local system of references); the Physics subpackage for this is Tetrads
withTetrads
Setting lowercaselatin letters to represent tetrad indices
Defined as tetrad tensors see ?Physics,tetrads,𝔢a,μ,ηa,b,γa,b,c,λa,b,c
Defined as spacetime tensors representing the NP null vectors of the tetrad formalism see ?Physics,tetrads,lμ,nμ,mμ,m&conjugate0;μ
IsTetrad,NullTetrad,OrthonormalTetrad,PetrovType,SegreType,TransformTetrad,e_,eta_,gamma_,l_,lambda_,m_,mb_,n_
To query about the definition of any of these tensors, of the Tetrads or Physics packages, index them with the keyword definition
l_definition
lμ⁢l⁢μ⁢μ=0,lμ⁢n⁢μ⁢μ=−1,lμ⁢m⁢μ⁢μ=0,lμ⁢m&conjugate0;⁢μ⁢μ=0,gμ,ν=−lμ⁢nν−lν⁢nμ+mμ⁢m&conjugate0;ν+mν⁢m&conjugate0;μ
gamma_definition
γa,b,c=𝒟ν⁡𝔢a,μ⁢𝔢bμbμ⁢𝔢cνcν
This is the orthonormal tetrad 𝔢a,μ related to the metric set in set by the package (you can change these tetrads in different ways using the TransformTetrad command)
e_
𝔢a,μ=sinh⁡x⁢k0000k0000k0000sin⁡θ⁢k
You can check these components directly using the definition. In this case, the right-hand side is the (orthonormal) tetrad metric
e_definition
𝔢a,μ⁢𝔢bμbμ=ηa,b
eta_
This shows that, for the tetrad components given by (13), the definition (14) holds
TensorArray
−1=−10=00=00=00=01=10=00=00=00=01=10=00=00=00=01=1
One frequently works with a different signature and null tetrads; set that, and everything related automatically
Setuptetradmetric = null
tetradmetric=1,2=−1,3,4=1
𝔢a,μ=2⁢sinh⁡x⁢k22⁢k2002⁢sinh⁡x⁢k2−2⁢k200002⁢k2I2⁢2⁢sin⁡θ⁢k002⁢k2−I2⁢2⁢sin⁡θ⁢k
Verifying (17) using the definition (14)
0=0−1=−10=00=0−1=−10=00=00=00=00=00=01=10=00=01=10=0
You can also change the signature and everything gets automatically recomputed as well, from the components of the tensors to the definition of Weyl scalars. You can query the value of the signature using Setup
Setupsignature
signature=- + + +
Set the signature to be + - - - , compare the components of the null tetrad metric with the components of and verify the tetrad 𝔢a,μ
Setupsignature = `+---`
signature=+ - - -
𝔢a,μ=−I2⁢2⁢sinh⁡x⁢k−I2⁢2⁢k00−I2⁢2⁢sinh⁡x⁢kI2⁢2⁢k0000−I2⁢2⁢k2⁢sin⁡θ⁢k200−I2⁢2⁢k−2⁢sin⁡θ⁢k2
0=01=10=00=01=10=00=00=00=00=00=0−1=−10=00=0−1=−10=0
The related 16 Riemann invariants
r0=0,r1=1k4,r2=0,r3=14⁢k8,w1=0,w2=0,m1=0,m2=0,m3=0,m4=0,m5=0
The ability to query rapidly, set things in one go, change everything and have all the quantities automatically adjusted are at the realm of the design of the General Relativity functionality of the Physics package, resulting in a rather flexible computational environment.
These are the metrics by Kaigorodov; next are those published in 1962
g_Kaigorodov
12,34,1=Authors=Kaigorodov (1962),Cahen (1964),Siklos (1981),Ozsvath (1987),PrimaryDescription=Einstein,SecondaryDescription=Homogeneous,Comments=All metrics with _epsilon <> 0 are equivalent to the cases _epsilon = +1, -1, _epsilon = 0 is anti-deSitter space
12,35,1=Authors=Kaigorodov (1962),Cahen (1964),Siklos (1981),Ozsvath (1987),PrimaryDescription=Einstein,SecondaryDescription=Homogeneous,SimpleTransitive
38,2,1=Authors=Kaigorodov (1962),PrimaryDescription=Hypersurface-Homogeneous,SecondaryDescription=AlgebraicallySpecial
38,3,1=Authors=Kaigorodov (1962),PrimaryDescription=Hypersurface-Homogeneous,SecondaryDescription=AlgebraicallySpecial,Comments=Note, there is a typo (the x^4 term) in this metric in Stephani which we were not able to verify from the original. This metric may not be correct!
g_`1962`
12,13,1=Authors=Ozsvath, Schucking (1962),PrimaryDescription=Vacuum,SecondaryDescription=Homogeneous,PlaneWave,Comments=geodesically complete, no curvature singularities
12,14,1=Authors=Petrov (1962),PrimaryDescription=Vacuum,SecondaryDescription=Homogeneous,SimpleTransitive
18,2,1=Authors=Ehlers and Kundt (1962),PrimaryDescription=vacuum solutions ,SecondaryDescription=degenerate static vacuum solution,Comments=Table 18.2, case AI, _Psi2=-_b/(2*r^3)
18,2,2=Authors=Ehlers and Kundt (1962),PrimaryDescription=vacuum solutions ,SecondaryDescription=degenerate static vacuum solution,Comments=Table 18.2, case AII, _Psi2=_b*(1/(2*z^3))
18,2,3=Authors=Ehlers and Kundt (1962),PrimaryDescription=vacuum solutions ,SecondaryDescription=degenerate static vacuum solution,Comments=Table 18.2, case AIII, _Psi2=1/(2*z^3)
18,2,4=Authors=Ehlers and Kundt (1962),PrimaryDescription=vacuum solutions ,SecondaryDescription=degenerate static vacuum solution,Comments=Table 18.2, case BI, _Psi2=-_b*(1/(2*r^3))
18,2,5=Authors=Ehlers and Kundt (1962),PrimaryDescription=vacuum solutions ,SecondaryDescription=degenerate static vacuum solution,Comments=Table 18.2, case BII, _Psi2=_b*(1/(2*z^3))
18,2,6=Authors=Ehlers and Kundt (1962),PrimaryDescription=vacuum solutions ,SecondaryDescription=degenerate static vacuum solution,Comments=Table 18.2, case BIII, _Psi2=1/(2*z^3)
18,2,7=Authors=Ehlers and Kundt (1962),PrimaryDescription=vacuum solutions ,SecondaryDescription=degenerate static vacuum solution,Comments=Table 18.2, case C, _Psi2=(1/2)*(x+y)^3, + case
18,2,8=Authors=Ehlers and Kundt (1962),PrimaryDescription=vacuum solutions ,SecondaryDescription=degenerate static vacuum solution,Comments=Table 18.2, case C, _Psi2=-(1/2)*(x+y)^3, - case
22,59,1=Authors=Misra and and Radhakrishna (1962),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=Cylindrically-Symmetric,Comments=Null EM field, _W=Int((diff(_Theta(u),u)^2+diff(_eta(u),u)^2),u), u=(t-rho)/sqrt(2), the non-null case is with _Theta(t,rho)=_Psi(t,rho)*cos(_alpha), _Theta(t,rho)=_Psi(t,rho)*sin(_alpha)
26,21,1=Authors=Newman and Tamburino (1962),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=Non-Aligned, Non-Null,Comments=Spherical class
26,22,1=Authors=Newman and Tamburino (1962),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=Non-Aligned, Non-Null,Comments=Cyllindrical class
26,23,1=Authors=Newman and Tamburino (1962),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=Non-Aligned, Non-Null,Comments=Cyllindrical class
27,27,1=Authors=Robinson and Trautman 1962, Debney et al. 1969, Talbot 1969, Robinson et al. 1969a, Lind 1974 ,PrimaryDescription=Generic,SecondaryDescription=geodesic, shearfree and diverging null congruence,Comments=rho:=1/(-(r+r0+I*Sigma)): Sigma:=-2*I*P^2*((diff(L,zb)-L*diff(L,u))-(diff(Lb,z)-L*diff(Lb,u))):
27,37,1=Authors=Robinson and Trautman (1962),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=Non-Aligned, non-null,Comments=admits geodesic, shearfree, twistfree null congruence, _rho=-1/r=_rho_b
28,8,1=Authors=Robinson-Trautman (1962),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=diff(2*_P(u,z,zb)^2*diff(2*_P(u,z,zb)^2*ln(_P(u,z,zb)),z,zb),z,zb)+12*_m(u)*diff(ln(_P(u,z,zb)),u)-4*diff(_m(u),u)=0
28,12,1=Authors=Robinson-Trautman (1962),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=2*P(u, z, zb)^2*diff(ln(P(u, z, zb)),z,zb)=_K(u), _K(u)=0,+1,-1
28,16,1=Authors=Robinson-Trautman (1962),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=The coordinate zeta is changed to xi,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,26,1=Authors=Robinson, Trautman (1962),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=One can use the diffeo r -> -r and u -> -u to make the assumption r > 0,The case _m = 0 is Stephani, [28, 16,1],The metric is type D at points where r = 3*_m/(xi1+xi2) and type II on either side of this hypersurface. For convenience, it is assumed that 3*_m - r*(xi1 + xi2) > 0,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,26,2=Authors=Robinson, Trautman (1962),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=One can use the diffeo r -> -r and u -> -u to make the assumption r > 0,The case _m = 0 is Stephani, [28, 16,1].,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,26,3=Authors=Robinson, Trautman (1962),PrimaryDescription=Vacuum,SecondaryDescription=RobinsonTrautman,Comments=One can use the diffeo r -> -r and u -> -u to make the assumption r > 0,The case _m = 0 is Stephani, [28, 16,1].,AlternativeOrthonormalTetrad1 and AlternativeNullTetrad1 are adapted to the shear-free null geodesic congruence (Robinson-Trautman tetrads)
28,37,1=Authors=Robinson-Trautman (1962),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=RobinsonTrautman,Comments=diff(2*P^2*diff(2*P(u,z,zb)^2*ln(P(u,z,zb)),z,zb),z,zb)+12*_m(u)*diff(ln(P(u,z,zb)),u)-4*diff(_m(u),u)=4*_kappa0*P^2*_h(u,z,zb)*_hb(u,z,zb), _Q(u,z,zb)*diff(_Qb(u,z,zb),u)-_Qb(u,z,zb)*diff(_Q(u,z,zb),u)=2*_P(u,z,zb)^2*(_hb(u,z,zb)*diff(_Qb(u,z,zb),zb)-_h(u,z,zb)*diff(_Q(u,z,zb),b)), diff(_h(u,z,zb),z)=diff(_Qb(u,z,zb)/2/_P(u,z,zb)^2,u), diff(_m(u,z,zb),z)=_kappa0*_hb(u,z,zb)*_Qb(u,z,zb)
28,43,1=Authors=Robinson, Trautman (1962),PrimaryDescription=EinsteinMaxwell,SecondaryDescription=PureRadiation,RobinsonTrautman,Comments=h1(u) is the conjugate of h(u)
31,40,1=Authors= Petrov (1962),PrimaryDescription=Vacuum,SecondaryDescription=Non-diverging
35,19,1=Authors=Ehlers and Kundt 1962,PrimaryDescription=Generic,SecondaryDescription=PPWave,Comments=constant null bivector, _k[mu]=-D_[mu](u)
The search can also be done visually, by properties. The following is the only solution in the database that is a Pure Radiation solution, of Petrov Type "D", Plebanski-Petrov Type "O" and that has Isometry Dimension equal to 1:
DifferentialGeometry:-Library:-MetricSearch
Set the solution, and everything related to work with it, in one go
g_28,74,1
Systems of spacetime Coordinates are: X=u,η,r,y
Default differentiation variables for d_, D_ and dAlembertian are: X=u,η,r,y
The Frolov and Khlebnikov (1975) metric in coordinates u,η,r,y
Parameters: κ0,m⁡u,b,d
Comments: Wⅈth m(u) = constant, thⅇ mⅇtrⅈc ⅈs Rⅈccⅈ flat anⅆ bⅇcomⅇs 28.24 ⅈn Stⅇphanⅈ.
Resetting the signature of spacetime from "+ - - -" to `- + + +` in order to match the signature in the database of metrics:
The related Riemann invariants:
r0=0,r1=0,r2=0,r3=0,w1=6⁢m⁡u2r6,w2=−6⁢m⁡u3r9,m1=0,m2=0,m3=0,m4=0,m5=0
To list all the triads of numbers associated to the solutions digitized in the database (each triad indicates the Chapter and equation number with which each of these metrics appear in the book), enter
DifferentialGeometry:-Library:-RetrieveStephani,1
8,33,1,8,34,1,12,6,1,12,7,1,12,8,1,12,8,2,12,8,3,12,8,4,12,8,5,12,8,6,12,8,7,12,8,8,12,9,1,12,9,2,12,9,3,12,9,4,12,9,5,12,9,6,12,9,7,12,12,1,12,12,2,12,12,3,12,12,4,12,13,1,12,14,1,12,16,1,12,18,1,12,19,1,12,21,1,12,23,1,12,23,2,12,23,3,12,24.1,1,12,24.2,1,12,24.3,1,12,26,1,12,27,1,12,28,1,12,29,1,12,30,1,12,31,1,12,32,1,12,34,1,12,35,1,12,36,1,12,37,1,12,37,2,12,37,3,12,37,4,12,37,5,12,37,6,12,37,7,12,38,1,12,38,2,12,38,3,12,38,4,12,38,5,13,1,1,13,1,2,13,1,3,13,2,1,13,2,2,13,2,3,13,3,1,13,7,1,13,7,2,13,7,3,13,7,4,13,7,5,13,7,6,13,7,7,13,7,8,13,9,1,13,14,1,13,14,2,13,14,3,13,15,1,13,15,2,13,15,3,13,15,4,13,17,1,13,19,1,13,22,1,13,31,1,13,32,1,13,35,1,13,46,1,13,48,1,13,49,1,13,49,2,13,51,1,13,53,1,13,55,1,13,56,1,13,57,1,13,58,1,13,59,1,13,59,2,13,60,1,13,60,2,13,60,3,13,60,4,13,60,5,13,60,6,13,60,7,13,60,8,13,61,1,13,61,2,13,62,1,13,62,2,13,62,4,13,62,6,13,63,1,13,63,2,13,63,3,13,63,4,13,64,1,13,64,2,13,64,3,13,64,4,13,65,1,13,67,1,13,67,2,13,69,1,13,71,1,13,72,1,13,73,1,13,74,1,13,74,2,13,74,3,13,76,1,13,77,1,13,77,2,13,79,1,13,79,2,13,80,1,13,81,1,13,83,1,13,84,1,13,84,2,13,84,3,13,85,1,13,85,2,13,86,1,13,87,1,14,6.1,1,14,6.2,1,14,6.3,1,14,7,1,14,8.1,1,14,8.2,1,14,8.3,1,14,9.1,1,14,9.2,1,14,10,1,14,10,2,14,12,1,14,12,2,14,12,3,14,14,1,14,14,2,14,15,1,14,15.1,2,14,15.2,2,14,15.3,2,14,16,1,14,16,2,14,17,1,14,18,1,14,18,2,14,19,1,14,20,1,14,21,1,14,21,2,14,21,3,14,22,1,14,23,1,14,24,1,14,25,1,14,26,1,14,26,2,14,26,3,14,26,4,14,27,1,14,28,1,14,28,2,14,28,3,14,29,1,14,30,1,14,31,1,14,32,1,14,33,1,14,35,1,14,37,1,14,38,1,14,38,2,14,38,3,14,39,1,14,39,2,14,39,3,14,39,4,14,39,5,14,39,6,14,40,1,14,41,1,14,42,1,14,46,1,15,3,1,15,3,2,15,4,1,15,4,2,15,4,3,15,9,1,15,10,1,15,12,1,15,12,2,15,12,3,15,12,4,15,12,5,15,12,6,15,17,1,15,17,2,15,17,3,15,17,4,15,18,1,15,19,1,15,19,2,15,20,1,15,21,1,15,21,2,15,22,1,15,23,1,15,23,2,15,24,1,15,24,2,15,25,1,15,25,2,15,26,1,15,26,2,15,27,1,15,27,2,15,27,3,15,27,4,15,27,5,15,27,6,15,27,7,15,27,8,15,28,1,15,29,1,15,30,1,15,31,1,15,32,1,15,34,1,15,34,2,15,34,3,15,43,1,15,43,2,15,43,3,15,50,1,15,50,2,15,50,3,15,50,4,15,50,5,15,50,6,15,62,1,15,62,2,15,62,3,15,63,1,15,63,2,15,63,3,15,65,1,15,65,2,15,66,1,15,66,2,15,66,3,15,75,1,15,75,2,15,75,3,15,77,1,15,77,2,15,77,3,15,78,1,15,79,1,15,81,1,15,81,2,15,81,3,15,82,1,15,82,2,15,82,3,15,83,1.1,15,83,1.2,15,83,2,15,83,3.1,15,83,3.2,15,83,4,15,84,1,15,85,1,15,85,2,15,85,3,15,86,1,15,86,2,15,86,3,15,87,1,15,87,2,15,87,3,15,87,4,15,87,5,15,88,1,15,89,1,15,90,1,16,1,1,16,1,2,16,1,3,16,1,4,16,1,5,16,1,6,16,1,7,16,1,8,16,1,9,16,1,10,16,1,11,16,1,12,16,1,13,16,1,14,16,1,15,16,1,16,16,1,17,16,1,18,16,1,19,16,1,20,16,1,21,16,1,22,16,1,23,16,1,24,16,1,25,16,1,26,16,1,27,16,14,1,16,14,2,16,14,3,16,14,4,16,14,5,16,14,6,16,14,7,16,14,8,16,14,9,16,14,10,16,14,11,16,14,12,16,14,13,16,14,14,16,14,15,16,14,16,16,14,17,16,14,18,16,14,19,16,14,20,16,18,1,16,19,1,16,20,1,16,22,1,16,24,1,16,24,2,16,43,1,16,45,1,16,45,2,16,46,1,16,47,1,16,50,1,16,51,1,16,54,1,16,61,1,16,63,1,16,66,1,16,66,2,16,66,3,16,67,1,16,71,1,16,72,1,16,73,1,16,74,1,16,75,1,16,76,1,16,77,1,16,77,2,16,77,3,16,78,1,17,4,1,17,4,2,17,5,1,17,9,1,17,14,1,17,15,1,17,15,2,17,16,1,17,20,1,17,22,1,17,23,1,17,24,1,17,24,2,17,26,1,17,27,1,17,27,2,17,28,1,17,28,2,17,29,1,17,29,2,17,30,1,17,31,1,18,2,1,18,2,2,18,2,3,18,2,4,18,2,5,18,2,6,18,2,7,18,2,8,18,48,1,18,48,2,18,49,1,18,50,1,18,64,1,18,64,2,18,64,3,18,65,1,18,66,1,18,67,1,18,71,1,18,75,1,19,17,1,19,17,2,19,21,1,20,3,1,20,4,1,20,5,1,20,7,1,20,8,1,20,9,1,20,10,1,20,11,1,20,12,1,20,13,1,20,15,1,20,16,1,20,17,1,20,20,1,20,21,1,20,23,1,20,27,1,20,28,1,20,29,1,20,32,1,20,34,1,20,36,1,20,38,1,20,38,2,20,38,3,20,44,1,20,46,1,20,54,1,20,57,1,20,57,2,21,1,1,21,1,2,21,1,3,21,4,1,21,5,1,21,6,1,21,7,1,21,10,1,21,10,2,21,11,1,21,16,1,21,17,1,21,17,2,21,20,1,21,22,1,21,22,2,21,24,1,21,28,1,21,30,1,21,30,2,21,30,3,21,31,1,21,35,1,21,41,1,21,52,1,21,57,1,21,58,1,21,59,1,21,60,1,21,61,1,21,61,2,21,61,3,21,61,4,21,61,5,21,70,1,21,71,1,21,72,1,21,73,1,21,74,1,21,74,2,21,74,3,21,74,4,22,3,1,22,4,1,22,4,2,22,5,1,22,6,1,22,6,2,22,7,1,22,8,1,22,8,2,22,8,3,22,8,4,22,8,5,22,11,1,22,12,1,22,13,1,22,14,1,22,15,1,22,16,1,22,17,1,22,18,1,22,18,2,22,19,1,22,22,1,22,23,1,22,24,1,22,27,1,22,28,1,22,29,1,22,34,1,22,34,2,22,34,3,22,34,4,22,34,5,22,47,1,22,48,1,22,49,1,22,50,1,22,51,1,22,52,1,22,53,1,22,59,1,22,63,1,22,64,1,22,67,1,22,67,2,22,70,1,23,1,1,23,2,1,23,2,2,23,2,3,23,2,4,23,3,1,23,4,1,23,6,1,23,6,2,23,7,1,23,8,1,23,8,2,23,8,3,23,11,1,23,12,1,23,13,1,23,13,2,23,13,3,23,13,4,23,14,1,23,14,2,23,14,3,23,14,4,23,14,5,23,14,6,23,14,7,23,15,1,23,15,2,23,15,3,23,15,4,23,15,5,23,16,1,23,16,2,23,16,3,23,16,4,23,16,5,23,17,1,23,17,2,23,18,1,23,18,2,23,19,1,23,19,2,23,20,1,23,20,2,23,20,3,23,21,1,23,21,2,23,22,1,23,23,1,23,24,1,23,25,1,23,26,1,23,26,2,23,26,3,23,26,4,23,26,5,23,26,6,23,26,7,23,26,8,23,30,1,23,31,1,23,32,1,23,32,2,23,32,3,23,33,1,23,33,2,23,33,3,23,34,1,23,35,1,23,36,1,23,36,2,23,36,3,23,37,1,23,38,1,23,38,2,23,38,3,23,39,1,23,39,2,23,40,1,23,40,2,23,40,3,23,40,4,23,41,1,23,41,2,23,41,3,23,42,1,23,43,1,23,44,1,23,44,2,23,45,1,23,46,1,23,47,1,23,48,1,23,48,2,23,49,1,23,50,1,23,51,1,23,52,1,23,53,1,24,2,1,24,2,2,24,2,3,24,2,4,24,2,5,24,2,6,24,2,7,24,2,8,24,2,9,24,2,10,24,2,11,24,21,1,24,22,1,24,28,1,24,35,1,24,37,1,24,37,2,24,37,3,24,37,4,24,37,5,24,37,6,24,37,7,24,37,8,24,37,9,24,38,1,24,38,2,24,40,1,24,46,1,24,46,2,24,46,3,24,47,1,24,51,1,25,2,1,25,5,1,25,5,2,25,6,1,25,9&a