![](https://img-blog.csdnimg.cn/20201014180756927.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数学
文章平均质量分 71
MatrixYg
暂无
展开
-
【BZOJ 4001】[TJOI2015]概率论
1.其实这是一个纯数学题。题目链接。额,网上很多人都说好猜出来,这个我是不敢苟同,至少我觉得我是猜不出来,打表也打不出来。因为毫无规律可言。2.其实这个题就是再说n个节点的二叉树叶子节点的期望是多少,ans=n*(n+1)/(2*(2*n-1)),n是节点个数。首先n个节点的二叉树种类是C_n个,C就是卡特兰数,这个题其实算一下所有的叶子节点有多少个就行了.但是很难算。考虑DP[i][j]代表...原创 2019-11-06 23:49:16 · 227 阅读 · 0 评论 -
【一些常用的递推通项公式的求法】
1.前言:在很多题目里会涉及到递推公式的使用,然而大多数递推公式都是有很强的前后依赖性,第n项的求解依赖于n前边的几项。比如大名鼎鼎的fib数列。但是很多时候,我们要O(1)公式。本文简单的谈一下一些常见的递推公式的求解。2.一阶线性递推。形如: 这种序列,显然是可以构造一个...原创 2019-10-17 20:52:15 · 2457 阅读 · 0 评论 -
【luogu 3868】中国剩余定理模板题
1.题目链接。考虑这样一个问题,存在一个同余方程组,其中模mi两两互质,求出满足这个方程组的最小的非负整数。 中国剩余定理就是干这件事的,其实也没什么神奇的,就是通过不断合并解得到最终答案。 #include<iostream&...原创 2019-06-16 16:44:46 · 324 阅读 · 0 评论 -
【HDU 5885】XM Reserves
1.题目链接。题目大意:一个n*md的矩阵,每个格点有一个p[i[j],让你选取一个位置,对于所有在以这个点为中心,半径为r的圆内的点的Σp(i,j)/(d+1)最大, d为点到圆心的距离.2.我们可以对于每个点统计它能对那些格点造成的贡献,然后我们遍历每一个格点就能找到答案了我们从坐标变换的角度来思考,假如一个点的坐标为(x,y),圆心跟它的坐标差为(dx,dy)那么圆心的坐标就是(x+dx...原创 2019-06-10 12:02:27 · 212 阅读 · 0 评论 -
【HDU 4609】3-idiots
1.题目链接。题目大意:给出n根棍子,现在从中任意选三种,能够组成三角形的概率。2.其实就是说在这些棍子能够组成多少种三角形。我们预处理a+b,然后枚举第三条边c.预处理a+b采用FFT加速即可。然后统计答案的时候,注意把不合法的减了。#include<bits/stdc++.h>using namespace std;#pragma warning(disable:49...原创 2019-06-09 22:05:51 · 166 阅读 · 0 评论 -
【BZOJ 1013】球形空间产生器sphere
1.题目链接。题目大意:给出一个n维球面上的n+1个点的坐标,求解这个球心的坐标。2.首先,如果解存在,一定是唯一确定的。数据保证解存在,所以不用管无解的情况。设球心是(x1,x2,x3,x4...xn)。根据球的性质我们可以得到n+1个方程,然后逐项相减,就可以得到一个n项的线性方程组。为了方便,就拿第一个方程作为基准,每个方程都减去它。做出增广矩阵后Gauss校园完事。#includ...原创 2019-06-19 15:37:02 · 246 阅读 · 0 评论 -
【HDU 5446】Unknown Treasure
1.题目链接。题目大意:已知C(n,m)%p[i]=a[i],求C(n,m)%M.其中M=p1*p2*..pn,p都是素数。2.Lucas定理+中国剩余定理裸题了。首先使用Lucas定理求出每个a[i],然后就可以得到一个模方程组,采用中国剩余定理合并即可。#include<bits/stdc++.h>using namespace std;typedef long lo...原创 2019-06-19 14:51:16 · 293 阅读 · 0 评论 -
【HDU 5917】Instability
1.题目链接。题目大意:给定一个无向图,求这个图有多少个点的集合,每个集合至少包含一个关于这个图的独立集或者团。并且集合的元素个数要大于等于3.2。其实这个题目就是拉姆齐定理的应用,根据拉姆齐定理,6个人中一定有三个人互相认识或者互相不认识。这个认识和不认识就是我们这道题目中的边,然后更具这个题目我们可以算出来大于等于六个人一定满足的,总共的选择方案是(C(n,0)+C(n,1)....C(n...原创 2019-05-23 23:35:17 · 237 阅读 · 0 评论 -
【HDU 5916】Harmonic Value Description
1.题目链接。是一个构造的题目,其实题目中已经有了一定的提示:2*k<=n.首先我们知道对于自然数1234...n这个排列一定是权值最小的,因为他们两两互质,gcd(ai,ai+1)=1.不可能比他更小,它的权值和是n-1。那么第二小是多少?第二小最少是n。那么有没有可能是n??其实很简单,我们知道gcd的几个性质:gcd(k,2*k)=k,所以我们每次只需要调换这个k的位置,就会让这个序列...原创 2019-05-22 21:15:38 · 160 阅读 · 0 评论 -
【焦作网赛】Participate in E-sports
1.题目链接。题意简洁:判断n和C(n,2)这两个数是不是都是平方数。2.由于n很大,但是理论上来说log的算法也是过不掉的啊,只有O(1)的算法才能过。然后尝试着用JAVA写了一发LOG的算法,是可以过的。import java.util.Scanner;import java.math.BigInteger;class Main{ public static boolean ...原创 2019-05-06 22:18:49 · 211 阅读 · 0 评论 -
【HDU 6298】Maximum Multiple
1.题目链接。找一下规律,其实看看看i,j,k的性质,发现最优的比例一般都是1:1:1或者1:1:2.#include<bits/stdc++.h>using namespace std;#define ll long long#pragma warning(disable:4996)int main(){ int T; scanf("%d", &T); ...原创 2019-05-14 10:56:45 · 129 阅读 · 0 评论 -
【HDU 6507】 Store The Matrix
1.题目链接。题目大意:把一个矩阵分解成若干个矩阵的乘积。一个m*n的矩阵需要m*n个元素,但是分解之后,是所有小矩阵元素个数之和。问最少需要多少个元素才能把这个矩阵表示出来。2.看到这个题有点懵逼,难道不是固定的mn个???,然后看了一眼样例,其实还是比较良心的,看到了样例种矩阵的秩为1,突然发现(猜想)如果一个矩阵式满秩的,是没办法减少元素数量的。但是如果不是满秩的,那些零行似乎就可以在其...原创 2019-05-04 23:58:17 · 331 阅读 · 0 评论 -
【HDU5974】思维
1.题目链接。题目大意:x+y=a,lcm(x,y)=b.求x,y.不存在就输出:No Solution.2.分析:首先, 这里我们得到了一个类似一元二次方程的东西,如果gcd(x,y)是一个定值,那么这题就解决了。基于这个猜想,我们发现:gcd(x,y)=gcd(a,b)....原创 2019-04-19 17:01:41 · 218 阅读 · 0 评论 -
关于gcd的一些证明思路。
我们都知道gcd(n,i)=gcd(n,n-i)对于所有的0<=i<=n都是成立的。那么我们来证明一下这个结论: 首先假设gcd(n,i)=k,gcd(n,n-i)=m,那么: 所以我们可以得出:gcd(n,i)=gcd(n,n-i)。同理,我们也可以得出:gcd(n,i)=gcd(n,n+i)其实我们用这个思路可以证明很多关于gcd的东西。...原创 2019-04-19 00:19:01 · 400 阅读 · 0 评论 -
【POJ 2429】大数分解+二进制枚举
1.题目链接。题目大意:给出两个数的gcd和lcm,求这两个数和最小是多少?2.分析:。那么我们可以简单的推导一下: 这里n是一个整数,并且: 说明n一定是由两个互质的数乘起来的。那么我们枚举n的因子,然后找一下最小...原创 2019-04-18 23:24:27 · 362 阅读 · 0 评论 -
【POJ 1284】Primitive Roots
1.题目链接。求一个奇素数原根的数量。结论题,如果一个数存在原根,那么数量是phi(phi(n))个。奇素数一定是存在原根的。#include<iostream>#include<stdio.h>#include<cstdio>#include<cstring>using namespace std;const int N = 755...原创 2019-06-11 19:09:46 · 157 阅读 · 0 评论 -
【HDU 2619】Love you Ten thousand years
1.题目链接。原根的一个性。当a的阶等于phi(m)的时候,说明a是m的一个原根。首先这里原根的存在性是不用判断的,因为奇素数一定存在原根。当a是原根时,a^0,a^1,a^2....a^(odr(a))构成了一个关于m的简化剩余系。简单点就是说这个集合和{1,2,...odr(a)}这个集合相等。又因为m时质数,所以phi(m)=m-1。所以这个得到的简化剩余系就是{1,2,...m-1},只要...原创 2019-06-12 12:12:12 · 317 阅读 · 0 评论 -
【HDU 4992】Primitive Roots
1.题目链接。原根的基本求法。根据下面的性质求: 流程:判断是不是存在原根?---->求出最小----->从最小得到全部。#include<bits/stdc++.h>using namespace std;typedef long long LL;#pragma warn...原创 2019-06-12 13:00:07 · 275 阅读 · 0 评论 -
【BBP 算法】HDU 6217
1.最近看到有学弟在说BBP算法,然后就写一下这个东东。2.BBP其实是三个人,然后这个算法是干嘛的嘞,其实是快速的求某些无理数的第n为数字或者第n位以后的一串数字。比如求PI,log2,log3等等。下面来说明一下它是怎么求PI的在16进制下的第n位数字的。(什么,为什么是16进制?,看完就明白了)3.首先给出一个PI的级数求和公式: ...原创 2019-10-07 22:15:45 · 1527 阅读 · 0 评论 -
【十进制快速幂】
1.国庆的牛客没有打,于是乎坐在车上觉得还挺无聊,就看了看题目其中有一个这样的签到题,矩阵快速幂加了一点点的限制,让n十分大,其实这是快速幂的另一种使用,叫做十进制快速幂。牛客多校dreamoon出过一次。原理也很简单了,思考一下快速幂其实就是把这个数搞成二进制,然后按位考虑复杂度就是这个数二进制位数的个数,那么对于一个长度1e6的十进制数,转成二进制不会超过3e6,所以我们做类似的操...原创 2019-10-07 16:37:24 · 374 阅读 · 0 评论 -
集合与代数系统
集合代数集合(1)二元关系(2)关系的性质(3)等价关系和等价类(4)偏序关系代数系统(1)群(2)环(3)格开篇废话。首先,我不是从一个专业的数学研究者的角度来谈论这个话题,只是从一个学生的角度来写下自己的看法,所以文中会有很多的不严谨的地方,多多包涵,也欢迎指出。集合集合的基本运算:交并补就不说了,提几个不常见的:集合的差:两个集合的差是一个新的集合,这个集合中的元素是在A中出现,但是...原创 2019-08-16 22:17:37 · 2215 阅读 · 0 评论 -
【BZOJ 2142 】礼物
1.题目链接。首先物品总量大于等于需求就可满足条件,否则才是Impossible。然后就是简单的推理了,高中组合数问题: 首先给第一个人分配,C(n,w1),第二个人:C(n-w1,w2)....第i个人:C(n-sum(wi),wi+1).根据分步乘法,答案就是这些数乘起来。...原创 2019-08-16 14:50:53 · 235 阅读 · 0 评论 -
【2019 牛客多校第九场 Quadratic equation 】
1.题面: 题意十分解明了。仔细观察这个式子,可以发现这个式子就是韦达定理。我们构造一个模意义下的二次函数,然后在这个模意义下解一下这个方程就好了。其实这就涉及到模意义下的一元二次方程怎么解的问题?和一般的解法是一样的,有解和无解的判断条件就是看一下delt是不是存在...原创 2019-08-15 20:36:49 · 219 阅读 · 0 评论 -
【HDU 6209】The Intersection
1.题目链接。给定一个k,求出一个分母小于100000的分数,最接近pow(k,2.0/3).2.这个题看起来像高精度?但是这题卡long double。这个题目其实考到了一些比较有趣的东西,SB树和分数二分。首先把整数部分求出来,然后小数部分采用分数二分求出一个满足条件最接近的解即可。#include<bits/stdc++.h>#define ld long doubl...原创 2019-07-28 16:23:32 · 341 阅读 · 0 评论 -
【BZOJ 4161】
1.题目链接。矩阵快速幂的多项式取模优化模板题。#include<bits/stdc++.h>using namespace std;const int maxn = 4000, mod = 1000000007;int a[maxn + 5], p[maxn + 5], ans[maxn + 5], num[maxn + 5];int h[maxn + 5], tmp[...原创 2019-07-21 11:12:50 · 268 阅读 · 0 评论 -
【luogu 4570】元素
1.题目链接。题意:给定一些物品,每种物品有两个属性,序号和魔力值,现在在这些物品中选出几个,要求选中的物品序号异或起来不能是0并且魔力值加起来最大,求出这个最大值。2.线性基的性质之一,一组元素的线性基不是唯一的,但是数量是固定的,所以按照魔力值排序,贪心的选择,如果当前元素能够被插入线性基,就对答案有贡献,否则不行。#include<bits/stdc++.h>usin...原创 2019-07-19 13:11:50 · 100 阅读 · 0 评论 -
【HDU 4305】Lightning
1.题目链接。题目大意:给定平面上n个点和一个距离r,任意两个点之间有边的条件是:两个点的欧几里得距离小于等于r并且这两个点之间的连线不再有其他的点。然后问有多少种方法使得所有的点都联通。2.显然,这个图的每个生成树都是满足条件的,所以这个题其实就是在问(假设已经连边连好了),给定的这个图有多少棵生成树??关于生成树的计数问题,似乎有个Cayley定理,他说一个n个点的完全图有n^(n-2)棵...原创 2019-07-10 16:33:11 · 166 阅读 · 0 评论 -
【hihcoder 1480】矩阵填数
1.题目链接。小Hi在玩一个游戏,他需要把1, 2, 3, ... NM填入一个N行M列的矩阵中,使得矩阵每一行从左到右、每一列从上到下都是递增的。 2.这是杨氏矩阵的一个性质,其实这个题目的本质就是在求1-nm这nm个数字,可以构造出多少个n*m杨氏矩阵。首先对于杨氏矩阵,n个数,可以构成的杨氏矩阵的数量,在不加形状限制的情况下:可以从递推公式得到。F[i]代表i个数构成的杨氏矩阵...原创 2019-07-14 15:49:08 · 360 阅读 · 0 评论 -
【HDU 5950】Recursive sequence
1.题目链接。题意:a[n]=a[n-1]+2*a[n-2]+pow(n,4).求an。2.十分裸的矩阵快速幂,构造个7*7的矩阵就完事了,在对n^4狗杂矩阵时其实很简单,就按照他们的二项展开式的系数填就行了。#include<bits/stdc++.h>using namespace std;const int N = 1e6 + 100;#define ll lo...原创 2019-07-14 13:48:58 · 156 阅读 · 0 评论 -
【POJ 2154】Color
1.题目链接。给一个正边形和n中颜色,用这n种颜色给这个正n边形染色,有多少种方案,这里的正n边形是可以旋转的,旋转得到为同一种方案被视为一种。2.十分裸的Polya问题。 Q是置换的个数,lamada(q)是每种置换下,可以被表示为循环的个数。m是m种颜色。直接对应到这个题就是:...原创 2019-07-10 12:40:06 · 210 阅读 · 0 评论 -
【HDU 5942】Just a Math Problem
1.题目链接。题目意思十分的明显,明摆着就是推导公式然后求值,实际上是一个很棒的数论题。2.首先我们可以看出来,f(i)代表能够整除i的素因子的种类。而g(i)=pow(2,f(i)).为什么是2的f(i)次方,这其实是这个题目的突破口,假设该数有k个因子: 其实就是把这k个素因子分为两个互不相交的集合,也就是k=P*...原创 2019-07-08 16:46:47 · 285 阅读 · 0 评论 -
【HUST网赛 B】Balls
1.题目链接。距离比赛结束还有1.5H开始看这道题目,当时已经分析出来了是一个组合问题,但是思路错了,最后找到一个可行的方法,验证发现需要整数分解,就gg,没有再仔细地思考下去。2.分析:其实这个题目就是在说这样一个事情: 这个式子就是在说构造一个这样的数组:长度为k,每个元素都大于等于1并且所有的元素...原创 2019-04-15 16:10:01 · 198 阅读 · 0 评论 -
【数论-一些关于gcd的公式的推导】
1.数论里面与gcd有关的考察点很多很多。下面就几个常用的式子来推导一下。(1) 证明: (2) 证明: (...原创 2019-04-02 21:28:31 · 952 阅读 · 0 评论 -
【HDU1714】Simpson积分的简单应用
1.题目链接。问题就是在求红色部分的面积。我们可以把面积分为两个部分,交点左边的是三角形,交点右边的是椭圆。我们把交点求出来之后,左边直接三角形的面积公式,右边Simpson积分。完美解决。但是这里注意一下坐标的范围,就是(x,y)的取值范围,因为x>=a的,所以斜率一定存在,并且很重要的是y有可能是负数,所以在计算三角形的面积的时候需要加上一个绝对值。代码如下:#include<...原创 2019-02-24 15:23:51 · 299 阅读 · 0 评论 -
【HDU1568】fib数列公式的使用
1.题目链接这个题目的意思就是输出fib数列每一项的前四位数字。这个题应该没有什么好的做法,只能用fib数列的通项公式了。2.fib的通项公式: 至于这个公式是怎么来的呢?三言两语也说不清。可能有很多别的推导办法吧,但是我所知道的时使用数列递推的特征方程,本质上就是矩阵的特征方程和特征根来解...原创 2019-02-18 21:07:24 · 616 阅读 · 0 评论 -
【HDU5017】模拟退火求距离最值
1.题目链接。题目大意就是给定一个椭球,在这个椭球面上找一个点使得这个点与原点距离最近。2.分析:emmmm,这个题怎么说嘞,使用模拟退火当然是很好写啦,但是我在考虑使用高数中的拉格朗日乘数法,也是可以求多元函数的最值的,当然了,这个计算量是很大的,求完偏导之后观察了一下就放弃了,感觉自己变懒了,计算量一大就不想动手了。好吧,懒也有懒得写法,那就直接模拟退火吧,注意系数设置要大于0.95,我设...原创 2019-02-16 23:15:29 · 253 阅读 · 0 评论 -
数论-增大GCD
1.题目链接。这个题目的意思是让你删去最少的数字,使得剩下的数据的CGD变大。输出删除的数字的个数2.分析:我们知道,要使GCD变大就是让这些数据的删除之后保留尽可能多的公因子,我们首先预处理数据,求出所有数据的gcd,然后把所有的数据除以gcd,这样得到的数据就是互质的,根据唯一分解定理,我们对每个数据进行质因数分解,统计每个素因子被访问的次数,其中那个访问最大的的次数就是我们需要保留的数量...原创 2019-01-30 21:49:27 · 247 阅读 · 0 评论 -
数论-唯一分解定理
1.唯一分解定理。这是一个很神奇的定理,具体是这样描述的:任何一个正整数可以被分解成若干个素数的乘积。数学符号的表示就是:a=(p1^k1)*(p2^k2)*....(pm^km),其中p1,p2,..pm是素数,k1,k2...km是正整数。需要说明的是,这种分解的形式的唯一的。2.唯一分解的实现:前边我们提到了如何再线性时间内进行素数打表,这里就用处很大了。我们可以用处理出来的素数来分解这...原创 2019-01-30 20:21:47 · 1838 阅读 · 0 评论 -
【HDU1411】四面体的体积公式
1.题目链接。题目大意:就是给出一个四面体的六条边,求出这个四面体的体积。2.这个,如果知道坐标是很好解决的,假设我们知道的是坐标: 体积就是混合积的六分之一。(什么?x,y,z是啥?不懂这个的可以去看一下高等数学:空间解析几何。x,y,z是其中任意一个点和其他三个点连接的向量的坐标...原创 2019-02-21 20:57:41 · 11474 阅读 · 0 评论 -
【HDU1724】自适应Simpson积分
1.题目链接。题目的意思很简单,就是求阴影部分的面积,图形是一个椭圆。2.思路:这个做法当然有很多种了,最暴力的就是直接积分了,用一下三角换元也是很简单的就搞定了。其实化简到最后可能涉及到反三角函数的化简,这个就量力而行了呗,有能力继续的,就可以继续(其实没必要,也不是做数学题),搞不明白反三角怎么化简得,直接就写公式上去。第二种方法就是Simpson积分法了,自适应得Simpson积分很简单...原创 2019-02-15 22:55:32 · 246 阅读 · 0 评论