给出一张n*m的图,有些点已被染色。求剩下未被染色的点能被1*2的小矩阵染的最大面积,输出最大染色方案的任意一组可行解。
二分图的行列式匹配。因为相邻的两点的i+j的奇偶性一定不同,可以根据这个性质进行建图,即X, Y部图分别对应奇偶,求最大匹配即可(奇偶性确实难想,但是这题的图不大,可以枚举每个点的 相邻点建图。好像是可以的= =)
#include<stdio.h>
#include<string.h>
#include<iostream>
#include<algorithm>
#include<vector>
using namespace std;
#define M 105
int mat[M][M];
vector<int> v[M];
int link[M], vis[M][M];
int n, m, k;
int dir[4][2] = {
1, 0, -1, 0, 0, 1, 0, -1
};
struct node{
int x, y;
}p[M][M];
void init()
{
for( int i = 1; i <= 105; i++ )
{
v[i].clear();
}
memset(mat, 0, sizeof(mat));
memset(p, 0, sizeof(p));
}
bool f( int x, int y )
{
if( x >= 1 && x <= n && y >= 1 && y <= m )
return 1;
return 0;
}
int find( int u )
{
int x = u / 105;
int y = u % 105;
for( int i = 0; i < 4; i++ )
{
int xx = x + dir[i][0];
int yy = y + dir[i][1];
if( !vis[xx][yy] && f(xx, yy) && !mat[xx][yy] )
{
vis[xx][yy] = 1;
if( p[xx][yy].x == 0 || find(p[xx][yy].x *105 + p[xx][yy].y) )
{
p[xx][yy].x = x;
p[xx][yy].y = y;
return 1;
}
}
}
return 0;
}
int match()
{
int x = 0;
for( int i = 1; i <= n; i++ )
{
for( int j = 1; j <= m; j++ )
{
memset(vis, 0, sizeof(vis));
if( (i+j) % 2 == 0 && !mat[i][j] )
x += find(i*105 + j);
}
}
return x;
}
int main()
{
int t = 0;
while(~scanf("%d%d", &n, &m) && (m || n))
{
if( t )
printf("\n");
t = 1;
init();
scanf("%d", &k);
int a, b;
while(k--)
{
scanf("%d%d", &a, &b);
mat[a][b] = 1;
}
int ans = match();
printf("%d\n", ans);
for( int i = 1; i <= n; i++ )
{
for ( int j = 1; j <= m; j++ )
{
if(p[i][j].x)
{
printf("(%d,%d)--(%d,%d)\n", i, j, p[i][j].x, p[i][j].y);
}
}
}
}
return 0;
}