#include <stdio.h>
#define N1 200000
#define N2 50000
#define LG 16 /* LG = ceil(log2(N2)) */
#define N_ (1 + N2 * (LG + 1))
#define INF 0x3f3f3f3f
int min(int a, int b) { return a < b ? a : b; }
int max(int a, int b) { return a > b ? a : b; }
int abs_(int a) { return a > 0 ? a : -a; }
unsigned int X = 12345;
int rand_() {
return (X *= 3) >> 1;
}
int *zz;
void sort(int *ii, int l, int r) {
while (l < r) {
int i = l, j = l, k = r, i_ = ii[l + rand_() % (r - l)], tmp;
while (j < k)
if (zz[ii[j]] == zz[i_])
j++;
else if (zz[ii[j]] < zz[i_]) {
tmp = ii[i], ii[i] = ii[j], ii[j] = tmp;
i++, j++;
} else {
k--;
tmp = ii[j], ii[j] = ii[k], ii[k] = tmp;
}
sort(ii, l, i);
l = k;
}
}
int ll[N_], rr[N_], kk[N_];
int update(int t, int l, int r, int i) {
static int _ = 1;
int t_ = _++;
kk[t_] = kk[t] + 1;
if (r - l > 1) {
int m = (l + r) / 2;
if (i < m)
ll[t_] = update(ll[t], l, m, i), rr[t_] = rr[t];
else
ll[t_] = ll[t], rr[t_] = update(rr[t], m, r, i);
}
return t_;
}
int query(int t, int l, int r, int ql, int qr) {
int m;
if (qr <= l || r <= ql || t == 0)
return 0;
if (ql <= l && r <= qr)
return kk[t];
m = (l + r) / 2;
return query(ll[t], l, m, ql, qr) + query(rr[t], m, r, ql, qr);
}
int xx[N1], yy[N1], xx_[N2], yy_[N2], ii[N2], jj[N2], tt[N2], n;
int t, i1, i2;
int check(int d_) {
if (t == 1) {
int i, xmn, xmx, ymn, ymx;
xmn = INF, xmx = -INF, ymn = INF, ymx = -INF;
for (i = 0; i < n; i++) {
xmn = min(xmn, xx[i]), xmx = max(xmx, xx[i]);
ymn = min(ymn, yy[i]), ymx = max(ymx, yy[i]);
}
for (i1 = 0; i1 < n; i1++)
if (xx[i1] - xmn <= d_ && xmx - xx[i1] <= d_ && yy[i1] - ymn <= d_ && ymx - yy[i1] <= d_)
return 1;
return 0;
} else {
static int iil[N2], iir[N2], iid[N2], iiu[N2], iil_[N2], iir_[N2], iid_[N2], iiu_[N2];
int i, i_, j, j_, l, r, d, u;
for (i = 0, j = 0; j < n; j++) {
while (i < n && xx[ii[j]] - xx[ii[i]] > d_)
i++;
iil[ii[j]] = i;
}
for (i = n - 1, j = n - 1; i >= 0; i--) {
while (j >= 0 && xx[ii[j]] - xx[ii[i]] > d_)
j--;
iir[ii[i]] = j;
}
for (i = 0, j = 0; j < n; j++) {
while (i < n && yy[jj[j]] - yy[jj[i]] > d_)
i++;
iid[jj[j]] = i;
}
for (i = n - 1, j = n - 1; i >= 0; i--) {
while (j >= 0 && yy[jj[j]] - yy[jj[i]] > d_)
j--;
iiu[jj[i]] = j;
}
for (i = 0; i < n; i++)
iil_[i] = 0, iir_[i] = n - 1, iid_[i] = 0, iiu_[i] = n - 1;
l = 0, r = n - 1, d = 0, u = n - 1;
for (i = 0, j = 0; j < n; j++) {
j_ = ii[j];
while (i < n && xx[j_] - xx[i_ = ii[i]] > d_) {
l = max(l, iil[i_]), r = min(r, iir[i_]);
d = max(d, iid[i_]), u = min(u, iiu[i_]);
i++;
}
iil_[j_] = max(iil_[j_], l), iir_[j_] = min(iir_[j_], r);
iid_[j_] = max(iid_[j_], d), iiu_[j_] = min(iiu_[j_], u);
}
l = 0, r = n - 1, d = 0, u = n - 1;
for (i = n - 1, j = n - 1; i >= 0; i--) {
i_ = ii[i];
while (j >= 0 && xx[j_ = ii[j]] - xx[i_] > d_) {
l = max(l, iil[j_]), r = min(r, iir[j_]);
d = max(d, iid[j_]), u = min(u, iiu[j_]);
j--;
}
iil_[i_] = max(iil_[i_], l), iir_[i_] = min(iir_[i_], r);
iid_[i_] = max(iid_[i_], d), iiu_[i_] = min(iiu_[i_], u);
}
l = 0, r = n - 1, d = 0, u = n - 1;
for (i = 0, j = 0; j < n; j++) {
j_ = jj[j];
while (i < n && yy[j_] - yy[i_ = jj[i]] > d_) {
l = max(l, iil[i_]), r = min(r, iir[i_]);
d = max(d, iid[i_]), u = min(u, iiu[i_]);
i++;
}
iil_[j_] = max(iil_[j_], l), iir_[j_] = min(iir_[j_], r);
iid_[j_] = max(iid_[j_], d), iiu_[j_] = min(iiu_[j_], u);
}
l = 0, r = n - 1, d = 0, u = n - 1;
for (i = n - 1, j = n - 1; i >= 0; i--) {
i_ = jj[i];
while (j >= 0 && yy[j_ = jj[j]] - yy[i_] > d_) {
l = max(l, iil[j_]), r = min(r, iir[j_]);
d = max(d, iid[j_]), u = min(u, iiu[j_]);
j--;
}
iil_[i_] = max(iil_[i_], l), iir_[i_] = min(iir_[i_], r);
iid_[i_] = max(iid_[i_], d), iiu_[i_] = min(iiu_[i_], u);
}
for (i1 = 0; i1 < n; i1++) {
int l = iil_[i1], r = iir_[i1], d = iid_[i1], u = iiu_[i1];
if (l <= r && d <= u && query(tt[r], 0, n, d, u + 1) - (l == 0 ? 0 : query(tt[l - 1], 0, n, d, u + 1) > 0)) {
for (i2 = 0; i2 < n; i2++)
if (xx_[i2] >= l && xx_[i2] <= r && yy_[i2] >= d && yy_[i2] <= u)
break;
return 1;
}
}
return 0;
}
}
int main() {
int i, j, lower, upper;
scanf("%d%d", &t, &n);
for (i = 0; i < n; i++) {
int x, y;
scanf("%d%d", &x, &y);
xx[i] = x + y, yy[i] = x - y;
}
if (t == 2) {
for (i = 0; i < n; i++)
ii[i] = i;
zz = xx, sort(ii, 0, n);
for (i = 0; i < n; i++)
xx_[ii[i]] = i;
for (j = 0; j < n; j++)
jj[j] = j;
zz = yy, sort(jj, 0, n);
for (j = 0; j < n; j++)
yy_[jj[j]] = j;
for (i = 0; i < n; i++)
tt[i] = update(i == 0 ? 0 : tt[i - 1], 0, n, yy_[ii[i]]);
}
lower = -1, upper = INF;
while (upper - lower > 1) {
int d = (lower + upper) / 2;
if (check(d))
upper = d;
else
lower = d;
}
check(upper);
printf("%d\n", upper);
if (t == 1)
printf("%d\n", i1 + 1);
else
printf("%d %d\n", i1 + 1, i2 + 1);
return 0;
}
Compilation message
kosta.c: In function 'main':
kosta.c:173:2: warning: ignoring return value of 'scanf' declared with attribute 'warn_unused_result' [-Wunused-result]
173 | scanf("%d%d", &t, &n);
| ^~~~~~~~~~~~~~~~~~~~~
kosta.c:177:3: warning: ignoring return value of 'scanf' declared with attribute 'warn_unused_result' [-Wunused-result]
177 | scanf("%d%d", &x, &y);
| ^~~~~~~~~~~~~~~~~~~~~
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
204 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
204 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
43 ms |
1152 KB |
Output is correct |
2 |
Correct |
42 ms |
1228 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
68 ms |
1836 KB |
Output is correct |
2 |
Correct |
71 ms |
1836 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
0 ms |
332 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1 ms |
332 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
12 ms |
976 KB |
Output is correct |
2 |
Incorrect |
10 ms |
972 KB |
Output isn't correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
11 ms |
972 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
172 ms |
10404 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Incorrect |
200 ms |
11832 KB |
Output isn't correct |
2 |
Halted |
0 ms |
0 KB |
- |