玲珑杯-1112 （倍增+二分）

DESCRIPTION

小 L 有一串 $Q$ 个数排成一排，从左往右是 $A_1$ 开始一直到 $A_Q$，它一直很喜欢这串数。

有一天，小 L 去高考报名了，小 J 和小 N 过来用这串数玩游戏。

一开始，小 N 手上有一个空序列，定义变量$C:=0$。小 J 从左往右依次把序列中的数取出，放到小 N 手上序列的最右边，如果放完之后，整个序列的混乱度超过了 $M$，小 N 就会把手上的所有数全部扔掉，然后令 $C:=C+1$。

定义一个长度为 $K$ 序列的混乱度 $S=\sum _{i=1}^K{B}_i\times V_i$，其中 $B_i$ 表示序列中第 $i$ 小的数，$V$ 为一个给定的序列。

小 J 和小 N 想知道，加入每个数之后，$C$ 是多少。

INPUT

第一行两个整数 $Q$， $M$接下来一行 $Q$个整数， 第 $i$个整数表示 $A_i$接下来一行 $Q$个整数， 第 $i$个整数表示 $V_i$

OUTPUT

一行 $Q$个整数， 第 $i$个整数表示加入第 $i$个数之后的 $C$， 相邻两个整数之间用空格隔开， 注意最后一个数后不要输出空格输出完请换行

SAMPLE INPUT

5 11 3 2 5 41 1 1 1 1

SAMPLE OUTPUT

0 1 2 3 4

HINT

$1\le Q\le 300000,0\le M\le {10}^{18},1\le A_i\le {10}^8,1\le V_i\le {10}^4$对于所有的 $1\le i\le Q-1,V_i\le V_{i+1}$

SOLUTION

http://www.ifrog.cc/acm/problem/1112

题解：倍增+二分

先说说自己在做这个题的时候卡在哪了吧。我想应该大部分就是卡在这么处理在添加的时候怎么确定第k小的数。想到这肯定去想用排序，但是排序又肯定会超时。

那么现在就来通过这个题来了解一下倍增二分的奥妙。那么根据这个题的题意，每次加一个数其实他的混乱度都是在增加的，然而添加又都是从左边开始一个一个的添加的。那么可想肯定有一个右边界使得在这个区域内的值得混乱度超过m。现在问题转化为了怎么确定右边界，当然最简单的办法就是一个一个的取寻找右边界，这个和上面讲个超时的做法是一样的，毫无改进。那么既然用到了查找寻值，肯定能想到用二分去查找呀。那么既然是二分就要用到求右边界的L和R,如果只是从左边界到n来二分排序的也会超时的。

现在问题就转化为了怎么来快速的确定一个右边界的大概区间然后用二分判断呢。这里就是用倍增的地方了。用倍增找一个混乱度大于m的r2边界，在找r2的时候，要记得更新r1边界。具体代码中也有解释了。

#include <iostream>
#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

int pow1[30];
ll a[300300],v[300300],temp[300300];
int ans[300300];
void init()
{
    pow1[0]=1;
    for(int i=1;i<25;i++)
    {
        pow1[i]=pow1[i-1]*2;
    }
}

bool judge(int l,int r,ll m)
{
    int i,j;
    for( i=0,j=l;j<=r;j++)
    {
        temp[i++]=a[j];
    }
    sort(temp,temp+i);///排序求这个区间的K值，看是否大于m
    ll sum=0;
    for(int j=0;j<i;j++)
    {
        sum+=temp[j]*v[j];
    }
    return sum>m;
}
int erfen(int r1,int r2,int l,ll m)
{
    while(r1<r2)
    {
        int mid=(r1+r2)/2;
        if(judge(l,mid,m))
        {
            r2=mid;
        }
        else
        {
            r1=mid+1;
        }
    }
    return r1;///确切的右端点
}
int cheak(int l,int n,ll m)
{
    int r1=l,r2=l;///寻找大于m的右端点一定在区间r1-r2内
    for(int i=0;;r2+=pow1[i],i++)///倍增法求大概区间
    {
        if(r2>n)r2=n;
        if(judge(l,r2,m))break;
        r1=r2+1;///注意r1的伟大作用
        if(r2==n)break;
    }
    return erfen(r1,r2,l,m);///然后利用二分在r1-r2范围内找出确切的第一个大于m的位置
}
void work(int n,ll m)
{
    for(int i=0;i<n;)///左端点
    {
        int j=cheak(i,n,m);///右端点
        ans[j]=ans[i==0?i:i-1]+1;///大于M的断点
        i=j+1;///左端点的更新
    }
    for(int i=1;i<n;i++)
    {
        ans[i]=max(ans[i],ans[i-1]);///求出所有的点的结果
    }
    return ;
}
int main()
{
    init();///初始化倍增值
    int n;
    ll m;
    scanf("%d%lld",&n,&m);
    for(int i=0;i<n;i++)scanf("%lld",&a[i]);
    for(int i=0;i<n;i++)scanf("%lld",&v[i]);
    memset(ans,0,sizeof(ans));///结果的存储

    work(n,m);

    for(int i=0;i<n;i++)
    {
        if(i==0)
        printf("%d",ans[i]);
        else printf(" %d",ans[i]);
    }
    printf("\n");
    return 0;
}