[NOIP 2024 模擬2]矩陣學說

maniubi發表於2024-09-12

[NOIP 2024 模擬2]矩陣學說

題意

給出 \(n\)\(m\) 列的矩陣,第 \(i\) 行第 \(j\) 列的元素為 \(a_{i,j}\),找出滿足以下條件的三元組 \((i, j, x)\) 的 數量:

  1. \(1 ≤ i ≤ n\), \(1 ≤ j \le m\), \(1 ≤ x ≤ \min(n − i + 1, m − j + 1)\)
  2. 矩陣的左上角 \((i, j)\) 到右下角 \((i + x − 1, j + x − 1)\) 恰好含 \(k\) 個不同的整數。

\(1\le a_{i,j} \le 100\) \(1\le n,m \le 2000\)

思路

考慮求出包含數的種類 \(\le k\) 的個數和 \(< k\) 的個數,相減得到答案。

列舉左上角,矩陣包含數的種類個數顯然單調不遞減,二分邊長即可。

快速求區間數的種類可以用二維 ST 表加 bitset。

時間複雜度:\(O(nm\log \min(n,m))\)

程式碼

#include <bits/stdc++.h>
#define ll long long
using namespace std;
const int N = 1505;
int n, m, k, a[N][N], lg[N];
bitset <100> st[N][N][12];
void init() {
	lg[1] = 0;
	for (int i = 2; i <= n; i ++) lg[i] = lg[i >> 1] + 1; 
	for (int i = 1; i <= n; i ++)
		for (int j = 1; j <= m; j ++)
			st[i][j][0][a[i][j] - 1] = 1;
	for (int l = 1; l <= 14; l ++) 
		for (int i = 1; i <= n && i + (1 << l) - 1 <= n; i ++) 
			for (int j = 1; j <= m && j + (1 << l) - 1 <= m; j ++) 
				st[i][j][l] = st[i][j][l - 1] | st[i + (1 << (l - 1))][j][l - 1] | st[i][j + (1 << (l - 1))][l - 1] | st[i + (1 << (l - 1))][j + (1 << (l - 1))][l - 1];
}
int query(int x, int y, int l) {
	int K = lg[l];
	int X = x + l - 1, Y = y + l - 1;
	return (st[x][y][K] | st[X - (1 << K) + 1][Y - (1 << K) + 1][K] | st[X - (1 << K) + 1][y][K] | st[x][Y - (1 << K) + 1][K]).count(); 
	return 0;
}
int calc(int num) {
	int res = 0;
	for (int i = 1; i <= n; i ++) {
		for (int j = 1; j <= m; j ++) {
			int l = 1, r = min(n - i + 1, m - j + 1), pos = 0;
			while (l <= r) {
				int mid = (l + r) >> 1;
				if (query(i, j, mid) <= num) pos = mid, l = mid + 1;
				else r = mid - 1;
			}
			res += pos;
		}
	}
	return res;
}
void solve() {
	cin >> n >> m >> k;
	for (int i = 1; i <= n; i ++) 
		for (int j = 1; j <= m; j ++) 
			cin >> a[i][j];
	init();
	cout << calc(k) - calc(k - 1) << "\n";
}
signed main() {
	freopen("mar.in", "r", stdin);
	freopen("mar.out", "w", stdout);
	int Case = 1;
//	cin >> Case;
	while (Case --)
		solve();
	return 0;
}
/*
2 3 4
1 2 3
4 5 6
*/