Problem Description
众所周知,度度熊喜欢的字符只有两个:B和D。
今天,它发明了一种用B和D组成字符串的规则:
S(1)=B
S(2)=BBD
S(3)=BBDBBDD
…
S(n)=S(n−1)+B+reverse(flip(S(n−1))
其中,reverse(s)指将字符串翻转,比如reverse(BBD)=DBB,flip(s)指将字符串中的B替换为D,D替换为B,比如flip(BBD)=DDB。
虽然度度熊平常只用它的电脑玩连连看,这丝毫不妨碍这台机器无与伦比的运算速度,目前它已经算出了S(21000)的内容,但度度熊毕竟只是只熊,一次读不完这么长的字符串。它现在想知道,这个字符串的第L位(从1开始)到第R位,含有的B的个数是多少?
今天,它发明了一种用B和D组成字符串的规则:
S(1)=B
S(2)=BBD
S(3)=BBDBBDD
…
S(n)=S(n−1)+B+reverse(flip(S(n−1))
其中,reverse(s)指将字符串翻转,比如reverse(BBD)=DBB,flip(s)指将字符串中的B替换为D,D替换为B,比如flip(BBD)=DDB。
虽然度度熊平常只用它的电脑玩连连看,这丝毫不妨碍这台机器无与伦比的运算速度,目前它已经算出了S(21000)的内容,但度度熊毕竟只是只熊,一次读不完这么长的字符串。它现在想知道,这个字符串的第L位(从1开始)到第R位,含有的B的个数是多少?
Input
第一行一个整数T,表示T(1≤T≤1000) 组数据。
每组数据包含两个数L和R(1≤L≤R≤1018) 。
每组数据包含两个数L和R(1≤L≤R≤1018) 。
Output
对于每组数据,输出S(21000)表示的字符串的第L位到第R位中B的个数。
Sample Input
3 1 3 1 7 4 8
Sample Output
2 4 3调了好久,今天状态不好,失败这题就是每次拆两半找规律的。#include<map> #include<set> #include<queue> #include<stack> #include<cmath> #include<cstdio> #include<bitset> #include<string> #include<vector> #include<cstring> #include<iostream> #include<algorithm> #include<functional> using namespace std; typedef __int64 LL; const int low(int x) { return x&-x; } const int INF = 0x7FFFFFFF; const int mod = 1e9 + 7; const int maxn = 2e5 + 10; int T; LL l, r; LL solve(LL x) { if (!x) return 0; LL i, ans = 0, k = 1; for (i = 1; i <= x; i = i << 1 | 1); for (i >>= 1; i; i >>= 1) { if (x == i) ans += i / 2 + 1, x -= i; if (x > i) { ans += i / 2 + 1 + k; if (x > i + i / 2 + 1) k = 0; else k = 1; x -= i + 1; } } return ans; } int main() { /*int dp[maxn]; int x[maxn]; x[1] = 1; for (int i = 1; i < 4096; i = i * 2 + 1) { x[i + 1] = 1; for (int j = i; j; j--) x[i + i + 2 - j] = x[j] ^ 1; } for (int i = 1; i <1024; i++) { dp[i] = dp[i - 1] + x[i]; if (solve(i) != dp[i]) printf("%d %d %d\n", i, dp[i], solve(i)); }*/ scanf("%d", &T); while (T--) { scanf("%I64d%I64d", &l, &r); LL x = solve(r) - solve(l - 1); cout << x << endl; } return 0; }