python：序列1~n的立方之和==序列1~n的和之平方

pip install sympy 或者  Anaconda 3 自带 sympy 点击 【Jupyter Notebook】

这是我最喜欢的代数恒等式之一

from IPython.display import Latex
Latex(r"$1^3+2^3+3^3+\cdots+n^3 = (1+2+3+\cdots+n)^2 $")
Latex(r"$\sum_{i=1}^n i^3 = (\sum_{i=1}^n i)^2 $")

# 求序列1~n的立方之和，序列1~n的和之平方

from sympy import *
i = symbols("i")
n = 100123
s3 = summation(i**3, (i,1,n))
print(s3)
s  = summation(i, (i,1,n))
print(s)
print(s*s)

编写  sum_cube_eq_square.py  如下

# -*- coding: utf-8 -*-
import sys

# 获取用户输入 正整数
n = int(input("input n= "))

if n <= 0:
    print('n>0, n is int.')
    sys.exit(1)

sc = 0
for i in range(1, n+1):
    sc += i**3
print('sum cube:', sc)

s = sum(range(1, n+1))
print('he:', s)
print('square:', s*s)

运行 python sum_cube_eq_square.py 

input n= 100123
sum cube: 25123728970920355876
he: 5012357626
square: 25123728970920355876

数学上如何证明这个代数恒等式？还没想明白。