SWERC 2011 / HDU 4196 Remoteland (数论&想法题)
2014-01-13 20:39
274 查看
Remoteland
Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 262144/131072 K (Java/Others)
http://acm.hdu.edu.cn/showproblem.php?pid=4196
思路:
显然取n!是最大的,但这不一定是一个完全平方数,需要把多余的部分除掉。
可以利用勒让德定理很快处理出n!所有素因子的指数,偶数保留,奇数减1,即可保证这些素因子最后乘积为完全平方数,也就是用n!除以所有指数为奇数的素因子的乘积。注意由于有取模运算,所以不能直接除。——那怎么办?
换个思路,先把n!除以<=n的所有质数,再乘上指数为偶数的质数!
(先乘再除转化为先除再乘!)
完整代码:
/*1718ms,90976KB*/ #include<cstdio> #include<cmath> const long long mod = 1000000007; const int mx = 10000001; const int sqrtmx = (int)sqrt((double)mx); long long fac[mx] = {1, 1}; bool vis[mx]; int prime[664580], cnt; void init() { for (int i = 2; i < mx; ++i) { fac[i] = fac[i - 1]; if (!vis[i]) { prime[cnt++] = i; if (i <= sqrtmx) for (int j = i * i; j < mx; j += i) vis[j] = true; } else fac[i] = (fac[i] * i) % mod; } } int main() { init(); int n; while (scanf("%d", &n), n) { long long res = fac ; for (int i = 0; i < cnt && prime[i] <= (n >> 1); ++i) { int r = 0, tmp = n; while (tmp) r += (tmp /= prime[i]); if ((r & 1) == 0) res = (res * prime[i]) % mod; } printf("%lld\n", res); } return 0; }
相关文章推荐
- 山寨腾讯“爱消除”游戏之框架总结
- POJ 2784 Buy or Build 最小生成树
- POJ 2784 Buy or Build 最小生成树
- IT与户外精神 推荐
- bh第一学期总结
- 2014加油
- 异步I/0基础
- telnet服务器源码安装
- URAL 1023 Buttons 博弈
- URAL 1023 Buttons 博弈
- iPhone学习第一天,iphoen基本控件的使用
- Octave对Windows操作系统的支持从3.6.1版本开始就已经全面升级了
- java 正则表达式(转)
- hpple使用几个注意点
- 机器学习&数据挖掘笔记_25(PGM练习九:HMM用于分类)
- 游戏开发随手记:cocos2d ccLayer响应触摸事件方法
- 套用今天的签名“最新的VS+最新的MATLAB+最新的WINDOWS=无尽的烦恼 ”,用了一天的时间去配置相关的环境,确实很囧。好在现在问题解决了,特发布过程如下,希望能给有需要的朋友带来帮助。
- 远程无法访问mysql Can't connect to MySQL server on "host" (111)
- 【leetcode】 Single Number
- 动态加载的 CTreeCtrl 消息响应