【HDU - 1452】 Happy 2004(因子和,逆元,快速幂)
题干:
Happy 2004
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 1863 Accepted Submission(s): 1361
Problem Description
Consider a positive integer X,and let S be the sum of all positive integer divisors of 2004^X. Your job is to determine S modulo 29 (the rest of the division of S by 29).
Take X = 1 for an example. The positive integer divisors of 2004^1 are 1, 2, 3, 4, 6, 12, 167, 334, 501, 668, 1002 and 2004. Therefore S = 4704 and S modulo 29 is equal to 6.
Input
The input consists of several test cases. Each test case contains a line with the integer X (1 <= X <= 10000000).
A test case of X = 0 indicates the end of input, and should not be processed.
Output
For each test case, in a separate line, please output the result of S modulo 29.
Sample Input
1 10000 0
Sample Output
6 10
Source
ACM暑期集训队练习赛(六)
题意:求2004^x的所有因子和对29取余的结果。
题解:因为所有数都可以被分解成a=p1^c1*p2^c2*...*pk^ck,并有约数和定理:sum=(p1^0+...+p1^c1)*(p2^0+...+p2^c1)*...*(pk^0+...+pk^ck)。每一项可以通过等比数列求和得到。2004=2^2*3*167,所以答案就是2^(2*n+1)*3^(n+1)*167^(n+1)/(2*166)。由于需要取余,所以不能直接除要求出2*166的逆元改为相乘。
AC代码:
#include<bits/stdc++.h>
#define debug cout<<"aaa"<<endl
#define mem(a,b) memset(a,b,sizeof(a))
#define LL long long
#define lson l,mid,root<<1
#define rson mid+1,r,root<<1|1
#define MIN_INT (-2147483647-1)
#define MAX_INT 2147483647
#define MAX_LL 9223372036854775807i64
#define MIN_LL (-9223372036854775807i64-1)
using namespace std;const int N = 100000 + 5;
const int mod = 29;int e_gcd(int a,int b,int &x,int &y) {if(b==0) {x=1,y=0;return a;}int q=e_gcd(b,a%b,y,x);y-=a/b*x;return q;}int quick(int a,int b){int ans=1;while(b){if(b&1){ans=(ans*a)%mod;}b>>=1;a=(a*a)%mod;}return ans;
}int main(){int n,x,y,a,b,c,ans;exgcd(166*2,29,x,y);while(~scanf("%d",&n)&&n){a=(quick(2,2*n+1)-1)%mod;b=(quick(3,n+1)-1)%mod;c=(quick(167,n+1)-1)%mod;ans=((a*b*c*x)%mod+mod)%mod;printf("%d\n",ans);}return 0;
}
【HDU - 1452】 Happy 2004(因子和,逆元,快速幂)相关推荐
- 数学--数论--HDU 4675 GCD of Sequence(莫比乌斯反演+卢卡斯定理求组合数+乘法逆元+快速幂取模)
先放知识点: 莫比乌斯反演 卢卡斯定理求组合数 乘法逆元 快速幂取模 GCD of Sequence Alice is playing a game with Bob. Alice shows N i ...
- hdu 5451 Best Solver 矩阵循环群+矩阵快速幂
http://acm.hdu.edu.cn/showproblem.php?pid=5451 题意:给定x 求解 思路: 由斐波那契数列的两种表示方法, 之后可以转化为 线性表示 F[n] = ...
- 2014 Super Training #7 F Power of Fibonacci --数学+逆元+快速幂
原题:ZOJ 3774 http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3774 --------------------- ...
- (组合数求模=乘法逆元+快速幂) Problem Makes Problem
题目: As I am fond of making easier problems, I discovered a problem. Actually, the problem is 'how ca ...
- hdu 5411 CRB and Puzzle(矩阵快速幂)
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5411 解题思路: 题目大意: 给定n个点 常数m 下面n行第i行第一个数字表示i点的出边数,后面给出 ...
- 【HDU - 5187】zhx's contest (快速幂+ 快速乘,模板)
题干: 2018百度之星复赛晋级名单出炉(增加20%晋级名额)~ zhx's contest Time Limit: 2000/1000 MS (Java/Others) Memory Limi ...
- 【HDU - 5015 】233 Matrix (矩阵快速幂)
题干: In our daily life we often use 233 to express our feelings. Actually, we may say 2333, 23333, or ...
- 【BZOJ-4522】密钥破解 数论 + 模拟 ( Pollard_Rho分解 + Exgcd求逆元 + 快速幂 + 快速乘)...
4522: [Cqoi2016]密钥破解 Time Limit: 10 Sec Memory Limit: 512 MB Submit: 290 Solved: 148 [Submit][Stat ...
- HDU 1757 A Simple Math Problem (矩阵快速幂)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1757 在吴神的帮助下才明白如何构造矩阵,还是好弱啊. 此处盗一张图 1 #include <io ...
- HDU 3306 Another kind of Fibonacci 矩阵快速幂
题目链接 因为S(N) , S(N) = A(0)^2 +A(1)^2+--+A(n)^2.所以构造的矩阵一定要维护A(n)^2 s[n-1]=s[n-2]+A[n-1]^2 A[n]=x*A[n-1 ...
最新文章
- 企业安全挑战 内网安全技术十大策略
- POJ1185:火炮(减少国家)
- 200 switching to ascii mode_2020年12月英语四级阅读200篇第106篇:白云之乡—新西兰_四级...
- 理解并实施:GLBP(ccna200-120新增考点)
- ubuntu 64 12.04 oracle,ubuntu server 12.04 x86_64 下安装oracle xe 11 x86_64
- python廖老师_Python3.5-20190518-廖老师-自我笔记-面向对象
- 网桥和交换机的工作原理及区别
- java开发业务流程图,什么是业务流程图?业务流程图如何绘制?
- 宽带波形测试软件,适用于5G时代的波形测试分析系统是怎样的? - 全文
- 【005】基于51单片机的多路热释电红外报警器proteus仿真与实物设计
- 狂神说Linux03:Vim使用及账号用户管理
- ERROR: Could not build wheels for cryptography which use PEP 517 and cannot be installed directly
- RSA算法实现(Python版)
- 计算机数字音乐分手,我的音乐史记之十四:那些关于分手的歌
- 特岗教师计算机专业面试题,2019特岗教师面试试题及参考答案
- 记忆力减退之----SP3232---STM32
- Duplicate keys detected “***“.This may cause can up date error,解决方案。
- 数字图像处理之雷登变换与图像重建(MATLAB)
- P1579 哥德巴赫猜想(升级版)
- cpu温度过高 ubuntu_开源技巧:Ubuntu中检查笔记本CPU温度
热门文章
- ZOJ -3203 三分+数学
- java添加信息_java – 向异常添加信息
- 个人博客网页设计_博客个人网页制作网页设计-1
- python网页爬虫循环获取_Python 爬虫第三篇(循环爬取多个网页)
- exception日志 php_PHP 错误与异常的日志记录
- 使用Jedis源码生成Jedis.jar
- jquery获取java对象的属性_jQuery - 获得内容和属性
- python接口自动化测试框架实战从设计到开发_Python接口自动化测试框架实战 从设计到开发...
- python游戏程序编码_python实现的生成随机迷宫算法核心代码分享(含游戏完整代码)...
- socket工具android,Android通过socket长连接实现推送