大概的思路就是用鏈表預處理一下,然後用倍增。
鏈表是因爲:我們每次增減的都是周圍的幾個數。
倍增是因爲:每一次之後的數都是固定的,就是惟一的數指向惟一的數。
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=1000005;
int n,a,b,l,r,j,i;
struct node
{
int h,id,l,r;
}d[maxn];
int f[maxn][35],fa[maxn][35],fb[maxn][35];
int na[maxn],nb[maxn],p[maxn];
int cmp(node x,node y)
{
return x.h<y.h;
}
int ok()
{
if(!l) return 0;
if(!r) return 1;
return d[j].h-d[l].h<=d[r].h-d[j].h;
}
int pd(int a,int b)
{
if(!a) return d[b].id;
if(!b) return d[a].id;
if(d[j].h-d[a].h<=d[b].h-d[j].h) return d[a].id;
return d[b].id;
}
void getab(ll x,int p)
{
a=b=0;
for(int i=19;i>=0;i--)
{
if(f[p][i]&&(ll)(a+b+fa[p][i]+fb[p][i])<=x)
{
a+=fa[p][i];b+=fb[p][i];p=f[p][i];
}
}
if(na[p]&&a+b+fa[p][0]<=x) a+=fa[p][0];
}
void st()
{
for(int i=1;i<=n;i++)
{
f[i][0]=nb[na[i]];
fa[i][0]=abs(d[p[i]].h-d[p[na[i]]].h);
fb[i][0]=abs(d[p[f[i][0]]].h-d[p[na[i]]].h);
}
for(int j=1;j<=19;j++)
{
for(int i=1;i<=n;i++)
{
f[i][j]=f[f[i][j-1]][j-1];
fa[i][j]=fa[i][j-1]+fa[f[i][j-1]][j-1];
fb[i][j]=fb[i][j-1]+fb[f[i][j-1]][j-1];
}
}
}
int main()
{
scanf("%d",&n);
for(int i=1;i<=n;i++) scanf("%d",&d[i].h);
for(int i=1;i<=n;i++) d[i].id=i;
sort(d+1,d+1+n,cmp);
for(int i=1;i<=n;i++) p[d[i].id]=i;
for(int i=1;i<=n;i++) d[i].l=i-1,d[i].r=i+1;
d[1].l=d[n].r=0;
for(int i=1;i<=n;i++)
{
j=p[i];l=d[j].l;r=d[j].r;
if(ok()) nb[i]=d[l].id,na[i]=pd(d[l].l,r);
else nb[i]=d[r].id,na[i]=pd(l,d[r].r);
if(l) d[l].r=r;
if(r) d[r].l=l;
}
st();
int x;scanf("%d",&x);
int m;scanf("%d",&m);
double minv=0x3f3f3f3f;int ans=0;
for(int i=1;i<=n;i++)
{
getab(x,i);
if(b&&1.0*a/b<minv)
{
minv=1.0*a/b;ans=i;
}
}
printf("%d\n",ans);
for(int i=1;i<=m;i++)
{
int j;ll x;scanf("%d%lld",&j,&x);getab(x,j);
printf("%d %d\n",a,b);
}
return 0;
}