求p^X的因子之和公式

最新推荐文章于 2021-11-01 10:34:56 发布
最新推荐文章于 2021-11-01 10:34:56 发布
阅读量258 收藏
点赞数
分类专栏： 公式
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/b_o_n_z_t_f/article/details/105003266
版权

如果 p 是素数 则其因子只有1和它本身，S(p^X) =1+p+p^ 2+…+p^X = (p^(X+1)-1)/(p-1)

霍雨浩——舞麟
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
基于改进粒子群算法结合惩罚因子的配电网无功优化的研究.pdf
09-28
通过公式（7）求取的配电网无功优化结果可以避免粒子进化方向单一性以及各粒子离散量引起误差，从而提高优化结果。 本文的贡献在于设计了一种基于改进粒子群算法结合惩罚因子的配电网无功优化方法，该方法可以避免...
整理自然数前n项p次方求和公式
weixin_34272308的博客
11-11 883
推广到p次方： 转载于:https://www.cnblogs.com/baalhuo/p/5762119.html
http://acm.nyist.net/JudgeOnline/problem.php?pid=420
pursuit的专栏
03-22 615
p次方求和 时间限制：1000 ms | 内存限制：65535 KB 难度：3 描述 一个很简单的问题，求1^p+2^p+3^p+……+n^p的和。 输入 第一行单独一个数字t表示测试数据组数。接下来会有t行数字，每行包括两个数字n,p， 输入保证0 输出 输出1^p+2^p+3^p+……+n^p对10003取余的结果，每个结果单独占一行。 样例输入 210 110 2 样例
n的x次方求和
nainiu的博客
04-01 1万+
:先求解 = 在求解 =........ ...... 求 = 最后将这n个式子相加得： 求出结果为 :需要在上面的基础下才能得出结果，因为式子中会出现。 同理 求结果过程中用到的公式： 扩展： ...
因子和的计算方法
DaDaguai001的博客
08-20 1万+
转发自 因子和的计算方法 2015年06月28日 10:45:22 阅读数：2163 储备知识: 一个数约数个数=所有素因子的次数+1的乘积 一个2000000000以内的数字不会有超过12个素因子 较小的数的指数一定大于等于较大的数的指数 因子和的计算方法 神马叫因子和...
LDA模型里Gibbs sampling后验概率详细推导过程
11-17
和贝叶斯公式可以得到 P(zawt,z_dt,ax,β,V) P(W,tz,c,β,甲)P(w,t,zcB, P(wtzd,a,阝,乎P(w, t, z-dil,β,v 根据 Graphical Model,wda,ta都是由za生成的,如果不考虑zd则无法考虑wd,td。从而得到 P(zdw,t,z_di,x,βv)...
一个算术函数与最大素因子函数的混合均值 (2010年)
05-17
在Smarandache函数S（n）及因子积数列{Pd（n）}的基础上，构造并研究了∑n≤x（S（Pf（n））-1/2d（n）P（n））2的一种均值分布性质，利用初等方法和素数定理研究了混合均值问题，给出了它的一个较强的渐进公式。
自适应局部邻域特征点提取和匹配的点云配准.docx
02-23
三维点云配准作为激光三维扫描和成像的关键技术之一，得到了广泛的研究和应用。 二、点云配准技术 点云配准技术多采用粗匹配和精配准结合的思想。在精配准中，应用最为广泛的是 Besl 等[5]提出的迭代最近点算法...
自平衡小车自平衡小车的原理硬件设计和源代码.zip
08-25
最直观的想法就是f(x) =Kp*(L-x)，简单的说就是你在左边我就向右推，你在右边我就向左推，这就是比例因子P； 现在考虑两种情况，同样是在x位置，小球静止和小球具有速度V这两种情况。很明显，如果V>0，我们只需要...
知识点 - 因数之和 因数个数公式
qq_41848675的博客
08-03 7768
知识点 - 因数之和 因数个数公式 解决问题类型： 问有几个因数，因数之和， 或者问某些特定约数之和，比如不能被大于4的平方数整除的约数之和（即质因数的次数都为1） 结论 若对 nnn 质因数分解得到 p1e1⋅p2e2⋯pkekp_1^{e_1} \cdot p_2^{e_2} \cdots p_k^{e_k}p1e1​​⋅p2e2​​⋯pkek​​则有 因数个数公式： d(n)=(e1+1)⋅...
python3练习题--求完数（所有因子之和等于本身）
热门推荐
PythonStory--记录一段苦逼变牛逼的往事
11-14 3万+
题目：求指定区域内的完数（10000为例），及所有因子之和恰好等于本身，如：6=1+2+3. 思路：第一步，将所有因子追加到一个列表中，第二步，将符合条件的数字打出来。 代码： #!/usr/bin/python3 l = [ ] for n in range (1,10000): for a in range (1,n): if n%a ==0:
因子和与因子数
weixin_52314777的博客
11-01 677
因子和 看这里 力扣1390. 四因数 暴力 class Solution { public: int sumFourDivisors(vector<int>& nums) { int res=0; for(auto t:nums){ vector<int>cnt; int j=t; for(int i=1;i<=j/i;i++){
python怎么求因子之和_P066 求n的所有因子之和 ★★
weixin_39805644的博客
12-30 3159
所属年份：2011.3请编写一个函数fun，它的功能是：计算并输出给定整数n的所有因子(不包括1与自身)之和。规定n的值不大于1000。例如，若主函数从键盘给n输入的值为856，则输出为sum＝763。#include int fun(int n){}int main(){int n,sum;FILE *out;printf("Input n:");scanf("%d", &n);s...
因子和，因子数，1到n的因子和，1到n的因子数
Angel-Yan的博客
04-26 3199
1 - 求n的因子和 因子和函数σ定义为整数n的所有正因子之和，记为σ(n)  它是一个积性函数 首先对n进行因子分解 (因子分解代码附后) n = p1^a1 * p2^a2 * ~~~ * px ^ ax σ(n) =((p1^(a1+1)-1)/(p1-1) * ((p2^(a2+1)-1)/(p2-1) * .... * ((pj^(aj+1)-1)/(pj-1)) = Π(j=1 -&g...
关于Windows 9x的vmm32问题解决方法
最新发布
06-15
关于Windows 9x的vmm32问题解决方法
Linux下用sdcc开发51单片机，该示例是中断处理程序.zip
06-15
该资源内项目源码是个人的课程设计、毕业设计，代码都测试ok，都是运行成功后才上传资源，答辩评审平均分达到96分，放心下载使用！ ## 项目备注 1、该资源内项目代码都经过测试运行成功，功能ok的情况下才上传的，请放心下载使用！ 2、本项目适合计算机相关专业(如计科、人工智能、通信工程、自动化、电子信息等)的在校学生、老师或者企业员工下载学习，也适合小白学习进阶，当然也可作为毕设项目、课程设计、作业、项目初期立项演示等。 3、如果基础还行，也可在此代码基础上进行修改，以实现其他功能，也可用于毕设、课设、作业等。 下载后请首先打开README.md文件（如有），仅供学习参考, 切勿用于商业用途。 该资源内项目源码是个人的课程设计，代码都测试ok，都是运行成功后才上传资源，答辩评审平均分达到96分，放心下载使用！ ## 项目备注 1、该资源内项目代码都经过测试运行成功，功能ok的情况下才上传的，请放心下载使用！ 2、本项目适合计算机相关专业(如计科、人工智能、通信工程、自动化、电子信息等)的在校学生、老师或者企业员工下载学习，也适合小白学习进阶，当然也可作为毕设项目、课程设计、作业、项目初期立项演示等。 3、如果基础还行，也可在此代码基础上进行修改，以实现其他功能，也可用于毕设、课设、作业等。 下载后请首先打开README.md文件（如有），仅供学习参考, 切勿用于商业用途。
施工作业人员配备与人员资格及职责分工表.docx
06-15
施工作业人员配备与人员资格及职责分工表.docx
求子数组最小公倍数中质因子数数目之和
03-27
首先，我们需要知道一个结论：一个数的质因子数目等于其所有质因子的指数之和加1。比如，$24=2^3\times3$，其质因子数目为$3+1=4$。 接下来，我们可以使用暴力枚举子数组的方法，对每个子数组求出其最小公倍数，然后计算其质因子数目之和即可。时间复杂度为$O(n^3)$，其中$n$为数组长度，由于子数组数量为$O(n^2)$，对每个子数组求最小公倍数的时间复杂度为$O(n)$。 代码如下： ```python def gcd(a, b): if b == 0: return a return gcd(b, a % b) def lcm(a, b): return a * b // gcd(a, b) def count_factors(num): count = 0 for i in range(2, int(num**0.5)+1): while num % i == 0: count += 1 num //= i if num > 1: count += 1 return count def subarray_lcm_factor_sum(arr): n = len(arr) ans = 0 for i in range(n): for j in range(i, n): subarr_lcm = arr[i] for k in range(i+1, j+1): subarr_lcm = lcm(subarr_lcm, arr[k]) ans += count_factors(subarr_lcm) return ans ``` 然而，上述方法的时间复杂度较高，不能通过本题。考虑优化。 观察题目的数据范围，发现数组中每个元素都小于$10^6$，因此一个子数组的最小公倍数也不可能超过$10^6$，这启示我们可以对每个质数预处理出其在$[1,10^6]$范围内的倍数中出现的次数。具体地，我们可以使用线性筛法预处理出所有小于等于$10^6$的质数，然后对每个质数$p$，预处理出其在$[1,10^6]$范围内的倍数中出现的次数$cnt_p$，即： $$cnt_p=\sum_{i=1}^{\lfloor 10^6/p\rfloor}\lfloor\frac{i}{p}\rfloor$$ 这里的$\lfloor x\rfloor$表示不超过$x$的最大整数。例如，当$p=2$时，有$cnt_2=\lfloor 10^6/2\rfloor+\lfloor 10^6/4\rfloor+\lfloor 10^6/8\rfloor+\cdots=500000+250000+125000+\cdots$。这个式子可以用等比数列求和公式简化为$cnt_2=10^6-1$。 预处理完成后，对于每个子数组，我们可以依次枚举$[1,10^6]$范围内的质数$p$，统计其在子数组的最小公倍数中出现的次数。具体地，设当前枚举到的质数为$p$，统计$p^{cnt_p}$在子数组的最小公倍数中出现的次数。这个数量等于子数组中所有元素中$p$的指数的最小值乘以$cnt_p$。 具体地，我们可以用一个长度为$10^6+1$的数组$factor$来存储每个数对应的最小质因子。我们可以先用线性筛法求出所有小于等于$10^6$的质数，并将它们的倍数都标记上最小质因子。然后，对于每个子数组，我们可以使用一个长度为$10^6+1$的数组$cnt$来统计每个质数的指数之和，然后根据上述公式计算质因子数目之和。 总时间复杂度为$O(n\log\log 10^6)$，其中$n$为数组长度。 代码如下： ```python def subarray_lcm_factor_sum(arr): n = len(arr) # 预处理每个数的最小质因子 factor = [0] * (10**6+1) primes = [] for i in range(2, 10**6+1): if factor[i] == 0: factor[i] = i primes.append(i) for p in primes: if p > factor[i] or i*p > 10**6: break factor[i*p] = p # 预处理每个质数在[1,10^6]范围内的倍数中出现的次数 cnt = [0] * (10**6+1) for p in primes: cnt[p] = 10**6 // p for i in range(2, cnt[p]+1): cnt[p] -= cnt[i*p] # 枚举子数组，统计质因子数目之和 ans = 0 for i in range(n): for j in range(i, n): subarr_lcm = arr[i] for k in range(i+1, j+1): subarr_lcm = lcm(subarr_lcm, arr[k]) for p in primes: if p > subarr_lcm: break exp_sum = float('inf') for x in arr[i:j+1]: exp_sum = min(exp_sum, get_exponent(x, p)) ans += exp_sum * cnt[p] return ans ``` 其中，$get\_exponent(x,p)$用于计算$x$在质因数分解中$p$的指数，代码如下： ```python def get_exponent(num, prime): exp = 0 while num % prime == 0: exp += 1 num //= prime return exp ```

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值