Description
“十字绣”就是在布上做针线活,交替地在布的两面穿线。
布是一个n*m的网格,线只能在网格的顶点处才能从布的一面穿到另一面。每一段线都覆盖一个单位网格的两条对角线之一,而在绣的过程中,一针中连续的两段线必须分处布的两面。给出布两面的图案,问最少需要几针才能绣出来?一针是指针不离开布的一次绣花过程。
Input
第1行两个数N和M(1<=N,M<=200)。
接下来N行每行M个数描述正面。
再接下来N行每行M个数描述反面。
每个格子用.(表示空),/(表示从右上角连到左下角),\(表示从左上角连到右下角)和X(表示连两条对角线)表示。
Output
一个数,最少要用的针数。
Sample Input
4 5
.....
.\...
..\..
.....
.....
....\
.\X..
.....
Sample Output
4
Source
Ural State University Problem Archive
【分析】
有一个很显然的结论,那就是各个连通块之间的针数互不影响。那么我们就可以单独考虑每个连通块,然后将答案累加。
为方便叙述,称正面图的边为正边,负面图的边为负边。
对于每个连通块,其每个节点都连了正边,或负边,而且刺绣的时候,正负一定是交替进行的。所以如果abs(正边数-负边数)=ki>0,那么表示以该点开始或结束的针数至少为ki。而如果某一针开始和结束在同一点,那么一定不会加在ki身上,所以我们修改刚才的叙述为:以该点开始或结束的针数一定为ki。对于每个连通块,我们记tempi=sigma(ki)。又由于每一针都有开始和结束,所以最终答案为ans=sigma(tempi/2)。
(有一点要注意,如果某个连通块的tempi为0,那么另其为1,因为该连通块只用了一针,而且开始和结束的点相同)。
【代码】
/*
ID:Ciocio
LANG:C++
DATE:2013-12-07
TASK:NKOJ 1105
*/
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <iostream>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
#define MAXM 800005
#define loc(x,y) ((x-1)*(M+1)+y)
int N,M,temp;
int Y[MAXM],Last[MAXM],Next[MAXM],degree[MAXM];
bool vis[MAXM];
int tot;
void _addedge(int a,int b){Y[++tot]=b;Next[tot]=Last[a];Last[a]=tot;}
void _read(char &x){for(x=getchar();x!='.'&&x!='/'&&x!='\\'&&x!='X';x=getchar());}
void _case1(int i,int j,int t) // for '/'
{
int x,y;
x=loc(i,j+1);y=loc(i+1,j);
_addedge(x,y);
_addedge(y,x);
degree[x]+=t;degree[y]+=t;
}
void _case2(int i,int j,int t) // for '\'
{
int x,y;
x=loc(i,j);y=loc(i+1,j+1);
_addedge(x,y);
_addedge(y,x);
degree[x]+=t;degree[y]+=t;
}
void _init()
{
scanf("%d%d",&N,&M);
char tt;
for(int i=1;i<=N;i++)
for(int j=1;j<=M;j++)
{
_read(tt);
if(tt=='/') _case1(i,j,1);
else if(tt=='\\') _case2(i,j,1);
else if(tt=='X') _case1(i,j,1),_case2(i,j,1);
}
for(int i=1;i<=N;i++)
for(int j=1;j<=M;j++)
{
_read(tt);
if(tt=='/') _case1(i,j,-1);
else if(tt=='\\') _case2(i,j,-1);
else if(tt=='X') _case1(i,j,-1),_case2(i,j,-1);
}
N=loc(N+1,M+1);
}
void _DFS(int x)
{
vis[x]=true;
temp+=abs(degree[x]);
for(int j=Last[x];j;j=Next[j])
if(!vis[Y[j]])
_DFS(Y[j]);
}
void _solve()
{
int ans=0;
for(int i=1;i<=N;i++)
if((!vis[i])&&Last[i])
{
temp=0;
_DFS(i);
ans+=temp?temp/2:1;
}
cout<<ans<<endl;
}
int main()
{
_init();
_solve();
return 0;
}