#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
#define fi first
#define se second
#define mp make_pair
#define fastIO ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
const int N = 2010;
const ll C = (ll)1e18;
struct frac{
ll A;
ll B;
frac operator+ (frac y){
return {A * y.B + B * y.A, B * y.B};
}
frac operator- (frac y){
return {A * y.B - B * y.A, B * y.B};
}
bool operator< (frac y){
if(A/B != y.A/y.B){
return A/B < y.A/y.B;
}
ll D = A/B;
ll temp0 = A - D * B;
ll temp1 = y.A - D * y.B;
return (temp0 * y.B - temp1 * B < 0);
}
};
ll V[N][N];
ll S[N][N];
int m;
frac calc(int id, frac P){
int las = P.A / P.B;
frac Y = {S[id][las], 1};
frac add = (P - (frac){las, 1});
add.A *= V[id][las + 1];
Y = Y + add;
return Y;
}
frac need[N];
bool use[N];
int pi[N];
int main(){
fastIO;
//freopen("in.txt", "r", stdin);
int n;
cin >> n >> m;
int vv;
for(int i = 1; i <= n; i ++ ){
need[i] = {0, n};
for(int j = 1; j <= m ; j ++ ){
cin >> vv;
V[i][j] = vv;
need[i].A += V[i][j];
S[i][j] = V[i][j] + S[i][j - 1];
}
pi[i] = -1;
}
frac las = {0, 1};
frac tis;
frac nx;
frac cur;
frac low;
int idx;
frac SS;
ll X;
int pos;
for(int it = 1; it < n; it ++ ){
low = {m, 1};
idx = -1;
for(int i = 1; i <= n; i ++ ){
if(!use[i]){
tis = calc(i, las);
int j = 0;
int pp;
for(int ln = 13; ln >= 0; ln -- ){
pp = (j + (1 << ln));
if(pp < m){
nx = (frac){S[i][pp], 1ll};
if(tis < nx) {
SS = nx - tis;
SS = SS - need[i];
if(SS.A < 0) j = pp;
}
else{
j = pp;
}
}
}
j ++ ;
nx = (frac){S[i][j], 1ll};
SS = nx - tis;
SS = SS - need[i];
SS.B /= n;
X = SS.A / SS.B;
tis = (frac){j, 1} - (frac){X, n * 1ll * V[i][j]};
if(tis < low){
low = tis;
idx = i;
}
}
}
pi[it] = idx;
use[idx]=true;
las = low;
cout << (long long)low.A << " " << (long long)low.B << "\n";
}
for(int i = 1; i <= n; i ++ ){
if(!use[i]){
pi[n]=i;
}
}
//return 0;
for(int i = 1; i <= n; i ++ ){
cout << pi[i] << " ";
}
cout << "\n";
return 0;
}
Compilation message
naan.cpp: In function 'int main()':
naan.cpp:72:10: warning: unused variable 'cur' [-Wunused-variable]
72 | frac cur;
| ^~~
naan.cpp:77:9: warning: unused variable 'pos' [-Wunused-variable]
77 | int pos;
| ^~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
340 KB |
Output is correct |
6 |
Correct |
0 ms |
340 KB |
Output is correct |
7 |
Correct |
0 ms |
340 KB |
Output is correct |
8 |
Correct |
1 ms |
340 KB |
Output is correct |
9 |
Correct |
1 ms |
340 KB |
Output is correct |
10 |
Correct |
1 ms |
340 KB |
Output is correct |
11 |
Correct |
1 ms |
340 KB |
Output is correct |
12 |
Correct |
1 ms |
340 KB |
Output is correct |
13 |
Correct |
1 ms |
340 KB |
Output is correct |
14 |
Correct |
1 ms |
340 KB |
Output is correct |
15 |
Correct |
1 ms |
340 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
348 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
468 KB |
Output is correct |
4 |
Correct |
1 ms |
468 KB |
Output is correct |
5 |
Correct |
1 ms |
468 KB |
Output is correct |
6 |
Correct |
1 ms |
340 KB |
Output is correct |
7 |
Correct |
1 ms |
340 KB |
Output is correct |
8 |
Correct |
1 ms |
340 KB |
Output is correct |
9 |
Correct |
1 ms |
468 KB |
Output is correct |
10 |
Correct |
1 ms |
468 KB |
Output is correct |
11 |
Correct |
1 ms |
468 KB |
Output is correct |
12 |
Correct |
0 ms |
340 KB |
Output is correct |
13 |
Correct |
1 ms |
352 KB |
Output is correct |
14 |
Correct |
1 ms |
468 KB |
Output is correct |
15 |
Correct |
1 ms |
468 KB |
Output is correct |
16 |
Correct |
1 ms |
468 KB |
Output is correct |
17 |
Correct |
1 ms |
468 KB |
Output is correct |
18 |
Correct |
1 ms |
468 KB |
Output is correct |
19 |
Correct |
1 ms |
468 KB |
Output is correct |
20 |
Correct |
1 ms |
468 KB |
Output is correct |
21 |
Correct |
1 ms |
468 KB |
Output is correct |
22 |
Correct |
1 ms |
468 KB |
Output is correct |
23 |
Correct |
0 ms |
340 KB |
Output is correct |
24 |
Correct |
1 ms |
468 KB |
Output is correct |
25 |
Correct |
1 ms |
468 KB |
Output is correct |
26 |
Correct |
1 ms |
340 KB |
Output is correct |
27 |
Correct |
1 ms |
468 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
340 KB |
Output is correct |
2 |
Correct |
1 ms |
340 KB |
Output is correct |
3 |
Correct |
1 ms |
340 KB |
Output is correct |
4 |
Correct |
1 ms |
340 KB |
Output is correct |
5 |
Correct |
1 ms |
340 KB |
Output is correct |
6 |
Correct |
0 ms |
340 KB |
Output is correct |
7 |
Correct |
0 ms |
340 KB |
Output is correct |
8 |
Correct |
1 ms |
340 KB |
Output is correct |
9 |
Correct |
1 ms |
340 KB |
Output is correct |
10 |
Correct |
1 ms |
340 KB |
Output is correct |
11 |
Correct |
1 ms |
340 KB |
Output is correct |
12 |
Correct |
1 ms |
340 KB |
Output is correct |
13 |
Correct |
1 ms |
340 KB |
Output is correct |
14 |
Correct |
1 ms |
340 KB |
Output is correct |
15 |
Correct |
1 ms |
340 KB |
Output is correct |
16 |
Correct |
1 ms |
348 KB |
Output is correct |
17 |
Correct |
1 ms |
340 KB |
Output is correct |
18 |
Correct |
1 ms |
468 KB |
Output is correct |
19 |
Correct |
1 ms |
468 KB |
Output is correct |
20 |
Correct |
1 ms |
468 KB |
Output is correct |
21 |
Correct |
1 ms |
340 KB |
Output is correct |
22 |
Correct |
1 ms |
340 KB |
Output is correct |
23 |
Correct |
1 ms |
340 KB |
Output is correct |
24 |
Correct |
1 ms |
468 KB |
Output is correct |
25 |
Correct |
1 ms |
468 KB |
Output is correct |
26 |
Correct |
1 ms |
468 KB |
Output is correct |
27 |
Correct |
0 ms |
340 KB |
Output is correct |
28 |
Correct |
1 ms |
352 KB |
Output is correct |
29 |
Correct |
1 ms |
468 KB |
Output is correct |
30 |
Correct |
1 ms |
468 KB |
Output is correct |
31 |
Correct |
1 ms |
468 KB |
Output is correct |
32 |
Correct |
1 ms |
468 KB |
Output is correct |
33 |
Correct |
1 ms |
468 KB |
Output is correct |
34 |
Correct |
1 ms |
468 KB |
Output is correct |
35 |
Correct |
1 ms |
468 KB |
Output is correct |
36 |
Correct |
1 ms |
468 KB |
Output is correct |
37 |
Correct |
1 ms |
468 KB |
Output is correct |
38 |
Correct |
0 ms |
340 KB |
Output is correct |
39 |
Correct |
1 ms |
468 KB |
Output is correct |
40 |
Correct |
1 ms |
468 KB |
Output is correct |
41 |
Correct |
1 ms |
340 KB |
Output is correct |
42 |
Correct |
1 ms |
468 KB |
Output is correct |
43 |
Correct |
41 ms |
8876 KB |
Output is correct |
44 |
Correct |
273 ms |
38176 KB |
Output is correct |
45 |
Correct |
128 ms |
22028 KB |
Output is correct |
46 |
Correct |
15 ms |
4052 KB |
Output is correct |
47 |
Correct |
180 ms |
30180 KB |
Output is correct |
48 |
Correct |
275 ms |
28908 KB |
Output is correct |
49 |
Correct |
78 ms |
15208 KB |
Output is correct |
50 |
Correct |
469 ms |
56132 KB |
Output is correct |
51 |
Correct |
179 ms |
33156 KB |
Output is correct |
52 |
Correct |
410 ms |
54272 KB |
Output is correct |
53 |
Correct |
359 ms |
49376 KB |
Output is correct |
54 |
Correct |
1 ms |
468 KB |
Output is correct |
55 |
Correct |
120 ms |
14052 KB |
Output is correct |
56 |
Correct |
223 ms |
41036 KB |
Output is correct |
57 |
Correct |
190 ms |
36380 KB |
Output is correct |
58 |
Correct |
301 ms |
44956 KB |
Output is correct |
59 |
Correct |
198 ms |
35504 KB |
Output is correct |
60 |
Correct |
593 ms |
63420 KB |
Output is correct |
61 |
Correct |
605 ms |
63232 KB |
Output is correct |
62 |
Correct |
590 ms |
63080 KB |
Output is correct |
63 |
Correct |
631 ms |
63392 KB |
Output is correct |
64 |
Correct |
598 ms |
63304 KB |
Output is correct |
65 |
Correct |
493 ms |
63292 KB |
Output is correct |
66 |
Correct |
508 ms |
63324 KB |
Output is correct |
67 |
Correct |
484 ms |
63308 KB |
Output is correct |
68 |
Correct |
239 ms |
37196 KB |
Output is correct |
69 |
Correct |
237 ms |
39396 KB |
Output is correct |
70 |
Correct |
286 ms |
40884 KB |
Output is correct |
71 |
Correct |
362 ms |
50412 KB |
Output is correct |