高中物理总要留下点什么才好……于是有了这题(其实搞的原题)。
想了想,这题还是回归学物理的原点好了。初中物理的第一课就是声吧。
声波(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;
}