之前在学习计数问题的时候也在网上找了很多关于行列式的资料 但是发现很多地方都只介绍2\3阶的情况 一些论文介绍的方法又看不懂 然后就一直耽搁着 今天恰好出到这样的题目 发现标算的代码简介明了 还挺开心的function det(a:arr):int64;var i,j,k:longint...
分类：
其他好文 时间：
2015-04-20 20:29:39
阅读次数：
142
题目大意：给定一张有向无环图，有恰好k个无入度的点和k个无出度的点，对于一个边集如果这个边集恰好形成了从每个无入度的点到每个无出度的点的k条不相交的路径，那么这k对点就会对答案有一个贡献；如果对应关系如果是一个奇排列，对答案的贡献为-1，否则为+1。求所有贡献的和
首先不考虑路径是否相交
令f[i][j]为从第i个无入度的点走到第j个无出度的点的方案数，那么这个矩阵的行列式的值就是答案
那么...
分类：
其他好文 时间：
2015-04-14 19:42:43
阅读次数：
413
最近几天学习高等代数老师说要写个程序算行列式的结果，闲来无事就简单写了一下。不多说了，上代码 1 using System; 2 using System.Collections.Generic; 3 using System.Linq; 4 using System.Text; 5 6 name....
线性代数课程，无论你从行列式入手还是直接从矩阵入手，从一开始就充斥着莫名其妙。比如说，在全国一般工科院系教学中应用最广泛的同济线性代数教材（现在到了第四版），一上来就介绍逆序数这个“前无古人，后无来者”的古怪概念，然后用逆序数给出行列式的一个极不直观的定义，接着是一些简直犯傻的行列式性质和习题——把...
分类：
其他好文 时间：
2015-04-12 20:43:28
阅读次数：
126
矩阵的基本概念
假设 aij∈Ra_{ij} \in \mathfrak{R}, 其中 i=1,2,...,mi=1,2,...,m; j=1,2,...,nj=1,2,...,n. 我们定义如下的行列式：
A=??????a11a21?am1a12a22?am2???a1na2n?amn??????\begin{equation*}
\mathbf{A}=
\begin{bmatrix}...
分类：
其他好文 时间：
2015-04-12 12:09:49
阅读次数：
172
点击打开链接 习题 行列式
高斯消元问题类型:
用LCM 保持整型
1. 基本的高斯消元,裸模板 HDU3359
2. 开关问题,用^操作代替 -, 求x[i]时候一样用* poj 1222 1830 1753
3. 枚举自由变元, return -1 是因为出现[0,0,0,0,a]这种情况,return 0 是唯一解,否则是有自由变元
4. 取模方程 (a1*x1+...
分类：
其他好文 时间：
2015-03-03 15:15:51
阅读次数：
220
步骤如下：
1，做beginPoint和endPoint的差vectorDiff。
2，对vectorDiff和该三角面片的法向量normal做叉积cross。
3，beginPoint和cross相加得到vectorAdd。
4，做三个向量(1,1,z),(1,0,z),(0,1,z).这三个向量的第三个元素都是z，要保证线性无关（三个线性无关的向量确定一个平面）。
5，对这三个向量以及vectorAdd做四点行列式，如果行列式的值小于0，反转，否则不反转。...
分类：
其他好文 时间：
2015-03-02 22:29:25
阅读次数：
313
var a:array[0..208,0..208] of int64; mo:int64; i,j,n,temp:longint;function guess:int64;var i,j,k,t:longint; ans:int64;begin ans:=1; for...
分类：
其他好文 时间：
2015-03-01 11:44:05
阅读次数：
105
这道题第一眼是生成树计数，n是100，是可以用O(n^3)的求基尔霍夫矩阵的n-1阶的子矩阵的行列式求解的，但是题目中并没有说取模之类的话，就不好办了。用高精度？有分数出现。用辗转相除的思想，让它不出现分数。但过程中会出现负数，高精度处理负数太麻烦。用Python打表？好吧，Python还不熟，写不...
分类：
其他好文 时间：
2015-02-15 23:00:56
阅读次数：
224
2.6 .1 符号矩阵分析
常用的矩阵分析指令如下：
det(A)
行列式A
diag(A)
取对角元素构成向量，或根据向量构成对角矩阵
[V,D]=eig(A)
特征值分解，使 AV=VD （没看懂）
expm(A)
矩阵指数e^A
inv(A)
A的逆矩阵
[V,J]=jordan(A)
Jordan分解，使AV=VJ
poly(A)
矩阵的特征多项式
rank...
分类：
其他好文 时间：
2015-02-08 16:56:48
阅读次数：
201
