Mysterium-CSDN博客

原创 prolog : Grammar and Parsing

今天意外发现，prolog原来有parsing macro，这样用它来做文法分析就更方便了，写起来很直观，就跟过去学的计算理论里面的写法一模一样。譬如有文法 S->aST. S->e. T->aSb. T->c.借助parsing macro 只需要写如下代码：s

2013-02-09 14:50:57 489

原创 Near Duplicate Image Detection

O. Chum, J. Philbin, A. Zisserman, Near Duplicate Image Detection: min-Hash and tf-idf Weighting, British Machine Vision Conference, 2008.这篇论文主要是提供了一种修改后的minhash方法去做图像的邻近重复检测。我对论文中提出的weighted hi

2013-01-28 12:22:14 2114

原创 Prolog 模拟 Turing Machine

既然写了PDA Parser，干脆再写个TM Parser。TM Parser 最麻烦的地方在于对整个TM的定义，就是An Example下面的部分。因此我只选了一个正则的简单语言来作为例子。%%2013-01-28%% Turing Machine Parser%% Author : Mysterium%%% Turing machine parser%%% sp

2013-01-28 10:58:41 555

原创 Prolog 最长公共子序列 LCS 问题

Prolog 解决这一类递归问题实在是极为简洁。%% LCSmax(A,B,M) :- A > B,M is A.max(A,B,M) :- A =lcs([],[],0).lcs([],_,0).lcs(_,[],0).lcs(X,X,1).lcs([A|R1],[A|R2],N) :-lcs(R1,R2,N1),N is N1+1

2013-01-27 20:41:02 521

原创 Prolog 模拟 PDA

在http://www.csupomona.edu/~jrfisher/www/prolog_tutorial/contents.html上有一个prolog的DFA parser，于是通读理解以后，模仿着写了一个PDA parser。PDA(pushdown automata)是接受上下文无关语言的装置，这里我使用了一个例子，其生成的语言为{a,b}上的回文。%2013-01-27

2013-01-27 18:41:03 729

原创 Lenstra因子分解椭圆曲线方法

这几天一直在研究ECM因子分解算法。这个方法主要基于椭圆曲线上的加法运算，一种群运算。椭圆曲线是代数几何里面非常重要的曲线，虽然我对代数几何了解甚少，但是椭圆曲线的大名早已如雷贯耳。从纯数学的费马大定理的证明到密码学的椭圆曲线加密体制……椭圆曲线的作用远不止此。下图是一个椭圆曲线的图像。可以证明椭圆曲线是一个加法群，我们可以定义其上的加法：

2012-08-30 17:45:17 4500

翻译 PRIMES is in P

PRIMES is in P的主要内容翻译，这是以前作为作业的时候翻译的，略去了一些次要内容。放在这里普及一下AKS。

2012-08-28 14:55:47 6976

原创整数因子分解的Pollard-rho方法

Pollard-rho方法该方法由Pollard于1975年在论文A monte carlo method for factorization中提出，论文可以见http://www.cs.cmu.edu/~avrim/451f11/lectures/lect1122_Pollard.pdf这个算法来源于一个事实：

2012-08-28 14:34:57 1669

原创整数因子分解的Pollard p-1方法

整数因子分解最直观的方法当数“试除法”，数论中的Mertens定理告诉我们76%的奇数都有小于100的素因子，因此对于大多数整数，“试除法”已经足够，但是对于特殊的数，特别是素因子普遍较大的时候，“试除法"的效率便明显不足。和素数检验类似，目前几乎所有实用的分解方法都是概率性的算法, 目标是找到能计算x 的算法, 使得(x,N) > 1 的概率较大(而最大公因子可以很快地计算).

2012-08-27 13:29:39 10382 2

原创 Solovay-Strassen素数检验

现在有效的素数检验算法基本都基于概率性的检测，虽然2002诞生了AKS算法，但是在实际应用中，依然依赖概率性检测算法。Solovay, Strassen于1974年最早提出了概率性的素数检验方法，我们现在称为Solovay-Strassen检测代码：（其中expmod是幂取模（见rabin-miller的代码），jacobi是jacobi符号计算函数）(defun so

2012-08-26 13:15:52 5420

原创 Jacobi符号计算的一种实现

Jacobi符号是Legendre符号的推广，在素性检验以及因子分解等诸多方面极为有用，所以先写一个备着，写其他算法时方便直接调用。(defun jacobi (n m) (let ((e 0) (a1 0)(a n) (s 0)(b m) (n1 0)) (cond ((= a 0) 0) ((= a 1) 1)

2012-08-26 12:15:19 2326

原创 Rabin-Miller概率素数检验算法

一个应用非常普遍的素数检测算法：我的实现：(defun expmod (base exp modulo) (if (= exp 0) 1 (if (evenp exp) (mod (expt (expmod base (/ exp 2) modulo) 2) modulo) (mod (* base (expmod base (-

2012-08-25 17:26:32 763

原创 (二)垂直几何、度量几何与常用几何(1)

我们先回顾一下Desargues（平面）几何的定义：以关联公理、平行公理、无限公理、Desargues公理为基础的几何。而且存在惟一的Desargues数系N与Desargues平面一一对应，称为几何附属的Desargues数系。这一数系一般不满足乘法交换公理（即a*b不等于b*a)，为了使Desargues数系变成一个数域，我们必须引进其他的公理使得乘法交

2012-08-23 15:02:35 1026

原创练手的程序

文章太少了，于是贴点东西上来好了。刚开始学习lisp，写了一些极其简单的程序来练手。;;;欧拉函数的计算;;;递归触发函数(defun Euler (n) (e n 2));;;这个版本遇素数会溢出(defun euler-recursive (n b) (if (= n 1) 1 (if (= (mod n b) 0)(if (= (mod (/ n b) b)

2012-08-22 18:38:32 366

原创 (一)Desargues几何与Desargues数系(2)

这一章的东西主要还是关于几何基础的，因此我的笔记较为简略。Desaruges数系定义：假设一个集合N，其中元素之间有两种二元运算加法与乘法，满足以下三类公理，则称N为一个Desargues数系，N中的元素称为数。三类公理为：关联公理计算公理

2012-08-22 18:28:07 558

原创 (一)Desargues几何与Desargues数系(1)

近来一直在看一些计算理论的资料同时pick up了一些common lisp， SRTP搁浅很久了，自从上次在论文里遇到看不懂的方法，就一直没弄明白论文里面的计算方法，没有任何推进，不明白又没地方可以寻求帮助，只好把它搁在一边了，况且SRTP的内容不是很理论很数学，也就没了什么兴趣。昨天突然心血来潮，打算将在硬盘里沉睡许久的吴文俊《机器证明》那本书看完，于是决定每天看一点，并写下读书笔记，以便鞭策

2012-08-22 12:08:32 684

Mysterium的专栏

转载从对偶问题到KKT条件