![](https://img-blog.csdnimg.cn/20201014180756927.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数论
ljpadamaa
这个作者很懒,什么都没留下…
展开
-
poj 2689计算素数间距离
利用两次素数筛选,第一次筛选 1--500000,第二次利用第一次的结果筛选出l--r间的素数,注意如果l是1,,则筛选不出,应单独判断。最后穷举最大最小距离 #include #include using namespace std; const int N=50001; bool prime[N+1]; long long p[N+1],pri[1000005]; int w; void g原创 2013-08-08 10:15:04 · 255 阅读 · 0 评论 -
hdu 4686 构造矩阵乘法 待更新
这题给我坑坏了 构造矩阵乘法,时刻要注意%mod 还有n==0时输出0 #include #include #include using namespace std; const long long maxn=10; long long n,m,k; const long long mod=1000000007; typedef struct { long long m[maxn][原创 2013-08-20 20:15:22 · 280 阅读 · 0 评论 -
poj 2649 素因子分解 ------ n!的素因子
n!的素因子分解中的素数p的个数为 n/p+n/p^2+n/p^3+n/p^4+...... 原理:n/p 计算出了1~n中含有素因子p的个数 n/p^2计算出了1~n中含有因子p^2的个数 ......... 求和就是n!中素数p的个数 注意:n或m为0时,都不能整除原创 2013-08-09 14:43:30 · 331 阅读 · 0 评论 -
poj 3233 Matrix Power Series 矩阵快速幂 二分
计算 s(k)= a+a^2+a^3+.....a^k a^k 用矩阵的快速幂 因为k比较大,不能直接循环相加,因吃应用二分的方法 s(k)=s(k/2)+s(k/2)*a^k/2 (k为奇,偶略有不同) #include #include using namespace std; const int maxn=40; int n,m,k; typedef struct { i原创 2013-08-14 15:12:54 · 266 阅读 · 0 评论