poj 3061 Subsequence

#include <iostream>
#include <algorithm>
#include <cmath>
#include <stdio.h>
#include <cstring>
#include <string>
#include <cstdlib>
#include <queue>
#include <stack>
#include <set>
#include <vector>
#include <map>
#include <list>
#include <iomanip>
 #include <fstream>
using namespace std;
typedef long long ll;
const int maxn=1000009;
ll a[maxn];
int main()
{
    //初始化前缀和数组a[maxn]
    memset(a,0,sizeof(a));
    int cas;
    scanf("%d",&cas);    //读入测试样例组数 
        ll n,sum,ans,s,pos,l;
    while(cas--)
    {
        ans=100005;    sum=0;
        scanf("%I64d%I64d",&n,&s);
        for(int i=1;i<=n;++i)
        {
            scanf("%I64d",&a[i]);
            sum+=a[i];
            a[i]+=a[i-1];//一边读一边计算前缀和 
        }
        if(sum<s){//如果所有的数相加都没有要求的数字高 那么直接跳过这个样例 
            printf("0\n");
            continue;
        }
    
        
        for(int i=1;i<=n;++i)
        {
            l=i-1;sum=0;pos=1;//倍增的思想，每次向右移动 pow(pos,2)的样子 
            while(pos)//当它移动到步长为0时 
            {
                while(l+pos>n) pos>>=1;//因为是倍增，所以要担心它是否越界超出数组的实际长度，之前用if，还不够！！因为它减半一次以后仍然有可能超出范围必须用while 
                if(sum+a[l+pos]-a[l]<s)//之前的和sum加上之后区间的和如果仍然小于要求的数 
                {
                    sum+=a[l+pos]-a[l];//继续加上这段区间的和 
                    l+=pos;//位置倍增，原来的右区间现在变成左区间 
                    pos<<=1;//步长左移，相当于乘2 
                }                                  
                else{//发现sum加上之后区间的和大于等于要求的数，符合题目要求 
                    ans=min(ans,l+pos-i+1);//那么答案就是原来答案与这整段区间的长度的较小值 
                    pos>>=1;//因为步长太大了，所以pos步长需要右移一位，相当于除以2 
                }
            }
            //我一开始还在想，它假如pos除以2，回到上一步，上一步pos又乘2，不是在两者之间反复横跳吗...然后用草稿一步一步算以后发现
            //它只是左端点不动，右边慢慢压缩pos 
        }
        if(ans>100000)//发现答案长度超过数组实际长度，这不科学 
            printf("0\n");                                                                       
        else{
            printf("%d\n",ans);
        }
    }

    return 0;
}