HDU 5358 First One 数学+尺取法

多校的题，摆明了数学题，但是没想出来，蠢爆了，之前算了半天的s[i][j]的和，其实是积。其实比赛的时候我连log(s[i][j])+1是s[i][j]的位数都没看出来，说出来都丢人。

知道了这个之后，就枚举二进制数的每一位，因为元素都是非负数，所以sum数组是非降的，这里用到了尺取法，之前也是听说过，应该是做过吧，不太记得了。

因为[2^k-1,2^k)的位数是k，枚举时，固定左端点，在sum数组找到最小的大于等于2^k-1，最大的小于2^k的点，这中间的点和左端点的s[i][j]就对于当前的k满足条件了，就把这些答案加到答案中，复杂度就是O(nlogn)

其实我的代码并没有AC，一直是超时，但是我做了一些极限的数据在本机上也是秒出，当然是我程序的原因，不过我现在还不知道怎么回事，花了这么多时间了，就这样吧。

哈哈原来是编译器的问题，之前交的C++，刚交了发G++就A了，

 #include <iostream>
 #include <cstdio>
 #include <fstream>
 #include <algorithm>
 #include <cmath>
 #include <deque>
 #include <vector>
 #include <queue>
 #include <string>
 #include <cstring>
 #include <map>
 #include <stack>
 #include <set>
 #define LL long long
 #define INF 0x3f3f3f3f
 #define OPEN_FILE
 #define MAXN 100005
 using namespace std;
 LL sum[MAXN], p[MAXN];
 LL ans;
 int main()
 {
 #ifdef OPEN_FILE
     freopen("in.txt", "r", stdin);
     //freopen("out.txt", "w", stdout);
 #endif // OPEN_FILE
     int T;
     scanf("%d", &T);
     p[] = ;
     p[] = ;
     for(int i = ; i <= ; i++){
         p[i] = p[i - ] * ;
     }
     for(int cas = ; cas <= T; cas++){
         int n;
         scanf("%d", &n);
         memset(sum, , sizeof(sum));
         int x;
         for(int i = ; i <= n; i++){
             scanf("%d", &x);
             sum[i] = sum[i - ] + x;
         }
         ans = ;
         for(int i = ; i <= ; i++){
             LL left = , right = ;
             for(int j = ; j <= n; j++){
                 left = max((LL)j, left);
                 right = max((LL)j, right);
                 while(left <= n && sum[left] - sum[j - ] < p[i]){
                     left++;
                 }
                 //int right = left;
                 while(right <= n && sum[right] - sum[j - ] < p[i + ]){
                     right++;
                 }
                 right--;
                 if(right > n){
                     right = n;
                 }
                 if(left <= right){
                     ans += ((((left + right)*(right - left + )) / ) + (right - left + ) * j)  * (LL)(i + );
                 }
             }
             if(p[i + ] > sum[n]){
                 break;
             }
         }
         printf("%I64d\n", ans);
     }
 }