| 24小時熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 604 | 回復(fù): 4 | ||
shuizl木蟲 (正式寫手)
|
[求助]
1stopt 的運行結(jié)果之輸出結(jié)果全部為0 求助
|
|
1stopt 的運行結(jié)果之輸出結(jié)果全部為0 源代碼: Title "AC2O"; Parameters k_1[0,1e18],k_2[0,1e18],k_3[0,1e18],k_4[0,1e18],k_5[0,1e18],k_6[0,1e18],E1[2e4,2e5],E2[2e4,2e5],E3[2e4,2e5],E4[2e4,2e5],E5[2e4,2e5],E6[2e4,2e5]; Variable TT,PDMEin,CCH3Iin,CRhin,PCOin,CLiIin; Variable C1out[Output],C2out[Output]; StartProgram [VB]; Dim iter As Integer Dim RK As Integer Dim j As Integer Dim RK_i As Integer Dim C_1 As Double Dim C_2 As Double Dim k_1 As Double Dim k_2 As Double Dim k_3 As Double Dim k_4 As Double Dim k_5 As Double Dim k_6 As Double Dim E1 As Double Dim E2 As Double Dim E3 As Double Dim E4 As Double Dim E5 As Double Dim E6 As Double Dim ti As Double Dim T As Double Dim TT As Double Dim PDMEin As Double Dim CCH3Iin As Double Dim CRhin As Double Dim PCOin As Double Dim CLiIin As Double Dim C1out As Double Dim C2out As Double Dim PDMEi As Double Dim CCH3I As Double Dim CRh As Double Dim PCO As Double Dim CLiI As Double Dim ee1 As Double Dim ee2 As Double Dim ee3 As Double Dim ee4 As Double Dim ee5 As Double Dim x_x1 As Double Dim x_x2 As Double Dim y_y1 As Double Dim y_y2 As Double Dim C1 As Double Dim C2 As Double Dim K1 As Double Dim K2 As Double Dim K3 As Double Dim K4 As Double Dim K5 As Double Dim K6 As Double Dim d_c1 As Double Dim d_c2 As Double Dim c As Double For RK=0 to Datalength-1 T=TT(iter) PDMEi=PDMEin(iter) CCH3I=CCH3Iin(iter) CRh=CRhin(iter) PCO=PCOin(iter) CLiI=CLiIin(iter) C_1=0; C_2=0; ti=6/20; For RK_i=1 to 20 ee1=0.5 ee2=0.5 ee3=1 ee4=1 ee5=0.5 x_x1=C_1 x_x2=C_2 y_y1=C_1 y_y2=C_2 For j=1 To 4 C1=x_x1 C2=x_x2 K1=k_1*exp(-E1/T/8.314) K2=k_2*exp(-E2/T/8.314) K3=k_3*exp(-E3/T/8.314) K4=k_4*exp(-E4/T/8.314) K5=k_5*exp(-E5/T/8.314) K6=k_6*exp(-E6/T/8.314) d_c1=K1*(PDMEi-8.314*T*(C1+C2))^2*CCH3I^2*CRh*PCO/PDMEi^3-K2*(PDMEi-8.314*T*(C1+C2))*CCH3I*CRh*C1/PDMEi^2 d_c2=K3*(PDMEi-8.314*T*(C1+C2))*CRh*PCO*C2/PDMEi/C1-K4*CRh*C2+K5*CCH3I*C1^2*PDMEi/(PDMEi-8.314*T*(C1+C2))-K6*CCH3I*C1*C2*CRh*PDMEi/CLiI/(PDMEi-8.314*T*(C1+C2)) c=d_c1*ti if j=1 then x_x1=C_1+ee1*c y_y1=y_y1+ee2*c/3 elseif j=2 then x_x1=C_1+ee2*c y_y1=y_y1+ee3*c/3 elseif j=3 then x_x1=C_1+ee3*c y_y1=y_y1+ee4*c/3 elseif j=4 then x_x1=C_1+ee4*c y_y1=y_y1+ee5*c/3 end if c=d_c2*ti if j=1 then x_x2=C_2+ee1*c y_y2=y_y2+ee2*c/3 elseif j=2 then x_x2=C_2+ee2*c y_y2=y_y2+ee3*c/3 elseif j=3 then x_x2=C_2+ee3*c y_y2=y_y2+ee4*c/3 elseif j=4 then x_x2=C_2+ee4*c y_y2=y_y2+ee5*c/3 end if Next C_1=y_y1 C_2=y_y2 Next C1out(iter)=y_y1 C2out(iter)=y_y2 Next ENDPROGRAM; Data; 453.00 150649.68 126.76 0.41 3049350.32 37.31 19.32 8.80 453.00 150649.68 126.76 0.83 3049350.32 37.31 13.00 19.28 453.00 150649.68 126.76 1.24 3049350.32 37.31 9.20 26.32 453.00 150649.68 126.76 1.66 3049350.32 37.31 8.72 29.52 453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82 453.00 150649.68 126.76 2.48 3049350.32 37.31 6.95 30.41 453.00 150649.68 126.76 2.90 3049350.32 37.31 5.59 26.17 453.00 150649.68 42.25 2.07 3049350.32 37.31 6.85 13.43 453.00 150649.68 70.42 2.07 3049350.32 37.31 8.57 20.71 453.00 150649.68 84.51 2.07 3049350.32 37.31 10.04 26.91 453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82 453.00 150649.68 169.01 2.07 3049350.32 37.31 6.47 29.89 453.00 150649.68 211.27 2.07 3049350.32 37.31 5.34 29.84 453.00 150649.68 253.52 2.07 3049350.32 37.31 4.40 29.96 453.00 150649.68 126.76 2.07 3049350.32 7.46 12.08 25.16 453.00 150649.68 126.76 2.07 3049350.32 14.93 10.42 27.91 453.00 150649.68 126.76 2.07 3049350.32 22.39 8.96 29.84 453.00 150649.68 126.76 2.07 3049350.32 29.85 8.45 30.95 453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82 453.00 150649.68 126.76 2.07 3049350.32 44.78 7.75 32.12 453.00 150649.68 126.76 2.07 3049350.32 52.24 7.17 32.55 453.00 75324.84 126.76 2.07 3124675.16 37.31 16.12 3.88 453.00 112987.26 126.76 2.07 3087012.74 37.31 12.25 17.75 453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82 453.00 188312.10 126.76 2.07 3011687.90 37.31 8.73 41.09 453.00 225974.52 126.76 2.07 2974025.48 37.31 9.48 50.40 453.00 263636.94 126.76 2.07 2936363.06 37.31 9.44 52.55 453.00 301299.36 126.76 2.07 2898700.64 37.31 7.60 43.20 453.00 150649.68 126.76 2.07 2349350.32 37.31 12.85 13.63 453.00 150649.68 126.76 2.07 2599350.32 37.31 12.80 19.71 453.00 150649.68 126.76 2.07 2849350.32 37.31 11.47 26.70 453.00 150649.68 126.76 2.07 3099350.32 37.31 7.93 31.95 453.00 150649.68 126.76 2.07 3349350.32 37.31 4.74 35.13 453.00 150649.68 126.76 2.07 3599350.32 37.31 4.02 35.82 453.00 150649.68 126.76 2.07 3849350.32 37.31 2.99 36.91 433.00 143998.48 126.76 2.07 3056001.52 37.31 13.52 18.40 438.00 145661.28 126.76 2.07 3104338.72 37.31 12.55 20.47 443.00 147324.08 126.76 2.07 3052675.92 37.31 11.55 22.17 448.00 148986.88 126.76 2.07 3101013.12 37.31 9.59 28.23 453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82 458.00 152312.48 126.76 2.07 3097687.52 37.31 6.04 31.06 463.00 153975.28 126.76 2.07 3046024.72 37.31 4.40 27.24 運行結(jié)果: "AC2O" 迭代數(shù): 16 計算用時(時:分:秒:微秒): 00:00:00:322 優(yōu)化算法: 準(zhǔn)牛頓法(BFGS) + 通用全局優(yōu)化法 計算結(jié)束原因: 達到收斂判定標(biāo)準(zhǔn) 均方差(RMSE): 22.1009039084701 殘差平方和(SSE): 41029.7961 相關(guān)系數(shù)(R): 0 相關(guān)系數(shù)之平方(R^2): 0 決定系數(shù)(DC): 8.8666753259632E-16 F統(tǒng)計(F-Statistic): 0 參數(shù) 最佳估算 -------------------- ------------- k_1 1.5926178987138E17 k_2 5.34790054196492E17 k_3 5.79585216706619E17 k_4 5.18475014017895E17 k_5 5.73410185752437E17 k_6 1.99592327931896E17 E1 76320.0160069391 E2 34522.3072776571 E3 0.659441070165485 E4 1.32150601642206 E5 67320.5888597295 E6 20693.1685144082 ====== 輸出結(jié)果 ===== No 實測C1out 計算C1out 實測C2out 計算C2out 1 19.32 0 8.8 0 2 13 0 19.28 0 3 9.2 0 26.32 0 4 8.72 0 29.52 0 5 8.06 0 31.82 0 6 6.95 0 30.41 0 7 5.59 0 26.17 0 8 6.85 0 13.43 0 9 8.57 0 20.71 0 10 10.04 0 26.91 0 11 8.06 0 31.82 0 12 6.47 0 29.89 0 13 5.34 0 29.84 0 14 4.4 0 29.96 0 15 12.08 0 25.16 0 16 10.42 0 27.91 0 17 8.96 0 29.84 0 18 8.45 0 30.95 0 19 8.06 0 31.82 0 20 7.75 0 32.12 0 21 7.17 0 32.55 0 22 16.12 0 3.88 0 23 12.25 0 17.75 0 24 8.06 0 31.82 0 25 8.73 0 41.09 0 26 9.48 0 50.4 0 27 9.44 0 52.55 0 28 7.6 0 43.2 0 29 12.85 0 13.63 0 30 12.8 0 19.71 0 31 11.47 0 26.7 0 32 7.93 0 31.95 0 33 4.74 0 35.13 0 34 4.02 0 35.82 0 35 2.99 0 36.91 0 36 13.52 0 18.4 0 37 12.55 0 20.47 0 38 11.55 0 22.17 0 39 9.59 0 28.23 0 40 8.06 0 31.82 0 41 6.04 0 31.06 0 42 4.4 0 27.24 0 |

鐵桿木蟲 (職業(yè)作家)
木蟲 (正式寫手)
|
Title "AC2O"; Parameters k_1[0,1e180],k_2[0,1e180],k_3[0,1e180],k_4[0,1e180],k_5[0,1e180],k_6[0,1e180],E1[2e4,2e5],E2[2e4,2e5],E3[2e4,2e5],E4[2e4,2e5],E5[2e4,2e5],E6[2e4,2e5]; Variable TT,PDMEin,CCH3Iin,CRhin,PCOin,CLiIin; Variable C1out[Output],C2out[Output]; StartProgram [VB]; Dim iter As Integer Dim RK As Integer Dim j As Integer Dim RK_i As Integer Dim C_1 As Double Dim C_2 As Double Dim k_1 As Double Dim k_2 As Double Dim k_3 As Double Dim k_4 As Double Dim k_5 As Double Dim k_6 As Double Dim E1 As Double Dim E2 As Double Dim E3 As Double Dim E4 As Double Dim E5 As Double Dim E6 As Double Dim ti As Double Dim T As Double Dim TT As Double Dim PDMEin As Double Dim CCH3Iin As Double Dim CRhin As Double Dim PCOin As Double Dim CLiIin As Double Dim C1out As Double Dim C2out As Double Dim PDMEi As Double Dim CCH3I As Double Dim CRh As Double Dim PCO As Double Dim CLiI As Double Dim ee1 As Double Dim ee2 As Double Dim ee3 As Double Dim ee4 As Double Dim ee5 As Double Dim x_x1 As Double Dim x_x2 As Double Dim y_y1 As Double Dim y_y2 As Double Dim C1 As Double Dim C2 As Double Dim K1 As Double Dim K2 As Double Dim K3 As Double Dim K4 As Double Dim K5 As Double Dim K6 As Double Dim d_c1 As Double Dim d_c2 As Double Dim c As Double For RK=0 to Datalength-1 T=TT(iter) PDMEi=PDMEin(iter) CCH3I=CCH3Iin(iter) CRh=CRhin(iter) PCO=PCOin(iter) CLiI=CLiIin(iter) C_1=0; C_2=0; ti=6/20; For RK_i=1 to 20 ee1=0.5 ee2=0.5 ee3=1 ee4=1 ee5=0.5 x_x1=C_1 x_x2=C_2 y_y1=C_1 y_y2=C_2 For j=1 To 4 C1=x_x1 C2=x_x2 K1=k_1*exp(-E1/T/8.314) K2=k_2*exp(-E2/T/8.314) K3=k_3*exp(-E3/T/8.314) K4=k_4*exp(-E4/T/8.314) K5=k_5*exp(-E5/T/8.314) K6=k_6*exp(-E6/T/8.314) d_c1=K1*(PDMEi-8.314*T*(C1+C2))^2*CCH3I^2*CRh*PCO/PDMEi^3-K2*(PDMEi-8.314*T*(C1+C2))*CCH3I*CRh*C1/PDMEi^2 c=d_c1*ti if j=1 then x_x1=C_1+ee1*c y_y1=y_y1+ee2*c/3 elseif j=2 then x_x1=C_1+ee2*c y_y1=y_y1+ee3*c/3 elseif j=3 then x_x1=C_1+ee3*c y_y1=y_y1+ee4*c/3 elseif j=4 then x_x1=C_1+ee4*c y_y1=y_y1+ee5*c/3 end if C1=x_x1 d_c2=K3*(PDMEi-8.314*T*(C1+C2))*CRh*PCO*C2/PDMEi/C1-K4*CRh*C2+K5*CCH3I*C1^2*PDMEi/(PDMEi-8.314*T*(C1+C2))-K6*CCH3I*C1*C2*CRh*PDMEi/CLiI/(PDMEi-8.314*T*(C1+C2)) c=d_c2*ti if j=1 then x_x2=C_2+ee1*c y_y2=y_y2+ee2*c/3 elseif j=2 then x_x2=C_2+ee2*c y_y2=y_y2+ee3*c/3 elseif j=3 then x_x2=C_2+ee3*c y_y2=y_y2+ee4*c/3 elseif j=4 then x_x2=C_2+ee4*c y_y2=y_y2+ee5*c/3 end if Next C_1=y_y1 C_2=y_y2 Next C1out(iter)=y_y1 C2out(iter)=y_y2 Next ENDPROGRAM; Data; 453.00 150649.68 126.76 0.41 3049350.32 37.31 19.32 8.80 453.00 150649.68 126.76 0.83 3049350.32 37.31 13.00 19.28 453.00 150649.68 126.76 1.24 3049350.32 37.31 9.20 26.32 453.00 150649.68 126.76 1.66 3049350.32 37.31 8.72 29.52 453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82 453.00 150649.68 126.76 2.48 3049350.32 37.31 6.95 30.41 453.00 150649.68 126.76 2.90 3049350.32 37.31 5.59 26.17 453.00 150649.68 42.25 2.07 3049350.32 37.31 6.85 13.43 453.00 150649.68 70.42 2.07 3049350.32 37.31 8.57 20.71 453.00 150649.68 84.51 2.07 3049350.32 37.31 10.04 26.91 453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82 453.00 150649.68 169.01 2.07 3049350.32 37.31 6.47 29.89 453.00 150649.68 211.27 2.07 3049350.32 37.31 5.34 29.84 453.00 150649.68 253.52 2.07 3049350.32 37.31 4.40 29.96 453.00 150649.68 126.76 2.07 3049350.32 7.46 12.08 25.16 453.00 150649.68 126.76 2.07 3049350.32 14.93 10.42 27.91 453.00 150649.68 126.76 2.07 3049350.32 22.39 8.96 29.84 453.00 150649.68 126.76 2.07 3049350.32 29.85 8.45 30.95 453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82 453.00 150649.68 126.76 2.07 3049350.32 44.78 7.75 32.12 453.00 150649.68 126.76 2.07 3049350.32 52.24 7.17 32.55 453.00 75324.84 126.76 2.07 3124675.16 37.31 16.12 3.88 453.00 112987.26 126.76 2.07 3087012.74 37.31 12.25 17.75 453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82 453.00 188312.10 126.76 2.07 3011687.90 37.31 8.73 41.09 453.00 225974.52 126.76 2.07 2974025.48 37.31 9.48 50.40 453.00 263636.94 126.76 2.07 2936363.06 37.31 9.44 52.55 453.00 301299.36 126.76 2.07 2898700.64 37.31 7.60 43.20 453.00 150649.68 126.76 2.07 2349350.32 37.31 12.85 13.63 453.00 150649.68 126.76 2.07 2599350.32 37.31 12.80 19.71 453.00 150649.68 126.76 2.07 2849350.32 37.31 11.47 26.70 453.00 150649.68 126.76 2.07 3099350.32 37.31 7.93 31.95 453.00 150649.68 126.76 2.07 3349350.32 37.31 4.74 35.13 453.00 150649.68 126.76 2.07 3599350.32 37.31 4.02 35.82 453.00 150649.68 126.76 2.07 3849350.32 37.31 2.99 36.91 433.00 143998.48 126.76 2.07 3056001.52 37.31 13.52 18.40 438.00 145661.28 126.76 2.07 3104338.72 37.31 12.55 20.47 443.00 147324.08 126.76 2.07 3052675.92 37.31 11.55 22.17 448.00 148986.88 126.76 2.07 3101013.12 37.31 9.59 28.23 453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82 458.00 152312.48 126.76 2.07 3097687.52 37.31 6.04 31.06 463.00 153975.28 126.76 2.07 3046024.72 37.31 4.40 27.24 |

木蟲 (正式寫手)

鐵桿木蟲 (職業(yè)作家)
|
改動如下,可以正常運行,效果不好,不知模型邏輯構(gòu)成方面是否正確,自己慢慢改動運算: Root of Mean Square Error (RMSE):11.5499416460334 Sum of Square Error:11205.6967702493 Correlation Coef. (R): -0.025193150477488 R-Square: 0.000634694830981356 Determination Coef. (DC): -3.6524379110701 F-Statistic: 0.379799714772977 Parameter Best Estimate -------------------- ------------- k_1 0.0443887227013897 k_2 1.48988141904544 k_3 0.15806362988123 k_4 0.551727467459332 k_5 1.27971229562335 k_6 3.90080527521901 e1 34064.6206926927 e2 33558.1459707385 e3 46339.7668465127 e4 32936.9747130519 e5 33108.0993428242 e6 32934.5560026088 ====== Output Results ====== File: Data File-1 No Target c1out Calculated c1out Target c2out Calculated c2out 1 19.32 3.72446132159712 8.8 1.45941820895285 2 13 6.53507135688917 19.28 5.43948599766212 3 9.2 8.3844585821214 26.32 10.9146624801039 4 8.72 9.57375357447746 29.52 17.4176359596032 5 8.06 10.2830886011808 31.82 24.8297659261248 6 6.95 10.8024412001823 30.41 27.2516499967275 7 5.59 11.3555309700269 26.17 39.278137463255 8 6.85 0.481409991498034 13.43 34.6582295922687 9 8.57 2.98065943110798 20.71 34.7680122723834 10 10.04 5.92518851451692 26.91 18.8384744557624 11 8.06 10.2830886011808 31.82 24.8297659261248 12 6.47 11.9449557638363 29.89 26.0025570402153 13 5.34 39.9686546702136 29.84 23.268772170038 14 4.4 25.6463940788978 29.96 36.9515358726961 15 12.08 11.2082725915566 25.16 12.5511124159842 16 10.42 10.6780692280548 27.91 17.9381468602351 17 8.96 10.4668996961607 29.84 20.9568431588697 18 8.45 10.3536541774117 30.95 23.0860941405943 19 8.06 10.2830886011808 31.82 24.8297659261248 20 7.75 10.2348920489271 32.12 26.4876454219494 21 7.17 10.2000642722205 32.55 28.5810076025036 22 16.12 10.2439751831954 3.88 6.01419821549014 23 12.25 9.95652354737421 17.75 15.1345857359885 24 8.06 10.2830886011808 31.82 24.8297659261248 25 8.73 10.1175326870607 41.09 20.3228972672854 26 9.48 9.44413185248538 50.4 19.1642620585202 27 9.44 8.58943190637761 52.55 19.8006045683615 28 7.6 7.69302617652497 43.2 22.6514804670537 29 12.85 9.45991413360167 13.63 15.8282959515173 30 12.8 9.81660871068178 19.71 18.6483271853985 31 11.47 10.0999291047028 26.7 21.769488046076 32 7.93 10.3237446946982 31.95 25.8185269212792 33 4.74 10.532154045967 35.13 35.2973893006523 34 4.02 10.7334779879877 35.82 44.776189116222 35 2.99 32.9322212693294 36.91 83.6274114464888 36 13.52 9.63412367023005 18.4 8.29143322070249 37 12.55 10.0548474744031 20.47 11.3109586906801 38 11.55 10.2064105002571 22.17 14.4464803808188 39 9.59 10.3619096378366 28.23 19.3200897771154 40 8.06 10.2830886011808 31.82 24.8297659261248 41 6.04 10.3746296215376 31.06 36.6337375772535 42 4.4 10.0909841087201 27.24 30.7266741146468 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +9 | 天天好運來上岸?/a> 2026-03-24 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 求調(diào)劑一志愿海大,0703化學(xué)學(xué)碩304分,有大創(chuàng)項目,四級已過 +7 | 幸運哩哩 2026-03-22 | 11/550 |
|
|
[考研] 071000生物學(xué)求調(diào)劑,初試成績343 +7 | 小小甜面團 2026-03-25 | 7/350 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 085600,材料與化工321分求調(diào)劑 +9 | 大饞小子 2026-03-28 | 9/450 |
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-26 | 3/150 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[基金申請] 請教下大家 2026年國家基金申請是雙盲審嗎? +3 | lishucheng1 2026-03-22 | 5/250 |
|