8rfuz的博客

原创 Dijkstra算法

Dijkstra算法 能够实现求得，图中某一源点s分别到其它所有点的最小距离上面的主要内容就是：1.Dijkstra算法的框架2.Dijkstra算法的运行时间效率分析：V记为顶点数，E即为边数a.     最朴素的情形，Dijkstra算法时间复杂度为O(V^2+E)b.     若算法中使用最小二叉堆实现，则时间复杂度为O(VlogV+E

原创 堆

原创 算法导论摘录

原创 floyd算法原理

floyd算法的本质是一个动态规划的过程。  下面  前k个点即：1，2,...,k状态转移方程：f[k][i][j]=min(f[k-1][i][j],f[k-1][i][k]+f[k-1][k][j])f[k][i][j]表示只经过前k个点(包括k)，从i到j的最小值。 我们熟悉的用二维数组的写法实际上是对空间的一种压缩。f[k][i][j]表示i到

原创 stl里自带的堆heap

最近要使用堆，懒得写了，所以想直接使用stl里面封装好的堆。结果发现功能好少，尤其少了一个最重要的adjust_heap函数来维护堆的性质。只找到了push_heap,pop_heap,make_heap，is_heap这些操作，虽然似乎够用，但没有adjust_heap函数来使用的话，真心很不方便。只好自己去写了...

转载 优化算法（模拟退火算法 与 遗传算法）

大白话解析模拟退火算法遗传算法

原创 Bezier曲线的几何连续性

计算机图形学链接：网络课程

原创 插值多项式

高阶等差数列问题其实跟 数值计算方法里的 插值问题差不多。还有一些其它的比如差商的许多公式：差商差分之间关系

原创 高阶等差数列的一些性质

原创 Bezier曲线绘制方法

de Casteljau算法能简单快速地求出某个t值的曲线值，复杂度是O(n^2)的乘法与加法。而如果是绘制一条曲线的话，如果要用m个等间距的点来逼近这条曲线的话，复杂度是O(m*n^2)，有乘法。感觉不应该用这个方法，太慢了。比如用高阶差分的方法，同样用m个等间距的点来逼近曲线的值，复杂度只要O(m*n)，而且只有加减法。(差了一个n，实际n取比较小，不过再小也是肉，n=3时

原创 1441 士兵的数字游戏

看懂题目后就知道，其实这道题就是求区间[1,a]内的所有数的质因子个数之和。方法利用线性筛素数的方法，可以以O(n)的复杂度求出所有素数。同时对程序稍作变动，就可以在求素数的同时求质因子个数，并累加#include//using namespace std;const int MaxPrimNum = 183072, MaxSize = 5000001; //, M

原创 poisson方程组矩阵

列出公式，然后离散化。最后对每个像素都满足公式，所以能得到稀疏矩阵A。矩阵A为弱对角占优矩阵。比如若其中对角线为4，则该行其它元素有4个-1，其它均为0。由于有边界条件，边界值已知，所以对于公式中的fq若是边界值，则可以移到右边去。所以最终矩阵A满足：某一些行为强对角占优，其它行为弱对角占优。则矩阵A很容易看出为对称矩阵。   并且，是非奇异矩阵。大致证明思路：

原创 Painting Fence题解

codeforces里题目链接：http://codeforces.com/problemset/problem/448/C要翻墙...这道题一看题目，应该一开始就想到DP，然后复杂度O(n2)，这个平方复杂度我就没去细想怎么实现了。然后就是想有没有更好的方法，复杂度更低。很容易能想到，如果某一列没有竖刷，那么必然被横刷，那么一定要从低处往高处都要横刷。

原创 1437 迈克步

51nod题目链接：http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1437有n只熊。他们站成一排队伍，从左到右依次1到n编号。第i只熊的高度是ai。一组熊指的队伍中连续的一个子段。组的大小就是熊的数目。而组的力量就是这一组熊中最小的高度。迈克想知道对于所有的组大小为x（1 ≤ x ≤

原创 一道有趣题目

求图中阴影面积。只要想到方法就很简单了。 很容易可以得到中间的阴影形状是个4面对称，旋转90对称的图形。并且每个弧边是30度角对应的弧。 所以阴影面积可以切割为4个弧面积+中间四边形面积，而且很容易得到中间四边形是个正方形...    进一步，可以把里面所有分割的图形面积都算出来。

原创 Color2Gray: Salience-Preserving Color Removal实现

上面的都是在实现论文过程中自己加的一些优化，有一部分优化暂时还没实现，还是很想实现的，不过暂时没时间，所以留待以后有时间再实现。

原创 Transferring Color to Greyscale Images

可以看出来，实现上有问题，主要是，真正的实现的话，需要对一堆二元值进行最近点查询。为了实现简单和效率，尝试加权后对一元值进行快排后在附近查询。不过似乎程序中有些bug。不过也没时间修改了，以后有时间再修改bug。

原创 Color Transfer between Images实现

原创 PoissonImageEdit实验报告

实验报告百度云链接： http://pan.baidu.com/s/1mgxnH9Y

原创 ImageWarping实 验 报 告 （IDW算法与RBF算法）

实验报告百度云 ： http://pan.baidu.com/s/1pJlCwsV

原创 1436 方程的解数

51nod网站http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1436这道题的方法很简单，关键在于这道题有很多角度都能切入思考，能不能快速想到合适的解法。方法步骤：将k化为二进制，则对应bit位置若为0，则a1,a2,...,an相应的bit位设为b1,b2,...,bn

转载 矩阵乘法的优化

/* for (i = 0; i < 2; ++i) for (j = 0; j < 2;++j) { result[i][j] = 0; for (k = 0; k < 2;++k) { result[i][j] = (result[i][j] + (A[i][k] * B[k][j])%m)%m; } }*/ result[0][0] = resu

原创 机器学习中的数学(5)-强大的矩阵奇异值分解(SVD)及其应用

机器学习中的数学(5)-强大的矩阵奇异值分解(SVD)及其应用http://www.cnblogs.com/LeftNotEasy/archive/2011/01/19/svd-and-applications.htmlhttp://www.cnblogs.com/LeftNotEasy/archive/2011/01/19/svd-and-applications.html

原创 摘录

学习shader之前必须知道的东西之计算机图形学http://game.ceeger.com/forum/read.php?tid=10970http://game.ceeger.com/forum/read.php?tid=10970

原创 QSlider控件使用记录

QLineEdit控件与QSlider控件： lineEdit = new QLineEdit("100"); //新建一个水平方向的滑动条QSlider控件 slider = new QSlider(Qt::Horizontal); //设置滑动条控件的最小值 slider->setMinimum(0); //设置滑动条控件的最大值 slider->