![](https://img-blog.csdnimg.cn/20201014180756913.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
--------数学知识--------
文章平均质量分 65
Faithfully__xly
这个作者很懒,什么都没留下…
展开
-
0818-求欧拉函数的两种方法
【欧拉函数】在数论,对正整数n,欧拉函数是小于或等于n的正整数中与n互质的数的数目(φ(1)=1)。【三性质】如果i是质数,那么phi[i]=i-1; 如果a mod b==0&&b是质数,那么phi[a*b]=phi[a]*b 如果a mod b!=0&&b是质数,(也就是说a,b互质),那么phi[a*b]=phi[a]*phi[b]【证明】性质一:i为质数,显然1~i-1中...原创 2018-08-18 20:56:29 · 216 阅读 · 0 评论 -
数论 质因数分解 - 樱花(BZOJ 2721)
题意输入n,求有多少对正整数x,y,满足Analysis设n!=z,y=z+d1/x+1/y=1/z1/x+1/(z+d)=1/z(x+z+d)/(xz+dx)=1/zz(x+z+d)=xz+dxz^2+dz=dxx=z^2/d+z发现就是求z^2的约数个数( – By hzwer)Code#include<bits/stdc++.h>#define P ...原创 2018-10-31 07:55:28 · 331 阅读 · 0 评论 -
数论 - GREEN友好的地址对amicbl(USACO_2001_FALL_GREEN)
Analysis手敲暴力过了30分暴力思路很好想:就是枚举L~H中的每一个点,然后算出其因子之和,再判断一下。算因子之和,就是O(n)枚举 1~n-1 中所有数然后看是否为其约数显然时间复杂度就卡在这里了我们该怎么想办法来优化一下呢?既然从数推因子之和会爆时间,那么我们换个思路,枚举因子,看他能被多少个数包含因为因子最多也只有H/2个。然后再枚举这个因子在L~H中的倍数,将其ans...原创 2018-10-31 08:02:38 · 154 阅读 · 0 评论 -
【提高组NOIP2017】小凯的疑惑(附证明)
小凯的疑惑描述小凯手中有两种面值的金币,两种面值均为正整数且彼此互素。 每种金币小凯都有 无数个。 在不找零的情况下, 仅凭这两种金币,有些物品他是无法准确支付的。 现在小 凯想知道在无法准确支付的物品中,最贵的价值是多少金币? 注意:输入数据保证存在 小凯无法准确支付的商品输入输入数据仅一行, 包含两个正整数 a 和 b, 它们之间用一个空格隔开, 表示小凯手中金币的面值输出输出文件...原创 2018-10-31 13:49:39 · 4784 阅读 · 0 评论 -
高斯消元 - 求解线性方程组
距离NOIP只有9天了可是我连高斯消元都还不会,嘤嘤%>_<% 快学习快学习所以这到底是个什么东东呢???就是我们平时解一次方程组的方法,只不过用程序实现罢了(系数是具体存在矩阵(数组)里的)其思路大致不发生改变,就是挨个消元,消到最后就可以回代了举个栗子:现在我们的系数长这样[xyzval2111621−1−2217]\begin{bmatrix} x &amp...原创 2018-10-31 19:57:22 · 250 阅读 · 0 评论 -
高斯消元 - [JSOI2008]球形空间产生器sphere
传送门Analysis第一次写高斯消元,式子倒是很快就推出来了然而,高斯消元板子写挂了……注意最后一项啊!!!!随便推一下:设球心的坐标为(x1,x2,x3,....xn)(x_1,x_2,x_3,....x_n)(x1,x2,x3,....xn)(从二维往高维推)就会发现:设圆心(x1,y1)(x_1,y_1)(x1,y1),给定的点(a,b)(a,b)到圆心的距离为...原创 2018-10-31 20:50:29 · 149 阅读 · 0 评论 -
容斥原理+dfs剪枝 - ost数(WOJ 2592)
cost数描述“给你一个有n个正整数的数列{an}。一个正整数x若满足在数列{an}中存在一个正整数ai,使x≡17(mod ai),那么x就是一个‘cost数’。请问1到m的正整数中,有多少个‘cost数’?”输入第一行两个正整数n和m,意义见问题描述。第二行n个正整数,分别为数列{an}中的n个数。输出输出一个整数,表示1到m中“cost数”的个数。样例输入3 10018 ...原创 2018-11-01 18:13:21 · 159 阅读 · 0 评论 -
NOIP2016提高组DAY2T1 - 组合数问题
描述组合数表示的是从n个物品中选出m个物品的方案数。举个例子,从(1,2,3) 三个物品中选择两个物品可以有(1,2),(1,3),(2,3)这三种选择方法。根据组合数的定 义,我们可以给出计算组合数的一般公式:其中n! = 1 × 2 × · · · × n小葱想知道如果给定n,m和k,对于所有的0 <= i <= n,0 <= j <= min(i,m)有多少对...原创 2018-11-01 20:34:15 · 515 阅读 · 0 评论 -
广义欧拉定理+打表 - Last Digits(POJ2720)
传送门Analysis好啊。。。。看起来蛮递归降幂,蛮容易hhhh……时间过不去原来要打表还要记录(多组数据)因为指数与模数不一定满足互质,所以在降幂的时候要分类讨论(指数是否大于φ(mod))由于直接讨论不知道怎么做所以……1e7以下的打表算出来剩下的讨论即可Code#include<cstdio>#include<algorithm>#in...原创 2019-08-23 16:45:53 · 222 阅读 · 0 评论 -
容斥原理 - 分糖(SOJ 747)
分糖题目描述有 N 个(相同的)糖果,M 个(不同的)小朋友。M 和 N 满足:1≤M≤N≤100000(105)。要求:1.每个小朋友都至少有一个糖果。2.不存在正整数 X(X>=2),使得每个小朋友的糖果数都是 X 的倍数。3.糖果不能剩余。求分糖方法总数。答案模 1000000007(109+7)输入格式第一行为数据组数:T<=100000。接下来 N 行,每...原创 2019-08-23 21:17:46 · 349 阅读 · 0 评论 -
第二类斯特林数 - Rank(HDU 2643)
传送门Analysis枚举总共的名次 i ,将 n 个人分到 i 个非空集合中 <–裸的第二类斯特林数然后对于每一种方案还需要乘以 i 的全排列ans=Σi=1nS2(n,i)∗i!ans=\Sigma_{i=1}^{n}S2(n,i)*i!ans=Σi=1nS2(n,i)∗i!Code#include<bits/stdc++.h>#define re reg...原创 2019-08-20 14:50:30 · 143 阅读 · 0 评论 -
第一类斯特林数+组合数 - Pole Arrangement(UVA 1638)
传送门Analysis大多数题解都提供的是dp做法这里推广一个不知哪个神犇想出来的神仙做法O(1)计算可还行?除了最高的那根杆,左边还能看到 l - 1 根,右边还能看到 r - 1 根。把这些杆子和他们挡住的杆子看做一个圆排列,那么共有[l+r−2n−1][^{n-1}_{l+r-2}][l+r−2n−1]种分配方案。左边要选l−1l-1l−1根故答案为ans=[l+r...原创 2019-08-20 19:32:25 · 156 阅读 · 1 评论 -
组合数学の学习笔记
组合数常用公式Cnm=n!(n−m)!∗m!C_n^m=\frac{n!}{(n-m)!*m!}Cnm=(n−m)!∗m!n!Cnm=Cnn−mC_n^m=C_n^{n-m}Cnm=Cnn−mCnm=Cn−1m−1+Cn−1mC_n^m=C_{n-1}^{m-1}+C_{n-1}^{m}Cnm=Cn−1m−1+Cn−1mCn0+Cn1+Cn2+Cn3...+Cnn=2nC_n...原创 2019-08-20 19:33:35 · 590 阅读 · 0 评论 -
第二类斯特林数 - Rhyme Schemes(POJ 1671)
传送门第一次应用Stirling,并没有看出来不过看了题解后,直呼太妙Analysis题意:n行诗存在的韵律组合数。思路:第二类Stirling数是把n个元素放入k个等价类的方法数目(每个等价类都不能为空)。注意到此题的'行数'>='韵律数',可以把'行数'抽象为n个元素,'韵律数'抽象为k个等价类.所求即为:把n行放入k个‘韵律’中的方法数( 其中:1=< k...原创 2019-08-22 09:38:01 · 181 阅读 · 0 评论 -
第二类斯特林数 - Push Botton Lock(POJ 3088)
传送门Analysis手动分析一波样例就出来了枚举参与组合的数 i (因为题目中说了可以不用完所有的数),再枚举将其分为的集合数 j (这就是第二类Stirling数的标准长相,因为对于每一个集合中的数我们不考虑顺序),那么最终的答案就是ans=(in)∗[ji]∗j!ans=(_i^n)*[^i_j]*j!ans=(in)∗[ji]∗j!()–>组合数[]–>第二类斯...原创 2019-08-22 10:12:27 · 230 阅读 · 0 评论 -
快速计算n!中素因子的个数
Learning我们来一个样例说明一下:1 2 3 4 5 6 7 8 我们求得在8!中2的个数 1 1 1 首先我们先计算出2的倍数的个数:8/2=4 1 1 其次我们计算出4的倍数的个数: 8/4=2 1 最后我们解出第三层的2的个数: 8/8=1我们把4+2+1=7,所以一共7个2出现了。即:cnt(x)=[n/(x1)]+[n/(x2)]+[n...原创 2019-08-22 14:32:09 · 707 阅读 · 0 评论 -
1025 -康托展开&逆康托展开 - Cow Line[USACO11FEB]
传送门分析本来以为好难好难结果除了不知道为什么,还是很简单的嘛[笑哭]我们就记结论,并理解代码就好了康托展开先上定义:康托展开是一个全排列到一个自然数的双射,常用于构建hash表时的空间压缩。设有n个数(1,2,3,4,…,n),可以有组成不同(n!种)的排列组合,康托展开表示的就是 当前排列组合在n个不同元素的全排列中的名次。再讲作用:对于一个具体的1~n的排列P,康托...原创 2018-10-25 16:09:47 · 210 阅读 · 0 评论 -
1018 - 数学期望 - 百事世界杯之旅(luogu 1291)
题目描述 输出如果输出是整数,那么直接输出一个整数如果有分数,分数必须不可再约。 分3行打印 第一行是分子,前面有整数部分的位数+1个空格第二行是整数部分+1个空格+分母位数个、‘-’第三行是分母,前面有整数部分的位数+1个空格 标签SHTSC2002 牢骚这个输出很扯啊……是真的在搞笑 分析假设现在已经有k个球星的名字,那么要使球星的 ...原创 2018-10-18 21:15:44 · 149 阅读 · 0 评论 -
1018 - 数论之扩展欧几里得 - 吃蛋糕
吃蛋糕描述Beny 想要用蛋糕填饱肚子。Beny 一共想吃体积为 c 的蛋糕,他发现有两种蛋糕可以吃,一 种体积为 a,一种体积为 b,但两种蛋糕各有特色。Beny 想知道他一共有多少种不同吃法, 使得他恰好可以填饱肚子。输入第一行一个 t接下来 t 行,每行三个正整数 a,b,输出对于每个 a,b,c,输出一个整数表示有几种不同吃法样例输入 132 3 43...原创 2018-10-18 20:32:20 · 266 阅读 · 0 评论 -
0830-扩展欧几里得算法+例题
phew~终于看懂了,以前一直以为很高深很高深的算法,结果还是很简单嘛-->参考资料<--证明什么的大家都写的很好啊,蒟蒻就不再bibi了,在这里解释一下代码吧,我看了好多博客都只讲了思路和证明,像宝宝这种代码能力不强的就只好自己想想想想想void exgcd(int a,int b,int &x,int &y){ if(b==0) { x=1;...原创 2018-08-30 20:59:38 · 7202 阅读 · 4 评论 -
【2018/08/30】T2-矩阵快速幂+扩展欧几里得-k-斐波那契(SDOJ 3739)
k-斐波那契题目描述:输入格式:一行两个整数 n,P输出格式:从小到大输出可能的 k,若不存在,输出 None样例数据:输入5 5输出2备注:【样例解释】 f[0] = 2f[1] = 2f[2] = 4f[3] = 6 mod 5 = 1f[4] = 5 mod 5 = 0f[5] = 1【数据范围】30% 的数据保...原创 2018-08-30 22:26:28 · 187 阅读 · 0 评论 -
【2018/08/29】T1-线性筛素数(水题)-pmatrix(SDOJ 3727)
哥德巴赫矩阵 【问题描述】 根据哥德巴赫猜想(每个不小于 6 的偶数都可以表示为两个奇素数之和),定义 哥德巴赫矩阵 A 如下:对于正整数对(i,j),若 i+j 为偶数且 i,j 均为奇素数,则Ai,j = 1, 否则Ai,j = 0。现在有若干询问(x1,y1,x2,y2),你需要回答下列式子的值 【输入】 第一行一个整数 m 接下来 m 行,每行四个整...原创 2018-08-29 13:06:56 · 181 阅读 · 0 评论 -
0829-全排列函数-POJ 1256
传送门 大致题意给定字典序: 'A'<'a'<'B'<'b'<...<'Z'<'z'.让你对输入的一串字母,按照上面定义的大小关系,从小到大输出所有的全排列情况 分析普遍做法: dfs但今天get了一个超棒的求全排列的函数,虽然复杂度和手写的 dfs 差不多,但一下子代码就简单起来了,多棒啊!-->参考<--...原创 2018-08-29 16:50:55 · 150 阅读 · 0 评论 -
0901-证明欧拉函数phi的积性
【写在前面】一个数若既与 m 互质又与 n 互质,那么他便和 m*n 互质完全剩余系:一个整数的集合,对 m 取模后,余数遍历了 0; 1; 2; 3; ……m欧拉函数:phi(n)-->不超过n且与n互质的整数的个数特别的:phi(1)=1现在我们需要证明phi(m*n ) = phi (m )*phi(n)【m,n互质】这个矩阵里列举了从1到m*n的所...原创 2018-09-01 17:19:35 · 965 阅读 · 0 评论 -
0901-Miller_Rabin素数测试算法+例题
看了好久终于把这个Miller_Rabin搞懂了,觉得自己棒棒哒~~~最后是在下面那篇博客里搞懂的,这里推荐给大家-->参考<--【写在前面】费马定理 and 二次探测<证明来源>然后费马定理是一个必要条件,也就是说素数一定满足这个定理,但满足这个定理的不一定是素数,比如说Carmichael数(我没研究过,有兴趣的同学自己百度吧,反正这种数就是反例)。...原创 2018-09-01 21:53:07 · 1342 阅读 · 6 评论 -
0902-求解逆元的三种方法(附证明)
【逆元】若gcd (a,m) = 1,且 ab = 1 (mod m),则称 b 为 a 模 m 意义下的逆元。显然,模意义下的除法可以用乘 逆元来代替。【三种方法求解】方法一:利用费马小定理,当 p 为质数,有 a^p−1=1 (mod p)很容易发现a * a^p-2 = 1(mod p),那么 a 在模 p 意义下的逆元为 a^p-2然后用快速幂计算即可板子in...原创 2018-09-02 15:04:07 · 1767 阅读 · 0 评论 -
0902-区间素数筛-POJ 2689 Prime Distance
传送门大致题意给定 L, R,保证 L, R在int范围以内, R − L ≤ 10^6 ,问 [L, R] 中距离最近的两个相邻质数,和距离最远的两个相邻质数。区间内若没有相邻的质数则输出:There are no adjacent primes.分析这个范围告诉我们任何一种我们熟知的求素数的方法都会TLE(except Miller_Rabin),那我们就来两次素数筛先筛出...原创 2018-09-02 16:54:21 · 131 阅读 · 0 评论 -
0902-数论之组合数-POJ 3421 X-factor Chains
【写在前面】明天是2018年9月3号,是博主开启高中生涯的第一天。进入高中的自己请一定要加倍努力,平衡好文化课和竞赛之间的关系,尽自己的努力不留遗憾就好了,你知道的竞赛是高风险的一项任务,但同时也是高回报的。为什么不去试一试呢?基础再薄弱只要肯练习,就总会补上去的,亲爱的自己,你常常会因为自己不如别人那般总能想到正解而沮丧,懊恼,甚至怀疑人生,可真的没有必要。你一步一步来吧,不管别人怎么说,...原创 2018-09-02 22:08:58 · 124 阅读 · 0 评论 -
0907-欧拉函数-POJ 3090
传送门大致题意给出范围为 (0, 0) 到 (n,n) 的整点,你站在原点处,问有多少个整点可见。分析哈,我喜欢简简单单的题显然一个点(x,y)若不能被看见,则它一定被(x/g,y/g)这个点遮挡了,所以对于没有被遮挡的点,其gcd(x,y)=1。那么我们用欧拉函数去搞,由于欧拉函数只能处理小于等于x的质数,我们就在答案处理的时候多加一次即可。那么这样就完了 吗?不...原创 2018-09-07 13:59:57 · 144 阅读 · 0 评论 -
【2018/10/04】T2 - 找规律 - 排队
排队描述在成都某中学有m个男生与n个女生排队,这个学校的女生比较古怪,从某个位置(包含这个位置)开始往前数,男生的数量超过了女生的数量,女生会感觉不安全,于是会大叫起来,为了构建和谐校园,安排队伍时应该避免这样的情况。请你计算出不会引发尖叫的排队方案的概率。(排队方案不同定义:当且仅当某个某个位置人不一样,如男生A、男生B ,与男生B、男生A ,2个排列是不同方案)输入第一行1个整...原创 2018-10-04 15:25:08 · 177 阅读 · 0 评论 -
【2018/10/05】T1 - 分解质因数 - 阶乘
我好菜啊……好菜啊……好菜啊……阶乘 (fact)【问题描述】 有 n个正整数 a[i],设它们乘积为 p,你可以给 p乘上一个正整数 q,使 p*q刚好为正整数m的阶乘,求m的最小值。 【输入】 共两行。 第一行一个正整数n。 第二行n个正整数a[i]。 【输出】 共一行 一个正整数m。 【输入样例】 1 6 【输出样例】 3 样例解释: 当p=6,q=1时,p*...原创 2018-10-05 17:15:45 · 225 阅读 · 0 评论 -
1012 - 曼哈顿距离&切比雪夫距离
什么是切比雪夫距离?什么是曼哈顿距离?傻傻分不清,没关系,看:曼哈顿距离设平面空间内存在两点,它们的坐标为(x1,y1),(x2,y2)则dis=|x1−x2|+|y1−y2|即两点横纵坐标差之和切比雪夫距离设平面空间内存在两点,它们的坐标为(x1,y1),(x2,y2)则dis=max(|x1−x2|,|y1−y2|)即两点横纵坐标差的最大值比如这...原创 2018-10-12 16:01:18 · 512 阅读 · 0 评论 -
【2018/10/12】校内模拟赛
题一题二题三 分析(题一) 一道模拟题但模拟也是有很多坑的……写出一道题后,一定要自己试很多很多很多组不同的数据比如这道题,几乎我试的每一组数据都找出了我程序的bug,好险~ 分析(题二)好题!!!!应该是和松鼠聚会这道题本质一样的首先我们需要推导出任意两个点(x,y) (a,b)之间能够到达彼此所需要走的最短步数,由于可以斜着走,那么我们肯定会...原创 2018-10-12 20:20:05 · 194 阅读 · 0 评论 -
1014 - 结论题&dp? - Mr. Young's Picture Permutations(POJ 2279)
传送门 题意给出n行,每行有人数限制num[i],并且num[i]>=num[i+1],总人数暂且称为tot=∑num[i],把1~tot这些数字填入矩阵,使得矩阵满足每行单调递增,每列单调递增,求满足要求的矩阵数目 分析一开始是打着做dp的名号,发现了这道题在《算法竞赛》这本书上,是这样讲的线性dp但实际上,由于这道题空间限制比较严,我们dp的话会爆...原创 2018-10-14 19:16:53 · 107 阅读 · 0 评论 -
组合数奇偶性的判断(附证明)
Learning方法一:计算一下,然后看它的奇偶性;但是这个时间以及数据范围上都不允许;方法二对于给定C(n,m),检查n!中2因子的个数与m!和(n-m)!中2因子个数和的关系,假设n!中2因子个数为a,m!中2因子个数为b,(n-m)!中2因子个数为c,则显然有a>=(b+c);并且当a==b+c时,一定为奇,否则为偶。方法三由方法2可以很容易(稍后给出证明)地看出,n!...原创 2019-08-22 15:12:15 · 2208 阅读 · 2 评论