[ABC276Ex] Construct a Matrix

没有题解，所以来写一篇。

Description

构造一个 $N\times N$ 的矩阵 $A$，其中 $A_{i,j}\in {0,1,2}$，要求同时满足 $Q$ 条限制。
每条限制形如：给定 $a,b,c,d,e$，要求 $A$ 满足 $\prod\limits_{i=a}^b\prod\limits_{j=c}^d A_{i,j} \equiv e\pmod 3$。

Solution

为贴合原题，设 $(a,b,c,d)$ 表示左上角为 $(a,c)$，右下角为 $(b,d)$ 的矩形。

根据连乘的性质，若 $e_i$ 不为 $0$，则 $(a_i,b_i,c_i,d_i)$ 中任意 $A_{i,j}$ 不能为 $0$。
而对于 $e_i=0$，要求 $(a,b,c,d)$ 中至少一个是 $0$。那么对于不被 $e_i\in\{1,2\}$ 限制的点，显然强制它们都为 $0$ 一定不劣。如果这种方式依然不能满足所有 $e_i=0$ 的限制，则无解。

现在我们已经确定哪些位置为 $0$，只需要考虑不为 $0$ 的位置与 $e_i\in\{1,2\}$ 的限制。而 $A_{i,j}=1$ 对矩阵积同样没有影响，矩阵积模 $3$ 的值只与 $A_{i,j}=2$ 的数量有关。

简单找规律（？）可以发现，$2$ 的奇数次幂 $\bmod 3$ 的值为 $2$，偶数次幂 $\bmod 3$ 的值为 $1$。
考虑把每个 $A_{i,j}$ 设为未知数解方程。
令 $x_{i,j}$ 的值 $0,1$ 对应 $A_{i,j}$ 的 ${1,2}$，那么会得到 $Q$ 组方程，形如：

$$x_{a,c}\oplus x_{a,c+1}\oplus\dots\oplus x_{b,d}=[e_i=1]$$

未知数个数为 $n^2$，直接高斯消元复杂度不能接受。
重定义 $x_{a,b}=\sum\limits_{i=1}^a\sum\limits_{j=1}^b A_{i,j} \pmod 2$，则每个方程可以写成：

$$x_{b,d}\oplus x_{a-1,c-1}\oplus x_{a-1,d}\oplus x_{b,c-1}=[e_i=1]$$

这样每次询问最多会产生 $4$ 个未知数，未知数个数降为 $4Q$。
发现系数只有 $0,1$，使用 bitset 优化，复杂度降为 $O(\frac{Q^3}{w})$，可以通过。

Code

通过研读部分提交记录才得以大致理解做法，代码可能有参考。

#include<bits/stdc++.h>
#define il inline
using namespace std;
il int read()
{
    int xr=0,F=1; char cr;
    while(cr=getchar(),cr<'0'||cr>'9') if(cr=='-') F=-1;
    while(cr>='0'&&cr<='9') 
        xr=(xr<<3)+(xr<<1)+(cr^48),cr=getchar();
    return xr*F;
}
const int N=2005;
int n,q;
struct node{
	int a,b,c,d,e;
};
vector<node> t;
bitset<(N<<2)> now,a[N<<2];
int id[N][N],s[N][N],tot;
il void add(int x,int y) 
{
	if(!x||!y) return;
	if(!id[x][y]) id[x][y]=++tot;
	now[id[x][y]]=1;
}
bool flag;
il void ins()
{
	for(int i=1;i<=tot;i++)
	{
		if(!now[i]) continue;
		if(!a[i][i]) {a[i]=now;return;}
		now^=a[i];
	}
	if(now[0]) {printf("No\n");flag=1;}
}
int ans[N][N];
int main()
{
	n=read(),q=read();
	for(int i=1;i<=q;i++)
	{
		int a=read(),b=read(),c=read(),d=read(),e=read();
		if(!e) t.push_back({a,b,c,d,e});
		else
		{
			s[a][c]++,s[b+1][d+1]++,s[b+1][c]--,s[a][d+1]--;
			now.reset();
			add(a-1,c-1),add(b,d),add(a-1,d),add(b,c-1);
			if(e==2) now[0]=1; ins();
			if(flag) return 0;
		}
	}
	//for(int i=1;i<=tot;i++,cout<<endl) for(int j=1;j<=n+1;j++) cout<<a[i][j]<<" ";
	
	for(int i=1;i<=n;i++)
		for(int j=1;j<=n;j++)
			s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+s[i][j];
	for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) s[i][j]=!s[i][j];
	for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) s[i][j]=s[i][j-1]+s[i-1][j]-s[i-1][j-1]+s[i][j];
	for(auto u:t)
		if(s[u.b][u.d]-s[u.b][u.c-1]-s[u.a-1][u.d]+s[u.a-1][u.c-1]==0) {printf("No\n");return 0;}
	
	
	now.reset();
	for(int i=tot;i;i--)
	{
		if(!a[i][i]) continue;
		now[i]=((now&a[i]).count()&1)^(a[i][0]);
	}
	for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) if(id[i][j]) ans[i][j]=now[id[i][j]];
	for(int i=n;i;i--) for(int j=n;j;j--) 
		ans[i][j]=ans[i][j]^ans[i-1][j]^ans[i][j-1]^ans[i-1][j-1];
	
	printf("Yes\n");
	for(int i=1;i<=n;i++,printf("\n"))
		for(int j=1;j<=n;j++)
		{
			if(s[i][j]-s[i-1][j]-s[i][j-1]+s[i-1][j-1]) printf("0 ");
			else if(ans[i][j]) printf("2 ");
			else printf("1 ");
		}
	return 0;
}