codeforces 215E 求有周期循环的二进制数的个数

101010属于周期循环的二进制，循环次数为3

现在给你十进制的区间，求出有周期循环的数的个数

思路：

    枚举循环的周期的长度，但注意周期长度所包含的因子长度不能重复计算。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
typedef long long LL;
int const N = 70;
LL dp[N],p2[N];
int bit[N],ln;
void pre()
{
     p2[0]=1;
     for(int i=1;i<=64;i++)p2[i]=(p2[i-1]<<1);
}
LL getsum1(LL n)
{
   LL ans=0;
   for(int i=1;i<=n/2;i++)
   {
       if(n%i==0)
       {
          dp[i]=p2[i-1];
          for(int j=1;j<=i/2;j++)if(i%j==0)dp[i]-=dp[j];
          ans+=dp[i];
       }
   }
   return ans;
}
LL getsum2(LL n)
{
   LL num=n;
   for(ln=0;num;bit[++ln]=num%2,num/=2);
   LL ans=0;
   for(int i=1;i<ln;i++)ans+=getsum1(i);
   for(int i=1;i<=ln/2;i++)
   {
       if(ln%i==0)
       {
          dp[i]=0LL;
          LL pre=bit[ln];
          for(int j=1;j<i;j++)
          {
              if(bit[ln-j])dp[i]+=p2[i-j-1];
              pre=pre*2+bit[ln-j];
          }
          LL d=pre;
          for(int j=1;j<ln/i;j++)pre=pre*p2[i]+d;
          dp[i]+=(pre<=n);
          for(int j=1;j<=i/2;j++)if(i%j==0)dp[i]-=dp[j];
          ans+=dp[i];
       }
   }
   return ans;
}
int main()
{
    LL l,r;
    pre();
    while(~scanf("%I64d %I64d",&l,&r))
    {
          l--;
          printf("%I64d\n",getsum2(r)-getsum2(l));
    }
    return 0;
}