高斯消元
alan_cty
蒟蒻一只
别打脸(⊙o⊙)哦
展开
-
求自然数幂和的各种方法(还有坑)
求∑i=1n1imodp\sum_{i=1}^{n}{i^k}\mod p高斯消元一个定理,k次方和一定可以由0~k-1次方和表示出来,设方程组接出来就好了。 O(k^3)倍增我们设f(n,k)=∑i=1nikf(n,k)=\sum_{i=1}^{n}i^k 怎么算呢? 我们采用分治思想。 如果n是奇数那么 f(n,k)=f(n−1,k)+nk原创 2016-03-18 19:00:45 · 1364 阅读 · 0 评论 -
【HNOI2013】游走
Description给出一张n个点,m条边的无向连通图。有一个人从点1开始随机游走,到点n结束。他每走过一条边就会得到其编号的分数。(可以重复走而重复得分)。现在让你安排每条边的编号,让他的得分期望值最小。求这个最小值。 n<=500Solution我们发现只要求出每条边的期望经过次数,然后从大到小排序,依次编号,就一定是最小值了。(即期望大的编号要小) 但是,每条边的期望经过次数要怎么求呢?原创 2016-03-21 21:02:02 · 1207 阅读 · 0 评论 -
高斯消元模板
高斯消元模板void gauss() { fo(i,1,n-1) { fo(j,i+1,n) if (abs(matrix[j].a[i])>abs(matrix[i].a[i])) swap(matrix[i],matrix[j]); fo(j,i+1,n) { db f=matrix[j].a[i]/matrix[原创 2016-03-21 21:03:46 · 668 阅读 · 0 评论 -
【NOI2017模拟6.3】子序列
Description n,q<=1e5Solution迟来的总结 比赛时只会O(n)Dp离线搞了60分 这个就是F[i]=2*F[i-1]-F[next[i]-1] 其中next[i]表示i前面第一个和i字符相同的位置正解的Dp长这样: 设s[i]=c,则F[i][c]=∑F[i-1][k],F[i][k]=F[i-1][k] 然后这样可以写成一个转移矩阵,并且是有逆矩阵的 所以我们原创 2017-08-04 22:31:26 · 364 阅读 · 0 评论 -
[51nod1323]完美平方
Description给出一个n*n的矩阵,求有多少种选择数的方案,使得每行每列均选择了奇数个数,并且选择的数的乘积为完全平方数 n<=20Solution完全平方数是什么,就是所有质因子的出现次数均为偶数 这样就变成了一堆奇偶的限制,异或方程组求自由元个数 高斯消元即可 感觉自己之前又打了假的高斯消元Code#include <cstdio>#include <cstring>#inc原创 2017-11-08 16:58:43 · 523 阅读 · 0 评论 -
图
Description这个世界是一个有向图,图中有n个点m条边且无重边无自环,每秒第i条边出现的概率是p[i]/100,一开始Samjia在1点,每一秒假设Samjia在点x上,那么Samjia要从存在的边中选一条来走,不可以不走,如果不存在可以走的边,那么Samjia就会掉出这个世界,假设Samjia绝顶聪明,问最后Samjia可以成功到达n的概率是多少。 nSolution首原创 2018-01-16 21:30:25 · 315 阅读 · 0 评论 -
[LOJ6360]复燃「恋之埋火」
Description古明地恋(koishi)和小石子(koishi)是好朋友。 旧地狱的空中散布着许多颗小石子。恋恋想找出一个位置,使得这个位置离最远的小石子的距离尽可能小。需要注意的是,这里的空间可能是高维空间。“在这幻想乡里,可不能被常理所束缚!”——东风谷早苗对于100%的数据,n≤20000,m≤5,0≤n\leq 20000,m\leq 5,0\leqn≤20000...原创 2018-08-17 21:58:26 · 540 阅读 · 0 评论