JZOJ2724. 【3.16XJ模拟题】圆

题目描述

二维坐标平面内有n个圆，第i个圆圆心在(Xi,Yi)，半径为Ri，权值Vi。任何两个圆都不会相交（也不会相切），但是圆与圆之间可能存在包含关系。当我们在一个圆里面的时候，我们必须经过它的边界一次，才能走出这个圆。

对于不同的两个圆A、B，如果可以从A到B且经过不超过K次边界，则称A、B是连通的。现在的问题是：对于所有的连通的一对圆(A,B)，权值差的绝对值最大有多少。即找到max{ |VA – VB|, 其中A B连通}。
若圆A包含圆B，那么圆A内的区域定义为，A挖去圆B的剩余部分。

对于20%的数据，n ≤ 1000
对于100%的数据，T ≤ 10, 2 ≤ n ≤ 50000, 2 ≤ K ≤ 100, 0 < Ri ≤ 200000000, -1000000000 ≤ Xi, Yi, Vi ≤ 1000000000，

分析

我们发现这些圆很微妙，他们不会相交，那么一个圆A包含另一个圆B，就可以看成圆A是B的祖先，而最小的包含B的就是它的父亲。那么整个图会转化成一棵树，而恰好，一个圆到另一个圆内，就是树中简单路径的长度，那么这样的话K这么小，设F[i][j]表示与i距离在j及以内的所有点的最小点权。那么转移先枚举j再枚举i及与其相连的点即可。答案就是max{v[i]-F[i][k]}了。
现在问题是如何转化？怎么找父亲？
暴力的话只有20分。
考虑扫描线，这个很关键。
考虑扫面线形如x=a。它跟一些圆肯定有交点嘛。那么这些交点代表了什么呢？
我们发现，假如一个圆B的两个交点都在另一个圆A的两个交点之内，那么肯定是包含关系。那么父亲肯定是所有圆A中两交点距离最小的一个。如图。

因为是对称的，我们只需要维护上交点或下交点即可。但是这个时候我们要分类讨论。
我们每一次扫到一个新的圆C，目标就是找到离他最近的上（下）交点X，然后这个时候有两种情况。
1，假如X的另一个交点Y与X在扫描线上形成的线段包含C与扫描线的切点，X所代表的圆D就是C父亲了。
2，假如不包含，那么他们形成并列关系。那么C的父亲就是D的父亲嘛。
怎么维护这些交点呢？

观察一下性质，假如移动了扫描线，交点的情况如何？
我们发现，实际上它们的相对顺序是不变的！（假如扫描线没有离开圆）
考虑使用平衡树维护交点点集。那么每次只需要树上二分查找新加入圆的切点的最近的交点是哪个，再判两种情况即可。离开了一个圆就直接删掉这个点。
这样NlogN解决了。

代码

定义struct和变量那一块打烂了。

#include<cstdio>
#include<algorithm>
#include<cmath>
#include<set>
using namespace std;
typedef long long ll;
typedef double db;
#define fo(i,j,k) for(i=j;i<=k;i++)
#define fd(i,j,k) for(i=j;i>=k;i--)
const int N=50005,mx=1e9+7;
struct vec
{
    int x,y;
    vec (int x0=0,int y0=0) {x=x0,y=y0;}
}c[N*2];
struct circle
{
    vec p;
    int r,v;
}a[N];
vec operator -(vec a,vec b)
{
    return vec(a.x-b.x,a.y-b.y);
}
ll operator *(vec a,vec b)
{
    return (ll)a.x*b.x+(ll)a.y*b.y;
}
ll Len2(vec a)
{
    return a*a;
}
int T,n,k,ans,i,j,T1,p,t1,now;
int first[N],next[N*2],b[N*2],tt,fa[N],f[N][105];
db pd(int i)
{
    return (db)a[i].p.x-sqrt((db)a[i].r*a[i].r-(db)(a[i].p.x-now)*(a[i].p.x-now));
}
struct cmp
{
    bool operator()(const int &x,const int &y)
    {
        return pd(x)>pd(y);
    }
};
multiset<int,cmp> tr;
multiset<int,cmp>::iterator dur;
void clear()
{
    ans=0;
    tr.clear();
    fo(i,1,n+1) first[i]=0;
    t1=tt=0;
}
bool cmp(vec i,vec j) 
{
    return a[i.x].p.x+i.y*a[i.x].r<a[j.x].p.x+j.y*a[j.x].r;
}

bool Include(int x,int y)//y include x
{
    return (ll)a[y].r*a[y].r>Len2(a[x].p-a[y].p);
}
void getfa()
{
    fo(i,1,t1)
    {
        now=a[c[i].x].p.x+c[i].y*a[c[i].x].r;
        if (c[i].y==-1) 
        {
            dur=tr.upper_bound(c[i].x);
            if (dur==tr.end())
                fa[c[i].x]=1;
            else
            if (Include(c[i].x,*dur))
                fa[c[i].x]=*dur;
            else fa[c[i].x]=fa[*dur];
            tr.insert(c[i].x);
        }
        else tr.erase(tr.find(c[i].x));
    }
}
void solve() 
{
    fo(j,1,k)
        fo(i,1,n+1)
        {
            f[i][j]=f[i][j-1];
            for(p=first[i];p;p=next[p])
                f[i][j]=min(f[i][j],f[b[p]][j-1]);
        }
    fo(i,2,n+1)
        ans=max(ans,a[i].v-f[i][k]);
}
void cr(int x,int y)
{
    tt++;
    b[tt]=y;
    next[tt]=first[x];
    first[x]=tt;
}
int main()
{
    freopen("t3.in","r",stdin);
    scanf("%d",&T);
    while (T--)
    {
        T1++;
        clear();
        scanf("%d %d",&n,&k);
        a[1].v=mx;
        f[1][0]=mx;
        fo(i,2,n+1)
        {
            scanf("%d %d %d %d",&a[i].p.x,&a[i].p.y,&a[i].r,&a[i].v);
            c[++t1].x=i;
            c[t1].y=1;
            c[++t1].x=i;
            c[t1].y=-1;
            f[i][0]=a[i].v;
        }
        sort(c+1,c+1+t1,cmp);
        getfa();
        fo(i,2,n+1)
        {
            cr(fa[i],i);
            cr(i,fa[i]);
        }
        solve();
        printf("Case %d: %d\n",T1,ans);
    }
    return 0;
}