Y sequence

Problem's Link: http://acm.hdu.edu.cn/showproblem.php?pid=5297


Mean:

有连续数列A={1,2,3,4,5,6,7,8 .....},将可以表示成a^b次方的数删除,a={1,2,3,4,5...},而2<=b<=r,删除后形成一个新的数列,求这个数列的第n项。

analyse:

很有趣的一道数论题。

对于给定的一个n,如果我们知道1~n中被删除的数字为k个,那么答案一定大于等于n+k,所以向后至少移动k个数。

但是n~n+k这一段中也可能含有被删除的数字,所以我们再求n~n+k这一段中被删除的数字的个数,假设为x1个,那么我们就得到了一个比x更精确的数字x1,即:我们要得到第n项,至少需要从n向后移动x1位。

但是移动x1位后同样还可能存在以上的问题,那么什么时候停止呢?答案是:当本次算的xi和上次算的xi相等时,就说明这个值是固定的了,也就是最终的精确值,代表第n项就是n+xi,也就是最终的答案。

有了这个理论基础后,我们来考虑如何求得1~n中有多少个数字能够表示成次幂形式。

解决这个问题的方法是反函数。幂的反函数是依然是幂函数,指数取导就行。

例:10以内能表示成a^2的数的个数为:pow(10,1/2) ; 100以内能表示成a^4的数的个数是:pow(100,1/4).

还有一个问题:在删除数的时候,a^6的数已经被a^2和a^3的数删过了,这样就造成了重复删除。

这儿就需要容斥原理来解决:只需要删质数次幂的数就行,而且还需要把质数的乘积(多删的)删除的数字加回来。

详情参考大牛博客:http://blog.csdn.net/firstlucker/article/details/46991885

Time complexity: O(N)

Source code: 

/*
* this code is made by crazyacking
* Verdict: Accepted
* Submission Date: 2015-07-22-15.56
* Time: 0MS
* Memory: 137KB
*/
#include <queue>
#include <cstdio>
#include <set>
#include <string>
#include <stack>
#include <cmath>
#include <climits>
#include <map>
#include <cstdlib>
#include <iostream>
#include <vector>
#include <algorithm>
#include <cstring>
#define LL long long
#define ULL unsigned long long
using namespace std;
long long t, n, r;
const int p[] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67};
vector<int> rongchi;
void get_rongchi()
{
rongchi.clear();
for( int i = 0; p[i] <= r; ++i )
{
int si = rongchi.size();
for( int j = 0; j < si; ++j )
{
if( abs( p[i]*rongchi[j] ) <= 63 )
rongchi.push_back( -p[i]*rongchi[j] );
}
rongchi.push_back( p[i] );
}
} long long cal( long long x )
{
if( x == 1 ) return 0;
long long ans = x;
for( int i = 0; i < rongchi.size(); ++i )
{
long long tmp = ( long long )( pow( x + 0.5 , 1.0 / abs( rongchi[i] ) ) ) - 1;
if( rongchi[i] < 0 ) ans += tmp;
else ans -= tmp;
}
return ans - 1;
}
long long solve()
{
get_rongchi();
long long ans = n;
while( true )
{
long long tmp = cal( ans );
if( tmp == n )break;
ans += n - tmp;
}
return ans;
} int main()
{
ios_base::sync_with_stdio( false );
cin.tie( 0 );
cin >> t;
while( t-- )
{
cin >> n >> r;
cout << solve() << endl;
}
return 0;
}
/* */

  

  

2015 Multi-University Training Contest 1 - 10010 Y sequence的更多相关文章

  1. 2015 Multi-University Training Contest 1 OO’s Sequence

    OO’s Sequence Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Others)T ...

  2. 2016 Multi-University Training Contest 5 Divide the Sequence

    Divide the Sequence 题意: 给你一个序列A,问你最多能够分成多少个连续子序列,使得每个子序列的所有前缀和均不小于0 题解: 这题是比赛时候的水题,但我比的时候也就做出这一题, = ...

  3. 2018 Multi-University Training Contest 1 - B Balanced Sequence (贪心)

    题意:对N个由(,)组成的字符串,求拼接后得到的最大的balance序列的长度.balance序列:空串/ A+B(A,B都是b序列)/ (+A+),A为b序列.此三种情况. 分析:在读入N每个字符串 ...

  4. 2015 Multi-University Training Contest 1 y sequence

    Y sequence Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)Total ...

  5. 2015 Multi-University Training Contest 8 hdu 5390 tree

    tree Time Limit: 8000ms Memory Limit: 262144KB This problem will be judged on HDU. Original ID: 5390 ...

  6. Root(hdu5777+扩展欧几里得+原根)2015 Multi-University Training Contest 7

    Root Time Limit: 30000/15000 MS (Java/Others)    Memory Limit: 262144/262144 K (Java/Others)Total Su ...

  7. hdu 5416 CRB and Tree(2015 Multi-University Training Contest 10)

    CRB and Tree                                                             Time Limit: 8000/4000 MS (J ...

  8. 2015 Multi-University Training Contest 10 hdu 5412 CRB and Queries

    CRB and Queries Time Limit: 12000/6000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Other ...

  9. 2015 Multi-University Training Contest 7 hdu 5378 Leader in Tree Land

    Leader in Tree Land Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Othe ...

随机推荐

  1. java.lang.OutOfMemoryError: PermGen space PermGen space & java.lang.OutOfMemoryError: Java heap space Heap siz

    java.lang.OutOfMemoryError: PermGen space PermGen space 由-XX:PermSize  -XX:MaxPermSize 引起 java.lang. ...

  2. DATE_FORMAT函数用法

    一.在oracle中,当想把字符串为‘2011-09-20 08:30:45’的格式转化为日期格式,我们可以使用oracle提供的to_date函数. sql语句为: SELECT to_date(' ...

  3. [SQL SERVER 2005]数据库差异备份及还原

    因为之前遇到还原差异备份,最开始遇到SQLServer报错:”无法还原日志备份或差异备份,因为没有文件可用于前滚“.查阅很多资料后,终于得到解决.收集整理成这篇随笔. 问题原因:出现这种错误绝大多数是 ...

  4. [原]unity5 AssetBundle打包

    本文unity版本5.1.3 一.现有的打包教程: 1.http://liweizhaolili.blog.163.com/blog/static/16230744201541410275298/ 阿 ...

  5. JVM 参数翻译汉化解释

    博客搬家,新地址:http://www.zicheng.net/article/38.htm Behavioral Options(行为参数) Option and Default Value Des ...

  6. 在spring,mybatis整合配置中走的弯路(1)

    在接触一个新东西,总免不了走一些弯路,也正是在这些弯路中,我们不断的成长. 从git上把之前写的代码扒下来,看看我在当初使用spring与mybatis中所走的弯路,路过的君子也可引以为戒. < ...

  7. GitHub使用简单记录

    根据<GotGitHub>[1]所做的一些整理 1. 浏览托管项目 在GitHub的页面中可以使用键盘快捷键 (1)按下问号(?)会在弹出窗口显示当前页面可用的快捷键. (2)在项目的代码 ...

  8. easyui plugin——etreegrid:CRUD Treegrid

    昨天写了一个ko+easyui的同样的实现,感觉写的太乱,用起来十分麻烦,于是今天照着edatagrid,写了一个etreegrid,这样再用ko绑定就方便多了. 使用很简单,$(tableId).e ...

  9. React 根据官方总结的规范

    1.语法上,根据生命周期方法执行的顺序编写代码 (1 生命周期方法[getDefaultProps, getInitialState, componentWillMount, componentDid ...

  10. 转战farbox

    好久没更新博客了,感觉在博客园这种技术为主的博客里面写太多东西写得太杂了,平时的阅读体验非常糟糕,我一直都是一个喜欢记录的人,以前经常有在笔记本上写日记的习惯,但是自从用了网上博客之后,好像就再也没有 ...