aboutsummaryrefslogtreecommitdiff
path: root/works/solutions/acwing/1221.cpp
diff options
context:
space:
mode:
authorcrupest <crupest@outlook.com>2021-03-09 16:57:08 +0800
committercrupest <crupest@outlook.com>2021-03-09 16:57:08 +0800
commite0abdf25655c78298937e8a530d56c3e32e804ae (patch)
tree0cdc59ec3dd0afa3b22d42af83374bb3624ad26a /works/solutions/acwing/1221.cpp
parentf08e0f78f20a17bc940a2e2c45949a8967b290db (diff)
downloadcrupest-e0abdf25655c78298937e8a530d56c3e32e804ae.tar.gz
crupest-e0abdf25655c78298937e8a530d56c3e32e804ae.tar.bz2
crupest-e0abdf25655c78298937e8a530d56c3e32e804ae.zip
import(solutions): Add problem 1221.
Diffstat (limited to 'works/solutions/acwing/1221.cpp')
-rw-r--r--works/solutions/acwing/1221.cpp45
1 files changed, 45 insertions, 0 deletions
diff --git a/works/solutions/acwing/1221.cpp b/works/solutions/acwing/1221.cpp
new file mode 100644
index 0000000..9ba4229
--- /dev/null
+++ b/works/solutions/acwing/1221.cpp
@@ -0,0 +1,45 @@
+#include <cmath>
+#include <iostream>
+
+const int M = 5000010;
+const int SM = std::sqrt(M / 2);
+
+int s[M];
+
+int main() {
+ int n;
+ std::cin >> n;
+
+ for (int i = 0; i <= SM; i++) {
+ const int i2 = i * i;
+ const int SMB = M - i2;
+ for (int j = i; j * j < SMB; j++) {
+ const int a = i2 + j * j;
+
+ if (s[a] == 0) {
+ s[a] = i + 1;
+ }
+ }
+ }
+
+ int sm = std::sqrt(n);
+
+ for (int a = 0; a <= sm; a++) {
+ int as = a * a;
+ int bsm = n - 3 * a * a;
+ for (int b = a; b * b <= bsm; b++) {
+ int bs = b * b + as;
+
+ int c = s[n - bs];
+
+ if (c != 0) {
+ c = c - 1;
+ std::cout << a << ' ' << b << ' ' << c << ' '
+ << std::sqrt(n - bs - c * c);
+ return 0;
+ }
+ }
+ }
+
+ return 0;
+}