博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
nyoj-小明的密钥(362)--数论
阅读量:4286 次
发布时间:2019-05-27

本文共 1073 字,大约阅读时间需要 3 分钟。

题目链接:

题意:求N的因子的约数的个数的立方和,

可以把A分解素因子A = p1^a1 * p2^a2....pk^ak

  所以A 的因子的个数为(1+a1)*(1+a2)....*(1+ak)

1~p1^a1 的因子的约数的个数的立方和为s1=1^3 + 2^3 +......+(1+a1)^3 = (a1+1)^2*(a1+2)^2/4

.

.

1~pk^ak 的因子的约数的个数的立方和为sk=1^3 + 2^3 +......+(1+ak)^3 = (ak+1)^2*(ak+2)^2/4

所以 A的因子的约数的个数的立方和为s1*s2*....*sk;

如果A = 72 = 2^3 * 3^2

A的因子个数为(1+3)(1+2)=12    有1 2 3 4 6 8 9 12 18 24 36 72

#include
#include
typedef long long LL;int main(){ LL a,b,p = 1; while(scanf("%lld%lld",&a,&b)!=EOF) { LL k,ret,sum = 1 ; for(LL i = 2 ; i <= sqrt(a);i ++) { k = 0 ; if(a % i == 0) { while(a % i == 0) { k ++ ; a /= i; } ret = ((k*b + 1)*(k*b + 2) / 2) % 10007 ; ret *= ret ; sum = (sum*ret) % 10007; } } if(a > 1) { ret = ((1+b)*(b+2) / 2) % 10007 ; ret *= ret ; sum = (sum*ret) % 10007; } printf("Case %lld: %lld\n",p++,sum); }}

转载地址:http://zysgi.baihongyu.com/

你可能感兴趣的文章
Android系统组件之Activity的生命周期
查看>>
JavaScript利用URL向后台传入中文参数乱码问题解决之道!
查看>>
网页布局框架之Iframe
查看>>
Tomcat发布多个项目抛出异常
查看>>
Hibernate坑爹之实体类与数据库字段
查看>>
SpringMvc前端提交多个对象,后台接收
查看>>
关于Request
查看>>
AJAX——核心XMLHttpRequest对象
查看>>
Java处理js输入特殊字符(如“+、@、¥”)
查看>>
Mysql注释
查看>>
MySQL(root用户)密码重置
查看>>
grant授权
查看>>
MySQL创建用户与授权方法
查看>>
MySql数据类型
查看>>
MySql简单sql使用
查看>>
"未能加载文件或程序集“MySql.Data, Version=6.9.3.0”或它的某一个依赖项。
查看>>
CodeFirst for MySql
查看>>
Code Frist for Mysql 实例
查看>>
Visual Studio 开源控件扩展 NuGet 常用命令及常用组件
查看>>
mysql局域网访问设置
查看>>