GDOI2016模拟8.19物理

高中物理总要留下点什么才好……于是有了这题（其实搞的原题）。

想了想，这题还是回归学物理的原点好了。初中物理的第一课就是声吧。

声波（Sound Wave或Acoustic Wave）是声音的传播形式。声波是一种机械波，由物体（声源）振动产生，声波传播的空间就称为声场。在气体和液体介质中传播时是一种纵波，但在固体介质中传播时可能混有横波。人耳可以听到的声波的频率一般在20赫兹至20000赫兹之间。

不同动物和不同的人能听到的声波频率都不一样。这里就不上图了。

现在，sg在和妹子玩一个很好玩的游戏。

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

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

这题我们可以考虑题目给的限制，首先对于一个小弟的方案，它是由前面距离l、频率与xi，yi有交集的小弟的最小值+1得到。

那好，我们可以先里散掉距离，由于有频率的限制想得会比较复杂，我们由简入难。

现在假设没有频率的限制，也就是所有小弟只要距离小于等于l，那么都可以到达，这样我们显然可以用单调队列来维护最优方案，每次对于小弟i，若队首与它的距离大于l，则删除（对后面没有影响），反则方案为队首方案+1，再将它插入队尾，若方案i比队尾的要小或等于，则队尾对后面的方案的贡献肯定不会比i大，去除。

那现在有了频率的限制，我们可以把它变成没限制：

用线段树来维护频率范围，每个节点开两个vector来存单调队列（队列的更新和前面的差不多），一个存频率全覆盖该区间的方案，一个存频率没有覆盖完的方案，（若不分开，那么得打标记，而且标记每次不能叠加，要下传，那么空间会很大）这样，我们可以沿小弟i的频率限制往线段树的叶子走，路径上，我们算全覆盖的，到了限制能覆盖整个区间时，我们算没有覆盖完的，记得取最小值。

若能找到方案，ans就是这个了，反之ans=-1。

找到方案后我们要将i加入线段树里，也是分全覆盖和没有覆盖完的来更新。

这是我第一次用vector，感觉还比较好用，没有过多细节，还是比较好打的。

贴代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#define N 250001
#include<vector>
using namespace std;
int n,L,sum,d;
int b[N],a[N][3],ans[N],c[N+N];
struct node{
    vector<int>all,part;
    int allf,partf;
}f[N*8];
void init(){
    scanf("%d %d",&n,&L);
    n--;
    c[0]=n+n;
    for (int i=1;i<=n;i++)
        scanf("%d %d %d",&a[i][0],&a[i][1],&a[i][2]),c[i]=a[i][0],c[i+n]=a[i][1],b[i]=i;
}
bool cmp(int x,int y){
    return a[x][2]<a[y][2];
}
void pre(){
    static int x;
    sort(c+1,c+c[0]+1);
    sort(b+1,b+b[0]+1,cmp);
    x=c[0],c[0]=1;
    for (int i=2;i<=x;i++)
        if (c[c[0]]!=c[i])c[++c[0]]=c[i];
}
void ins(int l,int r,int s,int ll,int rr,int v){
    static int ss;
    if (c[l]>rr||c[r]<ll)return;
    if (ll<=c[l]&&c[r]<=rr){
        ss=f[s].all.size();
        while (ss-1>=f[s].allf&&ans[f[s].all[ss-1]]>=ans[v])ss--,f[s].all.pop_back();
        f[s].all.push_back(v);
        return;
    }
    ss=f[s].part.size();
    while (ss-1>=f[s].partf&&ans[f[s].part[ss-1]]>=ans[v])ss--,f[s].part.pop_back();
    f[s].part.push_back(v);
    ins(l,(l+r)/2,s+s,ll,rr,v),ins((l+r)/2+1,r,s+s+1,ll,rr,v);
}
void find(int l,int r,int s,int ll,int rr,int v){
    static int ss;
    if (c[l]>rr||c[r]<ll)return;
    ss=f[s].all.size();
    while (f[s].allf<ss&&a[f[s].all[f[s].allf]][2]<v)f[s].allf++;
    if (f[s].allf!=ss){
        if (d==-1||ans[f[s].all[f[s].allf]]<d)
            d=ans[f[s].all[f[s].allf]];
    }
    if (ll<=c[l]&&c[r]<=rr){
        ss=f[s].part.size();
        while (f[s].partf<ss&&a[f[s].part[f[s].partf]][2]<v)f[s].partf++;
        if (f[s].partf!=ss){
            if (d==-1||ans[f[s].part[f[s].partf]]<d)
                d=ans[f[s].part[f[s].partf]];
        }
        return;
    }
    find(l,(l+r)/2,s+s,ll,rr,v),find((l+r)/2+1,r,s+s+1,ll,rr,v);
}
void work(){
    static int l,r;
    ins(1,c[0],1,0,c[c[0]],0);
    for (int i=1;i<=n;i++){
        d=-1;
        find(1,c[0],1,a[b[i]][0],a[b[i]][1],a[b[i]][2]-L);
        if (d!=-1){
            ans[b[i]]=d+1;
            ins(1,c[0],1,a[b[i]][0],a[b[i]][1],b[i]);
        }else ans[b[i]]=-1;
    }
}
void write(){
    for (int i=1;i<=n;i++)
        printf("%d\n",ans[i]);
}
int main(){
    freopen("3.in","r",stdin);
    freopen("3.out","w",stdout);
    init();
    pre();
    work();
    write();
    return 0;
}