![](https://img-blog.csdnimg.cn/20201014180756780.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数学
overcastt
这个作者很懒,什么都没留下…
展开
-
欧拉公式
前置技能:复数复数 x=a+bix=a+bix=a+bi ,其中 aaa 和 bbb 是实数,iii 是虚数单位。在复数 a+bia+bia+bi 中,a=Re(x)a=Re(x)a=Re(x) 称为实部,b=Im(x)b=Im(x)b=Im(x) 称为虚部。 计算上 i0=1,i1=i,i2=−1,i3=−i,i4=1,i5=i...i0=1,i1=i,i2=−1,i3=−i,i...原创 2018-07-30 15:04:40 · 3740 阅读 · 0 评论 -
bzoj2705(洛谷P2303): [SDOI2012]Longge的问题 (欧拉函数)
题意求 ∑i=1ngcd(i,n)∑i=1ngcd(i,n)\sum\limits_{i=1}^ngcd(i,n)题解这题…真的要写一下了… 其实吧,结果就是 ∑i|ni∗φ(ni)∑i|ni∗φ(ni)\sum\limits_{i|n}i*\varphi(\frac{n}{i}) 然后!怎么求 varphivarphivarphi 这里给出一个超级厉害的做法!线筛+杜教...原创 2018-08-02 16:11:04 · 150 阅读 · 0 评论 -
bzoj 3240 (洛谷P1397): [Noi2013]矩阵游戏 (欧拉定理、数学)
题意F(i,j)=⎧⎩⎨1F(i,j)=a⋅F(i,j−1)+bF(i,1)=c⋅F(i−1,m)+di=1,j=1j≠1i≠1F(i,j)={1i=1,j=1F(i,j)=a⋅F(i,j−1)+bj≠1F(i,1)=c⋅F(i−1,m)+di≠1F(i,j)=\begin{cases}1& i=1,j=1\\F(i,j)=a\cdot F(i,j-1)+b& j\neq1\\F(i,...原创 2018-08-01 23:34:21 · 166 阅读 · 0 评论 -
LOJ #525. 「LibreOJ β Round #4」多项式
LOJ #525. 「LibreOJ β Round #4」多项式题意给定一个正整数 kkk,你需要寻找一个系数均为 000 到 k−1k−1k−1 之间的非零多项式 f(x)f(x)f(x),满足对于任意整数 xxx 均有 f(x)≡0(modf(x)≡0(modf(x)\equiv 0 (mod k)k)k) 题解我们知道欧拉定理是:若 a 与 n 互质,则 aϕ(...原创 2018-08-01 18:09:48 · 395 阅读 · 0 评论 -
bzoj2005(洛谷P1447): [Noi2010]能量采集
题意(1,1)到(n,m)每个点,与(0,0)点连线上的整点个数若为 kkk ,其权值即为 2k+12k+12k+1,求权值和。题解画个样例我们能发现,整点个数也就是 gcd(i,j)−1gcd(i,j)−1gcd(i,j)-1那么结果就是 ans=∑i=1n∑j=1m2(gcd(i,j)−1)+1=∑i=1n∑j=1m(2gcd(i,j)−1)=2∑i=1n∑j=1mgc...原创 2018-08-01 16:54:14 · 159 阅读 · 0 评论 -
bzoj4659 Lcm (莫比乌斯反演)
题意求 ∑i=1n∑j=1mlcm(i,j)(gcd(i,j)无平方质因子)∑i=1n∑j=1mlcm(i,j)(gcd(i,j)无平方质因子)\sum\limits_{i=1}^n\sum\limits_{j=1}^m lcm(i,j) (gcd(i,j) 无平方质因子)题解无平方质因子,即 μ2=1μ2=1\mu^2=1 因此式子可以转化为 ∑i=1n∑j=1mlcm(i...原创 2018-08-01 15:18:19 · 502 阅读 · 0 评论 -
bzoj2154 (洛谷P1829) : Crash的数字表格 (莫比乌斯反演)
bzoj2154 (洛谷P1829) : Crash的数字表格题意求 ∑i=1n∑j=1mlcm(i,j)∑i=1n∑j=1mlcm(i,j)\sum\limits_{i=1}^{n}\sum\limits_{j=1}^{m} lcm(i,j)题解显然 lcm(i,j)=ijgcd(i,j)lcm(i,j)=ijgcd(i,j)lcm(i,j)=\frac{ij}{gcd(...原创 2018-08-01 10:45:08 · 247 阅读 · 0 评论 -
莫比乌斯反演之随便写写...
莫比乌斯真的…难推…有时间..再推吧【捂脸】莫比乌斯反演呢,其实说的就是g(n)=∑d|nμ(nd)f(d)=∑d|nμ(d)f(nd)g(n)=∑d|nμ(nd)f(d)=∑d|nμ(d)f(nd)g(n)=\sum\limits_{d|n} {\mu(\frac{n}{d})f(d)}=\sum\limits_{d|n}\mu(d)f(\frac{n}{d}) 变下型就是f(k)=∑x=...原创 2018-05-03 13:23:32 · 212 阅读 · 1 评论 -
bzoj1923(洛谷P2447): [Sdoi2010]外星千足虫 (高斯消元解异或方程组)
题意n只虫子,单数只脚的是外星的虫子,双数只脚的是地球的虫子 m种探测,每种探测选了一些虫子(选了的为1,没选的为0),告诉你这些虫子总和是单数还是双数 问每只虫子是地球的还是外星的题解显然是解异或方程组裸题…代码#include <cstdio>#include <bitset>#include <algorithm&...原创 2018-08-04 08:32:39 · 180 阅读 · 0 评论 -
bzoj2466: [中山市选2009]树 (树形dp/高斯消元+深搜)
题意给定一棵树。每个点被按以后,与它直接相连的点会发生变换(亮变不亮,不亮变亮),每个点最初都是熄灭的,问最少多少次可以使所有点都变暗题解显然每个点最多按一次。 我们把相连的两个点连边,解异或方程组 这样所有的主元就是必须要按的 然后取枚举自由元选不选,深搜就可以了代码#include <cstdio>#include <cstri...原创 2018-08-04 08:31:59 · 247 阅读 · 0 评论 -
bzoj1013: [JSOI2008]球形空间产生器sphere (高斯消元)
题意已知 n+1n+1n+1 个点的坐标,求球心坐标题解我们可以通过 2−n2−n2-n 个点,与第一个点列出方程组 即 (xi,1−a1)2+(xi,2−a2)2+⋯+(xi,n−an)2=(x1,1−a1)2+(x1,2−a2)2+⋯+(x1,n−an)2(xi,1−a1)2+(xi,2−a2)2+⋯+(xi,n−an)2=(x1,1−a1)2+(x1,2−a2)2+⋯+...原创 2018-08-04 08:31:19 · 127 阅读 · 0 评论 -
bzoj3503(洛谷P3164): [Cqoi2014]和谐矩阵 (高斯消元解异或方程组)
题意有一个 n∗mn∗mn*m 的矩阵,每个格子可以填 111 或 000 现要求每个格子上下左右自己(五个格子)中有偶数个 111,请输出一种合理的解题解我们可以对每个格子向其周围五个点连边 然后又是解异或方程组解的裸题了 我们给自由元设为 1 ,这样就不会出现全零了代码#include <cstdio>#include <bit...原创 2018-08-04 08:30:46 · 181 阅读 · 0 评论 -
SPOJ104 Highways (矩阵树定理模板)
题目链接:https://vjudge.net/problem/SPOJ-HIGH 一个有 nnn 座城市的组成国家,其中一些城市之间可以修建 高速公路。现在,需要有选择的修建一些高速公路,从而组成一个交通网络。求有多少种方案,使得任意两座城市之间恰好只有一条路径题解显然是求生成树个数这里我们说一下矩阵树定理,用于计算生成树个数,即 无向图的生成树个数=度数矩阵-邻接...原创 2018-08-04 08:29:43 · 189 阅读 · 0 评论 -
bzoj4031(洛谷P4111) [HEOI2015]小Z的房间(矩阵树定理)
题解连一些边,然后矩阵树定理就可以了但是注意这里的模数不是质数! 因此不能用逆元啥的,需要辗转相除!!!代码#include <cstdio>#include <algorithm>using namespace std;#define ll long long#define N 90#define mod 1000000000ch...原创 2018-08-04 08:28:55 · 147 阅读 · 0 评论 -
泰勒展开
泰勒展开,本质上就是为了在某个点附近,用多项式函数去近似其他函数泰勒展开式 f(x)=f(x0)+f′(x0)1!⋅(x−x0)+f′′(x0)2!⋅(x−x0)2+...+f(n)(x0)n!⋅(x−x0)n+Rn(x)f(x)=f(x0)+f′(x0)1!⋅(x−x0)+f″(x0)2!⋅(x−x0)2+...+f(n)(x0)n!⋅(x−x0)n+Rn(x)f(x)=f(x_0)+\fr...原创 2018-07-30 15:04:10 · 6989 阅读 · 0 评论 -
bzoj3629(洛谷P4397): [JLOI2014]聪明的燕姿(数学+深搜)
题意求约数个数和为 nnn 的数有哪几个题解首先我们来看一个定理:约数和定理(下文会用到)对于任意一个大于1的正整数N可以分解正整数:x=Pa11Pa22…Pannx=P1a1P2a2…Pnanx=P_1^{a_1} P_2^{a_2}…P_n^{a_n} ,则由约数个数定理可知 xxx 的正约数有 (a1+1)(a2+1)(a3+1)…(an+1)(a1+1)(a2+1)...原创 2018-08-02 20:12:13 · 229 阅读 · 0 评论