进化序列 jzoj3518 倍增

Description

---

Abathur采集了一系列Primal Zerg 的基因样本，这些基因构成了一个完整的进化链。为了方便，我们用A0,A1…An-1 这n 个正整数描述它们。

一个基因Ax 可以进化为序列中在它之后的基因Ay。这个进化的复杂度，等于Ax | Ax+1…| Ay的值，其中| 是二进制或运算。

Abathur 认为复杂度小于M 的进化的被认为是温和的。它希望计算出温和的进化的对数。

Input

---

第一行包含两个整数n,m。

接下来一行包含A0,A1…An-1 这n 个正整数，描述这n 个基因。

Output

---

第一行包含一个整数，表示温和的进化的对数。

Data Constraint

---

对于30% 的数据，1 <= n <=1000。

对于100% 的数据，1 <= n<= 100000，0 <= m <= 2^30，1<= Ai<= 2^30。

Solution

---

读题还是要仔细点，一开始把或(or)运算看成是异或(xor)了

已知$A\ or\ B = C$，那么不难得到$C\geq A$且$C\geq B$

嗯，也就是说我们由一个数字向右作或运算时所得到的结果是不减的。那么就考虑用二分＋区间查询或者是倍增。这里的区间查询用线段树比较好写

就酱

Code

---

#include <stdio.h>
#include <math.h>
#define rep(i, st, ed) for (int i = st; i <= ed; i += 1)
#define drp(i, st, ed) for (int i = st; i >= ed; i -= 1)
#define ll long long
#define N 1 << 17
int f[N][18], t[N];
int main(void){
    freopen("evolve.in", "r", stdin);
    freopen("evolve.out", "w", stdout);
    int n, m;
    scanf("%d%d", &n, &m);
    rep(i, 1, n){
        scanf("%d", &t[i]);
        f[i][0] = t[i];
    }
    int lim = log2(n);
    rep(j, 1, lim){
        rep(i, 1, n){
            f[i][j] = f[i][j - 1] | f[i + (1 << (j - 1))][j - 1];
        }
    }
    ll ans = 0;
    rep(i, 1, n){
        int now = i;
        int tot = 0;
        drp(j, lim, 0){
            if ((now + (1 << j) <= n + 1) && ((tot | f[now][j]) <= m)){
                tot |= f[now][j];
                now += 1 << j;
            }
        }
        // printf("%d, %d\n", i, now);
        if (now != i + 1){
            ans += now - i - 1;
            // printf("%d, %d\n", i, now - 1);
        }
    }
    printf("%lld\n", ans);
    return 0;
}