第一个方程mma能解 In[12]:= Solve[a==2207-1/a,a] ` 2 2207 + 987 Sqrt[5] `Out[12]= {{a -> ------------------}, {a -> ------------------}} ` 2207 + 987 Sqrt[5] 2

用mathematica解不出那个8次方程 Solve[((a+b*Sqrt[c]) /d) ^8 == (2207+987*Sqrt[5] )/2, {a,b,c,d},Integers]

wns3的通项公式能用于w(n,6,4) >>> for i in range(0,7):set_a=set();sortlst_a=[];x=[];xflag=[];pa=[];lenx=0;cnt=0;print(i+4,w3(4+i,6,3)) ... 4 165348 5 507708 6 1214076 <-从这行开始匹配 7 2478846 8 4537896 9 7668570 10 12189684 In[8]:= f[n_]:=wns3[n,6] In[13]:= f/@Range[4,10] Out[13]= {165360, 507714, 1214076, 2478846, 4537896, 7668570, 12189684}

如果扩大c=10**(23.135*3)，U(c,152)

若使c=sqrt(c),连乘积的个数要少1半 c=10**23.135 cnt=0 t=10**5 for i in range(t): if U(c,49)<1:cnt+=1 print(float(cnt)/t)