QOJ.ac

QOJ

IDProblemSubmitterResultTimeMemoryLanguageFile sizeSubmit timeJudge time
#878410#9971. 新本格魔法少女りすかzlt0 12831ms25360kbC++144.9kb2025-02-01 15:11:592025-02-01 15:12:00

Judging History

This is the latest submission verdict.

  • [2025-02-01 15:12:00]
  • Judged
  • Verdict: 0
  • Time: 12831ms
  • Memory: 25360kb
  • [2025-02-01 15:11:59]
  • Submitted

answer

// Problem: P11365 [Ynoi2024] 新本格魔法少女りすか
// Contest: Luogu
// URL: https://www.luogu.com.cn/problem/P11365
// Memory Limit: 128 MB
// Time Limit: 15000 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 = 500100;
const int B = 400;

int n, m, K, a[maxn], bel[maxn], L[maxn], R[maxn], b[maxn], c[maxn];
ll ans[maxn];
pii d[maxn], e[maxn];

namespace BIT {
	int c[maxn];
	ull b[maxn];
	
	inline void update(int x, int d) {
		b[x >> 6] |= (1ULL << (x & 63));
		for (int i = (x >> 6) + 1; i <= (n >> 6) + 1; i += (i & (-i))) {
			c[i] += d;
		}
	}
	
	inline int query(int x) {
		int res = __builtin_popcountll(b[x >> 6] << (63 - (x & 63)));
		for (int i = (x >> 6); i; i -= (i & (-i))) {
			res += c[i];
		}
		return res;
	}
	
	inline void clear(int x) {
		b[x >> 6] = 0;
		for (int i = (x >> 6) + 1; i <= (n >> 6) + 1; i += (i & (-i))) {
			int &x = c[i];
			if (x) {
				x = 0;
			} else {
				break;
			}
		}
	}
}

void solve() {
	n = read();
	m = read();
	for (int i = 1; i <= n; ++i) {
		a[i] = read();
		bel[i] = (i - 1) / B + 1;
		if (!L[bel[i]]) {
			L[bel[i]] = i;
		}
		R[bel[i]] = i;
	}
	for (int i = 1, k, l, r; i <= m; ++i) {
		k = read();
		e[i].fst = K + 1;
		while (k--) {
			l = read();
			r = read();
			d[++K] = mkp(l, r);
			if (bel[l] == bel[r]) {
				for (int j = l; j <= r; ++j) {
					ans[i] += BIT::query(a[j]);
				}
				for (int j = l; j <= r; ++j) {
					BIT::update(a[j], 1);
				}
			} else {
				for (int j = l; j <= R[bel[l]]; ++j) {
					ans[i] += BIT::query(a[j]);
				}
				for (int j = L[bel[r]]; j <= r; ++j) {
					ans[i] += BIT::query(a[j]);
				}
				for (int j = l; j <= R[bel[l]]; ++j) {
					BIT::update(a[j], 1);
				}
				for (int j = L[bel[r]]; j <= r; ++j) {
					BIT::update(a[j], 1);
				}
			}
		}
		e[i].scd = K;
		for (int j = e[i].fst; j <= e[i].scd; ++j) {
			pii p = d[j];
			int l = p.fst, r = p.scd;
			if (bel[l] == bel[r]) {
				for (int j = l; j <= r; ++j) {
					BIT::clear(a[j]);
				}
			} else {
				for (int j = l; j <= R[bel[l]]; ++j) {
					BIT::clear(a[j]);
				}
				for (int j = L[bel[r]]; j <= r; ++j) {
					BIT::clear(a[j]);
				}
			}
		}
	}
	for (int k = 1; k <= bel[n]; ++k) {
		mems(b, 0);
		for (int i = L[k]; i <= R[k]; ++i) {
			++b[a[i]];
		}
		for (int i = 1; i <= n; ++i) {
			b[i] += b[i - 1];
		}
		for (int i = 1; i <= n; ++i) {
			if (i < L[k]) {
				c[i] = c[i - 1] + R[k] - L[k] + 1 - b[a[i]];
			} else if (i <= R[k]) {
				c[i] = c[i - 1];
			} else {
				c[i] = c[i - 1] + b[a[i] - 1];
			}
		}
		for (int i = 1; i <= m; ++i) {
			bool fl = 0;
			ll s = 0;
			for (int j = e[i].fst; j <= e[i].scd; ++j) {
				pii p = d[j];
				int l = p.fst, r = p.scd;
				int bl = (l - 1) / B + 1, br = (r - 1) / B + 1;
				if (bl < k && k < br) {
					fl = 1;
				} else if (r < L[k]) {
					s += c[r] - c[l - 1];
				} else if (l > R[k]) {
					if (!fl) {
						break;
					}
					if (bl == br) {
						s += c[r] - c[l - 1];
					} else {
						s += c[min(bl * B, n)] - c[l - 1] + c[r] - c[br * (B - 1)];
					}
				}
			}
			if (fl) {
				ans[i] += s;
			}
		}
	}
	for (int i = 1; i <= m; ++i) {
		writeln(ans[i]);
	}
}

int main() {
	int T = 1;
	// scanf("%d", &T);
	while (T--) {
		solve();
	}
	return 0;
}

Details

Tip: Click on the bar to expand more detailed information

Subtask #1:

score: 0
Wrong Answer

Test #1:

score: 0
Wrong Answer
time: 12831ms
memory: 22560kb

input:

500000 50174
453952 363686 491616 32825 57465 422945 471561 73291 421205 416554 23295 133266 242225 330448 25039 340064 28713 465664 162059 323880 28978 273728 101338 161413 294941 214275 63951 267981 294251 202822 253124 272510 3268 37918 139343 385983 111577 311901 487665 482827 347449 416029 3065...

output:

37714946698
35224187569
44561035410
34502447258
44351966662
29641299969
42679630560
29616881397
34903331914
31440038101
45691749710
43176054605
25150426865
48392970986
40301032979
32344912915
32513693783
32128416760
45175935038
46804307531
35315057819
31968319960
34254687251
30951250043
44206067265
...

result:

wrong answer 1st numbers differ - expected: '37140482224', found: '37714946698'

Subtask #2:

score: 0
Wrong Answer

Test #6:

score: 10
Accepted
time: 2263ms
memory: 25360kb

input:

500000 5
157360 289139 98034 293691 150262 268366 36782 147093 365410 444658 343224 375392 278298 357620 389673 167019 7747 119244 102126 83512 3649 459230 197365 245259 38071 249539 34617 213697 292553 389625 395778 280152 280038 239519 301475 232272 145919 370004 422791 271143 488283 185166 351026...

output:

50666226791
50440151159
50719090228
50631079493
50747050985

result:

ok 5 number(s): "50666226791 50440151159 50719090228 50631079493 50747050985"

Test #7:

score: 0
Wrong Answer
time: 2228ms
memory: 21548kb

input:

500000 5
70752 248917 41910 13653 177839 45858 54229 174152 10090 332231 437916 391622 432270 53875 446421 91921 461870 243336 363086 249844 338371 495447 423857 350363 365324 255747 170681 435791 177298 199582 240768 449011 302158 378518 233809 267343 362784 187864 114604 322031 255693 54447 406202...

output:

56270227296
59467754167
55928542923
57169004523
62373267184

result:

wrong answer 1st numbers differ - expected: '48041849427', found: '56270227296'

Subtask #3:

score: 0
Skipped

Dependency #1:

0%