【一句话题意】给一个有弧形的图,问(x,y)所在的区域的面积大小。
【分析】由于x和y都是整型,极大地限制了图的可能性,所以直接按点广搜并且稍微分类讨论一下就可以了。所有询问都可以放在离线上做。
【code】
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define pi 3.14159265
using namespace std;
const int maxn=110;
int n,m,q;
int g[maxn][maxn];
char s[maxn];
double block[maxn*maxn<<2];
int clr[maxn<<1][maxn<<1],tot;
int used[maxn<<1][maxn<<1];
struct node{
int x,y;
};
void bfs(int px,int py){
clr[px][py]=++tot;
block[tot]=0;
queue<node>q;
q.push((node){px,py});
while(!q.empty()){
node u=q.front();q.pop();
int x=u.x>>1,y=u.y>>1;
if(x>0&&y>0&&!used[u.x][u.y]){
if(g[x][y]) block[tot]+=4-pi/2,clr[u.x-1][u.y-1]=tot,clr[u.x-2][u.y-2]=tot,used[u.x][u.y]=1,used[u.x-1][u.y-1]=1,q.push((node){u.x-2,u.y-2});
else block[tot]+=pi/4,used[u.x][u.y]=1;
}
if(x>0&&y<m&&!used[u.x][u.y+1]){
if(g[x][y+1]) block[tot]+=pi/4,used[u.x][u.y+1]=1;
else block[tot]+=4-pi/2,clr[u.x-1][u.y+1]=tot,clr[u.x-2][u.y+2]=tot,used[u.x][u.y+1]=1,used[u.x-1][u.y+2]=1,q.push((node){u.x-2,u.y+2});
}
if(x<n&&y>0&&!used[u.x+1][u.y]){
if(g[x+1][y]) block[tot]+=pi/4,used[u.x+1][u.y]=1;
else block[tot]+=4-pi/2,clr[u.x+1][u.y-1]=tot,clr[u.x+2][u.y-2]=tot,used[u.x+1][u.y]=1,used[u.x+2][u.y-1]=1,q.push((node){u.x+2,u.y-2});
}
if(x<n&&y<m&&!used[u.x+1][u.y+1]){
if(g[x+1][y+1]) block[tot]+=4-pi/2,clr[u.x+1][u.y+1]=tot,clr[u.x+2][u.y+2]=tot,used[u.x+1][u.y+1]=1,used[u.x+2][u.y+2]=1,q.push((node){u.x+2,u.y+2});
else block[tot]+=pi/4,used[u.x+1][u.y+1]=1;
}
}
}
int main(){
freopen("flooding.in","r",stdin);
freopen("flooding.out","w",stdout);
cin>>n>>m;
for(int i=1;i<=n;i++){
scanf("%s",s+1);
for(int j=1;j<=m;j++)
g[i][j]=s[j]-'0';
}
for(int i=0;i<=2*n;i+=2)
for(int j=0;j<=2*m;j+=2)
if(!clr[i][j])
bfs(i,j);
cin>>q;
while(q--){
int x,y;
// scanf("%d%d",&x,&y);
cin>>x>>y;//学校oj对scanf极不友好,所以不是快读就是cin
if((x&1)^(y&1)) printf("%.4lf\n",0);
else printf("%.4lf\n",block[clr[x][y]]);
}
return 0;
}