[数位DP] 好数

对于一个正整数X,如果把X化成二进制数后,如果X的二进制数至少有三个连续的1或者至少有3个连续的0(不能有前导0),那么X就是“好数”。例如8就是“好数”,因为8对应的二进制数是1000,有三个连续的0。整数15也是“好数”,因为15对应的二进制数是1111,也有三个连续的1。整数27就不是“好数”,因为27对应的二进制数是11011,既没有连续的三个1也没有连续三个0。

现在给出两个整数Low、UP,求Low和UP范围内有多少个“好数”。


#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define cl(x) memset(x,0,sizeof(x))
using namespace std;
typedef long long ll;

inline char nc(){
  static char buf[100000],*p1=buf,*p2=buf;
  if (p1==p2) { p2=(p1=buf)+fread(buf,1,100000,stdin); if (p1==p2) return EOF; }
  return *p1++;
}

inline void read(int &x){
  char c=nc(),b=1;
  for (;!(c>='0' && c<='9');c=nc()) if (c=='-') b=-1;
  for (x=0;c>='0' && c<='9';x=x*10+c-'0',c=nc()); x*=b;
}

const int N=35;

int w[N],len;
int f[N][2][4];
int f1[N];
int ans[N];

inline int Solve(ll n){
  cl(f); cl(f1); cl(ans);
  if (n<=0) return 0;
  n++; len=0; int last,cnt;
  while (n) w[++len]=n&1,n>>=1; reverse(w+1,w+len+1);
  f1[1]=1;
  last=1; cnt=1;
  int flag=0;
  for (int i=2;i<=len;i++){
    if (w[i]==1){
      if (!flag){
	if (last==1)
	  f[i][0][1]++;
	else{
	  if (cnt==1)
	    f[i][0][2]++;
	  else
	    ans[i]++;
	}
      }else
	ans[i]++;
    }
    if (w[i]==last) cnt++; else last=w[i],cnt=1;
    if (cnt==3) flag=1;
    ans[i]+=ans[i-1]*2;
    //d==0
    f[i][0][1]+=f[i-1][1][1]+f[i-1][1][2];
    f[i][0][2]+=f[i-1][0][1];
    f1[i]+=f1[i-1];
    ans[i]+=f[i-1][0][2];
    //d==1
    f[i][1][1]+=f[i-1][0][1]+f[i-1][0][2]+f1[i-1];
    f[i][1][2]+=f[i-1][1][1];
    ans[i]+=f[i-1][1][2];
  }
  return ans[len];
}

int main(){
  int l,r;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  read(l); read(r);
  printf("%d\n",Solve(r)-Solve(l-1));
  return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值