各种数学
文章平均质量分 71
iYUNDI
iYUNDI
展开
-
[hoj 1076]Ordered Fractions[Farey序列]
生成分母不大于N的Farey序列. 这是递归求法. 用扩展欧几里得定理可以求紧邻着的下一项, 公式 K1*L2 - K2*L1 = 1 但是目前还没有学习它... #include #include #include #include using namespace std; int n; typedef struct node { int up,down原创 2013-08-11 17:56:20 · 913 阅读 · 0 评论 -
[hoj 2576 2577]Simple Computing & II
基本思想都是用容斥原理。 hoj 2576 给出一组数x1...xn,问从1到m中能有多少个数能够整除这组数中的至少一个数。 hoj 2577 给出一组数x1...xn,问从1到m中能有多少个数能够整除这组数中的唯一的数。 2576 x1 + x2 + x3 - x1*x2 - x1*x3 - x2*x3 + x1*x2*x3; 即奇加偶减。 2677 x1 + x2 + x3转载 2013-08-11 19:48:23 · 1108 阅读 · 0 评论 -
[poj 3090]Visible Lattice Point[欧拉函数]
找出N*N范围内可见格点的个数. 只考虑下半三角形区域,可以从可见格点的生成过程发现如下规律: 若横纵坐标c,r均从0开始标号,则 (c,r)为可见格点 <=>r与c互质 证明: 若r与c有公因子1<b<min(r,c),则(c/b, r/b)在线段(0, 0)(c, r)上,则(c, r)不是可见格点.(充分性) 若r与c互质,显然线段上不存在整点,则(c, r)不是可见格点.(必要性)原创 2013-08-12 11:57:43 · 1302 阅读 · 0 评论 -
[poj 2115]C Looooops[扩展欧几里德][模线性方程]
题意: for (variable = A; variable != B; variable += C) statement; 给出A,B,C和k(k表示变量是在k位机下的无符号整数),判断循环次数,不能终止输出"FOREVER". 思路: 需要求解 (A + x * C) % mod = B 变形之后即 C * x + mod * y = B -原创 2013-08-13 12:09:45 · 1158 阅读 · 0 评论 -
Miller_Rabin素数测试[Fermat小定理][二次探测定理][同余式][Wilson定理]
Miller和Rabin两个人的工作让Fermat素性测试迈出了革命性的一步,建立了Miller-Rabin素性测试算法。新的测试基于下面的定理: 如果p是素数,x是小于p的正整数,且,那么要么x=1,要么x=p-1。 这是显然的,因为相当于p能整除,也即p能整除(x+1)(x-1)。 由于p是素数,那么只可能是x-1能被p整除(此时x=1) 或 x+1能被p整除(此时x=p-1)。原创 2013-08-14 19:06:09 · 4407 阅读 · 1 评论 -
[HDU 4112]Break the Chocolate[YY推公式]
题意: 一块N*M*K的巧克力, 徒手掰一次掰一层, 用刀切可重叠. 问得到1*1*1的块最少的次数. 思路: 推公式... 首先考虑徒手掰: 对于一维的, n-1. 二维的, (m-1) + m*(n-1) = m*n - 1. 三维的, (k - 1) + k*(m*n - 1) = k*m*n - 1. 然后是用刀切: 对于一维的, lg n. 对于二维的, lg n原创 2013-10-01 14:30:59 · 934 阅读 · 0 评论