牛客 215E 黄魔法师 题解

Description

给出 n,kn, k,求一个长度为 nn 的数组 aa, 满足有恰好 kk 对数对 (i,j)(1i<jn)(i, j) (1 \leq i < j \leq n) 满足 ai+aja_i + a_j 为完全平方数。如果不存在,输出 1-1

link

Solution

显然如果 k>(n2)k>\binom{n}{2} 就一定无解。

构造时会发现肯定要尽量弄成相同的然后进行微调,那么设 mm 为最大的数满足 (m2)k\binom{m}{2}\leq kr=k(m2)r=k-\binom{m}{2}

这个时候直接选 mm22 会发现多出来的 rr 很难微调,因为这个时候要是调整加数的话最小就是增加 mm 了。

这时可以把 mm 个数拆成 rrAAmrm-rBB,然后找 11CC 使得 A+B,2A,2B,A+CA+B,2A,2B,A+C 均为完全平方数并且其他的都不是完全平方数,剩下多的 nm1n-m-1DD 只要随便找一个数使得加出来不是完全平方数即可。

经枚举 A=2,B=98,C=7,D=1A=2,B=98,C=7,D=1 可以满足条件。

时间复杂度:O(n)O(n)

Code

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
#include <bits/stdc++.h>

#define int int64_t

const int kMaxN = 1e5 + 5;

int n, k, m, r;

bool check(int x) {
int y = sqrtl(x);
for (; y * y > x; --y) {}
for (; (y + 1) * (y + 1) <= x; ++y) {}
return x == y * y;
}

void dickdreamer() {
std::cin >> n >> k;
if (k > n * (n - 1) / 2) return void(std::cout << "-1\n");
for (m = 1; m * (m + 1) / 2 <= k; ++m) {}
if (n == m) {
for (int i = 1; i <= n; ++i) std::cout << "2 ";
return;
}
r = k - m * (m - 1) / 2;
for (int i = 1; i <= r; ++i) std::cout << "2 ";
for (int i = 1; i <= m - r; ++i) std::cout << "98 ";
std::cout << "7 ";
for (int i = 1; i <= n - m - 1; ++i) std::cout << "1 ";
}

int32_t main() {
#ifdef ORZXKR
freopen("in.txt", "r", stdin);
freopen("out.txt", "w", stdout);
#endif
std::ios::sync_with_stdio(0), std::cin.tie(0), std::cout.tie(0);
int T = 1;
// std::cin >> T;
while (T--) dickdreamer();
// std::cerr << 1.0 * clock() / CLOCKS_PER_SEC << "s\n";
return 0;
}

牛客 215E 黄魔法师 题解
https://sobaliuziao.github.io/2024/05/04/post/53ed080d.html
作者
Egg_laying_master
发布于
2024年5月4日
许可协议