一刀不二

原创 Meta-Learning with Latent Embedding Optimization

LEO 思想MAML 虽然可以通过少量梯度下降就能找到适应新任务的最优参数，然而对extreme low-data regimes在高维参数空间上操作时，还是过于不便。而LEO 则通过学习关于模型参数的a data-dependent latent generative representation，并在这种低维的latent space 中执行基于梯度下降，可以一定程度上绕过这些痛点。LEO 算法能够捕获数据中的不确定性，将基于梯度的自适应优化过程，从模型参数的高维空间中分离出来，并通过实验...

原创 过去的三年

白驹过隙，三年过去了，又学了一些东西，可是很久没写博客，因为懒。今天突然想吐个槽。然，语文不好，思维跳跃，想到什么写什么。代数系统：《抽象代数》-->《泛代数》-->《范畴论》。这些学科最好都学一下，它们不仅在“密码学”，“形式化”中有很多应用，在神经网络中也有不少应用，比如将 Backpropagation 抽象为 functor；比如研究 Symmetry Net...

原创 多维空间树 kd-tree

#lang scheme( define ( make-node point axis left-child right-child )   ( define ( dispatch msg . args )      ( cond          [ ( eq? msg 'point ) point ]         [ ( eq? msg 'axis ) axis

原创 Emacs 降魔系列之制作 Python IDE

个人的武装配置比较搓。希望大神指正下。我们先去这里收集主人公闯关道具 —— 昊天塔 emacs for windows：http://ftp.gnu.org/gnu/emacs/windows/我使用的是冥界版 24.3。释放 emacs 的灵压，将其解压 emacs 在你需要的目录下，在解压后生成的 bin 目录下，点击 addpm.exe，会在 “C:\Users\Ad

原创 百世经纶一页书

大学 4 年的计算机专业的生活快结束了，看完了大部分的书，还有一小部分没看完，记之。看了比较杂，有些也不是很好分类。>>>>>>>>>>>>>>>>>>>>>>>道格拉斯 《图论导引》王树禾《图论及其应用》《有向图理论，算法及其应用》《欧拉图理论》《图论与代数结构》《复杂网络理论及其应用》《网络群体与市场》>>>>>>>>>>>>>>>>

原创 高度平衡树 -- AVL 树

原本打算用 Cpp，Ada，Python 各重写一份，觉得这种行为比较2.因为还是 Scheme 的表达最优雅.#lang scheme( define nil '() )( define ( root tree )( car tree ) )( define ( left-tree tree )( cadr tree ) )( define ( right

原创 二分图最大匹配 -- 匈牙利算法

Algorithm.( Augmenting Path Algorithm )Input:    An X-Y bigraph G, a matching M in G,    and the set U of M-unsaturated vertices in X.        Idea:    Explore M-alternating paths form

原创 DFA最小化 -- Hopcroft算法 Python实现

问了 30 个技术群，也问了无数的前辈，真是各种不礼貌，吃了无数闭门羹，还是自己看着有点眉目了还有 wiki 的伪代码看了总觉得奇怪，于是看了同一页面其他语言翻译过来的伪代码，发现葡萄牙语和俄罗斯语那里的 if 判断都还缺少一个条件国内的资料比较少，这几份学习资料不错，比我稀里糊涂的思路要好，分享下：http://www.liafa.univ-paris-diderot.f

原创 矩阵树定理速证

凯莱公式:spanning_trees_num( G ) = spanning_trees_num( G - e ) + spanning_trees_num( G · e )矩阵树定理：G 对应的拉普拉斯矩阵（度矩阵 - 邻接矩阵）L( G )  删除任意一行一列得到的行列式的值det( L*( G ) )即生成树的个数，即spanning_trees_num( G

原创 并查集小记

实现不相交集合数据结构的几种方式：单纯的快速查找：若id相同则在一个集合中，下图中，（ 2， 3， 4， 9 ）为一集合， 3 和 6 不在一个集合中合并集合时，需逐个比较将两个集合的 id 统一，慢缺陷：合并慢==================================================2.单纯的快速合并：

原创 有向图的强连通算法 -- tarjan算法

在有向图 G 中，若两个顶点相互可达，则称两个顶点强连通(strongly connected)。如果有向图G的每两个顶点都强连通，称G是一个强连通图。非强连通图有向图的极大强连通子图，称为强连通分量(strongly connected components)。比如下面第一幅图( a, b, e ), ( d, c, h ), ( f, g ) 分别为三个 SCC。

原创 python 装饰器与AOP

无高见1.缓存from functools import wrapslineseq = '==' * 20def memo( func ): cache = {} @wraps( func ) def wrapper( *args ): result = cache.get( args ) if resul

原创 python中lambda的另类使用

Lambda挺强大，有兴趣的人看下关于lambda的理论，就清楚邱奇编码的实现了。带if/else：( lambda x, y: x if x 科里化：( lambda x: ( lambda y: ( lambda z: x + y + z )( 1 ) )( 2 ) )( 3 )递归：func = lambda n: 1 if n == 0 e

原创 图的一因子

图特定理：对于任意点集S ∈V( G )，满足o( G - S ) ≤ | S |，G则存在一因子。（ o代表分量中点数为奇数的分量个数 ）证明：必要性：很显然，若是G存在一因子，G - S中的奇分量的个数必然要小于S的个数，不然G - S中的某个奇分量中的点得不到匹配。充分性：反证，假设G满足任意S，o( G - S ) ≤ | S |，但是G不存在一因子。但

原创 蚁群算法求解旅行商问题 python

这组数据，刚开始只运行到2879m，于450m差的比较远，后来发现是路径距离算错了，再后来变成了1789m，发现是初始化点没选择，改了下，最后运行结果是440m。python有点慢，需要运行一段时间。import randomimport copyimport timeimport sys( ALPHA, BETA, RHO, Q ) = ( 1.0, 2.0, 0.

原创 压入与重标记算法（预流推进算法）

最大流算法之压入与重标记算法（预流推进算法）今天是平安夜，明天是圣诞节，明天也是我的生日，到目前都是直接贴代码，懒动手，今天好好写一篇博客。                                             — —Pandora_Scheme1.算法发现者:      Goldberg && Tarjan( 87 );            不了解

原创 Van Kampen 定理

我们知道任何一个闭曲面去掉一点后可以强形变收缩为曲面上的一个圆束。因此求圆束的基本群大有用途。

原创 为什么定义链群时需要用定向单形而不是非定向的

链群的一个重要应用是定义边界算子，边界算子将一个单形映射到其边界的链。为了确保边界算子的定义是自洽的，则需要单形的定向。例如，一个定向的三角形 σ=[a,b,c] 的边界是其三条边的线性组合，每条边的系数取决于 σ 的定向。如果没有定向，边界算子的定义将变得模糊不清。这时候是说 [a,b] ∈ boundary(σ) 还是 [b,a] ∈ boundary(σ)?[b,a]+[a,b] 是等于 2[a,b] 还是等于 0？

原创 同调群的计算实例

仅当 ( i = j = k ) 时，此条件才满足。中没有 3-单形，我们有。中任何单形的边界，因此。根据同构定理，我们得到。

原创 持续同调 —— 同调类的消亡

晚上想明白了，从算法角度，一个单纯形的添加，若导致其 boundary class 消亡，则它必然是变成 trivial 的，但是这过程中必然也伴随着其他同调类的 merge，这些同调类的合并必然是新类和旧类的合并。第二定义更通用，它其实包含了第一种（因为 trivial class 可以视为 older class），但是从算法角度，可以只用第一种定义，因为注意力只需要要放在 boundary class 上。

原创 Blender 怎么导出透明背景图片

原创 blender 制作莫比乌斯带

创建 Curve -> Cycle在 Edit 模式下，选择：选中两个点，按 delete 删除 Segment如下选中：选中最上面的点，然后按 E 将它拖到右边的点上。按 R 旋转 90 度。依次调整参数：回到 Object 模式下，调整参数：

转载 APMO 1998 不等式

原创 sympy 矩阵求导零碎

1. Frobenius 范数对矩阵求导。

原创 christoffel 符号的直观解释

看的 bilibili 的“冰侍狐”的讲解，就不整理了，直接截图。

原创 python 目录巴拉巴拉

《End-to-End Differentiable Proving》

原创 blender 导入动画到 unity

看自由小径教程，自己记的乱七八糟的笔记。1.骨架设置：物体模式下，添加头部骨架。改名head。（object mode 下的骨架透视：选中骨架，在 viewport display 中选中 in front。）编辑模式下，选中头骨，按 E 拉出 3 条腿骨，删除与头部相接的腿骨。复制旋转出另外两条腿。改名 leg.up.L，leg.down.L 等等。（姿态模式下，观察骨架运动关系。）编辑模式下，依次选中 3 条 up 腿，最后选中 head，Ctrl + P 保持偏移量（出现虚线.

原创 有点迷糊的 B 树 delete 操作

1.如果键 在节点 中，且 是叶子节点，则从 中删除键 。2. 如果键 在节点 中，且 是内部节点，请执行以下操作：a. 如果节点 中位于 之前的子节点 至少有 个键，则在以为根的子树中找到 的前驱值。将中的 替换为​​​​​​,并递归地删除 。b. 对称地，如果节点 中位于 后面的子节点 至少有 个键，则在以 为根的子树中找到 的后继值 。将 中的 替换为 ，并递归地删除 。c. 否则，如果 和 都只有 键，那么...

原创 Standard ML 无穷数据

datatype 'a seq = Nil| Cons of 'a * ( unit -> 'a seq );exception Empty;fun head( Cons( x, xf ) ) = x  | head( Nil ) = raise Empty;fun tail( Cons( x, xf ) ) = xf()  | tail( Nil ) = rais

转载 归结原理

归结原理程序设计：http://www.docin.com/p-1293107814.html

转载 sat-solver

python 实现的 sat-solver，用了三种算法，DPLL，Genetic，HillClimbing ：https://github.com/DRTooley/PythonSatSolver。

原创 计算群论 Schreier-Sim algorithm

计算群论 Schreier-Sim algorithm：http://www.m8j.net/data/List/Files-118/Documentation.pdf

原创 Kunth 的 LR parsing 论文

http://www.dc.uba.ar/materias/tl/2010/c2/files/Knuth-1965-TranslationofLanguages.pdf

原创 Beautifulsoup 小用

用 beautifulsoup 爬了下伯克利大学 programming languages  and compilers 的课件import reimport requestsfrom bs4 import BeautifulSoupr = requests.get( "http://inst.eecs.berkeley.edu/~cs164/fa11/lectures/in

原创 win7 下 bochs 运行 Hello World OS

安装 bochs，dd for windows，nasm用 bochs 中 bximage.exe 创建新的 img 文件编写 Hello.asm org 07c00h ; 告诉编译器程序加载到 7c00处 mov ax, cs mov ds, ax mov es, ax call DispStr ; 调用显示字符

原创 CrackMe 之Afkayas 破解

rtcAnsiValueBstr 是取输入字符串的第一个字母的 ascii 码注册码是 len( 输入 ) * 0x17CFB + ascii( 输入的第一个字符 )

原创 CrackMe 之 Acid Burn 破解

注册码是蓝线之上就生成好的，但是，是在红色断点处 call 之后才载入到数据区域的。注册码算法是取输入的第一个字节乘以 0x29 再翻倍后，头部接 CW-，尾部接 -CRACKED。

原创 汇编语言 冒泡排序

.section .datamsg:.string "value is %d\n"arr:.int 1, 10, 5, 2, 20, 55, 7, 124, 77, 55.section .text.globl mainmain: pushl %ebpmovl %esp, %ebpsubl $32, %espmovl $10, -4(%ebp

原创 python rmq

import mathdef init( arr ): arr_len = len( arr ) exp = int( math.log( arr_len, 2 ) ) DP = [ [ 0 for col in range( exp + 1 ) ] for row in range( arr_len + 1 ) ] for index, item in

原创 用 ply 制作计算器

import ply.leximport ply.yacctokens = ( "NAME", "NUMBER" )literals = [ "=", "+", "-", "*", "/", "(", ")" ]t_NAME = r"[a-zA-Z_][a-zA-Z0-9_]*"def t_NUMBER( t ): r"\d+" t.value = int( t.