51nod - 1310 Chandrima and XOR

有这样一个小到大排列的无穷序列S：1, 2, 4, 5, 8......，其中任何一个数转为2进制不包括2个连续的1。给出一个长度为N的正整数数组A，A1, A2......An记录的是下标（下标从1开始）。求S[A1] Xor S[A2] Xor S[A3] ..... Xor S[An]的结果（Xor 为异或运算），由于该数很大，输出Mod 1000000007的结果。

 

例如：A = {1, 2, 3}，对应S[1] = 1, S[2] = 2, S[3] = 4。1 Xor 2 Xor 4 = 7。

Input

第1行：1个数N，表示数组A的长度(1 <= N <= 50000)。
第2 - N + 1行：每行一个数，对应数组A的元素A[i]（1 <= A[i] <= 10^18)。

Output

输出一个数，S[A1] Xor S[A2] Xor S[A3] ..... Xor S[An]的结果Mod 1000000007。

Input示例

3
1
2
3

Output示例

7

思路：

以每一个2的整数次方作为分割点，把这个数列分割成很多块。 有1个数在[1,2)，1个数在[2,4)，2个数在[4,8)，3个数在[8,16)…它是按斐波那契数列的和递增。

#include <iostream>
using namespace std;

typedef long long ll;
const int MAXF = 88;
const int MOD = 1e9 + 7;

long long Fib[MAXF] = {1, 1, 2};
long long Pow[MAXF] = {1, 2, 4};

int XOR[MAXF];
int maxPos = 0;

void init()
{
    for (int i = 3; i < MAXF; i++)
    {
        Fib[i] = Fib[i - 1] + Fib[i - 2];
        Pow[i] = Pow[i - 1] * 2 % MOD;
    }
}

ll BToD()
{
    ll res = 0;
    for (int i = 0; i < maxPos; i++)
    {
        res = (res + XOR[i] * Pow[i]) % MOD;
    }
    return res;
}

int main(int argc, const char * argv[])
{
    init();

    int N;
    cin >> N;

    ll s;
    for (int i = 0; i < N; i++)
    {
        cin >> s;
        while (s != 0)
        {
            int l = 1, r = MAXF;
            while (l < r - 1)
            {
				int mid = l + (r - l) / 2;
				if (Fib[mid] > s)
				{
					r = mid;
				}
				else
				{
					l = mid;
				}
            }
			if (Fib[r] <= s)
			{
				l = r;
			}
            XOR[l - 1] ^= 1;
            maxPos = max(maxPos, l);
            s -= Fib[l];
        }
    }

    cout << BToD() << '\n';

    return 0;
}