矩阵A的零空间就Ax=0的解的集合。
零空间的求法：对矩阵A进行消元求得主变量和自由变量；给自由变量赋值得到特解；对特解进行线性组合得到零空间。
假设矩阵如下：
对矩阵A进行高斯消元得到上三角矩阵U，继续化简得到最简矩阵R：
由于方程Ax=0的右侧是零向量，所以只对矩阵A进行消元不会影响解,因此不需要增广矩阵，所以有：
从上面的高斯消元的结果可以看出，矩阵...
分类：
其他好文 时间：
2014-10-13 11:27:46
阅读次数：
158
题目大意：给定一个半径为为r的圆x^2+y^2=r^2，求圆上多少个点的坐标为整数
卡了很久的一道题。。。我之前用了两个公式，理论上可以O(√n)出解，可惜这两个公式并不能涵盖所有勾股数。。。
于是去找了下题解，发现这样一种方法：
x^2+y^2=r^2
化简为 y^2=(r-x)(r+x)
我们令d=gcd(r-x,r+x)
则(r-x)/d与(r+x)/d一定互质，二者相乘为完全平...
分类：
其他好文 时间：
2014-10-09 02:15:37
阅读次数：
159
题目链接：HDU 2454 Degree Sequence of Graph G
题意：给出N个点的度（简单图），问能能否画出个图，（其实就是给出一个串非负的序列是否有对应的图存在）
没见过这个定理 题意真的难懂。
havel定理就是一个给出一串非负的序列，存在一个无向图使得图中各点的度与此序列一一对应，则称此序列可图化。简单图的话就是，可简单图化。
可简单图化的判定（Havel定理）：把...
分类：
其他好文 时间：
2014-10-06 12:15:30
阅读次数：
295
/*
dp[i][j]表示(i,j)到(R,C)需要消耗的能量
则：
dp[i][j]=p1[i][j]*dp[i][j]+p2[i][j]*dp[i][j+1]+p3[i][j]*dp[i+1][j]+2;///+2 转移到下一个能量要消耗2
化简得：
dp[i][j]=((p2[i][j]*dp[i][j+1])+(p3[i][j]*dp[i+1][j])+2)/(1-p1[i][j]);
*...
分类：
其他好文 时间：
2014-10-06 11:47:00
阅读次数：
213
简单的指数型生成函数题，就是按照泰勒级数化简而已。#include #include #define MOD 100#define LL __int64using namespace std;int quick(int a,LL b){ int ans=1; while(b){ if(b&1) .....
分类：
其他好文 时间：
2014-10-05 18:05:48
阅读次数：
186
原文:项目笔记---CSharp图片处理项目笔记---CSharp图片处理 最近由于项目上需要对图片进行二值化处理，就学习了相关的图片处理上的知识，从开始的二值化的意义到动态阀值检测二值化等等，并用C#得以应用，学到了很多的知识和大家分享下我个人的经验，希望对大家有帮助。 二值化
二值化简而言之是对...
分类：
其他好文 时间：
2014-09-29 01:38:07
阅读次数：
737
项目笔记---CSharp图片处理最近由于项目上需要对图片进行二值化处理，就学习了相关的图片处理上的知识，从开始的二值化的意义到动态阀值检测二值化等等，并用C#得以应用，学到了很多的知识和大家分享下我个人的经验，希望对大家有帮助。二值化二值化简而言之是对一副彩色图片进行0/1运算，最终显示一副黑白相...
分类：
其他好文 时间：
2014-09-28 22:25:35
阅读次数：
533
其实求的这个数的式子化简一下，就是C(N,M).....#include #include #include #define LL __int64using namespace std;LL N,M;LL myc(LL n,LL r){ LL sum=1; for(LL i=1;i<=r;i++) ...
分类：
其他好文 时间：
2014-09-27 12:59:39
阅读次数：
142
生成函数是说，构造这么一个多项式函数g(x)，使得x的n次方系数为f(n)。 如：序列{0,1，2，3，4，5...n}的生成函数为：$f(x)=0+x+2x^2+3x^3+4x^4+...+nx^n$生成函数最绝妙的是，某些生成函数可以化简为一个很简单的函数。也就是说，不一定每个生成函数都是用一长...
分类：
其他好文 时间：
2014-09-25 14:19:39
阅读次数：
202
题目大意给定一个函数找出满足条件 等于 k的最小的xm,k,d已知 其中 m，k 很大需要使用高精度存储思路：对 函数f(m)进行化简 ，令t=ceil( log(d,m) )可以得到 f(m)=d ^ t * ( a[ m / (d^t)] ) + d ^ (t-1) * ( b[ m/( d^(...
分类：
其他好文 时间：
2014-09-24 19:31:47
阅读次数：
354
