ADV-2-矩阵翻转

问题描述

Ciel有一个N*N的矩阵，每个格子里都有一个整数。

N是一个奇数，设X = (N+1)/2。Ciel每次都可以做这样的一次操作：他从矩阵选出一个X*X的子矩阵，并将这个子矩阵中的所有整数都乘以-1。

现在问你经过一些操作之后，矩阵中所有数的和最大可以为多少。

输入格式

第一行为一个正整数N。

接下来N行每行有N个整数，表示初始矩阵中的数字。每个数的绝对值不超过1000。

输出格式

输出一个整数，表示操作后矩阵中所有数之和的最大值。

样例输入

3
-1 -1 1
-1 1 -1
1 -1 -1

样例输出

9

数据规模与约定

1 <= N <= 33，且N为奇数。

PS：

锦囊1

使用贪心。

锦囊2

令dp[i][j]表示矩阵中(i, j)位置的数是否要乘以-1。

注意dp[i][j]^dp[i][x]^dp[i][j+x] == 0,

dp[i][j]^dp[x][j]^dp[i+x][j] == 0 ;

所以只要确定了dp[x][j], dp[i][x], 就能将dp[i][j] 与 dp[i+x][j], dp[i][j+x]绑定起来了。

首先我们可以枚举dp[x][j]--可以通过枚举使用了哪些以(x, j)为左上角的矩形，那么这样，我们可以同时将dp[x][j+x]确定下来。

然后我们依次dp[j][x](j<x)--我们可以发现，对于每一个dp[j][x],其影响的只有dp[j][t](t<x)与dp[j][t+x]的关系以及dp[j+x][x]，所以任两个dp[j][x]是线性无关的。

那么我们就可以一个个枚举dp[j][x]是否为1，然后运用贪心。

//-----C----------------
#include<stdio.h>
int x[33][33],ans,N;
void fun1(int n)
{
	int i,j,lin=0,aa,bb;
	for(j=0;j<N;j++)
		lin+=x[n-1][j];
	for(i=0;i<n-1;i++)
	{
		aa=-1000000000;
		bb=x[i][n-1]+x[i+n][n-1];
		for(j=0;j<n-1;j++)
			bb+=abs(x[i][j]+x[i+n][j]+x[i][j+n]+x[i+n][j+n]);
		aa=aa>bb?aa:bb;
		bb=-x[i][n-1]-x[i+n][n-1];
		for(j=0;j<n-1;j++)
			bb+=abs(-x[i][j]-x[i+n][j]+x[i][j+n]+x[i+n][j+n]);
		aa=aa>bb?aa:bb;
		lin+=aa;
	}
	ans=ans>lin?ans:lin;
}
void fun(int n)
{
	int i,j,k;
	for(k=0;k<(1<<n-1);k++)
	{
		for(i=0;i<n-1;i++)
			if((k&(1<<i))!=0)
				for(j=0;j<n;j++)
				{
					x[j][i]*=-1;
					x[j][i+n]*=-1;
				}
		fun1(n);
		for(i=0;i<n-1;i++)
			if((k&(1<<i))!=0)
				for(j=0;j<n;j++)
				{
					x[j][i]*=-1;
					x[j][i+n]*=-1;
				}
	}
}
int main(void)
{
	int i,j,k;
	scanf("%d",&N);
	for(i=0;i<N;i++)
		for(j=0;j<N;j++)
			scanf("%d",&x[i][j]);
	k=(N+1)/2;
	ans=-1000000000;
	fun(k);
	for(i=0;i<k;i++)
		for(j=0;j<k;j++)
			x[i][j]=-x[i][j];
	fun(k);
	printf("%d\n",ans);
	return 0;
}




//------C++-----------

#include<stdio.h>
#include<cmath>
using namespace std;
int x[33][33],ans,N;
void fun1(int n)
{
	int i,j,lin=0,aa,bb;
	for(j=0;j<N;j++)
		lin+=x[n-1][j];
	for(i=0;i<n-1;i++)
	{
		aa=-1000000000;
		bb=x[i][n-1]+x[i+n][n-1];
		for(j=0;j<n-1;j++)
			bb+=abs(x[i][j]+x[i+n][j]+x[i][j+n]+x[i+n][j+n]);
		aa=aa>bb?aa:bb;
		bb=-x[i][n-1]-x[i+n][n-1];
		for(j=0;j<n-1;j++)
			bb+=abs(-x[i][j]-x[i+n][j]+x[i][j+n]+x[i+n][j+n]);
		aa=aa>bb?aa:bb;
		lin+=aa;
	}
	ans=ans>lin?ans:lin;
}
void fun(int n)
{
	int i,j,k;
	for(k=0;k<(1<<n-1);k++)
	{
		for(i=0;i<n-1;i++)
			if((k&(1<<i))!=0)
				for(j=0;j<n;j++)
				{
					x[j][i]*=-1;
					x[j][i+n]*=-1;
				}
		fun1(n);
		for(i=0;i<n-1;i++)
			if((k&(1<<i))!=0)
				for(j=0;j<n;j++)
				{
					x[j][i]*=-1;
					x[j][i+n]*=-1;
				}
	}
}
int main(void)
{
	int i,j,k;
	scanf("%d",&N);
	for(i=0;i<N;i++)
		for(j=0;j<N;j++)
			scanf("%d",&x[i][j]);
	k=(N+1)/2;
	ans=-1000000000;
	fun(k);
	for(i=0;i<k;i++)
		for(j=0;j<k;j++)
			x[i][j]=-x[i][j];
	fun(k);
	printf("%d\n",ans);
	return 0;
}