![](https://img-blog.csdnimg.cn/20201014180756926.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论
lhfl911
这个作者很懒,什么都没留下…
展开
-
中国剩余定理(互质+非互质)
给定n组同余方程,求x。 x≡a1(modx\equiv a_1(mod m1) m_1) x≡a2(modx\equiv a_2(mod m2) m_2) x≡a3(modx\equiv a_3(mod m3) m_3) …… x≡an(modx\equiv a_n(mod mn) m_n)若mim_i两两互素,则x必定有解。令M=m1∗m2∗...∗mnM = m_1 * m_2 *原创 2016-04-25 17:36:00 · 1693 阅读 · 0 评论 -
hdu 5288 OO's sequence【因子分解】
题目大意: 给定长度为n(n<=10510^5)的整数序列(a<=10410^4),求出所有子区间中不被区间内其他数整除的数的个数和,结果取余109+710^9+7;思路: 一个个枚举区间肯定行不通,可以根据每个整数对结果的贡献求和; 对于第i个数aia_i,找出aia_i左右第一个被整除的边界l、r,则aia_i的贡献为(i-l)*(r-i); 利用标记数组可以在sqrt(a)内找出左右边原创 2016-05-11 20:16:51 · 307 阅读 · 0 评论 -
zoj 3903 Ant【推公式+逆元】
题目大意: t组数据,每组数据给一个n,表示长方体最长的一条边。 L表示从一顶点到其对顶点的表面最短路程。 计算所有可能的长方体的L2L^2和, 结果对1000000007取余。 t≤100,1≤n≤1014t\le100, 1≤n≤10^14注:(a, b, n)和(b, a, n)算重复。思路: 1≤a≤n1≤a≤n,a≤b≤na≤b≤n, sum += (a+b)2+n2(a+b)^原创 2016-04-17 16:36:45 · 416 阅读 · 0 评论 -
hdu 5673 robot【默慈金数or卡特兰数】
补充知识:卡特兰数 Catalan(n)=C2nnn+1Catalan(n) = \frac{C_n^{2n}}{n+1}Catalan(n)=4n−2n+1∗Catalan(n−1)Catalan(n) = \frac{4n-2}{n+1}*Catalan(n-1)默慈金数 Mn=∑n/2i=0C2in∗Catalan(i)M_n = \sum_{i=0}^{n/2}C_n^{2i}*C原创 2016-04-26 20:33:18 · 573 阅读 · 0 评论 -
hdu 5755 Gambler Bo【gauss】
思路: n*m个数变化,设置n*m个变量,每个位置对应5个变量发生改变,连立n*m个方程,高斯消元求解,注意取模。// hdu 5755#include <cstdio>#include <cstring>#include <algorithm>using namespace std;int a[950][950];int x[950];int gcd(int a, int b){原创 2016-07-29 13:23:05 · 343 阅读 · 0 评论 -
hdu2204 Eddy's 爱好 【容斥原理】
题目链接: hdu2204题意: 给一正整数n(1<=n<=1018)(1<= n <= 10^{18}), 求[1,n]中能表示成MK(K>1)M^K(K>1)的数的个数。第一次好好看容斥,感觉智商已经不够了,看了好多题解才明白。。。若MK<=nM^K <= n 满足, 则任意m <= M 都满足mK<=nm^K<= n ;满足MKM^K 的数中,若K不是素数,则MK=(Ma)pM^K =原创 2016-08-05 15:39:52 · 385 阅读 · 0 评论 -
hdu 5768 多校4 Lucky7【crt+容斥】
hdu5768题意: 求区间[l, r]中是7的倍数且不满足任一给定同余式的数的个数,(0<l<r<108)(0\lt l \lt r \lt 10^8)。分析: 满足不同同余式的集合不同且可能有重叠,用容斥原理解决,联立同余式求解用中国剩余定理,因为pi都是素数,可以用互质版求解**注意**crt中乘法会爆long long,写个快速乘。#include <cstdio>#include <原创 2016-07-30 16:05:33 · 584 阅读 · 0 评论 -
hdu 5667 Sequence【费马小定理+矩阵快速幂】
题意: Lcomyn 是个很厉害的选手,除了喜欢写17kb+的代码题,偶尔还会写数学题.他找到了一个数列: fn=⎧⎩⎨⎪⎪1,ab,abfcn−1fn−2,n=1n=2otherwisef_n=\left\{\begin{matrix}1 ,&n=1 \\ a^b,&n=2 \\ a^bf_{n-1}^cf_{n-2},&otherwise \end{matrix}\right.原创 2016-04-17 15:15:11 · 639 阅读 · 0 评论 -
hdu5902 GCD is Funny
传送门: hdu5902 题意: n个数每次选三个数删除,取其中两个数将gcd放回去两次,问最后剩的数可能是多少。一开始看是1001题,天真的以为两两gcd就是答案,被查的不忍直视。 假设某一次放回了x,x,下次操作选x,x,c,可以选择放回gcd(x, x)或者gcd(x,c)。所以答案为n的子集的gcd,因为第一次操作时只能舍弃一个数,因此答案不包括n的gcd,集合元素个数i, 2≤i<原创 2016-09-25 14:03:01 · 1138 阅读 · 0 评论