【JZOJ A组】着色

Description

Alice是一个奇怪的画家。她想对一副有N*N个像素点组成的画进行着色，N是2的幂（1，2，4，8，16等等）。每个像素点可以着成黑色或白色。
　　Alice着色方案不是唯一的，她采用以下不确定的规则：
　　•如果画作只有一个像素点，那可以直接着白色或黑色；
　　•否则，把画平均分成四块，然后进行以下操作：
　　(1) 选择一块全部着白色；
　　(2) 选择一块全部着黑色；
　　(3) 把剩下的两块当作是独立的画作并采用同样的方法进行着色。
　　对于每一幅画作，Alice心目中已经有一个蓝图，接下来请你帮她采用上述方法着色，要求选择跟心目中的蓝图差异最小的着色方案，当然要遵循上述的着色规则，两幅图的差异是指对应位置颜色不相同的像素点的个数。

Input

输入第一行包含整数N(1<=N<=512)，表示画作的尺寸为N*N，N保证是2的幂。
　　接下来N行每行包含N个0或1，描述Alice心目中的蓝图，0表示白色，1表示黑色。

Output

第一行输出最小的差异是多少。

Sample Input

输入1：
4
0001
0001
0011
1110

输入2：
4
1111
1111
1111
1111

输入3：
8
01010001
10100011
01010111
10101111
01010111
10100011
01010001
10100000

Sample Output

输出1：
1

输出2：
6

输出3：
16

Hint

【数据范围】
　　50%的数据N<=8

思路

DP

设f[i][j][k][0…2]表示以(i,j)为左下角的正方形，边长为2^k，染色状态为w(0…2)，与脑补图的不同点数

0代表全白，1代表全黑，2代表按要求自由染色

0和1很容易处理（把四个矩阵加起来）
2根据题意有12种情况。。。

具体看代码

代码

#include<bits/stdc++.h>
const int N=577;
using namespace std;
int f[N][N][13][3],mi[13],n;
int main()
{
	memset(f,127,sizeof(f));
	scanf("%d\n",&n);
	mi[0]=1;
	for(int i=1; i<=12; i++)
		mi[i]=mi[i-1]*2;
	for(int i=1; i<=n; i++) 
	{
		for(int j=1; j<=n; j++)
		{
			char ch;
			scanf("%c",&ch);
		    int t=ch-48;
		    f[i][j][0][0]=t;
		    f[i][j][0][1]=1-t;
		    f[i][j][0][2]=0;
		}
		scanf("\n");
	}
	for(int k=1; k<=log2(n); k++)
	{
		for(int i=1; i<=n-mi[k]+1; i++) for(int j=1; j<=n-mi[k]+1; j++)
		{
			f[i][j][k][0]=f[i][j][k-1][0]+f[i+mi[k-1]][j][k-1][0]+f[i][j+mi[k-1]][k-1][0]+f[i+mi[k-1]][j+mi[k-1]][k-1][0];
            f[i][j][k][1]=f[i][j][k-1][1]+f[i+mi[k-1]][j][k-1][1]+f[i][j+mi[k-1]][k-1][1]+f[i+mi[k-1]][j+mi[k-1]][k-1][1];

            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][0]+f[i+mi[k-1]][j][k-1][1]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][0]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][1]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][0]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][1]);

            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][1]+f[i+mi[k-1]][j][k-1][0]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][1]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][0]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][1]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][0]);

            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][0]+f[i][j+mi[k-1]][k-1][1]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);
            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][0]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][1]);
            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][1]+f[i][j+mi[k-1]][k-1][0]+f[i+mi[k-1]][j+mi[k-1]][k-1][2]);

            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][1]+f[i][j+mi[k-1]][k-1][2]+f[i+mi[k-1]][j+mi[k-1]][k-1][0]);
            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][0]+f[i+mi[k-1]][j+mi[k-1]][k-1][1]);
            f[i][j][k][2]=min(f[i][j][k][2],f[i][j][k-1][2]+f[i+mi[k-1]][j][k-1][2]+f[i][j+mi[k-1]][k-1][1]+f[i+mi[k-1]][j+mi[k-1]][k-1][0]);
		}
	}
	int pp=log2(n);
    printf("%d",f[1][1][pp][2]);
}