CCPC2020 - 秦皇岛 - G. Good Number (数学)

亚历克斯喜欢数字。

亚历克斯认为,正整数 x x x 是好数,当且仅当 ⌊ x k ⌋ \lfloor \sqrt[k]{x} \rfloor kx 整除 x x x

你能告诉他不超过 n n n 的正整数的个数吗?

输入

输入的第一行给出了测试用例的数量 T ( 1 ≤ T ≤ 10 ) T (1 \le T \le 10) T(1T10) 。接下来是 T T T 个测试用例。

对于每个测试用例,唯一的一行包含两个整数 n n n k ( 1 ≤ n , k ≤ 1 0 9 ) k (1 \le n,k \le 10^9) k(1n,k109)

输出

对于每个测试用例,输出一行包含 "Case #x: y"的内容,其中 x \texttt{x} x 是测试用例编号(从 1 1 1 开始), y \texttt{y} y 是答案。

Example
input

2
233 1
233 2

output

Case #1: 233
Case #2: 43

这道题需要找规律,如果打表把所有的能够符合要求的数打出来的话,就会发现他们有的是对应了连续几个 ⌊ x k ⌋ \lfloor \sqrt[k]{x} \rfloor kx
比如n为250,k为2的情况。
在这里插入图片描述
就比如 9 9 9 对应 3 3 3 12 12 12 15 15 15 也对应 3 3 3
这里这么看, 3 3 3 k k k 次方对应了 9 9 9 ,而 4 4 4 的k次方对应了 16 16 16
所以发现 [ n k , ( n + 1 ) k ) [ n^k , (n+1)^k) [nk,(n+1)k) 这个区间内开 k k k 次方,下取整后能够得到的数都是 n n n。所以就只需要求每一个区间内有多少能够整除 n n n 的数,我们看到 9 9 9 12 12 12 15 15 15,每个都相差了 3 3 3,所以对于区间 [ a , b ] [a,b] [a,b],能够被 n n n 整除的个数就是 ⌊ b n ⌋ − ⌊ a − 1 n ⌋ \lfloor\frac{b}{n}\rfloor - \lfloor\frac{a-1}{n}\rfloor nbna1

这么算完之后,如果遇见了 x k < = n < = ( x + 1 ) k x^k <= n <= (x+1)^k xk<=n<=(x+1)k 的情况,那么就在这一次取完答案之后跳出。

并且注意,因为题目给出要求 k k k 是小于等于 1 e 9 1e9 1e9,但是当 k k k 大于等于 31 31 31 的时候,整个 i n t int int 范围内的数都只能下取整到 1 1 1 了,所以这种情况直接特判。
还有 k = 1 k=1 k=1 的时候,直接输出 n n n

至于下取整,直接强转 i n t int int 即可。

#include<iostream>
#include<map>
#include<cmath>
using namespace std;
#pragma warning (disable:4996)
#define ll long long

ll qmi(ll a, ll k) {//需要快速幂,不然TLE
	ll res = 1;
	while (k) {
		if (k & 1)res = (ll)res * a;
		k >>= 1;
		a = (ll)a * a;
	}
	return res;
}

int main() {
	int T; cin >> T;
	for (int cases = 1; cases <= T; cases++) {
		int n, k; cin >> n >> k;

		if (k == 1 || k >= 31) {
			printf("Cases #%d: %d\n", cases, n);
			continue;
		}

		int res = 0;
		int i = 0;
		while (++i) {
			ll l = qmi(i, k), r = qmi(i + 1, k) - 1;
			
			if (l <= n && n <= r) {
				if (l == n)res++;
				else {
					res += (n - l) / i + 1;
				}
				break;
			}

			res += (r - l) / i + 1;
		}

		printf("Case #%d: %d\n", cases, res);
	}
	return 0;
}
  • 14
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值