D. Matrix Cascade Harbour.Space Scholarship Contest 2023-2024 (Div. 1 + Div. 2)

Problem - D - Codeforces

题目大意：有一个n*n仅由0和1构成的方阵，每次操作可以选择一个位置(i,j)，然后将这个位置和所有满足x>i且x-i>|y-j|的位置(x,y)上的数都取反，问最少几次操作能将整个方阵变成全0

2<=n<=3000

思路：可以发现，取反的区间即为以(i,j)为顶点，两条斜率为1和-1的直线所包络的三角形，且只有上面的层会影响下面的层，所以我们的最优策略肯定是从上往下遍历每一行，一次统计必须要取反的位置，然后逐层维护状态

考虑如何维护状态，我们要记录的就是每个位置当前该是0还是1，那么可以用异或和来维护，从第一行开始，如果(i,j)的位置进行了操作，那么sum[i,j]要异或1，对于下一行，就从(i-1,j-1)、(i-1,j)、(i-1,j+1)转移，sum[i][j]^=sum[i-1][j-1]^sum[i-1][j]^sum[i-1][j+1]，对于再下一行，如果依然这么转移的话，在中间三个位置会出现问题，解决办法就是从(i-1,j-1)、(i-1,j+1)转移，然后再异或一下(i-2,j)来消去多于的部分，也就是sum[i][j]^=sum[i-1][j-1]^sum[i-1][j+1]^sum[i-2][j]，同时我们要用op[i][j]记录哪些位置被操作过，sum[i][j]要异或op[i-1][j]，注意最左和最右列，因为他们不会有重合部分所以是不需要异或sum[i-2][j]的

//#include<__msvc_all_public_headers.hpp>
#include<bits/stdc++.h>
using namespace std;
const int N = 3e3 + 5;
typedef long long ll;
const int INF = 0x7fffffff;
int sum[N][N];
string s[N];
int op[N][N];
int n;
void init()
{
	for (int i = 0; i < n; i++)
	{
		for (int j = 0; j < n; j++)
		{
			sum[i][j] = op[i][j] = 0;
		}
	}
}
void solve()
{
	cin >> n;
	init();
	for (int i = 0; i < n; i++)
	{
		cin >> s[i];
	}
	int ans = 0;
	for (int i = 0; i < n; i++)
	{
		for (int j = 0; j < n; j++)
		{
			if (i == 0)
			{//第一行不需要转移

			}
			else if (i == 1)
			{
				if (j > 0)
				{
					sum[i][j] ^= sum[i - 1][j - 1];//从左上转移
				}
				sum[i][j] ^= sum[i - 1][j];//纵向转移
				if (j < n - 1)
				{
					sum[i][j] ^= sum[i - 1][j + 1];//从右上转移
				}
			}
			else
			{
				if (j > 0)
				{
					sum[i][j] ^= sum[i - 1][j - 1];/从左上转移
				}
				else
				{
					sum[i][j] ^= sum[i - 2][j];//左右两列不需要纵向转移
				}
				if (j < n - 1)
				{
					sum[i][j] ^= sum[i - 1][j + 1];
				}
				else
				{
					sum[i][j] ^= sum[i - 2][j];
				}
				sum[i][j] ^= sum[i - 2][j];//去掉多于的部分
				sum[i][j] ^= op[i - 1][j];//上一行操作过要异或一下，就像第二行那样
			}
			if (sum[i][j] ^ (s[i][j] - '0'))
			{//当前位置是1就必须要操作这里
				ans++;
				sum[i][j] ^= 1;//异或和再额外异或1
				op[i][j] = 1;//标记操作过的位置
			}
		}
	}
	cout << ans << endl;
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	int t;
	cin >> t;
	while (t--)
	{
		solve();
	}
	return 0;
}