高斯消元
文章平均质量分 74
BraketBN
这个作者很懒,什么都没留下…
展开
-
【BZOJ1013】[JSOI2008]球形空间产生器sphere【高斯消元】
将两个点的方程做差可以得到一个线性方程,就可以高斯消元了。 /* Footprints In The Blood Soaked Snow */ #include #include #include using namespace std; typedef double DB; const int maxn = 12; const DB eps = 1e-5; int n,原创 2016-03-17 15:41:58 · 408 阅读 · 0 评论 -
【SPOJ-BTTNS】Buttons【高斯消元】【异或方程组】
还是模板题... 一共n*m个方程。对于一个点,与它曼哈顿距离不小于k的点,方程系数为1,其余为0,。方程右边等于最后的状态。 直接消元就好了。 没找到好看的模板,就自己脑补了个,估计有点丑。 建方程组时不要枚举点,而是直接算出曼哈顿距离符合要求的点的坐标,直接填值。否则TLE。 注意n和m不要搞混了。 #include #include usin原创 2016-02-18 12:19:17 · 494 阅读 · 0 评论 -
【SPOJ-MIFF】Matrix inverse【高斯消元】
算是模板了吧。 经典方法: 左边写原矩阵,右边写单位矩阵,然后把左边高斯消元成单位矩阵,右边就是逆元。 #include #include #include using namespace std; const int maxn = 205; int n, m, p, a[maxn][maxn]; inline int qpow(int x, int n) {原创 2016-02-17 13:09:28 · 552 阅读 · 0 评论 -
【SPOJ-HIGH】Highways【高斯消元】【Matrix Tree定理】【行列式】
求生成树个数... n 注意开LL。 (又写了一发高斯消元) #include #include #include using namespace std; typedef long long LL; typedef long double LD; const int maxn = 15; const double eps = 1e-3;原创 2016-02-18 17:42:53 · 965 阅读 · 0 评论 -
【BZOJ1778】[Usaco2010 Hol]Dotp 驱逐猪猡【概率DP】【高斯消元】
【题目链接】 竟然被卡eps。 设xu为到达u点没有爆炸的概率。 方程为xu - (1 - p / q)∑xv / du[v] = [u = 1]。 高斯消元一发就好了。 答案要乘p / q。 /* Footprints In The Blood Soaked Snow */ #include #include #include using namespace st原创 2016-03-28 10:56:08 · 1072 阅读 · 0 评论 -
【BZOJ3143】[Hnoi2013]游走【高斯消元】【期望DP】【贪心】
【题目链接】 比较经典的模型。 先高斯消元求出经过每个点的期望次数,每条边经过的期望次数为两个端点经过的期望次数除以各自的度数。 然后按照经过次数贪心编号。 注意走到n点之后不会再继续走了。 /* Telekinetic Forest Guard */ #include #include #include #include using namespace std; typ原创 2016-05-24 10:44:46 · 536 阅读 · 0 评论 -
【Ural1041】Nikifor【拟阵】【线性无关】【高斯消元】【矩阵的秩】
【题目链接】 这题有毒... 论文 刘雨辰《对拟阵的初步研究》 里提到了这个题,于是来做一下。 这题做法是贪心,把每个向量的花费从小到大排序,然后能放就放。 主要是判断放入一个向量之后,是不是线性无关的。所以我们去求矩阵的秩,看是否增加了1。 和一般的高斯消元不同,一般的是用当前行去消其他行,而这个是用其他行消当前行...为了方便,记录一个vis[j]数组,表示第j列(第j列非0原创 2016-04-29 12:01:44 · 914 阅读 · 0 评论 -
【BZOJ3659】Which Dreamed It【有向图欧拉回路计数】【matrix tree定理】【BEST定理】【高斯消元】
【题目链接】 定理题... /* Think Thank Thunk */ #include #include #include using namespace std; typedef long long LL; const int maxn = 105, p = 1000003; int n, fact[200005], A[maxn][maxn], du[maxn];原创 2016-07-18 15:19:12 · 2383 阅读 · 0 评论