题意:
给定正整数n和r,定义Y数列为从正整数序列中删除所有能表示成a^b(2 ≤ b ≤ r)的数后的数列,求Y数列的第n个数是多少。例如n = 10, r = 3,则Y数列为2 3 5 6 7 10 11 12 13 14,第10个数是14。
思路:
首先想到二分,如果我们求出了对于mid的位置删去了多少个数,就可以通过二分来得到答案。实际上会卡一卡二分的常数,加一个优化就可以过了。这里用的是迭代的方法。
那么接下来的问题就是如何快速求1~n里面有多少个数能表示成某个数的幂的形式。
我们这样想,首先删去所有的平方数。有n^(1/2)个。再删去所有的三次方数。有n^(1/3)个 …… 这样枚举幂数b(从2开始),然后删去n^(1/b)个数,但是有重复的,比如某个数的6次方可能被平方数删去了一次,又被3次方数删去了一次,于是想到用容斥原理加回来。这样,只需要枚举幂为质数的b的值。枚举质数b的时候,当b超过r就不再枚举。然后再用容斥求最后结果即可。当然,用pow函数是可以的,不过要注意精度。另外容斥的时候注意先不要把1算上去,最后答案减去1(因为1无论如何都是要被删的)即可。
最后发现容斥出来的数的系数就是莫比乌斯系数,从而去学了一些关于容斥和莫比乌斯的内容。
// whn6325689
// Mr.Phoebe
// http://blog.csdn.net/u013007900
#include <algorithm>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <climits>
#include <complex>
#include <fstream>
#include <cassert>
#include <cstdio>
#include <bitset>
#include <vector>
#include <deque>
#include <queue>
#include <stack>
#include <ctime>
#include <set>
#include <map>
#include <cmath>
#include <functional>
#include <numeric>
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
#define eps 1e-9
#define PI acos(-1.0)
#define INF 0x3f3f3f3f
#define LLINF 1LL<<62
#define speed std::ios::sync_with_stdio(false);
typedef long long ll;
typedef unsigned long long ull;
typedef long double ld;
typedef pair<ll, ll> pll;
typedef complex<ld> point;
typedef pair<int, int> pii;
typedef pair<pii, int> piii;
typedef vector<int> vi;
#define CLR(x,y) memset(x,y,sizeof(x))
#define CPY(x,y) memcpy(x,y,sizeof(x))
#define clr(a,x,size) memset(a,x,sizeof(a[0])*(size))
#define cpy(a,x,size) memcpy(a,x,sizeof(a[0])*(size))
#define mp(x,y) make_pair(x,y)
#define pb(x) push_back(x)
#define lowbit(x) (x&(-x))
#define MID(x,y) (x+((y-x)>>1))
#define ls (idx<<1)
#define rs (idx<<1|1)
#define lson ls,l,mid
#define rson rs,mid+1,r
#define root 1,1,n
template<class T>
inline bool read(T &n)
{
T x = 0, tmp = 1;
char c = getchar();
while((c < '0' || c > '9') && c != '-' && c != EOF) c = getchar();
if(c == EOF) return false;
if(c == '-') c = getchar(), tmp = -1;
while(c >= '0' && c <= '9') x *= 10, x += (c - '0'),c = getchar();
n = x*tmp;
return true;
}
template <class T>
inline void write(T n)
{
if(n < 0)
{
putchar('-');
n = -n;
}
int len = 0,data[20];
while(n)
{
data[len++] = n%10;
n /= 10;
}
if(!len) data[len++] = 0;
while(len--) putchar(data[len]+48);
}
//-----------------------------------
const int mi[19]={-2,-3,-5,-7,-11,-13,-17,-19,-23,-29,-31,-37,-41,-43,-47,-53,-59,-61,-67};
long long n;
int r;
vector <int> mu;
void get_mu()
{
mu.clear();
for(int i=0; abs(mi[i])<=r; i++)
{
int temp=mu.size();
for(int j=0; j<temp; j++)
if(abs(mi[i]*mu[j]) <= 63)
mu.push_back(mi[i]*mu[j]);
mu.push_back(mi[i]);
}
}
ll cal(ll x)
{
if(x==1)
return 0;
ll ans=x;
for(int i=0; i<mu.size(); i++)
{
ll temp =ll(pow(x + 0.5, 1.0/abs(mu[i])))-1;
if(mu[i]<0)
ans-=temp;
else
ans+=temp;
}
return ans-1;
}
int main()
{
int T;
read(T);
while(T--)
{
read(n),read(r);
get_mu();
ll ans = n;
while(1)
{
ll temp=cal(ans);
if(temp==n)
break;
ans+=n-temp; //每次加上被删去的数的个数
}
write(ans),putchar('\n');
}
return 0;
}