![](https://img-blog.csdnimg.cn/20201014180756724.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数学
programmy
本人擅长Ai、Fw、Fl、Br、Ae、Pr、Id、PS等软件的安装与卸载,精通CSS、JavaScript、PHP、C、C++、C#、java、Ruby、Perl、Lisp、python、Objective-C、ActionScript等单词的拼写,熟悉windows、Linux、Mac、Android、IOS等系统的开关机。欢迎各大公司HR与本人联系。
展开
-
组合取模问题与Lucas定理
组合取模问题是求CmnC_n^m mod p的问题 CmnC_n^mmod p可以写成n!(n−m)!m!(modp){\frac{n!}{(n-m)!m!}}\pmod p这样的形式。一般做法结合逆元和费马小定理(ap−1≡1(modp){a^{p-1}}\equiv 1\pmod{p}),有: =>a∗ap−2≡1(modp)a*a^{p-2}\equiv 1\pmod p =>原创 2016-11-02 21:55:14 · 389 阅读 · 0 评论 -
【专题】线性同余方程
定义:a,b是整数,形如ax≡b(mod m)ax\equiv b(mod \ m),且x是未知整数的同余式称为一元线性同余方程。定理:a,b,m是整数且m>0,gcd(a,m)=dgcd(a,m)=d,如果d|b,则方程恰好有d个模m不同余的解,否则方程无解。由同余方程的定义式可得ax+my=bax+my=b,这个方程称为二元一次不定方程。解一元线性同余方程设d=gcd(a,m),由定理可知若不原创 2017-02-11 22:07:47 · 2245 阅读 · 0 评论 -
POJ 2115 H - C Looooops(一元线性同余方程)
POJ 2115题目大意给你四个数A B C kA \ B\ C \ k 问你for(int i=A;i!=B;i+=C)for(int \ i=A;i!=B;i+=C)语句循环的次数,死循环就输出FOREVER,所有的数字都是k−bit unsigned integerk-bit\ unsigned\ integer,mod 2kmod \ 2^k (也就是0≤A,B,C<2k,k≤320\l原创 2017-02-12 10:27:51 · 339 阅读 · 0 评论 -
POJ 3273 Monthly Expense(二分答案)
POJ 3272题目大意把N个数(aia_i)分成M份(1≤ai≤10000,1≤N≤100000,1≤M≤N)(1\le a_i\le 10000,1\le N \le 100000,1\le M\le N),每份得到一个和,问和的最大值的最小值是多少分析动态规划的方法是O(N3)O(N^3)的复杂度肯定超时。根据经验,问最大值中的最小值问题有一部分和二分有关,,对于这道题,我们可以对答案进行二分原创 2017-02-12 12:05:04 · 622 阅读 · 0 评论 -
POJ 3258 River Hopscotch(二分+贪心)
POJ 3258题目大意一条线段两个端点之间的距离是L,两端点之间分布着N个点,这N个点把线段分成了N+1份,现在让你最多去掉(第一次读错题想了很久不知道怎么做,remove是去掉不是移动,lll¬ω¬)M个点,问N+1份线段最小值的最大值是多少(1≤L≤109,0≤M≤N≤50000)(1\le L \le 10^9,0 \le M\le N\le 50000)。分析类似POJ 3273,也是用二原创 2017-02-12 16:27:41 · 380 阅读 · 0 评论 -
POJ 3122 Pie(二分水题)
POJ 3122题目大意有N个半径不同的饼,有F个人分,每个人分到的饼只能来自同一个饼并且每个人饼的大小必须相同,问每个人最大能得到多大的饼。分析二分水题代码#include<cstdio>#include<iostream>#include<cmath>#include<cstring>#include<cstdlib>#include<queue>#include<iomanip>原创 2017-02-13 11:03:10 · 499 阅读 · 0 评论 -
POJ 1905 Expanding Rods(推公式+二分求方程解析解)
POJ 1905题目大意已知圆的弦长和弧长求半径分析这道题直接解方程不好解,可以通过数值逼近的方法来求得方程的近似解,这道题可以用二分法来逼近答案。我一开始的做法是通过θ\theta和r之间的关系对θ\theta进行二分,WA,后面分析原因是可能是通过角度和边的关系来做精度达不到 ⎧⎩⎨⎪⎪⎪⎪sin(θ2)=L/2r2πr⋅θ2π=r∗2arcsin(L/2r)=L′L′=(1+n∗C)L\b原创 2017-02-12 22:02:31 · 299 阅读 · 0 评论 -
POJ 3096(map水题)
POJ 3096题目大意给你一个字符串,定义D-pair为在这个字符中间隔了D的两个字符组成的一对,如果这个字符串中对于相同的D,D-pair都是独特的话那么这个字符串就是surprising,让你判断一个字符串是否surprising。分析一个字符串最多长79,数据量不大,数据量大的话应该是要用hash的,为了方便就直接map来实现了,map更多是用来键值的查询,我的代码中通过count函数来判断原创 2017-02-13 12:13:01 · 405 阅读 · 0 评论 -
【数学建模】熵值法与多指标评价系统
一.熵的基本概念熵最初是热力学中的一个概念,后来延伸到信息领域。熵是对系统无序度的度量,系统越无序,熵越大。熵的统计学定义是: 某个宏观状态的微观状态数,取对数lnWlnW举个例子,抛10次硬币,正面朝上的次数是3是宏观状态,对应的微观状态数,也就是使得正面朝上的排列组合数为3。可以看出熵越大,该熵对应的状态出现的概率也越大,弄清楚定义之后也就不难理解为什么熵越大系统会越稳定了。二.熵值法在多原创 2017-01-26 22:22:19 · 33095 阅读 · 5 评论 -
DUTacm 1085 Water Problem(矩阵快速幂 找规律)
DUTacm 1085连接题目大意已知f(1),f(2)f(1),f(2)的值,且对于任意x>1x>1,有f(x+1)=f(x)+f(x−1)+sin(πx2)f(x+1)=f(x)+f(x-1)+sin(\frac{\pi x}{2}),求f(n)f(n)的值分析找规律之后可以发现这样的式子: f(n)=fib(n−2)∗f(1)+fib(n−1)∗f(2)+C(n)f(n)=fib(n-2)*原创 2017-03-09 10:20:54 · 780 阅读 · 0 评论 -
【专题】欧几里得算法、扩展欧几里得、乘法逆元
1.欧几里得用途 最大公因数和最小公倍数定理: gcd(a,b)=gcd(b,a%b)gcd(a,b)=gcd(b,a\%b) 证明: 我们令c=gcd(a,b)c=gcd(a,b) 令a=n∗ca=n*c , b=m∗cb=m*c a%b=a−k∗b=(n−m∗k)∗ca\%b=a-k*b=(n-m*k)*c,可知,c也是a%ba\%b的因子,现在只需证明c是b和原创 2016-10-31 23:29:02 · 2771 阅读 · 0 评论 -
【专题】因子和与因子个数
基本理论定义1:因子和函数 δ\delta定义为整数n的所有正因子之和,记δ(n)\delta(n)定义2:因子个数函数τ\tau 定义为正整数n的所有正因子个数,记为τ(n)\tau(n)定理1:如果ff是积性函数,那么ff的和函数F(n)=∑d|nf(d)F(n)=\sum_{d|n}f(d)也是积性函数推论:因子和函数δ\delta与因子个数函数τ\tau 是积性函数定理2:设p是一个素数,a原创 2017-02-08 21:45:04 · 720 阅读 · 0 评论 -
POJ 1845 Sumdiv(因子和 大数模 费马小定理)
POJ 1845题目大意 给你A和B(0≤A,B≤50000000)(0\le A,B\le 50000000),定义S为ABA^B的所有因子的和,让你输出S modmod 9901分析ABA^B可以写成(pa11pa22..pakk)B(p_1^{a_1}p_2^{a_2}..p_k^{a_k})^B因子的和S=∑(pa1b11pa2b22...pakbkk),1≤bi≤BS=\原创 2017-02-08 21:39:12 · 405 阅读 · 0 评论 -
hdu 5894 分位置(环上组合,16沈阳网络赛)
hdu5894题目大意 一个大小为 n的环,选 m 个位置涂黑,要求相邻两个黑点之间至少间隔 k个白点,问方案数。分析1. 思路一(较复杂) 令f(n,m)表示n个座位排成一排,选m个,相邻间隔不小于k个方案数。 将这环上n个座位编号从1到n. ①从1开始选的第一个>k :方案数为f(n-k,m) ②从1开始选的第一<=k :方案数为f(n-2k-1,m-1) 剩下原创 2016-11-03 10:04:21 · 526 阅读 · 0 评论 -
POJ 1722(gcd+思维)
POJ 1722链接题目大意 一次生日Party可能有p人或者q人参加,现准备有一个大蛋糕.问最少要将蛋糕切成多少块(每块大小不一定相等),才能使p人或者q人出席的任何一种情况,都能平均将蛋糕分食. 分析 我们可以这样思考:将一个蛋糕均匀切成p份可以满足p个人的情况,均匀切成q份可以满足q个人的情况。那现在这两种切法的刀印和在一起(也就是有最右边的图)一定可以满足这两种情况。问题就转化成了原创 2016-11-06 21:25:21 · 728 阅读 · 0 评论 -
概率论
我认为学习一门学科首要的是建立这门学科的知识框架,把握主干的知识脉络。基本概念1. 古典概型概率计算 古典概型中的事件A的概率为P(A)=A包含的基本事件数S中的基本事件数P(A)=\frac{A包含的基本事件数}{S中的基本事件数} 2. 条件概率 设A、B为两个事件,且P(B)>0,则A、B为两个事件,且P(B)>0,则 P(A|B)=P(AB)P(B)P(A|B)=\frac原创 2016-12-01 22:16:48 · 1728 阅读 · 0 评论 -
计算方法知识总结
计算方法author:AIDreamerblog:http://blog.csdn.net/mmy1996last modified on :2016/12/16一:线性方程组的直接解法基本方法 ⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪消去法{Gauss消去法:把原来的方程组化为等价的三角形方程组Gauss列主元消去法:在Gauss消去法的基础上的改进,用一原创 2016-12-23 18:36:05 · 9186 阅读 · 0 评论 -
POJ 1850 code(组合数学)
POJ 1850 code题目大意 将字母和单词(全部小写并且按照字典序递增)按照字典序编号,比如 a-1 b-2 … z-26 ab-27 … az-51 … 给你一个字母或单词,问你它的编号分析以树状的结构来分析这道题会比较直观。 按照树的节点依次编号,每个字母或单词就落在一个节点上。要求一个单词的编号,比较直接的想法是分别求出该原创 2017-02-05 22:21:33 · 368 阅读 · 0 评论 -
POJ 1942 Paths on a Grid(求组合数)
POJ 1942 Paths on a Grid题目大意 给定n,m(n,m≤232)(n,m \le 2^{32}) 求Cnn+mC_{n+m}^{n}分析裸的求组合数的题原公式是:Cmn=AmnAmm=n!m!(n−m)!=n∗(n−1)...∗(n−m+1)m!=n∗(n−1)∗...∗(m+1)(n−m)!C_n^m=\frac{A_n^m}{A_m^m}=\frac{n!}{m!(n原创 2017-02-06 20:30:45 · 323 阅读 · 0 评论 -
POJ 2635 The Embarrassed Cryptographer(千进制 高精度取模)
POJ 2635题目大意 给定一个K和L,K是由两个素数组成,如果这两个素数存在比L小(严格小于,因为这个原因WA了两次)的则输出最小的素数,都大于L输出“GOOD”(4≤K≤10100,2≤L≤106)(4\le K\le 10^{100},2\le L\le 10^{6})分析L比较小直接枚举L中的素数p,可以打表,然后高精度取模判断K是否是p的整数倍 大数求模可以用千进原创 2017-02-07 22:13:36 · 642 阅读 · 0 评论 -
POJ 3292 Semi-prime H-numbers(筛 打表)
POJ 3292题目大意定义H-number是满足模4余1的数系 H-number分为1,H-prime和H-composites. 能写成两个H-number的乘积的数称为H-composites.其余为H-prime 若一个H-number可以写成两个H-prime的乘积则称为H-semi-primes. 给你一个H-number数h(h \le 1,000,001),问你在1到h之间(原创 2017-02-07 22:16:39 · 377 阅读 · 0 评论 -
【专题】欧拉函数
欧拉函数欧拉函数定义:欧拉函数ϕ(n)\phi(n)指不超过n且与n互素的正整数的个数定义及定理定义1(算数函数): 定义在所有正整数上的函数称为算数函数定义2(积性函数):算数函数ff如果满足对任意两个互素的正整数**n和m,均有f(mn)=f(m)f(n)f(mn)=f(m)f(n),就称为积性函数(或乘性函数).如果对任意的两个正整数**n和m,均有f(mn)=f(m)f(n)f(mn)=f原创 2017-02-08 12:29:33 · 532 阅读 · 0 评论 -
HDU 6143 Killer Names(组合数学 递推 17多校第八场)
题目大意 有m种字母,每种字母可以随便用,让你用这m种字母组成两个长为n的字符串,满足这两个字符串中不含有相同字母的方案数有多少种。分析先分类:组成第一个字符串的数目为ii,对ii进行分类 不难得出, ans=∑i=1min(n,m−1)Cim⋅f(n,i)⋅g(n,m−i)f(n,i)表示从由i种字母组成长为n的方案数,i种都要用到g(n,i)表示从由i种字母组成长为n的方案数,i种不一定原创 2017-08-20 11:42:13 · 356 阅读 · 0 评论