【HNOI2016模拟4.10】线性代数与逻辑

这题废了。。。
首先，题意一定要明白。
我们要构造一个y数组（一维，题目的）
然后再构造出我们要的X矩阵。（满足X_i，j=y_i^y_j）
由于题目要

也就是对于每一个位置（i，j），要使X[i,j]|(!A[i,j])

也就是说，如果A[i,j]=1，那么X[i,j]一定要等于1，否则随便

然后，我们要让

也就是说，y中0的个数和1的个数差最小。
答案就是

2*y中0的个数**y中1的个数（也就是n-y中0的个数）-A矩阵中1的个数

上标：

#include<cstdio>
#include<cstring>
#define max(x,y) x=x>y ? x:y
using namespace std;
int T,n,m,a[1010][1010],cnt=0,tot=0;
int fa[1010],tf[1010],sum[1010][2],z[1010][2];
int f[1010][1010],s;

inline int read()
{
	int x=0; char c=getchar();
	while (c<'0' || c>'9') c=getchar();
	while (c>='0' && c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
	return x;
}

int gf(int x)
{
	if (!fa[x]) return x;
	int ff=gf(fa[x]);
	tf[x]=tf[fa[x]]^tf[x];
	return fa[x]=ff;
}

int haha()
{
	memset(fa,0,sizeof(fa));
	memset(tf,0,sizeof(tf));
	memset(sum,0,sizeof(sum));
	memset(z,0,sizeof(z));
	for (int i=1;i<=n;i++) sum[i][0]=1;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=n;j++)
			if (a[i][j])
			{
				int x=gf(i),y=gf(j);
				if (x!=y)
				{
					if (tf[i]^tf[j]==0) tf[x]^=1;
					sum[y][0]+=sum[x][tf[x]];
					sum[y][1]+=sum[x][tf[x]^1];
					fa[x]=y;
				}
				else if (tf[i]==tf[j]) return -1;
			}
	for (int i=1;i<=n;i++)
		if (!fa[i]) z[++tot][0]=sum[i][0],z[tot][1]=sum[i][1];
	memset(f,-1,sizeof(f));
	f[0][0]=0;
	for (int i=1;i<=tot;i++)
		for (int j=0;j<=n;j++)
		{
			if (j>=z[i][0]) f[i][j]|=f[i-1][j-z[i][0]];
			if (j>=z[i][1]) f[i][j]|=f[i-1][j-z[i][1]];
		}
	for (int i=n/2;i>0;i--)
		if (f[tot][i]) return i*(n-i)*2-s;
	return -1;
}

int main()
{
	freopen("C.in","r",stdin);
//	freopen("C.out","w",stdout);
	T=read();
	while (T--)
	{
		n=read();s=0;
		for (int i=1;i<=n;i++)
			for (int j=1;j<=n;j++)
				s+=(a[i][j]=read());
		printf("%d
",haha());
	}
	return 0;
}