文章/答案/技术大牛

发布

社区首页 >问答首页 >如何在渐近中复制此Wolfram解题示例

问如何在渐近中复制此Wolfram解题示例
EN

Stack Overflow用户

提问于 2020-02-09 23:39:44

回答 1查看 44关注 0票数 0

我想使用渐近来重现用Wolfram语言获得的结果。

使用Wolframcloud，这个表达式

Solve[m^2+m*n==500 && m>n,{m,n},PositiveIntegers]

给出了我想要的结果：

{{m->20,n->5}}

我如何使用sympy重现它呢？

我试过了

import sympy as sp
m,n = sp.symbols('m n',integer=True)
sp.solve(m**2 + m*n - 500, m,n)

这给了我们

[(m, -m + 500/m)]

这是正确的，但不是特别有帮助。

请注意，此问题的灵感来自Project Euler问题9。

python

sympy

wolfram-language

回答 1

Stack Overflow用户

回答已采纳

发布于 2020-02-10 00:15:34

对于整数解，您应该使用丢番图解：

In [10]: m, n = symbols('m, n')                                                                                                                

In [11]: sols = diophantine(m**2 + m*n - 500, (m, n))                                                                                          

In [12]: sols                                                                                                                                  
Out[12]: 
{(-500, 499), (-250, 248), (-125, 121), (-100, 95), (-50, 40), (-25, 5), (-20, -5), (-10, -40), (-5, -95), (-4, -121), (-2, -248), (-1, -499)
, (1, 499), (2, 248), (4, 121), (5, 95), (10, 40), (20, 5), (25, -5), (50, -40), (100, -95), (125, -121), (250, -248), (500, -499)}

这给出了整数m，n的解决方案。您可以筛选出满足条件的解决方案：

In [13]: {(m, n) for m, n in sols if m > n and m > 0 and n > 0}                                                                                
Out[13]: {(20, 5)}

票数 3

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/60138293

复制

相似问题

问如何在渐近中复制此Wolfram解题示例
EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问如何在渐近中复制此Wolfram解题示例EN

回答 1

Stack Overflow用户

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问如何在渐近中复制此Wolfram解题示例
EN