![](https://img-blog.csdnimg.cn/20201014180756919.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
高斯消元
Coco_T_
东南大学本科毕业生,浙江大学在读研究生
展开
-
高斯消元(写(shui)题必备)
前言: 在一次学校hu测中, 遇到一道正解不用高斯消元,但是部分分需要的中档题 用舒老师的话说,只要是会高斯消元和树形dp 乱搞一下那道题就可以水到70 所以还是学习一下这个很有用算法:高斯消元简介 数学上,高斯消元法(或译:高斯消去法),是线性代数中的一个算法, 可用来为线性方程组求解,求出矩阵的秩,以及求出可逆方阵的逆矩阵。 当用于一个矩阵时,高斯消元法会产生出一个“行梯阵式”。原创 2017-10-07 08:58:52 · 1274 阅读 · 0 评论 -
SPOJ - DETER3 Find The Determinant III(行列式求值+gauss)
题目链接题目描述:求行列式的值(模一正数)分析: 我们要用gauss消元将行列式消成上三角,计算对角线乘积 因为模数不一定是质数,所以gauss消元时不能用费马小定理求逆元了 我们就要用辗转相除法的gauss 每次消元(第jj行)时,默认a[j][i]>a[now][i]a[j][i]>a[now][i](被消元的主元小) 记x=a[j][i]/a[now][i]x=...原创 2018-03-29 08:03:59 · 199 阅读 · 0 评论 -
UVa 10766 Organising the Organisation(矩阵树定理)
题目链接题意: 给出n,m,k,n,m,k,代表一家公司有nn个部门,有mm组关系,表示ii和jj不能直接联通,kk代表主管部门,询问有多少种分层方案题解: 首先题面看不大懂,直接找到前辈的题目描述Ctrl+C一下 其次,这道题的k可以忽略掉,所以ta的范围完全是吓唬人的那么这道题就是求nn个点的有编号无根树数量 有编号无根树的一种有力的解决方法就是Prufer编码 但是...原创 2018-03-29 06:57:17 · 253 阅读 · 0 评论 -
Matrix-Tree定理(sx之前填坑还来得及吗)
从入门到入土:矩阵树Matrix-Tree定理 参考blog在正式介绍Matrix_Tree定理之前,我们需要一些前置知识一些定义与定理对于一个无向图GGG,ta的生成树个数等于其基尔霍夫Kirchhhoff矩阵任何一个N−1N−1N-1阶主子式的行列式的绝对值所谓N−1N−1N-1阶主子式就是对于任意的一个r,将矩阵的第rrr行和第rrr列同时删去得到的新矩阵基尔霍夫Ki...原创 2018-03-28 20:53:39 · 719 阅读 · 0 评论 -
loli的毒瘤hu测 T2&&T3
T2.tree(矩阵树定理+gauss) 分析: 这道题的爆搜真的恶心。。。 可以放出来吓吓人:const ll p=1e9+7;const int N=100005;struct node{ int x,y;};node t[N][31];int n,m,Q[N],deep[N],tot=0;ll ans=0;map<ll,int> mp;...原创 2018-03-28 15:23:26 · 278 阅读 · 0 评论 -
hdu5755 Gambler Bo(gauss解同余方程)
题目链接题目描述: 一个n∗mn∗mn*m由0,1,2组成的矩阵,每次操作可以选取一个方格,使得ta加上2之后对3取模,周围的四个方格加上1后对3取模,在n∗mn∗mn*m操作次数内让整个矩阵变成0。输出一种方案。分析: 模数好啊,原数=逆元状压dp枚举每一行肯定不行了 然而想到了开灯问题,这道题可以用高斯消元!设直接操作(i,j)(i,j)(i,j)的次数为xi,jxi,j...原创 2018-03-26 21:35:10 · 188 阅读 · 0 评论 -
bzoj3143 [Hnoi2013]游走(概率+gauss)
题目链接题外话: 小Z,你的袜子找到了吗?分析: 题目要求使得小Z获得的总分的期望值最小 那么考虑我们知道那些算法可以达到这个目的: 二分(不可能),线性规划(有门,但是不等式呢),dp(有点靠谱)但是我们的目的是确定边的编号,dp之类的应该是没有办法的 期望=概率*权值想要期望小,实际上就是 “概率*权值” 尽量小 显然,我们希望概率越小的边权值越大(mmp又...原创 2018-03-06 07:28:18 · 285 阅读 · 0 评论 -
bzoj3270 博物馆(概率+gauss)
题目链接分析: 感觉这道题和聪聪可可有点像实际上我们还是设置状态:f(x,y)f(x,y)f(x,y) 表示Petya在xxx,Vasya在yyy的概率(准确人名良心) 设每个结点的度为deg[i]deg[i]deg[i] 从i到达与ta相邻的任意节点的概率out[i]out[i]out[i]就是1−p[i]deg[i]1−p[i]deg[i]1-p[i] \over deg[i]...原创 2018-03-05 20:51:21 · 507 阅读 · 0 评论 -
bzoj2337 [HNOI2011]XOR和路径(概率期望dp+gauss)
分析: 之前写过图上的线性基,所以一下子就想到了线性基我在想是不是要把所有的路径都求出来*概率? 这显然不是很现实。。。可以尝试着把二进制的每一位分开考虑(以下的讨论都是针对二进制的一位)所以如果路径上经过了奇数个1,那么 “2i2i2^i * 路径上经过奇数个1的概率” 就是这一位的贡献设计状态:f[x]f[x]f[x]为从xxx点到nnn点,XOR和为1的概率 f...原创 2018-03-07 11:07:25 · 248 阅读 · 0 评论 -
bzoj2707 [SDOI2012]走迷宫(tarjan缩点+拓扑+概率期望+gauss)
DescriptionMorenan被困在了一个迷宫里。迷宫可以视为N个点M条边的有向图,其中Morenan处于起点S,迷宫的终点设为T。可惜的是,Morenan非常的脑小,他只会从一个点出发随机沿着一条从该点出发的有向边,到达另一个点。这样,Morenan走的步数可能很长,也可能是无限,更可能到不了终点。若到不了终点,则步数视为无穷大。但你必须想方设法求出Morenan所走步数的期望值。...原创 2018-03-07 09:03:45 · 680 阅读 · 0 评论 -
UVa11542 - Square(gauss)
题目链接分析: 题目特意说到,不含有超过500的素因子 这就是在提示我们可以质因数分解 我们把这n个数分解成若干素数乘积的形式,并用一个向量表示如:4 6 10 15 这n=4个数中,涉及到的素数只有2,3,5这三个 分解后: 4=2^2*3^0*5^0 ——>{2,0,0} 6=2^1*3^1*5^0 ——>{1,1,0} 10=2^1*3^0*5^1 ——>{1,0,1}原创 2017-10-12 20:45:50 · 289 阅读 · 0 评论 -
poj1681 Painter's Problem(gauss+dfs判定)
DescriptionThere is a square wall which is made of n*n small square bricks. Some bricks are white while some bricks are yellow. Bob is a painter and he wants to paint all the bricks yellow. But there i原创 2017-10-07 17:43:00 · 258 阅读 · 0 评论 -
poj 2065 SETI(gauss--->≡)
题目链接 翻译链接分析: 体面看了好久 简单来说,就是有若干个同余方程 a0*k^0+a1*k^1+a2^k^2+…=x1 1<=k<=n,x1已知 求解aitip再用高斯消元同余方程的时候 注意需要使用abs 有一句话的作用我搞不大懂if (a[i][j]*a[j][i]<0) t=-t;看来只能记住了回带求值的过程在高斯消元解异或方程时是没有的//回带求值for (int i=原创 2017-10-08 17:04:38 · 188 阅读 · 0 评论 -
poj2947 Widget Factory(gauss--->≡)
DescriptionThe widget factory produces several different kinds of widgets. Each widget is carefully built by a skilled widgeteer. The time required to build a widget depends on its type: the simple原创 2017-10-08 14:31:16 · 209 阅读 · 0 评论 -
poj1830 开关问题(gauss)
poj上的中文题面,大家一定要亲眼见证一下分析: 这道题和poj1222大同小异 还是异或方程求解 唯一不同的就是,需要计算方案数 这就和那些不确定的元素有关了 每个不定元都有两种选择:开或关 所以方案数就是 2^(不定元个数)tip怎么确定有没有解呢for (int i=n;i>=1;i--) if (a[i][n+1]!=0&&a[i][i]==0) ret原创 2017-10-07 16:39:45 · 239 阅读 · 0 评论 -
poj1222 EXTENDED LIGHTS OUT(gauss)
题目翻译 题目链接分析: 说实话,要是单看这道题, 是不大可能往解方程的方面想的由于只有两个状态 开是1,关是0 按开关是1,不按是0 这样的对应关系符合XOR的特点 所以我们又可以列出若干异或方程由于我们求的是一个矩阵 每个灯的状态由ta的初始状态以及周围按钮的状态有关 所以得到的方程就是 [x,y]^[x-1,y]^[x,y-1]^[x+1,y]^[x,y+1]=当前灯的状态原创 2017-10-07 15:26:59 · 258 阅读 · 0 评论 -
bzoj1923 [Sdoi2010]外星千足虫(gauss)
Description Input 第一行是两个正整数 N, M。 接下来 M行,按顺序给出 Charles 这M次使用“点足机”的统计结果。每行 包含一个“01”串和一个数字,用一个空格隔开。“01”串按位依次表示每只虫 子是否被放入机器:如果第 i 个字符是“0”则代表编号为 i 的虫子未被放入,“1” 则代表已被放入。后面跟的数字是统计的昆虫足数 mod 2 的结果。 由于 NASA的原创 2017-10-07 14:30:33 · 334 阅读 · 0 评论 -
bzoj1013 [JSOI2008]球形空间产生器sphere(gauss)
Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体。现在,你被困在了这个n维球体中,你只知道球面上n+1个点的坐标,你需要以最快的速度确定这个n维球体的球心坐标,以便于摧毁这个球形空间产生器。Input 第一行是一个整数n(1<=N=10)。接下来的n+1行,每行有n个实数,表示球面上一点的n维坐标。每一个实数精确到小数点 后6位,且其绝对值都不超过20000。Ou原创 2017-10-07 10:55:30 · 269 阅读 · 0 评论 -
bzoj4031 [HEOI2015]小Z的房间(矩阵树定理)
题目链接分析: 从题目限制我们可以得到这种图中有哪些边可以存在 因为任意两个房间之间都只有一条通路,那么得到的一定是一棵树 因此题目转化成:图的生成树个数 矩阵树定理即可tip注意柱子是不能算结点的 注意取模,矩阵中不要有负数 一开始我把矩阵消成了对角线形,结果WA了,消成上三角性就A了 看来保险起见,求解行列式的时候还是消成上三角#include<cstdi...原创 2018-03-29 09:14:21 · 207 阅读 · 0 评论