[已解决]代码 怎么运算不起来? 望修正

caipiaofans

香川群子 发表于 2014-12-3 08:56
原问题帖的方程式为：
(1579-K1)*(1579-K1) + (5803-K2)*(5803-K2) + (1137-K3)*(1137-K3) = 217502634
...

s1,s2,s3未知

^2 是平方的意思

整体公式是这样
a3^2+ a2^2 + a1^2 + s1^2 + s2^2 + s3^2 - 2*(a3*s3 + a2*s2 + a1*s1) = 217502634

整理成这样:
(a1-s1)*(a1-s1) + (a2-s2)*(a2-s2) + (a3-s3)*(a3-s3) = 217502634

应该没错吧...


不然直接符合这个公式的运算也可以
a3^2+ a2^2 + a1^2 + s1^2 + s2^2 + s3^2 - 2*(a3*s3 + a2*s2 + a1*s1) = 217502634

caipiaofans

香川群子 发表于 2014-12-3 08:56
原问题帖的方程式为：
(1579-K1)*(1579-K1) + (5803-K2)*(5803-K2) + (1137-K3)*(1137-K3) = 217502634
...

老师好...


是不是 我的限定错误了?

一个小数 减  一个 大数 ,负数的平方 是正.
我感觉
是另外一种考虑...
不知道我的理解是否正确

香川群子

caipiaofans 发表于 2014-12-3 09:00
应该 是 可以解的..

我去……上了帅哥的当了……他的算法是错误的。

香川群子

用下面这个程序代码可以得到很多、很多解：

'要求例子:
(1579-k1)^2 + (5803-k2)^2 + (1137-k3)^2 = 217502634
或：
(x1)^2 + (x2)^2 + (x3)^2 = 217502634

1. 
   
2. Sub test()
   
3.     Dim ar&(1 To 65530, 1 To 3), a1&, a2&, a3&, s&, x1&, x2&, x3&, k&, tms#
   
4.     tms = Timer
   
5.    
   
6.     a1 = Range("a2"): a2 = Range("b2"): a3 = Range("c2"): s = Range("d2") '从单元格读取参数
   
7.     a1 = 9342: a2 = 9462: a3 = 7007: s = 217502634 '或直接赋值
   
8.    
   
9.     For x1 = 0 To Int(Sqr(s))
   
10.         For x2 = 0 To Int(Sqr(s - x1 ^ 2))
    
11.             x3 = Int(Sqr(s - x1 ^ 2 - x2 ^ 2))
    
12.             If x1 ^ 2 + x2 ^ 2 + x3 ^ 2 = s Then
    
13.                 k = k + 1: ar(k, 1) = x1: ar(k, 2) = x2: ar(k, 3) = x3
    
14.                 If k > 0 Then GoTo Ext '如果只需一组解则马上停止
    
15.             End If
    
16.         Next
    
17.     Next
    
18. Ext:
    
19.     MsgBox Format(Timer - tms, "0.000s ") & k
    
20.     If k Then Range("a5").CurrentRegion.Offset(1) = "": Range("a6").Resize(k, 3) = ar
    
21. End Sub
    

复制代码

当参数要求为：a1 = 9342: a2 = 9462: a3 = 7007: s = 217502634 时，
得到第一组解为：
x1 = 4: x2 = 10057: x3 = 10787

即：
k1 = 9342±4: k2 = 9462±10057: k3 = 7007±10787 这样就可以算一下子得到2*2*2=8组解了。

caipiaofans

香川群子 发表于 2014-12-3 12:09
用下面这个程序代码可以得到很多、很多解：

'要求例子:

嗯..谢谢老师..

我自己再 写个 计算公式 实现后面的 + - 法运算 也是可以的

caipiaofans

我测试 全部+的，得不到正确结果

4
10057
10787
10787

caipiaofans

测试有效...

刚才 计算有误...