V2EX = way to explore
V2EX 是一个关于分享和探索的地方
现在注册
已注册用户请  登录
推荐学习书目
Learn Python the Hard Way
Python Sites
PyPI - Python Package Index
http://diveintopython.org/toc/index.html
Pocoo
值得关注的项目
PyPy
Celery
Jinja2
Read the Docs
gevent
pyenv
virtualenv
Stackless Python
Beautiful Soup
结巴中文分词
Green Unicorn
Sentry
Shovel
Pyflakes
pytest
Python 编程
pep8 Checker
Styles
PEP 8
Google Python Style Guide
Code Style from The Hitchhiker's Guide
mingyun
V2EX  ›  Python

如何用 Python 解北大招生题

  •  
  •   mingyun · 2019-01-19 19:24:19 +08:00 · 1997 次点击
    这是一个创建于 2140 天前的主题,其中的信息可能已经有所发展或是发生改变。

    在 b 站看了个视频北大招生题

    北大招生题:a+b+c=1,a²+b²+c²=2,a³+b³+c³=3 求 abc 的值

    img

    结果 abc=1/6

    除了用数学公式怎么用 python 来解呢?搜了下有个 sympy 的库

    from sympy import *
    x=Symbol('x')
    y=Symbol('y')
    z=Symbol('z')
    solve([x+y+z-1,x**2+y**2+z**2-2,x**3+y**3+z**3-3],[x,y,z])#这个不行
    
    solve([x+y+z-1,x**2+y**2+z**2-2,x**3+y**3+z**3-3],[x*y*z])
    []
    
    5 条回复    2019-02-07 01:30:31 +08:00
    jiejiss
        1
    jiejiss  
       2019-01-19 20:58:50 +08:00 via Android
    可能需要先学线性代数,之后用线代那一套去解。
    qwertyegg
        2
    qwertyegg  
       2019-01-21 05:57:11 +08:00
    prolog
    necomancer
        3
    necomancer  
       2019-01-24 22:51:39 +08:00
    黑科技飘过

    Mathematica 11.2.0 Kernel for Linux x86 (64-bit)
    Copyright 1988-2017 Wolfram Research, Inc.

    In[6]:= Reduce[Exists[{x,y,z},x+y+z==1&&x^2+y^2+z^2==2&&x^3+y^3+z^3==3&&x*y*z==a],a] // ToString

    Out[6]= 1
    a == -
    6

    ……
    necomancer
        4
    necomancer  
       2019-01-24 22:52:30 +08:00
    空格被和谐掉了……呃(⊙﹏⊙)
    lovestudykid
        5
    lovestudykid  
       2019-02-07 01:30:31 +08:00
    DeleteDuplicates[ Table[a*b*c /.Solve[{a + b + c == 1, a^2 + b^2 + c^2 == 2,a^3 + b^3 + c^3 == 3}, {a, b, c}][[i]] // FullSimplify, {i,1,6}]]
    Sympy 大概可以用类似这种笨办法,mathematica 是黑科技
    关于   ·   帮助文档   ·   博客   ·   API   ·   FAQ   ·   实用小工具   ·   1310 人在线   最高记录 6679   ·     Select Language
    创意工作者们的社区
    World is powered by solitude
    VERSION: 3.9.8.5 · 19ms · UTC 23:27 · PVG 07:27 · LAX 15:27 · JFK 18:27
    Developed with CodeLauncher
    ♥ Do have faith in what you're doing.