tjut 3511

原创 2016年08月30日 08:34:18
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <set>
#include <algorithm>

using namespace std;

const int UP = 0;
const int DOWN = 1;
const int IN = 0;
const int OUT = 1;
const int N = 50005;

int Time;

struct circle {
    int x, y, r;
    int w;
    void read() {
        scanf("%d %d %d", &x, &y, &r);
        w = 0;
    }
    int getX(int flag) {
        if( flag == IN )     return x - r;
        else                return x + r;
    }
    double getY(int flag) {
        double ret = sqrt((double)r*r-(double)(Time-x)*(Time-x));
        if( flag == UP )     return (double)y + ret;
        else                return (double)y - ret;
    }
    
} cir[N];

struct event {
    int x, y, id;
    int flag;
    void get(int _x, int _y, int _id, int _flag) {
        x = _x;
        y = _y;
        id = _id;
        flag = _flag;
    }
    bool operator<(const event ev) const {
        return x < ev.x || x == ev.x && y > ev.y;
    }
} eve[N*2];

struct node {
    int id;
    int flag;
    node(){}
    node(int _id, int _flag) {
        id = _id;
        flag = _flag;
    }
    bool operator<(const node a) const {
        double y1 = cir[id].getY(flag);
        double y2 = cir[a.id].getY(a.flag);
        return y1 > y2 || y1 == y2 && flag < a.flag;
    }
};

int n, eveN;

set<node> line;
set<node>::iterator it, f, e, p;

inline int max(int a, int b) { return a > b ? a : b;}

void moveline() {
    line.clear();
    for(int i = 0; i < eveN; i ++) {
        Time = eve[i].x;
        if( eve[i].flag == OUT ) {
            line.erase(node(eve[i].id, UP));
            line.erase(node(eve[i].id, DOWN));
        } else {
            it = line.insert(node(eve[i].id, UP)).first;
            e = f = it;
            e ++;
            int id = it->id;
            if( it == line.begin() || e == line.end() ) {
                cir[id].w = 1;
            } else {
                f --;
                if( f->id == e->id ) {
                    cir[id].w = cir[f->id].w + 1;
                } else {
                    cir[id].w = max( cir[f->id].w, cir[e->id].w);
                }
            }
            line.insert(node(eve[i].id, DOWN));
        }
    }
}

int main() {
    while( scanf("%d", &n) != EOF ) {
        eveN = 0;        
        for(int i = 0; i < n; i ++) {
            cir[i].read();
            eve[eveN++].get(cir[i].getX(IN), cir[i].y, i, IN);
            eve[eveN++].get(cir[i].getX(OUT), cir[i].y, i, OUT);
            
        }
        sort(eve, eve + eveN);
        moveline();
        int ans = 0;
        for(int i = 0; i < n; i ++) {
            ans = max(ans, cir[i].w);
        }
        printf("%d\n", ans);
    }
}

版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

tjut 5412

#include using namespace std; const int Maxn=100020; int nq; int ty[Maxn*3],ql[Maxn*3],qr[Maxn*3],qk...

TL_IPC_HI3511_ESV1.6_V3.3.8_140731.tar.gz

  • 2014年10月30日 00:43
  • 2.11MB
  • 下载

【BZOJ3511】土地划分【最小割】

看了题解才会建图... 对于点x,从S向x连容量为va的边,从x向T连容量为vb的边。 对于边(u, v),从S向u和v分别连两条容量为ea / 2的边,从u和v向T分别连两条容量为eb ...

ZOJ3511 Cake Robbery,暴力

一个凸包,有多少个顶点,就有多少条边。 切的刀痕不想交,所以对于每一刀,只要枚举对于现有的每一个凸包,到底落在哪一个多边形上,然后在把这个多边形分成两份就可以了。 /***********...
  • neofung
  • neofung
  • 2012年07月11日 19:15
  • 759

BZOJ3511 土地划分

一眼最小割 考虑划为S集代表A国,划为T集代表B国 建图:S连每个点流量VA,每个点连T流量VB 对于每条边,两个端点点之间连双向边流量EC 新建一个点X,S连X流量EA,X连两个端点流量IN...

hdu 3511

  题目意思实际上就是找嵌套数最大的圆并输出其嵌套数。  一开始见时间有5000ms,就果断暴力,然后果断超时— — 之后用二分,终于不tle啦,但wa。。。TAT。上网找题解,发现有两种方法,一个是...

HI3511、3512SDK

  • 2012年03月09日 16:29
  • 8.43MB
  • 下载

BZOJ 3511: 土地划分

题意:Y国有N座城市,并且有M条双向公路将这些城市连接起来,并且任意两个城市至少有一条路径可以互达。 Y国的国王去世之后,他的两个儿子A和B都想成为新的国王,但他们都想让这个国家更加安定,不会用武...

Hi3511/Hi3512_硬件设计用户指南

  • 2012年07月05日 11:40
  • 2.31MB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:tjut 3511
举报原因:
原因补充:

(最多只允许输入30个字)