An Easy Math Problem
题目就是说任意选取 n n n的两个因子 p q pq pq,问 p q \frac{p}{q} qp的个数有多少,其中限制了 p ≤ q p \le q p≤q
显然当 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 p≤q的条件,对于每一个素因子,我们可以把它放在 p p p或 q q q中,可以同时不放入但不能同时放入,否则两者不互素,那么对于 x i x_i xi就有 2 ∗ y i + 1 2 * y_i + 1 2∗yi+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 ∏(2∗yi+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;
}