2023XCPC个人补题题解

An Easy Math Problem

题目就是说任意选取 n n n的两个因子 p q pq pq,问 p q \frac{p}{q} qp的个数有多少,其中限制了 p ≤ q p \le q pq

显然当 g c d ( p , q ) > 1 gcd(p,q) > 1 gcd(p,q)>1时,化为最简分数后一定会与其他互素的一组 p q pq pq相同,所以我们只需要考虑 p q pq pq互素的情况就好了。

那么如何考虑 p q pq pq互素呢?由于 p q pq pq均为 n n n的因子,那么 p q pq pq一定由 n n n的素因子组成,然后就可以考虑如何选数使得 p p p q q q互素。

由唯一分解定理,设 n = x 1 y 1 n={x_1}^{y_1} n=x1y1 x 2 y 2 {x_2}^{y_2} x2y2 ⋯ \cdots

先不考虑 p ≤ q p \le q pq的条件,对于每一个素因子,我们可以把它放在 p p p q q q中,可以同时不放入但不能同时放入,否则两者不互素,那么对于 x i x_i xi就有 2 ∗ y i + 1 2 * y_i + 1 2yi+1种方案数。除了 p = 1 p=1 p=1 q = 1 q=1 q=1的情况,其他的 p q \frac{p}{q} qp都是两两对应的,所以最终答案为 ∏ ( 2 ∗ y i + 1 ) / 2 + 1 \prod(2 * y_i + 1) / 2 + 1 (2yi+1)/2+1

时间复杂度 O ( q × n ) O(q \times \sqrt{n}) O(q×n )

#include <bits/stdc++.h>

using i64 = long long;

void solve() {
	i64 n;
	std::cin >> n;
	std::vector<std::pair<i64, int>> fac;
	for (i64 i = 2; i <= n / i; i++) {
		if (n % i == 0) {
			int cnt = 0;
			while (n % i == 0) {
				n /= i;
				cnt += 1;
			}
			fac.push_back({i, cnt});
		}
	}
	if (n != 1) {
		fac.push_back({n, 1});
	}

	int ans = 1;
	for (auto [_, y] : fac) {
		ans *= 2 * y + 1;
	}

	ans = ans / 2 + 1;
	std::cout << ans << "\n";
	return;
}

int main() {
	std::ios::sync_with_stdio(false);
	std::cin.tie(nullptr);

	int t;
	std::cin >> t;

	while (t--) {
		solve();
	}

	return 0;
}
Elimination Series Once More

淘汰赛制其实就是一个树形结构,仔细观察还能发现这是一个满二叉树, 2 n 2^{n} 2n个选手至多会打 n n n场比赛。

对于每一个位置,我们可以枚举他能够胜多少局,其能否取胜其实与他在这个子树中值的排名以及还有多少剩余的交换次数有关,因此我们只要按照选手的强弱值从小到大加入我们构造的这棵树就好了。

从这个题目中也可以学到对于一棵树,从子节点找父节点的方法:其实是通过编号,而不是具体的区间,顺便也复习一下树的编号吧。

#include <bits/stdc++.h>

using i64 = long long;

int main() {
	std::ios::sync_with_stdio(false);
	std::cin.tie(nullptr);

	int n, k;
	std::cin >> n >> k;

	const int N = 1 << n;
	std::vector<int> b(1 << (n + 1)), c(1 << (n + 1));
	for (int i = 1 << n; i < (1 << (n + 1)); i++) {
		int x;
		std::cin >> x;
		b[x] = i;
	}

	std::vector<int> ans(1 << n);
	for (int i = 1; i <= (1 << n); i++) {
		int sum = 0;
		for (int x = b[i], now = 1; x; x /= 2, now *= 2, sum += 1) {
			if (now * 2 > i) {
				break;
			}
			if (c[x / 2] + k + 1 < now * 2) {
				break;
			}
		}

		ans[b[i] - (1 << n)] = sum;
		for (int x = b[i]; x; x /= 2) {
			c[x] += 1;
		}
	}

	for (int i = 0; i < (1 << n); i++) {
		std::cout << ans[i] << " \n"[i == (1 << n) - 1];
	}

	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值