QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#879130#9973. 魔法少女网站第二部zltWA 2746ms223620kbC++147.4kb2025-02-01 21:28:462025-02-01 21:28:46

Judging History

This is the latest submission verdict.

  • [2025-02-01 21:28:46]
  • Judged
  • Verdict: WA
  • Time: 2746ms
  • Memory: 223620kb
  • [2025-02-01 21:28:46]
  • Submitted

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;
}

Details

Tip: Click on the bar to expand more detailed information

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