0%

因为写工程/写作业时频繁要(偷懒地)使用计算机辅助计算,之前一直用的是Mathematica,但是语法和界面一直不尽如人意。遂开始研究IPython+Numpy的组合,同时发现了超强的符号运算库Sympy.

温馨提醒:Ctrl+L可清空屏幕。

Numpy/Scipy

1
2
np.linalg
scipy.linalg.polar(np.array) # 极分解

Sympy

一些前置概念:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
from sympy import *
x, y, z, t = symbols('x y z t')
k, m, n = symbols('k m n', integer=True)
f, g, h = symbols('f g h', cls=Function)
init_printing()

Ms = MatrixSymbol('Ms', 3, 3) # 创建3×3的矩阵Symbol, not necessary

f.subs(x, VALUE) # 代入
f.subs({x: }) # 代入多个变量
f.subs({Ms: }).doit() # 矩阵运算是惰性的,需要手动
POLY(VALUE) # 计算多项式在VALUE处的取值
EXPR.evalf(n) # 输出n位小数

f.expand()
expand(f) # 展开
f.simplify()
simplify(f) # 化简
f.factor()
factor(f) # 因式分解
f.as_expr() # 将f作为表达式

solve(x-y+1, (x,y)) # 解方程组 x-y+1=0
limit(x**x,x,0) # 求极限
limit(1/x,x,oo,dir='+') # 右极限
Derivative(f,x,x,y,z)
f.diff(f,(x,2),y,z) # 计算导数
integrate() # 求积分
dsolve() # 解微分方程

gcdex() # 扩展欧几里得

A.exp()

A.jordan_form()
A.is_positive_semidefinite()
A.QRdecomposition()
A.LUdecomposition()
A.LUdecomposition_Simple() # 不满秩情况
A.LDLdecomposition(hermitian=True) # Return L,D; A = LDL^H
A.cholesky(hermitian=True)
A.bidiagonal_decomposition() # Return U,B,V^H; A = UBV^H
A.singular_value_decomposition()

A.eigenvals()
A.eigenvects()
A.left_eigenvects()
A.singular_values()

IPython

阅读全文 »

记录THU生活中的小天才瞬间。

  • 大风天的小天才停车法

  • 小天才的Windows推迟更新

  • 给你自己的颜值打个分,10分制我知道,2进制是吧,我打1.

  • 小天才赶DDL,灵感来自高代群:

    赶不上DDL了,怎么办?网络学堂先随便交个文件,然后写完了跟老师说,不好意思之前交错了,刚发现,重交一下。

  • 实验室生存指南:

    1
    2
    3
    4
    5
    solid
    promising
    reasonable
    it works
    make sense
  • 下周三再考虑吧。

  • 迷惑行为:

    IMG_2022-09-10T22-10-50.jfif

由于ZRQ太强啦,因此他的名场面可以单独列出来:

  • 对某个车子尝试开锁,结果直到车主过来,他才意识到这不是自己的车子

  • 把导师的名字记错了

  • 因此痛失4.0:

    image-20220416124633451

  • 投屏的王:

    1657893824888

更多张总语录请见这里

习题一

2.

假设G存在孤立节点,不妨设其编号为n,则1…n-1号点的度数最大为n-2,则m<=(n-1)(n-2)/2,这与m>(n-1)(n-2)/2矛盾,则原命题成立。

6.

证明9个人中如果不存在4人互相认识,则存在3人互相不认识。


反证法。假设不存在4人互相认识,也不存在3人互相不认识(*)。从9人中任选一人,称为A.

Lemma 1:A至少认识5人。

阅读全文 »

苦于写作业时经常分心/摸鱼,一不留神就如白驹过隙,光阴飞逝。从GZZ学长得到启发,遂决定尝试番茄工作法(Pomodoro)。

配置:

  • 浏览器:Strict Workflow

    实践番茄工作法:工作时屏蔽浪费时间的网站,休息时允许访问。

  • Andriod:iPomodoro

鉴于以上应用没有统计数据,遂开此贴以记录。自己为Strict Workflow插件添加了一些新功能,可以看这里.

2021

3.7: ❤️ × 15

3.8: ❤️ × 6

3.9: ❤️ × 7

阅读全文 »

Lovely i207M forever! Lovely OIer forever!

以下个人总结按照时间倒序排列,一字一句都是我OI生涯一点一滴积累起来的。真心希望能够对你有所帮助。

可能远古时期比较傻,请见谅

小技巧Trick-思路-个人总结

杂项个人总结

本来想把链接贴出来的,结果懒了,

请在“个人总结”分类下搜索以下关键字:

组合数学、数学、DP(动态规划)、贪心、数据结构、图论、字符串

阅读全文 »

李宏毅老师好评!

Pre

优化方法:

Stochastic Gradient Descent 随机梯度下降 Update for each example

image-20201018190125562

Normalize归一化:Make different features have the same scaling.

image-20201019201617440

一次函数逼近

阅读全文 »

From 洛谷日报

正则表达式(Regular Expression, 简写为regex)

选择

|:或运算,选择(最低优先级)。

数量限定

+:前面的字符至少出现一次;

?:最多出现一次;

*:出现任意次(包括0次)。

阅读全文 »