【CERC2015】【BZOJ4434】Ice Igloos

Description

一个坐落在北极冰冻湖面上的渔村因全球变暖而面临危险——湖面开始破碎。村庄内有n座球形冰屋，每座冰屋在湖面上占据了一块圆形区域。
一座冰屋可以视为平面座标系中的一个圆：其圆心座标为整数，其半径为小于1的一位浮点数。给出冰面断裂的座标，这些村民希望知道每次断裂会有多少个冰屋受到影响。
形式上，给出q个询问，每个询问包括一条由两个端点确定的线段，查询每条线段贯穿的冰屋数。如果线段与圆的内部有至少一个交点，则线段贯穿了这座冰屋。

Input

第一行包括一个整数n(1<=n<=100 000)——冰屋数。
接下来的n行每行包括三个数x，y和r——冰屋的座标和半径。x和y为整数且1<=x,y<=500。r是一位浮点数且0

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define eps 1e-6
#define SIZE 510
#define MAXN 100010
#define GET (ch>='0'&&ch<='9')
using namespace std;
int n,q;
double r[SIZE][SIZE];
double A,B,C,k,b;
inline void in(int &x)
{
    char ch=getchar();x=0;
    while (!GET)    ch=getchar();
    while (GET) x=x*10+ch-'0',ch=getchar();
}
inline void fin(double &x)
{
    char ch=getchar();x=0;double f=1;
    while (!GET)    ch=getchar();
    while (GET) x=x*10+ch-'0',ch=getchar();
    while (!GET)    ch=getchar();
    while (GET) f*=0.1,x=x+f*(ch-'0'),ch=getchar();
}
inline void init(int x1,int y1,int x2,int y2)
{
    A=(y2-y1);B=(x1-x2);
    double len=sqrt(A*A+B*B);
    A/=len;B/=len;C=-(A*x1+B*y1);
    k=(double)(y2-y1)/(x2-x1);b=y1-k*x1;
}
inline double f(int x)  {   return k*x+b;   }
inline double dis(int x,int y)  {   return fabs(A*x+B*y+C); }
inline int calc(int x,int y1,int y2)
{
    int ret=0;
    for (int i=y1;i<=y2;i++) ret+=(r[x][i]>0.0&&dis(x,i)<=r[x][i]);
    return ret;
}
int main()
{
    in(n);int x,y,x1,y1,x2,y2;double z;
    for (int i=1;i<=n;i++)   in(x),in(y),fin(z),r[x][y]=z;
    for (in(q);q;q--)
    {
        in(x1);in(y1);in(x2);in(y2);int ans=0;
        if (x1>x2||(x1==x2&&y1>y2))   swap(x1,x2),swap(y1,y2);
        if (x1==x2) for (int i=y1;i<=y2;i++) ans+=(r[x1][i]>0.0);
        else    if (y1==y2) for (int i=x1;i<=x2;i++) ans+=(r[i][y1]>0.0);
        else
        {
            init(x1,y1,x2,y2);
            if (y2>y1)
            {
                int ty1=y1,ty2=ceil(f(x1+1)-eps);
                ans+=calc(x1,ty1,ty2);
                for (int i=x1+1;i<x2;i++)
                {
                    int tmp=ty2;ty2=ceil(f(i+1)-eps);
                    ans+=calc(i,ty1,ty2);ty1=tmp-1;
                }
                ans+=calc(x2,ty1,ty2);
            }
            else
            {
                int ty1=y1,ty2=floor(f(x1+1)+eps);
                ans+=calc(x1,ty2,ty1);
                for (int i=x1+1;i<x2;i++)
                {
                    int tmp=ty2;ty2=floor(f(i+1)+eps);
                    ans+=calc(i,ty2,ty1);ty1=tmp+1;
                }
                ans+=calc(x2,ty2,ty1);
            }
        }
        printf("%d\n",ans);
    }
}