此题100分的解法就是先预处理出每个点的下一个点之后倍增就好了。其实并没有太大难度。
pbihao用双向链表写过了此题。在本地上我treap狂操他,but在rqnoj上,我依靠反复提交才A掉此题(最后一组数据为999ms)。treap好看脸啊。
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 100005
#define inf (1LL<<60)
#define ll long long
#define pli pair<long long,int>
#define pll pair<long long,long long>
#define FI "drive"
int nxt[N][17],da[N][17],db[N][17],na[N],nb[N],H[N];
int ina; char inc,inb[1<<16],*ins=inb,*ine=inb; bool insign;
#define getc() ((ins==ine&&(ine=(ins=inb)+fread(inb,1,1<<16,stdin),ins==ine))?EOF:*ins++)
inline int geti() {
for(insign=false;(inc=getc())<'0'||inc>'9';insign|=inc=='-');
for(ina=inc-'0';(inc=getc())>='0'&&inc<='9';ina=(ina<<3)+(ina<<1)+inc-'0');
return insign?-ina:ina;
}
ll iabs_t;
#define iabs(a) ((iabs_t=(a))<0?-iabs_t:iabs_t)
inline void cal(int s,int x,int &t1,int &t2) {
t1=t2=0;
for(int i=14;~i;--i)
if(nxt[s][i]&&t1+da[s][i]+t2+db[s][i]<=x){
t1+=da[s][i],t2+=db[s][i];
s=nxt[s][i];
}
}
int totnode,ch[2][N],r[N],root,id[N]; ll v[N];
inline void rot(int &u,int d) {
int t=ch[d^1][u]; ch[d^1][u]=ch[d][t];
ch[d][t]=u; u=t;
}
void Ins(int &u,const ll&val,int idx) {
if(!u){
u=++totnode;v[u]=val,r[u]=rand(),id[u]=idx;
ch[0][u]=ch[1][u]=0; return;
}
int k=val<v[u]?0:1;
Ins(ch[k][u],val,idx);
if(r[ch[k][u]]>r[u]) rot(u,k^1);
}
inline int find(int val) {
int u=root;
while(u) {
if(v[u]==val) return u;
u=ch[v[u]<val][u];
}return 0;
}
int pred(int u,const ll&val) {
if(!u) return -1;
if(val<=v[u]) return pred(ch[0][u],val);
int t=pred(ch[1][u],val); if(t<0) t=u;
return t;
}
int succ(int u,const ll&val) {
if(!u) return -1;
if(val>=v[u]) return succ(ch[1][u],val);
int t=succ(ch[0][u],val); if(t<0) t=u;
return t;
}
int main() {
freopen(FI".in","r",stdin); freopen(FI".out","w",stdout);
int n,i,j,s,x,m; pll t[5]; root=totnode=0;
for(n=geti(),i=1;i<=n;++i) H[i]=geti(); srand(time(NULL));
Ins(root,inf,0); Ins(root,-inf,0);
for(i=n;i;--i) { Ins(root,H[i],i);
t[1].second=v[pred(root,H[i])],t[2].second=v[succ(root,H[i])];
if(t[1].second>-inf) t[3].second=v[pred(root,t[1].second)];
else t[3].second=-inf;
if(t[2].second<inf) t[4].second=v[succ(root,t[2].second)];
else t[4].second=inf;
for(j=1;j<5;++j) t[j].first=iabs(t[j].second-H[i]);
sort(t+1,t+5);
na[i]=id[find(t[2].second)],nb[i]=id[find(t[1].second)];
(na[i])?da[i][0]=t[2].first:1;
(nb[na[i]])?db[i][1]=iabs(H[na[i]]-H[nb[na[i]]]),da[i][1]=t[2].first:1;
nxt[i][0]=na[i],nxt[i][1]=nb[na[i]];
for(j=2;j<15;++j)
if(nxt[nxt[i][j-1]][j-1]) {
nxt[i][j]=nxt[nxt[i][j-1]][j-1];
da[i][j]=da[i][j-1]+da[nxt[i][j-1]][j-1];
db[i][j]=db[i][j-1]+db[nxt[i][j-1]][j-1];
}else break;
}
int t1,t2; x=geti(); double mi=inf,_t;
for(i=1;i<=n;++i) {
cal(i,x,t1,t2);
_t=t2?(double)t1/t2:1e60;
if(_t<mi||(_t==mi&&H[j]<H[i])) j=i,mi=_t;
}printf("%d
",j);
for(m=geti();m;--m) {
s=geti(),x=geti();
cal(s,x,t1,t2);
printf("%d %d
",t1,t2);
}return 0;
}
update:今天写了一份用指针的,快了好多。
#include <ctime>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
using namespace std;
#define inf (1LL<<60)
#define ll long long
#define pli pair<long long,int>
#define pll pair<long long,long long>
char INB[1<<16],*INS=INB,*INE=INB;
#define getc() (INS==INE&&(INE=(INS=INB)+fread(INB,1,1<<16,stdin),INS==INE)?EOF:*INS++)
inline int geti() {
register int a; register char c,f=0;
while(c=getc(),c<'0')f|=c=='-';a=c-'0';
while(c=getc(),'-'<c)a=(a<<3)+(a<<1)+c-'0';
return f?-a:a;
}
ll ibs_t;
#define ibs(a) ((ibs_t=(a))<0?-ibs_t:ibs_t)
/*template*/
#define N 100005
int nxt[17][N],da[17][N],db[17][N],na[N],nb[N],H[N];
inline void cal(int s,int x,int&t1,int&t2) {
t1=t2=0;
for(register int i=16;~i;--i)
if(nxt[i][s]&&t1+da[i][s]+t2+db[i][s]<=x) {
t1+=da[i][s];
t2+=db[i][s];
s=nxt[i][s];
}
}
struct node{
int r,id; ll v;
node*ch[2];
}CD[N],*cd=CD,*root=NULL;
inline void rot(node*(&u),const int&d) {
node*t=u->ch[d^1];u->ch[d^1]=t->ch[d];
t->ch[d]=u; u=t;
}
void ins(node*(&u),const ll&val,const int&idx) {
if(!u) {
u=cd++; u->v=val,u->id=idx,u->r=rand();
u->ch[0]=u->ch[1]=NULL; return;
}
int k=val<u->v?0:1;
ins(u->ch[k],val,idx);
if(u->ch[k]->r>u->r) rot(u,k^1);
}
inline node*find(const ll&val) {
node*u=root;
while(u) {
if(u->v==val) return u;
u=u->ch[u->v<val];
}return NULL;
}
inline node*pred(const ll&val) {
node*u=root,*r=NULL;
while(u){
if(u->v<val) r=u,u=u->ch[1];
else u=u->ch[0];
}
return r;
}
inline node*succ(const ll&val) {
node*u=root,*r=NULL;
while(u) {
if(val<u->v) r=u,u=u->ch[0];
else u=u->ch[1];
}
return r;
}
int main() {
register int n,i,j,s,x,m; pll t[5];
for(n=geti(),i=1;i<=n;++i) H[i]=geti();
srand(time(NULL));
ins(root,inf,0); ins(root,-inf,0);
for(i=n;i;--i) {
ins(root,H[i],i);
t[1].second=pred(H[i])->v,t[2].second=succ(H[i])->v;
if(t[1].second>-inf) t[3].second=pred(t[1].second)->v;
else t[3].second=-inf;
if(t[2].second<inf) t[4].second=succ(t[2].second)->v;
else t[4].second=inf;
for(j=1;j<5;++j) t[j].first=ibs(t[j].second-H[i]);
sort(t+1,t+5);
if(nxt[0][i]=na[i]=find(t[2].second)->id) da[0][i]=t[2].first;
nb[i]=find(t[1].second)->id;
if(nxt[1][i]=nb[na[i]]) db[1][i]=ibs(H[na[i]]-H[nb[na[i]]]),da[1][i]=t[2].first;
for(j=2;j^17;++j)
if(nxt[j][i]=nxt[j-1][nxt[j-1][i]]) {
da[j][i]=da[j-1][i]+da[j-1][nxt[j-1][i]];
db[j][i]=db[j-1][i]+db[j-1][nxt[j-1][i]];
}else break;
}
int t1,t2; x=geti(); double mi=inf,_t;
for(i=1;i<=n;++i) {
cal(i,x,t1,t2);
_t=t2?(double)t1/t2:1e60;
if(_t<mi||(_t==mi&&H[j]<H[i])) j=i,mi=_t;
}printf("%d
",j);
for(m=geti();m;--m) {
s=geti(),x=geti();
cal(s,x,t1,t2);
printf("%d %d
",t1,t2);
}return 0;
}