数学定理
Marcus-Bao
这个作者很懒,什么都没留下…
展开
-
整数因子的分解
感谢! 一、相关概念 素数:又叫质数,就是只能被1和自己整除的整数。例如:2,3,5,7…….合数:与素数相对,即除了能被1和自己整除还能被其他整数整除的整数。例如:4,6,8,9……注意:1就不是素数,也不是合数。定理:任意一个正整数都能分解成若干个素数乘积的形式即:证明:存在性:唯一性:二、用c#实现整数的素因子转载 2017-03-24 00:41:31 · 2033 阅读 · 0 评论 -
hdu 1573 X问题 CRT不互质解法
Problem Description求在小于等于N的正整数中有多少个X满足:X mod a[0] = b[0], X mod a[1] = b[1], X mod a[2] = b[2], …, X mod a[i] = b[i], … (0 Input输入数据的第一行为一个正整数T,表示有T组测试数据。每组测试数据的第一行为两个正整数N,M (0 Outpu原创 2017-10-10 19:27:48 · 766 阅读 · 0 评论 -
HDU - 4596 Yet another end of the world exgcd性质
题目链接 题意:给定若干组x,y,z,问是否能找到一个整数d,使得至少两组x,y,z满足y≤d%x≤z。 思路:首先可以想到直接找两组的就可以了. 假设该整数为W, 有: Wmod(x1)=c1 W mod (x_1) = c_1 Wmod(x2)=c2 W mod (x_2) = c_2 可以得到 W=k1∗x1+c1 W = k_1*x_1 + c_1原创 2017-10-02 17:10:08 · 323 阅读 · 0 评论 -
HYSBZ - 1951 古代猪文 数论全集
题意:题目描述 猪王国的文明源远流长,博大精深。 iPig在大肥猪学校图书馆中查阅资料,得知远古时期猪文文字总个数为N。当然,一种语言如果字数很多,字典也相应会很大。当时的猪王国国王考虑到如果修一本字典,规模有可能远远超过康熙字典,花费的猪力、物力将难以估量。故考虑再三没有进行这一项劳猪伤财之举。当然,猪王国的文字后来随着历史变迁逐渐进行了简化,去掉了一些不常用的字。 iPig打算研究古时某个朝代的原创 2017-10-09 23:44:28 · 475 阅读 · 0 评论 -
CodeForces - 724C Ray Tracing 扩展欧几里得+思维
点击打开链接题意:n*m的矩形内有k个点,四周有围墙围起来。从(0,0)45度发射小球,速度为2√每次遇到墙都正常反弹,直到射到顶点被吸收。问每个点第一次被经过的时刻。思路:我们试着将这个反射的过程看成穿过,那么需要将整个矩阵展开(即变成一条直线).即所有穿过的都是关于 x==2*k*n,或者y == 2*kk*m进行对称的,根据轴对称性计算坐标进而可以原创 2017-09-30 22:31:10 · 337 阅读 · 0 评论 -
hdu 5446 Unknown Treasure lucas+CRT || 扩展lucas
题目链接题意:求C(n,m)%(∏pi)。pi小于105,m,n,以及答案都是1018。思路:裸地lucas定理,分别对每个素数取模然后CRT合成就行。我用的扩展lucas过得..#include using namespace std;typedef long long ll;const int N = 1e5;ll prim[N原创 2017-09-20 17:04:24 · 321 阅读 · 0 评论 -
hdu 1370 中国剩余定理+exgcd
题目链接题意:一个人有体力,感情,智商三个周期,周期分别为23天,28天,33天,然后告诉你某一次三者达到顶峰的是哪一天,问下一次三者达到顶峰的日子。思路:暴力可过,但是还是想用中国剩余定理来解,中国剩余定理:其中Mi为Ni 在%ni 下的逆元. PS : 最基础的中国剩余定理,也就是上面我们提到的,要求%原创 2017-08-27 14:56:10 · 329 阅读 · 0 评论 -
HDU - 5514 Frogs 欧拉函数||容斥定理
题目链接 题意:有 n 个青蛙,第 i 个青蛙每次只能够跳 ai 步,现在有 m 个石头围成一圈,编号为 0 到 m−1,现在青蛙可以围着这个石头组成的圆跳无限次,每跳一次就会占领这个石头,可以无限占领,现在问你的是这 n 个青蛙占领的石头的编号的总和是多少。思路: 先说第一种方法: 我们可以发现对于每个ai,他所能经过的石头为 k*gcd(m,ai).但是我们发现比如第一个样例 2 1原创 2017-09-28 00:14:22 · 503 阅读 · 0 评论 -
HDU - 5768 Lucky7 CRT +容斥
题目链接题意:对于给出的区间[x, y]找出有多少个符合要求的数:能被7整除.给出不超过15组(pi, ai),其中pi为质数;要求找出的数x满足 x % pi != ai;思路: 按照数位dp那一顿怼啊,后来发现好像不大对啊。我们知道如果给你很多(pi,ai)让你找出最小的x%pi = ai。也就是解同余方程,我们肯定就知道是CRT原创 2017-09-01 23:55:33 · 357 阅读 · 0 评论 -
2017 ACM-ICPC 亚洲区(西安赛区)网络赛 B. Coin(构造二项式定理 or dp矩阵优化)
题目链接 题意:一枚硬币,投一次正面朝上概率为q/p, 求投k次,正面朝上的次数为偶数次的概率X/Y,输出(X∗Y−1)mod(1e9+7)(X∗Y^-1)mod(1e9+7) 思路: 观察发现所有可能的次数,是一个二项式定理展开式, 所以我们构造一个 ((p/q)x+(p−q)/p)k((p/q)x+(p-q)/p)^k 可以发现我们首先将x赋值为1,那么全部为正的,然后再将x赋原创 2017-09-16 20:42:06 · 622 阅读 · 0 评论 -
2017 ACM-ICPC 亚洲区(西安赛区)网络赛 F. Trig Function (切比雪夫多项式)
题意:思路:手算了几项,然后把n为6的那一行输到OEIS,就搜到了切比雪夫多项式。切比雪夫多项式前几项为:The triangle a(n,m) begins:n\m 0 1 2 3 4 5 6 7 8 9 100:原创 2017-09-16 20:07:49 · 539 阅读 · 0 评论 -
hdu4704 && FZU1759 && HYSBZ3884 &&POJ3090 欧拉函数 || 欧拉降幂 ||费马小定理
hdu 4704思路: 题目本身很简单,就是讲N划分为K个数的和有几种划分(1.....K)对应每一个K有 C(N-1,K-1)中划分法,一共就是 2^(N-1)。但是发现N最多为 10^100000 哪怕快速幂也会爆.方法1: 这也是我自己的一个做法吧,因为取模了,可以用欧拉降幂公式.PS : 1e9+7 是大素数,所以他的欧拉函数值为1e9原创 2017-08-27 17:23:17 · 619 阅读 · 0 评论 -
HDU - 5728 PowMod 欧拉函数(降幂)+数学推倒
题目链接题意:k=∑mi=1φ(i∗n) mod 1000000007k=∑i=1mφ(i∗n) mod 1000000007 nn is a square-free number. φφ is the Euler's totient function. find: ans=kkkk...k mod pans=kkkk...k mod p T原创 2017-08-28 23:55:58 · 767 阅读 · 0 评论 -
POJ 1061 青蛙的约会 扩展欧几里得
题目链接思路:这个题目可以转化为: x + km = y + kn (% L)k(m-n)= y - x (% L)k(m-n) + t*L = y -x 进而转化到扩展欧几里得.ll exgcd(ll a,ll b,ll &x,ll &y) { if(b == 0) { x = 1;原创 2017-08-27 14:53:19 · 376 阅读 · 0 评论 -
基础数学公式 长期更新
1.取模运算(a + b) % p = (a % p + b % p) % p (1)(a - b) % p = (a % p - b % p) % p (2)(a * b) % p = (a % p * b % p) % p (3)a ^ b % p = ((a % p)^b) % p (4)a/b%p=a%(b*p)/b%p;原创 2017-06-02 20:36:56 · 398 阅读 · 0 评论 -
数学定理及部分代码(长期更新)
1.卡特兰数前几项为:1, 1, 2, 5, 14, 42, 132,catalan数满足递推式[1] :h(n)= h(0)*h(n-1)+h(1)*h(n-2) + ... + h(n-1)*h(0) (n>=2)h(n)=h(n-1)*(4*n-2)/(n+1);h(n)=C(2n,n)/(n+1) (n=0,1,2,...)h(n)=c原创 2017-04-26 14:35:57 · 599 阅读 · 0 评论 -
默慈金数
感谢大佬~!今天,我来讲一种比较特殊的数,可能很多人都没有听过这种数,它叫默慈金数。但事实是它早就已经进入了ACM竞赛中了。好了,接下来让我们一起来认识它,并会讲述一些它的重要应用。 在百度百科上,是这样定义默慈金数的:一个给定的数的默慈金数是在一个圆上的个点间,画出彼此不相交弦的全部方法的总数。比如为4时,方法数为9,如下图转载 2017-04-26 14:38:55 · 586 阅读 · 0 评论 -
UVA - 11754 Code Feat CRT (big技巧枚举 small暴力)
题目链接 题意:给定一个c个x, y1,y2,y3..yk形式,前s小的答案满足s % x在集合y1, y2, y3 … yk中 (c <= 9 ,k <= 100,S<= 10) 思路:考虑搜索枚举每个集合的yi,那么复杂度就是所有k的乘积,所以这个方法当k很小的时候我们可以dfs枚举每个集合选一个yi的组合,然后CRT暴力求解输出S个正整数解。 然而当k乘积很大的时候我们怎么解决?考原创 2017-10-10 23:09:33 · 322 阅读 · 0 评论