Mathematica的RootReduce语句能不能优化

hejoseph

这段Mathematica代码最后一句是找一个多面体体积所满足的方程，运行很长时间都没结果，能不能优化？

1. b=Root[65536 #1^10-270336 #1^9+387328 #1^8-207872 #1^7-3328 #1^6+68352 #1^5-47008 #1^4-10560 #1^3+12976 #1^2+4464 #1+345&,4]
   
2. a=((2 b-1)^2-2)/((2 b-1)^2-4)
   
3. c=(8388608000 b^21-79586918400 b^20+291095183360 b^19-470908928000 b^18+201526083584 b^17+159530680320 b^16+493590937600 b^15-1485855981568 b^14+815036301312 b^13+601976651776 b^12-504076419072 b^11-148765050880 b^10-9481250816 b^9+102792349184 b^8+95526461440 b^7-36947707648 b^6-42856999168 b^5-301373312 b^4+7127721504 b^3+2029968900 b^2+162185892 b-2537283)/(15099494400 b^20-146360238080 b^19+548751278080 b^18-886463004672 b^17+109206568960 b^16+1666276392960 b^15-2054359089152 b^14+299851579392 b^13+788277886976 b^12-202246488064 b^11-41992994816 b^10-193101914112 b^9-30723567616 b^8+149127761920 b^7+25154799616 b^6-40151213056 b^5-11388126720 b^4+3264281152 b^3+1620897568 b^2+164732976 b-1566984)
   
4. RootReduce[(1+a) Sqrt[1-a^2]+(1+a) (2 b-1)/3 Sqrt[1-a^2]+(4 a+3)/12 Sqrt[3-4 a^2]+(1+a)/3 Sqrt[2 (1+a) (1-2 a)]+(1+2 a)/6 Sqrt[2 (1-2 a-2 a^2)]]

复制代码

mathematica

什么注释都没有，这样的代码谁愿意看？

hejoseph

不愿意看你回复干什么？Mathematica里面的命令都很清楚了，我只想能不能在现有的命令基础上优化。

nyy

通常我都是先计算高精度数值解，然后RootApproximant来得到解析解，
如果解析解能写成代数数的话，通常得到的答案都是正确的。

有时候软件太笨了，化简不过来！