JZOJ 3518. 【NOIP2013模擬11.6A組】進化序列(evolve)

3518. 【NOIP2013模擬11.6A組】進化序列(evolve)

(File IO): input:evolve.in output:evolve.out
Time Limits: 1000 ms Memory Limits: 262144 KB

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

第一行包含一個整數，表示溫和的進化的對數。

Sample Input

4 6
1 3 5 1

Sample Output

2

Data Constraint

對於30% 的數據，1 <= n <=1000。
對於100% 的數據，1 <= n<= 100000，0 <= m <= 2^30，1<= Ai<= 2^30。

題解

兩種解法：
一種是類似RMQ的倍增算法
另一種是……（我也不知道，有點像單調隊列……反正是隊列）

我用的是第二種，雖然不知道叫什麼算法，但是也講講

用隊列que表示當前選擇
a [ i ] 表示第i位的1的個數
num表示當前進化複雜度

如果當前值x，x|num>m 就把隊首丟掉

代碼

#include<cstdio>
#include<queue>
#include<cmath>
#define lowbit(a) ((a)&-(a))
#define qu(q) ((long)log2(lowbit(q)))
#define N 32
using namespace std;
queue<long>que;
long a[N];
int main()
{   long n,m,i,q,num,x,ans=0;
    freopen("evolve.in","r",stdin);
    freopen("evolve.out","w",stdout);
    scanf("%ld%ld",&n,&m);
    num=0;
    for(i=1;i<=n;i++){
        scanf("%ld",&x);
        while((num|x)>m){
            for(q=que.front();q;q^=lowbit(q)){
                a[qu(q)]--;
                if(!a[qu(q)])
                    num^=lowbit(q);
            }
            que.pop();
        }
        num|=x;
        que.push(x);
        for(q=x;q;q^=lowbit(q))
            a[qu(q)]++;
        ans+=que.size()-1;
    }
    printf("%ld\n",ans);
    return 0;
}