题目链接:
BD String
Time Limit: 2000/1000 MS (Java/Others)
Memory Limit: 65536/65536 K (Java/Others)
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≤10^18) 。
每组数据包含两个数L和R(1≤L≤R≤10^18) 。
Output
对于每组数据,输出
S(21000)表示的字符串的第L位到第R位中B的个数。
Sample Input
3
1 3
1 7
4 8
Sample Output
2
4
3
题意:
思路:
L[i],R[i]分别表示node==i时这个串的左右端点,dp[i][0]表示node==i时B的个数,dp[i][1]表示D的个数;
给一个区间[l,r]进行分治,mid=(L[i]+R[i])/2,如果r<mid容易想,当l>mid时这时就可以把求node==i-1时的[l,r]的B的
个数转化成求node==i-1时的[l',r']的D的个数,在转化成求r'-l'-[l',r']里面的B的个数,递归就可以了,
当l<mid<r的时候就是把前面两个整合一下,详细看代码ba ;
AC代码:
#include <bits/stdc++.h> /* #include <iostream> #include <queue> #include <cmath> #include <map> #include <cstring> #include <algorithm> #include <cstdio> */ using namespace std; #define Riep(n) for(int i=1;i<=n;i++) #define Riop(n) for(int i=0;i<n;i++) #define Rjep(n) for(int j=1;j<=n;j++) #define Rjop(n) for(int j=0;j<n;j++) #define mst(ss,b) memset(ss,b,sizeof(ss)); typedef long long LL; const LL mod=1e9+7; const double PI=acos(-1.0); const int inf=0x3f3f3f3f; const int N=1e5+25; int n; LL L[70],R[70],dp[70][2]; int fun() { L[1]=1; R[1]=1; dp[1][0]=1; dp[1][1]=0; for(int i=2;i<=65;i++) { L[i]=1; R[i]=2*R[i-1]+1; dp[i][0]=dp[i-1][0]+1+dp[i-1][1]; dp[i][1]=dp[i-1][0]+dp[i-1][1]; } } LL dfs(LL l, LL r,int node) { if(l>r||node<1)return 0; if(l==L[node]&&r==R[node])return dp[node][0]; LL mid=(L[node]+R[node])>>1; if(r<=mid) { if(r==mid)return 1+dfs(l,r-1,node-1); else return dfs(l,r,node-1); } else { if(l>=mid) { if(l==mid)return r-l+1-dfs(R[node]-r+1,R[node]-l,node-1); else return r-l+1-dfs(R[node]-r+1,R[node]-l+1,node-1); } else return dfs(l,mid-1,node-1)+1+r-mid-dfs(R[node]-r+1,R[node]-mid,node-1); } } int main() { int t; scanf("%d",&t); fun(); while(t--) { LL l,r; scanf("%I64d%I64d",&l,&r); LL ans=dfs(l,r,62); printf("%I64d\n",ans); } return 0; }