#include <bits/stdc++.h>
struct Frac {
__int128_t n, d;
Frac(__int128_t _n, __int128_t _d) {
n = _n, d = _d;
__int128_t g = std::__gcd(n, d); n /= g, d /= g;
if (d < 0) n *= -1, d *= -1;
}
Frac(__int128_t _n) : Frac(_n, 1) {}
Frac() : Frac(0) {}
friend Frac abs(Frac F) { return Frac(abs(F.n), F.d); }
friend bool operator < (const Frac& l, const Frac& r) { return l.n * r.d < r.n * l.d; }
friend bool operator <= (const Frac& l, const Frac& r) { return l.n * r.d <= r.n * l.d; }
friend bool operator > (const Frac& l, const Frac& r) { return l.n * r.d > r.n * l.d; }
friend bool operator >= (const Frac& l, const Frac& r) { return l.n * r.d >= r.n * l.d; }
friend bool operator == (const Frac& l, const Frac& r) { return l.n == r.n && l.d == r.d; }
friend bool operator != (const Frac& l, const Frac& r) { return !(l == r); }
Frac operator - () const { return Frac(-n, d); }
friend Frac operator + (const Frac& l, const Frac& r) { return Frac(l.n * r.d + r.n * l.d, l.d * r.d); }
friend Frac operator - (const Frac& l, const Frac& r) { return Frac(l.n * r.d - r.n * l.d, l.d * r.d); }
friend Frac operator * (const Frac& l, const Frac& r) { return Frac(l.n * r.n, l.d * r.d); }
friend Frac operator * (const Frac& l, int r) { return l * Frac(r, 1); }
friend Frac operator * (int r, const Frac& l) { return l * r; }
friend Frac operator / (const Frac& l, const Frac& r) { return l * Frac(r.d, r.n); }
friend Frac operator / (const Frac& l, const int& r) { return l / Frac(r, 1); }
friend Frac operator / (const int& l, const Frac& r) { return Frac(l, 1) / r; }
friend Frac& operator += (Frac& l, const Frac& r) { return l = l + r; }
friend Frac& operator -= (Frac& l, const Frac& r) { return l = l - r; }
template <class T> friend Frac& operator *= (Frac& l, const T& r) { return l = l * r; }
template <class T> friend Frac& operator /= (Frac& l, const T& r) { return l = l / r; }
};
inline char gc() { // like getchar()
static char buf[1 << 16];
static size_t bc, be;
if (bc >= be) {
buf[0] = 0, bc = 0;
be = fread(buf, 1, sizeof(buf), stdin);
}
return buf[bc++]; // returns 0 on EOF
}
int read_int() {
int a, c;
while ((a = gc()) < 40);
if (a == '-') return -read_int();
while ((c = gc()) >= 48) a = a * 10 + c - 480;
return a - 48;
}
int main() {
using namespace std;
typedef __int128_t ll;
ios_base::sync_with_stdio(0);
int N, L;
N = read_int();
L = read_int();
vector<vector<ll>> V(N, vector<ll>(L));
for (int i = 0; i < N; i++)
for (int j = 0; j < L; j++) {
int x; x = read_int();
V[i][j] = x;
}
vector<vector<Frac>> split(N, vector<Frac>(N + 1));
for (int i = 0; i < N; i++) {
Frac need = 0;
for (int j = 0; j < L; j++)
need += V[i][j];
need /= N;
split[i][0] = 0;
Frac have = 0;
int it = 0;
Frac at = 0;
for (int j = 1; j <= N; j++) {
while (have < need) {
Frac add = (1 - at) * V[i][it];
if (have + add <= need) {
have += add;
at = 0;
it++;
} else {
Frac go = (need - have) / V[i][it];
have = need;
at += go;
}
}
split[i][j] = it + at;
have = 0;
}
}
const int INF = 1e9;
vector<Frac> splits(N);
vector<int> ord(N);
vector<bool> used(N);
for (int i = 1; i <= N - 1; i++) {
pair<Frac, int> best = {INF, INF};
for (int j = 0; j < N; j++)
if (!used[j])
best = min(best, {split[j][i], j});
splits[i - 1] = best.first;
ord[i - 1] = best.second;
used[best.second] = true;
}
for (int i = 0; i < N; i++) {
if (!used[i])
ord.back() = i;
}
for (int i = 1; i < N - 1; i++)
splits[i] = max(splits[i], splits[i - 1]);
for (int i = 0; i < N - 1; i++) {
assert(splits[i - 1] <= splits[i]);
}
for (int i = 0; i < N - 1; i++)
cout << (long long) splits[i].n << " " << (long long) splits[i].d << '\n';
for (int i = 0; i < N; i++)
cout << ord[i] + 1 << " ";
cout << '\n';
return 0;
}
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
364 KB |
Output is correct |
2 |
Correct |
1 ms |
364 KB |
Output is correct |
3 |
Correct |
1 ms |
492 KB |
Output is correct |
4 |
Correct |
1 ms |
492 KB |
Output is correct |
5 |
Correct |
1 ms |
492 KB |
Output is correct |
6 |
Correct |
1 ms |
364 KB |
Output is correct |
7 |
Correct |
1 ms |
364 KB |
Output is correct |
8 |
Correct |
1 ms |
364 KB |
Output is correct |
9 |
Correct |
1 ms |
364 KB |
Output is correct |
10 |
Correct |
1 ms |
492 KB |
Output is correct |
11 |
Correct |
1 ms |
492 KB |
Output is correct |
12 |
Correct |
1 ms |
492 KB |
Output is correct |
13 |
Correct |
1 ms |
492 KB |
Output is correct |
14 |
Correct |
1 ms |
492 KB |
Output is correct |
15 |
Correct |
1 ms |
492 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
364 KB |
Output is correct |
2 |
Correct |
1 ms |
492 KB |
Output is correct |
3 |
Correct |
1 ms |
492 KB |
Output is correct |
4 |
Correct |
2 ms |
492 KB |
Output is correct |
5 |
Correct |
1 ms |
492 KB |
Output is correct |
6 |
Correct |
1 ms |
492 KB |
Output is correct |
7 |
Correct |
1 ms |
364 KB |
Output is correct |
8 |
Correct |
1 ms |
364 KB |
Output is correct |
9 |
Correct |
2 ms |
492 KB |
Output is correct |
10 |
Correct |
2 ms |
492 KB |
Output is correct |
11 |
Correct |
1 ms |
492 KB |
Output is correct |
12 |
Correct |
1 ms |
364 KB |
Output is correct |
13 |
Correct |
1 ms |
492 KB |
Output is correct |
14 |
Correct |
2 ms |
492 KB |
Output is correct |
15 |
Correct |
1 ms |
492 KB |
Output is correct |
16 |
Correct |
2 ms |
492 KB |
Output is correct |
17 |
Correct |
2 ms |
492 KB |
Output is correct |
18 |
Correct |
2 ms |
620 KB |
Output is correct |
19 |
Correct |
2 ms |
620 KB |
Output is correct |
20 |
Correct |
2 ms |
492 KB |
Output is correct |
21 |
Correct |
2 ms |
620 KB |
Output is correct |
22 |
Correct |
2 ms |
620 KB |
Output is correct |
23 |
Correct |
1 ms |
364 KB |
Output is correct |
24 |
Correct |
1 ms |
492 KB |
Output is correct |
25 |
Correct |
1 ms |
492 KB |
Output is correct |
26 |
Correct |
1 ms |
492 KB |
Output is correct |
27 |
Correct |
2 ms |
492 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
364 KB |
Output is correct |
2 |
Correct |
1 ms |
364 KB |
Output is correct |
3 |
Correct |
1 ms |
492 KB |
Output is correct |
4 |
Correct |
1 ms |
492 KB |
Output is correct |
5 |
Correct |
1 ms |
492 KB |
Output is correct |
6 |
Correct |
1 ms |
364 KB |
Output is correct |
7 |
Correct |
1 ms |
364 KB |
Output is correct |
8 |
Correct |
1 ms |
364 KB |
Output is correct |
9 |
Correct |
1 ms |
364 KB |
Output is correct |
10 |
Correct |
1 ms |
492 KB |
Output is correct |
11 |
Correct |
1 ms |
492 KB |
Output is correct |
12 |
Correct |
1 ms |
492 KB |
Output is correct |
13 |
Correct |
1 ms |
492 KB |
Output is correct |
14 |
Correct |
1 ms |
492 KB |
Output is correct |
15 |
Correct |
1 ms |
492 KB |
Output is correct |
16 |
Correct |
1 ms |
364 KB |
Output is correct |
17 |
Correct |
1 ms |
492 KB |
Output is correct |
18 |
Correct |
1 ms |
492 KB |
Output is correct |
19 |
Correct |
2 ms |
492 KB |
Output is correct |
20 |
Correct |
1 ms |
492 KB |
Output is correct |
21 |
Correct |
1 ms |
492 KB |
Output is correct |
22 |
Correct |
1 ms |
364 KB |
Output is correct |
23 |
Correct |
1 ms |
364 KB |
Output is correct |
24 |
Correct |
2 ms |
492 KB |
Output is correct |
25 |
Correct |
2 ms |
492 KB |
Output is correct |
26 |
Correct |
1 ms |
492 KB |
Output is correct |
27 |
Correct |
1 ms |
364 KB |
Output is correct |
28 |
Correct |
1 ms |
492 KB |
Output is correct |
29 |
Correct |
2 ms |
492 KB |
Output is correct |
30 |
Correct |
1 ms |
492 KB |
Output is correct |
31 |
Correct |
2 ms |
492 KB |
Output is correct |
32 |
Correct |
2 ms |
492 KB |
Output is correct |
33 |
Correct |
2 ms |
620 KB |
Output is correct |
34 |
Correct |
2 ms |
620 KB |
Output is correct |
35 |
Correct |
2 ms |
492 KB |
Output is correct |
36 |
Correct |
2 ms |
620 KB |
Output is correct |
37 |
Correct |
2 ms |
620 KB |
Output is correct |
38 |
Correct |
1 ms |
364 KB |
Output is correct |
39 |
Correct |
1 ms |
492 KB |
Output is correct |
40 |
Correct |
1 ms |
492 KB |
Output is correct |
41 |
Correct |
1 ms |
492 KB |
Output is correct |
42 |
Correct |
2 ms |
492 KB |
Output is correct |
43 |
Correct |
295 ms |
14572 KB |
Output is correct |
44 |
Correct |
1762 ms |
85996 KB |
Output is correct |
45 |
Correct |
554 ms |
35564 KB |
Output is correct |
46 |
Correct |
38 ms |
3948 KB |
Output is correct |
47 |
Correct |
1022 ms |
52972 KB |
Output is correct |
48 |
Correct |
2766 ms |
124140 KB |
Output is correct |
49 |
Correct |
568 ms |
31980 KB |
Output is correct |
50 |
Correct |
3291 ms |
162412 KB |
Output is correct |
51 |
Correct |
1095 ms |
61676 KB |
Output is correct |
52 |
Correct |
2685 ms |
143468 KB |
Output is correct |
53 |
Correct |
2492 ms |
125164 KB |
Output is correct |
54 |
Correct |
1 ms |
492 KB |
Output is correct |
55 |
Correct |
1300 ms |
64876 KB |
Output is correct |
56 |
Correct |
1523 ms |
87312 KB |
Output is correct |
57 |
Correct |
1236 ms |
72684 KB |
Output is correct |
58 |
Correct |
2265 ms |
122228 KB |
Output is correct |
59 |
Correct |
1252 ms |
73580 KB |
Output is correct |
60 |
Correct |
2942 ms |
188524 KB |
Output is correct |
61 |
Correct |
3052 ms |
188356 KB |
Output is correct |
62 |
Correct |
3601 ms |
187500 KB |
Output is correct |
63 |
Correct |
3421 ms |
188652 KB |
Output is correct |
64 |
Correct |
2964 ms |
188652 KB |
Output is correct |
65 |
Correct |
3634 ms |
188524 KB |
Output is correct |
66 |
Correct |
3471 ms |
189036 KB |
Output is correct |
67 |
Correct |
3000 ms |
189028 KB |
Output is correct |
68 |
Correct |
1392 ms |
86240 KB |
Output is correct |
69 |
Correct |
1248 ms |
83052 KB |
Output is correct |
70 |
Correct |
1781 ms |
106604 KB |
Output is correct |
71 |
Correct |
2061 ms |
128080 KB |
Output is correct |