bzoj2844: albus就是要第一个出场

Description

已知一个长度为n的正整数序列A（下标从1开始），令 S = { x | 1 <= x <= n }, S 的幂集2^S定义为S 所有子集构成的集合。

定义映射 f : 2^S -> Z

f(空集) = 0
f(T) = XOR A[t] , 对于一切t属于T

现在albus把2^S中每个集合的f值计算出来，从小到大排成一行，记为序列B（下标从1开始）。给定一个数，那么这个数在序列B中第1次出现时的下标是多少呢？

Input

第一行一个数n, 为序列A的长度。
接下来一行n个数，为序列A，用空格隔开。
最后一个数Q，为给定的数.

Output

共一行，一个整数，为Q在序列B中第一次出现时的下标模10086的值.

Sample Input

3
1 2 3
1

Sample Output

3

【Hint】
样例解释：
N = 3, A = [1 2 3]
S = {1, 2, 3}
2^S = {空, {1}, {2}, {3}, {1, 2}, {1, 3}, {2, 3}, {1, 2, 3}}
f(空) = 0
f({1}) = 1
f({2}) = 2
f({3}) = 3
f({1, 2}) = 1 xor 2 = 3
f({1, 3}) = 1 xor 3 = 2
f({2, 3}) = 2 xor 3 = 1
f({1, 2, 3}) = 0
所以
B = [0, 0, 1, 1, 2, 2, 3, 3]

HINT

数据范围：

1 <= N <= 10,0000

其他所有输入均不超过10^9

见http://blog.csdn.net/popoqqq/article/details/39829237

code：

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cmath>
 4 #include<cstring>
 5 #include<algorithm>
 6 #define maxn 100005
 7 #define mod 10086
 8 using namespace std;
 9 char ch;
10 int n,m,q,a[maxn],ans;
11 bool ok;
12 void read(int &x){
13     for (ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if (ch=='-') ok=1;
14     for (x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());
15     if (ok) x=-x;
16 }
17 void gauss(){
18     int i,j,k;
19     for (i=30,k=1;i>=0;i--){
20         for (j=k;j<=n&&!(a[j]&(1<<i));j++);
21         if (j<=n){
22             swap(a[k],a[j]);
23             for (j=1;j<=n;j++) if (j!=k&&(a[j]&(1<<i))) a[j]^=a[k];
24             k++;
25         }
26     }
27     m=k-1;
28 }
29 int ksm(int a,int b){
30     int t=1;
31     for (;b;b>>=1){if (b&1) t=t*a%mod; a=a*a%mod;}
32     return t;
33 }
34 int main(){
35     read(n);
36     for (int i=1;i<=n;i++) read(a[i]);
37     gauss();
38     read(q),ans=0;
39     int now=0;
40     for (int i=1;i<=m;i++) if ((now^a[i])<q) ans+=ksm(2,m-i),ans%=mod,now^=a[i];
41     if (q) ans++,ans%=mod;
42     ans=ans*ksm(2,n-m)%mod;
43     ans++,ans%=mod;
44     printf("%d
",ans);
45     return 0;
46 }