物理( 线段树套单调队列 )

sg站在1号点,后面有他的n-1个小弟。现在妹子想告诉sg一些话,但是她不知道选谁。由于sg的小弟们构造奇特,他们能听到的声波频率也是不一样的。第i个人可以发出并接收声波频率为[xi,yi]的声音,并且每个人的声音最多只能传播不超过l的距离。每一次说话并被另一个人接收到的过程都要消耗1单位时间。

现在给出所有人能发出并接收的声波频率的范围和位置,sg的妹子想知道,每个小弟传话给sg要多久。由于sg很牛逼,他可以接收所有小弟发出的声波,当然前提要在距离之内。sg的位置为0。

Input

第一行:n,l,含义见题意。

第二行至第n行:xi,yi,li,[xi,yi]是第i个小弟能发出并接收的声波频率范围,li是第i个小弟的位置。

Output

第一行至第n-1行:第i行表示第i+1个小弟要传话给sg所消耗的时间。

对于暴力 n2 的dp显然f[i] = minSectioniSectionj!=nullANDposiposj<=L ( f[j] ) + 1

然后我们考虑维护min的部分;
考虑开一个线段树,以频率为下标。
然后呢? 怎么解决长度问题?总不可能在线段树每个节点塞一个线性表每次遍历扫一遍
显然 posi 是增函数 , 我们如果用单调队列代替线性表当然好多了。
怎么维护这个线段树呢?
每个节点开两个单调队列,cover[] && pass[]
对于Insert操作,每经过一个节点更新一下pass[],遇到整段配对的区间更新cover[]
对于Ask操作,每经过一个节点用cover[]来更新遇到整段配对的区间更新pass[],

时间虽然给了5s仍然有点紧,有些小优化,如
cover被更新的节点不用更新pass
不要用map离散化数据
只在Ask操作内删除队首非法节点

#include
#include
#include
#include
#include

using namespace std ;

#define N 250010
#define psb( x ) push_back( x ) 

int i , j , k , n , L ;

map< int , int > s ;

struct node {
    int x , y , l , ps ;
}z[N] ;

// We Store id 
struct SegmentTree_With_Queue {
    vector<int> cover , pass ;
    int cl , cr , pl , pr  ;
    // Remember to Initialize 
}tr[N*2*4] ;

int f[N] , demap[N*2] ;

struct For_Dis {
    int ps , va ;
    bool x ;
}a[N*2] ;

bool cmp( For_Dis a , For_Dis b ) {
    return a.va < b.va ;
}

bool _cmp( node a , node b ) {
    return a.l < b.l ;
}

void Ins( int _l , int _r , int l , int r , int va , int po ) {
    int m = l+r >> 1 ;
    if( _l==l && _r==r ) {
        while( tr[po].cl<=tr[po].cr && f[ tr[po].cover[ tr[po].cr ] ] > f[va] ) tr[po].cr-- ;
        if( ( tr[po].cr + 1 ) >= tr[po].cover.size() ) tr[po].cover.psb( 0 );
        tr[po].cover[ ++tr[po].cr ] = va ;
        return ;
    }
    while( tr[po].pl<=tr[po].pr && f[ tr[po].pass[ tr[po].pr ] ] > f[va] ) tr[po].pr -- ;
    if( ( tr[po].pr+1 ) >= tr[po].pass.size() ) tr[po].pass.psb( 0 ) ;
    tr[po].pass[ ++tr[po].pr ] = va ;
    if( _r<=m ) Ins( _l,_r,l,m,va,po*2 ) ;
        else if( _l>m ) Ins( _l,_r,m+1,r,va,po*2+1 ) ;
            else Ins( _l,m,l,m,va,po*2 ) , Ins( m+1,_r,m+1,r,va,po*2+1 );

}

int cur , ax = 0 ;

void Ask( int _l , int _r , int l , int r , int po  ) {
    int m = l + r >> 1 ;
    while( tr[po].cl<=tr[po].cr && -z[ tr[po].cover[ tr[po].cl ] ].l  + z[ cur ].l  > L ) tr[po].cl ++ ;
    if( tr[po].cl <= tr[po].cr ) if( ax==-1 || f[ax]>f[ tr[po].cover[ tr[po].cl ] ] ) ax = tr[po].cover[ tr[po].cl ] ;
    if( l==_l && r==_r ) {
        while( tr[po].pl<=tr[po].pr && -z[ tr[po].pass[ tr[po].pl ] ].l + z[ cur ].l  > L ) tr[po].pl ++ ;
        if( tr[po].pl <= tr[po].pr ) if( ax==-1 || f[ax]>f[ tr[po].pass[ tr[po].pl ] ] ) ax = tr[po].pass[ tr[po].pl ] ;
        return ;
    }     
    if( _r<=m ) Ask( _l , _r , l , m , po * 2 ) ;
        else if( _l>m ) Ask( _l , _r , m+1 , r , po * 2 + 1 ) ;
            else Ask( _l , m , l , m , po * 2 ) , Ask( m+1 , _r , m+1 , r , po * 2 + 1 ) ;
}

void BuildTree( int l , int r , int po ) {
    tr[po].cl = tr[po].pl = 1 , tr[po].cr = tr[po].pr = 0 ;
    tr[po].cover.psb( 0 ) , tr[po].pass.psb( 0 ) ;
    tr[po].cover.psb( 0 ) , tr[po].pass.psb( 0 ) ;
    if( l==r ) return ;
    BuildTree( l , l + r >> 1 , po << 1 ) ;
    BuildTree( ( r + l >> 1 ) + 1 , r , ( po << 1 ) + 1 ) ;
}

int ans[N] ;

int main() {
    scanf("%d%d",&n,&L ) ;
    --n  ;
    for( i=1 ; i<=n ; i++ ) {
        scanf("%d%d%d",&z[i].x,&z[i].y,&z[i].l ) ;
        a[i*2-1].va = z[i].x ;
        a[i*2].va = z[i].y ;
        a[i*2-1].ps = a[i*2].ps = i ;
        a[i*2-1].x = 1 ;
        z[i].ps = i ;
    }
    sort( a+1 , a+1+n*2 , cmp ) ;
    int disva = 0 ;
    a[0].va = -1 ;
    for( i=1 ; i<=n*2 ; i++ )   {
        if( a[i].va!=a[i-1].va ) ++ disva ;
        if( a[i].x ) z[ a[i].ps ].x = disva ;
            else z[ a[i].ps ].y = disva ;
    }
    sort( z+1 , z+1+n , _cmp ) ;
    BuildTree( 1 , disva , 1 ) ;
    tr[1].cover.psb( 0 ) ;
    tr[1].cl = tr[1].cr = 1 ;
    f[0] = 0 ;
    for( i=1 ; i<=n ; i++ ) {
        cur = i ;
        ax = -1 ;
        Ask( z[i].x , z[i].y , 1 , disva , 1 ) ;
        if( ax == -1 ) {
            ans[ z[i].ps ] = -1 ;
            continue ;
        }
        f[i] = f[ax] + 1 ;
        Ins( z[i].x , z[i].y , 1 , disva , i , 1 ) ;
        ans[ z[i].ps ] = f[i] ;
    }
    for( i=1 ; i<=n ; i++ ) printf("%d\n",ans[i] ) ;
}

DebugLogo
1 删队首不等式符号恰好写反
2 更新区间最小值判断语句符号写反

写比较长的布尔式要慢点写,图快反而得不偿失。

你可能感兴趣的:(物理( 线段树套单调队列 ))