2015 Multi-University Training Contest 1 - 10010 Y sequence
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相关推荐
- 2017 Multi-University Training Contest - Team 4 :Wavel Sequence
题目链接 题目翻译: 给出长度为n的序列a, 给出长度为m的序列b. 挑选一些下标f, 挑选一些下表g, 并且这两组下标都是升序排列的. 然后在选出Af == Bg,这些数,组成题目中所要求的波浪,求 ...
- HUST-2015 Multi-University Training Contest 9
2015 Multi-University Training Contest 9 solutions BY xudyh 1001.Expression 记dp_{l,r}dpl,r表示l,rl, ...
- HDU 6051 - If the starlight never fade | 2017 Multi-University Training Contest 2
/* HDU 6051 - If the starlight never fade [ 原根,欧拉函数 ] | 2017 Multi-University Training Contest 2 题意: ...
- 2018 Multi-University Training Contest 3 Problem F. Grab The Tree 【YY+BFS】
传送门:http://acm.hdu.edu.cn/showproblem.php?pid=6324 Problem F. Grab The Tree Time Limit: 2000/1000 MS ...
- 2015 HIAST Collegiate Programming Contest J
Polygons Intersection 题意:给2个凸多边形,求相交面积 思路:不会,套板子就是了 AC代码: #include "iostream" #include &qu ...
- 2015 UESTC Winter Training #10【Northeastern Europe 2009】
2015 UESTC Winter Training #10 Northeastern Europe 2009 最近集训都不在状态啊,嘛,上午一直在练车,比赛时也是刚吃过午饭,状态不好也难免,下次比赛 ...
- hdu 4925 Apple Tree--2014 Multi-University Training Contest 6
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4925 Apple Tree Time Limit: 2000/1000 MS (Java/Others ...
- 2015 German Collegiate Programming Contest (GCPC 15)
2015 German Collegiate Programming Contest (GCPC 15) B. Bounty Hunter II 给定一张DAG,求一种方案:用最少的路径将所有点覆盖. ...
- 2018 Multi-University Training Contest 4 Problem E. Matrix from Arrays 【打表+二维前缀和】
任意门:http://acm.hdu.edu.cn/showproblem.php?pid=6336 Problem E. Matrix from Arrays Time Limit: 4000/20 ...
最新文章
- 第二次作业——个人项目实战
- ABAP Debugging Script(调试器脚本)使用的一些实际例子
- ContextLoaderListener.contextInitialized - how is WebApplicationContext created
- C++之virtual 方法
- leetcode 27. 移除元素(双指针)
- 数据库连接 未将对象引用到实例
- 团伙(信息学奥赛一本通-T1385)
- 数据库-MySQL-Java数据库连接-JDBC
- QT每日一练day5:QLabel和按钮窗口打印功能
- Spring中的面向切面(AOP)
- Maven项目导出War包
- linux安装软件常用命令
- Android 动态权限申请
- Arduino实训-可移动机械臂(智能小车)
- MF,PMF算法比较
- 关于我转行嵌入式的那些事
- 关于DSP28335的CAN模块与上位机(PC)通信调试心得(1)
- 在SpringMVC项目中使用HIkariCP
- 【第01题】A + B | 基础输入输出,开启学习C语言打卡的序章
- OCP-1Z0-051 第120题 NATURAL JOIN自然连接