数论
SolarDomo
每天都被头像萌醒
展开
-
O(N) 的时间求 1 - N 关于 MOD 的逆元
新学的一个求逆元的方法:inv[i] = ( MOD - MOD / i ) * inv[MOD%i] % MOD证明:设t = MOD / i , k = MOD % i则有 t * i + k == 0 % MOD有 -t * i == k % MOD两边同时除以ik得到-t * inv[k] == inv[i] % MOD即inv[i] == -MOD / i * inv[MOD%i]即inv转载 2016-07-30 09:37:51 · 739 阅读 · 0 评论 -
POJ 1222 EXTENDED LIGHTS OUT 高斯消元
题目连接:http://poj.org/problem?id=1222 题意:一个矩阵上有灯,问怎么才能将矩阵上的所有的等都消掉我们假设Xi表示第i个按钮按还是不按,那么对于一个灯可以得到这样的一个方程 p表示对第i个灯的影响 L表示第i个等的状态 例如 对于第1个灯 p1 p2 p7 等于1 之后对这30个方程进行高斯消元既可看大佬代码写的版本://#include <bits/std原创 2016-11-02 19:56:12 · 316 阅读 · 0 评论 -
Vijos 1966 夜夜的旅游计划 高斯消元
题目链接:https://vijos.org/p/1966代码:#include <bits/stdc++.h>#define sf scanf#define pf printfusing namespace std;const int maxn = 200 + 5;const double eps = 1e-7;double adj[maxn][maxn];typedef double原创 2016-11-08 14:42:53 · 434 阅读 · 0 评论 -
BZOJ 3143 游走 高斯消元
题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3143中文题目。 F(v) 表示小Z在图上游走时,在v点走的次数 这样就可以根据图上的链接关系 构造出N个方程在求解既可代码:#include <cstdio>#include <cstring>#include <cmath>#include <algorithm>using n原创 2016-11-05 18:30:21 · 313 阅读 · 0 评论 -
HDU 5833 Zhu and 772002 高斯消元
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5833 题意:给你N个数,取其中K个数相乘,问有多少种选择方法似的乘机为一个完全平方数对每一个因子做一个方程,要求对方程组 高斯消元, 有唯一解方案数为0(因为所有的数都不选得到的一定是一个完全平方数) 方程有多解,假设自由变元有K个那么方案数为2的K次方 - 1(减去全都不选的方案)代码:#inc原创 2016-11-04 19:04:13 · 299 阅读 · 0 评论 -
POJ 1681 Painter's Problem 高斯消元 枚举自由变量
题目链接:http://poj.org/problem?id=1681和POJ 1753 一样代码:#include <iostream>#include <cstring>#include <algorithm>#include <cstdio>#define sf scanf#define pf printfusing namespace std;const int maxn =原创 2016-11-04 18:14:08 · 438 阅读 · 0 评论 -
POJ 1753 Flip Game 高斯消元 枚举自由变量
题目链接:http://poj.org/problem?id=1753构造号方程之后,先判断是否有解,之后在位运算枚举自由变量的值代码:#include <cstdio>#include <algorithm>#include <cstring>#define pf printf#define sf scanfusing namespace std;const int maxn = 20,原创 2016-11-04 18:12:09 · 446 阅读 · 0 评论 -
高斯消元法,高斯约旦消元法
1、高斯消元法的算法(设akk(k)不等于0)将非零阵A=(aij)m×n,经过行初等变化,变为上三角矩阵。步骤:当m>n 当m当m=n 例1A= 此方法常用于解线性方程组和矩阵的秩的计算。如例1中矩阵A的秩r(A)=3。2、高斯约旦消元法的算法转载 2016-11-04 15:57:20 · 1929 阅读 · 0 评论 -
C++大整数模板 BigInteger
网上找的模板#include <cstdio>#include <algorithm>#include <cstring>using namespace std;struct BigInteger{ int A[25]; enum{MOD = 10000}; BigInteger(){memset(A, 0, sizeof(A)); A[0]=1;} void转载 2016-08-13 20:56:00 · 2745 阅读 · 0 评论 -
POJ 1830 开关问题 高斯消元求自由变元
题目链接:http://poj.org/problem?id=1830代码:#include <algorithm>#include <cstdio>#include <cstring>#define sf scanf#define pf printfusing namespace std;const int maxn = 50;typedef int Matrix[maxn][max原创 2016-11-03 10:08:53 · 328 阅读 · 0 评论 -
POJ 2976 Dropping tests 01分数规划
题目链接:http://poj.org/problem?id=297601分数规划主要是理解F(L)函数 比如我们要求R的最大值的时候 ,如果F(L)能取到一个大于0的数,那么就可以得到在这个选择的情况下R的值大于L的值顺便提一下:如果和别人的对拍没有任何问题,但POJ还是WA的话,看看是不是输出用的%lf,POJ的G++输出用lf会wa,以上 二分法:#include <cstdio>#in原创 2016-10-11 09:17:36 · 246 阅读 · 0 评论 -
HDU 5690 All X 矩阵快速幂
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=5690 题意:问一个全部由X组成的M位数 对K取余之后是否等于c 构造一个公式 之后用矩阵快速幂 解即可! #include <bits/stdc++.h>#define sf scanf#define pf printfusing namespace std;typedef long lon原创 2016-10-10 14:15:31 · 334 阅读 · 0 评论 -
POJ 2778 DNA Sequence AC自动机 + 可达矩阵 + 矩阵快速幂
题目链接: http://poj.org/problem?id=2778代码:#include <cstdio>#include <iostream>#include <cstring>#include <queue>#define sf scanf#define pf printf//#define LOCALusing namespace std;typedef long lon原创 2016-08-13 13:19:23 · 387 阅读 · 0 评论 -
POJ 3070 Fibonacci 矩阵快速幂
题目链接:http://poj.org/problem?id=3070用矩阵求斐波那契数列的第N项 第一次用矩阵快速幂 和乘法快速幂一样 看模板在纸上模拟一下过程就清楚了 矩阵快速幂不建议用递归写 代码:#include <cstdio>#include <iostream>#include <cstring>#define sf scanf#define pf printfus原创 2016-08-13 11:14:43 · 298 阅读 · 0 评论 -
POJ 3233 二分求等比数列 矩阵快速幂
题目链接:http://poj.org/problem?id=3233 题意: 求 S =( A + A^2 + A^3 + … + A^k) mod m二分原理(转自AcDreamer大佬): (1)当时, (2)当时,那么有 (3)当时,那么有 代码:#include <cstdio>#include <iostream>#include <cstring>#define sf原创 2016-09-11 22:26:33 · 276 阅读 · 0 评论 -
CODRFORCES 27E 构造
题目链接:http://codeforces.com/problemset/problem/27/E代码:#include <cstdio>#include <iostream>#include <cstring>#define sf scanf#define pf printfusing namespace std;int prim[] = {2,3,5,7,11,13,17,19,2原创 2016-09-10 22:44:41 · 333 阅读 · 0 评论 -
POJ 1006 Biorhythms 中国剩余定理
题目链接:http://poj.org/problem?id=1006 直接套板子 代码://#include <bits/stdc++.h>#include <cstdio>#include <cstring>using namespace std;int A[10],B[10];//扩展欧几里得算法int Extended_Euclid(int a,int b,int &x,int原创 2016-12-07 22:45:45 · 429 阅读 · 0 评论