2020 CCPC 秦皇岛 G - Good Number（数学+分块）

传送门

---

题目大意

一个数$x$被定义为$Goodnumber$当且仅当$\lfloor \sqrt[k]{x}\rfloor$能整除$x$。对于给定的$n,k$，求出$[1,n]$中有多少$Goodnumber$。

解题思路

因为开$k$次方的算递减的非常快的，对于$int$范围内的数，显然$k\ge 32$时向下取整就已经变成$1$了，$1$一定整除所有数那么答案为$n$。还有一个特殊情况就是$k=1$时，每个数最后得到的都是它本身，那么答案也为$n$。

不难发现实际上结果是分块的，也就是对于某一$k$，$[p^k,(p+1{)}^k)$范围内的数开$k$次方得到的答案都为$p$。但是这个区间内整除$p$的数有多少个怎么算？

首先不难想到一个十分经典的数学问题$[1,n]$以内有多少数能被$x$整除——答案就是$\lfloor \frac{n}{x}\rfloor$，那么不难想到可以使用容斥去解决这个问题，也就是对于区间$[a,b]$求有多少个能被$x$整除的数：$\lfloor \frac{b}{x}\rfloor -\lfloor \frac{a-1}{x}\rfloor$。还有一种思路是因为左边界恰好被$p$整除，那么我们就单独计算左边界，然后求$(p^k+1,(p+1{)}^k)$有多少能被$p$整除，这时模$p$的剩余系是从$1$开始的，那么直接$\lfloor \frac{(p+1{)}^k-p^k}{p}\rfloor$就能得到正确结果——$1+\lfloor \frac{(p+1{)}^k-p^k}{p}\rfloor$。

最后就是边界问题，只需判断什么时候$p^k\le n\le (p+1{)}^k$，然后右界设为$n$，统计结束后$break$

扩展

前不久遇到了一个求$[1,n]$有多少平方数的问题，细想一下也就是$\lfloor \sqrt{n}\rfloor$。

感悟

现场赛被这题卡住实在是心态不够好

#include <bits/stdc++.h>
#include <unordered_map>
#include <unordered_set>

using namespace std;
#define fi first
#define se second
#define pb push_back
#define ins insert
#define Vector Point
#define ENDL "\n"
#define lowbit(x) (x&(-x))
#define mkp(x, y) make_pair(x,y)
#define mem(a, x) memset(a,x,sizeof a);
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef pair<double, double> pdd;
const double eps = 1e-8;
const double pi = acos(-1.0);
const int inf = 0x3f3f3f3f;
const double dinf = 1e300;
const ll INF = 1e18;
const int Mod = 1e9 + 7;
const int maxn = 2e5 + 10;

ll n, k;
int T;

ll qkp(ll x, ll n) {
    ll ans = 1;
    while (n) {
        if (n & 1) ans *= x;
        x *= x;
        n >>= 1;
    }
    return ans;
}

ll solve() {
    if (k == 1 || k >= 32) return n;
    ll ans = 0, l, r;
    for (int i = 1;; i++) {
        l = qkp(i, k), r = qkp(i + 1, k) - 1;
        //cout << l << " " << r << endl;
        if (l <= n && n <= r) {
            if (l == n) ans++;
            else {
                r = n;
                ans += 1 + (r - l) / i;
                //ans += r / i - (l - 1) / i;
            }
            break;
        }
        ans += 1 + (r - l) / i;
        //ans += r / i - (l - 1) / i;
    }
    return ans;
}

int main() {
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    ios_base::sync_with_stdio(0), cin.tie(0), cout.tie(0);
    int kase = 0;
    cin >> T;
    while (T--) {
        cin >> n >> k;
        cout << "Case #" << ++kase << ": ";
        cout << solve() << endl;
    }
    return 0;
}