SSL2787 2017年10月24日提高组 一个不成熟的想法（dfs）

2017年10月24日提高组 一个不成熟的想法

Description

Ymw因为十分厉害，刷题已经不能满足他了，他想去写一个小程序来实现一些功能，现在有一个2N*2M的矩形水池，被一些宽度不计的隔板所隔开，每个2*2的小矩形中都有两块弧形隔板，且只有可能是这两种图案中的一种：





现在你看到了ymw的水池的样子，ymw有Q个问题，每次想知道如果在(x,y)上倒了一同染料，则最终可以扩散多大的面积？

Input



Output



分析：显然可得染料在每个矩形中的扩散只有四种情况，即从4个顶点开始扩散，这样两种图案就有8种，dfs一下就好啦。

代码

#include <cstdio>
#include <cstring>
#include <string>
#include <iostream>
#define maxn 500
#define pi 3.141592653589793
using namespace std;

int a[maxn][maxn],n,m;
double ans;
bool f[maxn][maxn],v[maxn][maxn][5];

bool check(int x,int y,int tmp)
{
if (x<0||x>n-1||y<0||y>m-1) return false;
if (a[x][y]==0)
{
if (tmp==1||tmp==4)
if (v[x][y][tmp]) return false;
if (tmp==2||tmp==3)
if (v[x][y][2]||v[x][y][3]) return false;
}
else
{
if (tmp==2||tmp==3)
if (v[x][y][tmp]) return false;
if (tmp==1||tmp==4)
if (v[x][y][1]||v[x][y][4]) return false;
}
return true;
}

void dfs(int x,int y,int tmp)
{
if (!check(x,y,tmp)) return;
v[x][y][tmp]=true;
if (a[x][y]==0)
{
if (tmp==1)
{
ans+=pi/4;
return;
}
if (tmp==2)
{
ans+=4-pi/2;
dfs(x,y-1,4);
dfs(x+1,y,1);
dfs(x+1,y-1,2);
}
if (tmp==3)
{
ans+=4-pi/2;
dfs(x-1,y,4);
dfs(x,y+1,1);
dfs(x-1,y+1,3);
}
if (tmp==4)
{
ans+=pi/4;
return;
}
}
else
{
if (tmp==1)
{
ans+=4-pi/2;
dfs(x+1,y,2);
dfs(x,y+1,3);
dfs(x+1,y+1,1);
}
if (tmp==2)
{
ans+=pi/4;
return;
}
if (tmp==3)
{
ans+=pi/4;
return;
}
if (tmp==4)
{
ans+=4-pi/2;
dfs(x,y-1,2);
dfs(x-1,y,3);
dfs(x-1,y-1,4);
}
}
}

int main()
{
scanf("%d%d",&n,&m);
for (int ii=0;ii<=n-1;ii++)
for (int jj=0;jj<=m-1;jj++)
{
char ch;
cin>>ch;
a[ii][jj]=ch-'0';
int i=ii;
int j=jj;
if (ch=='0')
{
f[i*2][j*2+1]=true;
f[i*2+1][j*2]=true;
f[(i+1)*2][j*2+1]=true;
f[i*2+1][(j+1)*2]=true;
}
else
{
f[i*2][j*2+1]=true;
f[i*2+1][(j+1)*2]=true;
f[i*2+1][j*2]=true;
f[(i+1)*2][j*2+1]=true;
}
}
int q;
scanf("%d",&q);
for (int i=1;i<=q;i++)
{
int x,y;
scanf("%d%d",&x,&y);
memset(v,false,sizeof(v));
if (f[x][y]) printf("0.0000\n");
else
{
if (x%2==1&&y%2==1)
{
x/=2;
y/=2;
if (a[x][y]==1)
{
ans=4-pi/2;
v[x][y][1]=true;
dfs(x-1,y-1,4);
dfs(x,y-1,2);
dfs(x-1,y,3);
dfs(x+1,y,2);
dfs(x,y+1,3);
dfs(x+1,y+1,1);
}
else
{
ans=4-pi/2;
v[x][y][2]=true;
dfs(x-1,y,4);
dfs(x,y+1,1);
dfs(x-1,y+1,3);
dfs(x+1,y,1);
dfs(x+1,y-1,2);
dfs(x,y-1,4);
}
}
else
{
x/=2;
y/=2;
dfs(x,y,1);
dfs(x-1,y,3);
dfs(x-1,y-1,4);
dfs(x,y-1,2);
}
printf("%.4f\n",ans);
ans=0;
}
}
}