The Super Powers UVA 11752 分析分析 求无符号长整形以内的数满足至少可以用两种不同的次方来表示。比如64 = 2^6 = 8^2; 一个数的1次方不算数。
/**
题目:The Super Powers UVA 11752
链接:https://vjudge.net/contest/154246#problem/Y
题意:求无符号长整形以内的数满足至少可以用两种不同的次方来表示。比如64 = 2^6 = 8^2; 一个数的1次方不算数。
思路:
分析过程如下:
1 = 1^1 1^2 1^3
16 = 2^4 4^2
64 = 2^6 8^2
81 = 3^4 9^2
256 = 2^8 16^2
512 = 2^9 8^3 设Max为最大的可能获得的值。
Max = a^e1 or a^e2 or ... (2<=e1<e2) a <= Max开三次方根=[1,1e7] 前面我猜想一个满足条件的数是某一个素因子的合数次方。 (2*3)^4 = (4*9)^2 含有多个素数。实际上不止一个素因子。 那么底数在(Max开三次方根)内是任意可能的。 有一点可以肯定:一个数的合数次方一定是满足条件的数。 底数在[1,1e7]以内。然后合数范围在64以内。然后set存数并去重;时间复杂度过大。 避免使用set?尝试提前去重。 2^4, 2^6, 2^8, 2^9, 2^10, 2^12, 2^14, 2^15. 取2^8 = 4^4 重复了。 2^12 = 4^6 = 8^4 (保证幂是合数) 重复了。 而:2^4, 2^6, 2^9, 2^10, 2^14, 2^15不会再重复。因为幂都是两个素数相乘的结果。 如何通过现在的某些计算来预测将来的某些数不用再重复计算?
或者通过一些计算来判定当前的数已经算过了。 4^4, 4^6, 4^8, 4^9, 4^10, 4^12. 好像这些都是重复了的。那么。。。所有的底数为合数都会被重复。 只要计算底数为素数的就行了? (2*3)^4 = (4*9)^2 含有多个素数。从哪里得来? 再试试其他合数为底数的情况。 6^4, 6^6, 6^5, 6^8, 6^9, 6^10. 6^4没有重复过。并不是所有的底数为合数都会被重复。 2^12 = 8^4 可以得知如果以前出现过。那么底数应该可以开根号才行,即底数可以表示成一个数的次方,次方>1。 新结论:底数如果是某个数的次方,那么它的所有合数幂都重复过。其他情况均没有重复过。 问题好像重复了。都是某个数的次方。只不过题目为至少两个。 那么对每一个非次方数的底数a。它的贡献为a^e. (e为合数,且a^e<=2^64 - 1) 还是往这方面处理:如何通过现在的某些计算来预测将来的某些数不用再重复计算? 预处理64以内的合数的约数对,即:两个约数相乘等于本身为一对。这样通过约束对来判断是否那些底数会重复。 当两个约数至少一个>=4的合数时候,比如d1,d2;那么(x^d1)^d2 or (x^d2)^d1 都要去重。 解题方法:
预处理64以内的约数。以及每一个约数的约数对。
枚举底数a在[1,1e7]内。
对每一个底数a枚举它约数次方。并如果次方可以分解成两个约数,其中一个是>=4的约数。那么底数^另一个约数形成的数不能再作为计算的底数。 猜想:如果时间超限,适当降低[1,1e7]右边界的范围。 为了避免底数的合数次方超出unsigned long long 范围,所以用log对数来处理。 其实应该可以很快想到的,但是走了很多弯路。积累经验。多思考。
*/ #include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <set>
#include<vector>
#include <cmath>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int maxn = 1e5+;
const double eps = 1e-;
vector<int> v[];//存储约数对。
int isHe[];//isHe[i] == 1表示这是一个合数。
int he[], z;//存储合数。
ull ans[], cnt;//注意这里不能和maxn一直相同。因为ans数量可能更多。
//经测试最多67385个;
int flag[maxn];
ll f(ll a,ll b)
{
ll p = ;
while(b){
if(b&) p*=a;
a = a*a;
b >>= ;
}
return p;
}
void init()
{
memset(isHe, , sizeof isHe);
for(int i = ; i < ; i++){
if(isHe[i]==){
for(int j = i*i; j < ; j+=i){
isHe[j] = ;
}
}
}
z = ;
for(int i = ; i < ; i++){
if(isHe[i]){
he[z++] = i;
for(int j = ; j*j < ; j++){
if(i%j==){
if(he[j]||he[i/j]){
v[i].push_back(j);
v[i].push_back(i/j);
}
}
}
}
}
ull mas = (<<)-;
//cout<<"mas = "<<mas<<endl;
cnt = ;
for(int i = ; i < maxn; i++){
if(flag[i]) continue;
ull p = i, num = ;
for(int j = ; j < z; j++){
if(he[j]>log(mas)/log(i)-eps) break;///不加eps,会存在误差,比如2^64次方可能会算进去。然后p = 0;
while(num<he[j]){
num++;
p *= i;
}
if(num<he[j]) break;
ans[cnt++] = p;
int len = v[he[j]].size();
for(int k = ; k < len; k+=){
if(isHe[v[he[j]][k]]){
ull temp = f(i,v[he[j]][k+]);
if(temp>=maxn) continue;
flag[temp] = ;
}
if(isHe[v[he[j]][k+]]){
ull temp = f(i,v[he[j]][k]);
if(temp>=maxn) continue;
flag[temp] = ;
}
}
}
}
sort(ans,ans+cnt);
printf("1\n");
for(int i = ; i < cnt; i++){
printf("%llu\n",ans[i]);
}
}
int main()
{
init(); return ;
}
The Super Powers UVA 11752 分析分析 求无符号长整形以内的数满足至少可以用两种不同的次方来表示。比如64 = 2^6 = 8^2; 一个数的1次方不算数。的更多相关文章
- The Super Powers UVA - 11752
题目大意:将范围从1~pow(2,64)-1内的super power输出.super power的定义:一个数x至少存在两种x=pow(i,k),(k!=1). 题解: 注意数据范围2的64次方-1 ...
- The Super Powers UVA - 11752(合数幂)
题意: 求1~2^64-1之间所有的 至少是两个不同的正整数的幂的数 升序输出 一个数的合数次幂即为这样的数 找出1~2^64-1中所有数的合数次幂 用set存起来(既能防止重复 又能升序) 最后输 ...
- linux内核分析作业4:使用库函数API和C代码中嵌入汇编代码两种方式使用同一个系统调用
系统调用:库函数封装了系统调用,通过库函数和系统调用打交道 用户态:低级别执行状态,代码的掌控范围会受到限制. 内核态:高执行级别,代码可移植性特权指令,访问任意物理地址 为什么划分级别:如果全部特权 ...
- 【google面试题】求1到n的正数中1出现的次数的两种思路及其复杂度分析
问题描写叙述: 输入一个整数n,求从1到n这n个整数的十进制表示中1出现的次数.比如输入12,从1到12这些整数中包括1 的数字有1.10.11和12.1一共出现了5次. 这是一道广为流传的googl ...
- 求GCD(最大公约数)的两种方式
求GCD(最大公约数)的两种方式 这篇随笔讲解C++语言程序设计与应用中求GCD(最大公约数,下文使用GCD代替)的两种常用方式:更相减损法和辗转相除法,前提要求是具有小学数学的基本素养,知道GCD是 ...
- uva 11752 The Super Powers 素数+大数判断大小
题目链接: https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_proble ...
- python基础练习题(题目 求s=a+aa+aaa+aaaa+aa…a的值,其中a是一个数字。例如2+22+222+2222+22222(此时共有5个数相加),几个数相加由键盘控制)
day11 --------------------------------------------------------------- 实例018:复读机相加 题目 求s=a+aa+aaa+aaa ...
- 2014 Super Training #8 G Grouping --Tarjan求强连通分量
原题:ZOJ 3795 http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3795 题目大意:给定一个有向图,要求把点分为k个集 ...
- 创建TabHost的两种方式的简单分析
最近做了一个TabHost的界面,在做的过程中发现了一些问题,故和大家分享一下. 首先我的界面如下: 目前就我所知,创建TabHost有两种方式,第一种是继承TabActivity类,然后用getTa ...
随机推荐
- Asp.Net MVC part2 View、Controller详解
View详解Razor视图引擎简介HtmlHelper强类型页面 Razor视图引擎简介强大的@:表示使用C#代码,相当于aspx中的<%%>可以完成输出功能当遇到html标签时会认为C# ...
- SeaJS入门
Sea.js是一种模块加载工具 官网:http://seajs.org/docs/ 使用步骤: 1.下载sea.js 2.引入sea.js 3.加入配置代码 // seajs 的简单配置 seajs. ...
- expdp impdp 数据库导入导出命令详解
一.创建逻辑目录,该命令不会在操作系统创建真正的目录,最好以system等管理员创建.create directory dpdata1 as 'd:\test\dump'; 二.查看管理理员目录(同时 ...
- Kingdee Apusic 中间件有关资料
Kingdee Apusic 中间件有关资料: 1.官方网站:http://www.apusic.com 2.资料目录:http://www.apusic.com/dist 3.Apusic 8 资料 ...
- Plplot中line width 问题
Plplot延续了Pgplot的设计风格,线的宽度函数plwid(int width)只能用整型变量,不能精确控制线宽.用pscairo画出的曲线总是太粗,这是plplot一个很老的问题(issue) ...
- [ubuntu Setup] ubuntu 14.10 安装 JDK
from : http://www.cnblogs.com/plinx/archive/2013/06/01/3113106.html 1.到 Sun 的官网下载 http://www.oracle ...
- opencv实现camera模组的暗电流和lenshading补偿 .
目录(?)[-] 简介 基本原理 产生原因 校正补偿原理 具体实现 框架搭建 功能实现 暗电流 lenshading补偿 效果演示 图片处理 效果演示 简介 在接触过的qcom和mtk平台中,came ...
- 读取Properties配置文件, 四种方式都可以得到webroot/WEB-INF/classes这个路径
下面四种方式都可以得到webroot/WEB-INF/classes这个路径,有什么区别,哪种方式最好? String path = this.getClass().getResource(" ...
- PowerDesigner 将表的字段name属性设置到comment凝视
1.首先copy以下的vbs脚本.并将其另存为name2comment.vbs '*********************************************************** ...
- python——TypeError: 'str' does not support the buffer interface
import socket import sys port=51423 host="localhost" data=b"x"*10485760 #在字符串前加 ...