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;
        }
    }
}