QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#879130 | #9973. 魔法少女网站第二部 | zlt | WA | 2746ms | 223620kb | C++14 | 7.4kb | 2025-02-01 21:28:46 | 2025-02-01 21:28:46 |
Judging History
answer
// Problem: P11367 [Ynoi2024] 魔法少女网站第二部
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P11367
// Memory Limit: 512 MB
// Time Limit: 3000 ms
//
// Powered by CP Editor (https://cpeditor.org)
#include <bits/stdc++.h>
#define pb emplace_back
#define fst first
#define scd second
#define mkp make_pair
#define uint unsigned
#define mems(a, x) memset((a), (x), sizeof(a))
using namespace std;
typedef long long ll;
typedef double db;
typedef unsigned long long ull;
typedef long double ldb;
typedef pair<int, int> pii;
namespace IO {
const int maxn = 1 << 20;
char ibuf[maxn], *iS, *iT, obuf[maxn], *oS = obuf;
inline char gc() {
return (iS == iT ? iT = (iS = ibuf) + fread(ibuf, 1, maxn, stdin), (iS == iT ? EOF : *iS++) : *iS++);
}
template<typename T = int>
inline T read() {
char c = gc();
T x = 0;
bool f = 0;
while (c < '0' || c > '9') {
f |= (c == '-');
c = gc();
}
while (c >= '0' && c <= '9') {
x = (x << 1) + (x << 3) + (c ^ 48);
c = gc();
}
return f ? ~(x - 1) : x;
}
inline void flush() {
fwrite(obuf, 1, oS - obuf, stdout);
oS = obuf;
}
struct Flusher {
~Flusher() {
flush();
}
} AutoFlush;
inline void pc(char ch) {
if (oS == obuf + maxn) {
flush();
}
*oS++ = ch;
}
template<typename T>
inline void write(T x) {
static char stk[64], *tp = stk;
if (x < 0) {
x = ~(x - 1);
pc('-');
}
do {
*tp++ = x % 10;
x /= 10;
} while (x);
while (tp != stk) {
pc((*--tp) | 48);
}
}
template<typename T>
inline void writesp(T x) {
write(x);
pc(' ');
}
template<typename T>
inline void writeln(T x) {
write(x);
pc('\n');
}
}
using IO::read;
using IO::write;
using IO::pc;
using IO::writesp;
using IO::writeln;
const int maxn = 2000100;
const int logn = 22;
const int inf = 0x3f3f3f3f;
int n, m, a[maxn], b[maxn], f[logn][maxn], d[maxn];
pii p[maxn], q[maxn];
struct node {
int l, r;
} c[maxn];
ll ans[maxn];
inline int qmin(int l, int r) {
if (l > r) {
return inf;
}
int k = __lg(r - l + 1);
return min(f[k][l], f[k][r - (1 << k) + 1]);
}
inline int qmax(int l, int r) {
if (l > r) {
return 0;
}
int k = __lg(r - l + 1);
return max(f[k][l], f[k][r - (1 << k) + 1]);
}
struct List {
int hd[maxn * 3], len, val[maxn << 1], nxt[maxn << 1];
inline void add(int x, int y) {
val[++len] = y;
nxt[len] = hd[x];
hd[x] = len;
}
} G, T;
namespace BIT {
ll c[maxn];
inline void update(int x, int d) {
for (int i = x; i <= n; i += (i & (-i))) {
c[i] += d;
}
}
inline ll query(int x) {
ll res = 0;
for (int i = x; i; i -= (i & (-i))) {
res += c[i];
}
return res;
}
inline ll query(int l, int r) {
return query(r) - query(l - 1);
}
inline void clear(int x) {
for (int i = x; i <= n; i += (i & (-i))) {
c[i] = 0;
}
}
}
int fa[maxn];
int find(int x) {
return fa[x] == x ? x : fa[x] = find(fa[x]);
}
void update(int rt, int l, int r, int ql, int qr, int i) {
int mid = (l + r) >> 1;
if (qr <= mid) {
update(rt << 1, l, mid, ql, qr, i);
} else if (ql > mid) {
update(rt << 1 | 1, mid + 1, r, ql, qr, i);
} else {
T.add(rt, i);
}
}
void dfs(int rt, int l, int r) {
if (l == r) {
return;
}
int mid = (l + r) >> 1;
dfs(rt << 1, l, mid);
dfs(rt << 1 | 1, mid + 1, r);
inplace_merge(p + l, p + mid + 1, p + r + 1);
G.len = 0;
for (int i = l; i <= r; ++i) {
G.hd[i] = 0;
}
for (int i = l; i <= r; ++i) {
d[p[i].scd] = i - l + 1;
}
for (int i = l; i <= r; ++i) {
b[d[i]] = i;
}
for (int i = T.hd[rt]; i; i = T.nxt[i]) {
int j = T.val[i];
G.add(c[j].l, j);
G.add(c[j].r, j);
}
ll s = 0;
int mn = inf, mx = 0, tot = r - l + 1;
for (int i = 1; i <= tot; ++i) {
f[0][i] = (b[i] > mid ? b[i] : inf);
}
for (int j = 1; (1 << j) <= tot; ++j) {
for (int i = 1; i + (1 << j) - 1 <= tot; ++i) {
f[j][i] = min(f[j - 1][i], f[j - 1][i + (1 << (j - 1))]);
}
}
for (int i = 1; i <= tot; ++i) {
fa[i] = (b[i] <= mid ? i : i - 1);
}
fa[0] = 0;
for (int i = l; i <= mid; ++i) {
fa[d[i]] = q[i].fst = find(d[i] - 1);
}
for (int i = 1; i <= tot; ++i) {
fa[i] = (b[i] <= mid ? i + 1 : i);
}
fa[tot + 1] = tot + 1;
for (int i = l; i <= mid; ++i) {
fa[d[i]] = q[i].scd = find(d[i] + 1);
}
for (int i = mid; i >= l; --i) {
mn = min(mn, d[i]);
mx = max(mx, d[i]);
int j = q[i].fst, k = q[i].scd, p = inf;
if (j) {
s += abs(i - b[j]);
int t = qmin(j + 1, d[i] - 1);
p = min(p, t);
if (t <= r) {
BIT::update(t, (mid - max(i, b[j])) * 2);
}
}
if (k <= tot) {
s += abs(i - b[k]);
int t = qmin(d[i] + 1, k - 1);
p = min(p, t);
if (t <= r) {
BIT::update(t, (mid - max(i, b[k])) * 2);
}
}
if (j && k <= tot) {
s -= abs(b[j] - b[k]);
if (p <= r) {
BIT::update(p, -(mid - max(b[j], b[k])) * 2);
}
}
for (int _ = G.hd[i]; _; _ = G.nxt[_]) {
int j = G.val[_];
ans[j] += s + BIT::query(mid + 1, c[j].r);
if (qmin(mx + 1, tot) <= c[j].r) {
ans[j] += mid - b[mx];
}
if (qmin(1, mn - 1) <= c[j].r) {
ans[j] += mid - b[mn];
}
}
}
for (int i = mid + 1; i <= r; ++i) {
BIT::clear(i);
}
s = mx = 0;
mn = inf;
for (int i = 1; i <= tot; ++i) {
f[0][i] = (b[i] <= mid ? b[i] : 0);
}
for (int j = 1; (1 << j) <= tot; ++j) {
for (int i = 1; i + (1 << j) - 1 <= tot; ++i) {
f[j][i] = max(f[j - 1][i], f[j - 1][i + (1 << (j - 1))]);
}
}
for (int i = 1; i <= tot; ++i) {
fa[i] = (b[i] > mid ? i : i - 1);
}
fa[0] = 0;
for (int i = r; i > mid; --i) {
fa[d[i]] = q[i].fst = find(d[i] - 1);
}
for (int i = 1; i <= tot; ++i) {
fa[i] = (b[i] > mid ? i : i + 1);
}
fa[tot + 1] = tot + 1;
for (int i = r; i > mid; --i) {
fa[d[i]] = q[i].scd = find(d[i] + 1);
}
for (int i = mid + 1; i <= r; ++i) {
mn = min(mn, d[i]);
mx = max(mx, d[i]);
int j = q[i].fst, k = q[i].scd, p = 0;
if (j) {
s += abs(i - b[j]);
int t = qmax(j + 1, d[i] - 1);
p = max(p, t);
if (t >= l) {
BIT::update(t, (min(i, b[j]) - mid) * 2);
}
}
if (k <= tot) {
s += abs(i - b[k]);
int t = qmax(d[i] + 1, k - 1);
p = max(p, t);
if (t >= l) {
BIT::update(t, (min(i, b[k]) - mid) * 2);
}
}
if (j && k <= tot) {
s -= abs(b[j] - b[k]);
if (p >= l) {
BIT::update(p, -(min(b[j], b[k]) - mid) * 2);
}
}
for (int _ = G.hd[i]; _; _ = G.nxt[_]) {
int j = G.val[_];
ans[j] += s + BIT::query(c[j].l, mid);
if (qmax(mx + 1, tot) >= c[j].l) {
ans[j] += b[mx] - mid;
}
if (qmax(1, mn - 1) >= c[j].l) {
ans[j] += b[mn] - mid;
}
}
}
for (int i = l; i <= mid; ++i) {
BIT::clear(i);
}
}
void solve() {
mt19937 rnd(chrono::steady_clock::now().time_since_epoch().count());
n = m = 1000000;
for (int i = 1; i <= n; ++i) {
a[i] = i;
}
shuffle(a + 1, a + n + 1, rnd);
for (int i = 1; i <= n; ++i) {
p[i] = pii(a[i], i);
}
for (int i = 1; i <= m; ++i) {
c[i].l = rnd() % n + 1;
c[i].r = rnd() % n + 1;
if (c[i].l > c[i].r) {
swap(c[i].l, c[i].r);
}
if (c[i].l != c[i].r) {
update(1, 1, n, c[i].l, c[i].r, i);
}
}
dfs(1, 1, n);
}
int main() {
int T = 1;
// scanf("%d", &T);
while (T--) {
solve();
}
return 0;
}
详细
Test #1:
score: 0
Wrong Answer
time: 2746ms
memory: 223620kb
input:
1999999 1999998 8 9 6 4 1 7 5 13 3 16 2 20 14 10 12 17 21 11 27 29 19 23 15 18 24 26 22 31 28 37 36 25 38 41 34 43 35 30 33 45 51 32 52 42 50 44 39 40 58 56 48 49 46 59 47 63 60 57 54 55 66 53 65 62 67 61 72 75 68 69 71 79 77 64 85 82 74 73 70 87 76 81 78 88 91 90 84 97 80 89 93 99 83 98 96 86 95 10...
output:
result:
wrong answer Answer contains longer sequence [length = 1999998], but output contains 0 elements