hdu4588Count The Carries

链接

去年南京邀请赛的水题，当时找规律过的，看它长得很像数位dp，试了试用数位dp能不能过，d出每位上有多少个1，然后TLE了。。然后用规律优化了前4位，勉强过了。

附数位dp代码及找规律代码。

 #include <iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<queue>
 #include<stack>
 #include<cmath>
 #include<stdlib.h>
 #include<map>
 using namespace std;
 #define N 1000000
 #define LL __int64
 int dp[][][];
 int di[],o[],oo[];
 int s[N][],tg,pp[];
 bool f[N];
 int  dfs(int i,int e,bool z,int k)
 {
     if(i<)  return z;
     if(!e&&dp[i][k][z]!=-) return dp[i][k][z];
     int mk = e?di[i]:;
     int ans = ;
     for(int j = ; j <= mk ; j++)
     {
         int ct =  dfs(i-,e&&(j==mk),(z|((i==k)&&j)),k);
         ans+=ct;
     }
     return e?ans:dp[i][k][z] = ans;
 }
 int cc(int x,int k)
 {
     int p1 = (x-pp[k]+)/pp[k+];
     int qq = p1*pp[k];
     if((x-pp[k]+)%pp[k+]>pp[k]) qq += pp[k];
     else qq+=((x-pp[k]+)%pp[k+]);
     return qq;
 }
 //void cal(int x,int y)
 //{
 //    int g = 0,i;
 //    int tx = x;
 //    memset(o,0,sizeof(o));
 //    if(x<=0)
 //    o[0] = o[1] = 0;
 //    else{
 //        o[0] = (x+1)/2;
 //        o[1] = cc(x,1);
 //        if(x>=3)
 //        o[2] = cc(x,2);
 //        if(x>=7)
 //        o[3] = cc(x,3);
 //    }
 //    while(x>0)
 //    {
 //        di[g++] = x%2;
 //        x/=2;
 //    }
 //    for(i = 4 ; i <= g-1 ; i++)
 //    {
 //        o[i] = dfs(g-1,1,0,i);
 //    }
 //    int t = 0;
 //    LL ans = 0;
 //    if(y<=0)
 //    oo[0] = oo[1] = o[2] = o[3] = 0;
 //    else{
 //        oo[0] = ((y+1)/2);
 //        oo[1] = cc(y,1);
 //        if(y>=3)
 //        oo[2] = cc(y,2);
 //        if(y>=7)
 //        oo[3] = cc(y,3);
 //    }
 //    g = 0;
 //    while(y)
 //    {
 //        di[g++] = y%2;
 //        y/=2;
 //    }
 //    for(i = 4 ;i <= g-1 ; i++)
 //    {
 //        oo[i] = dfs(g-1,1,0,i);
 //    }
 //    for(i = 0 ; i <= g-1; i++)
 //    {
 //        t = (t+oo[i]-o[i])/2;
 //        ans+=t;
 //    }
 //    while(t)
 //    {
 //        t = (t)/2;
 //        ans+=t;
 //    }
 //    printf("%I64d\n",ans);
 //}
 LL cal(int a,int b)
 {
     LL ans = ;
     int t = ,s1,s2;
     if(a>=) s1 = (a+)/;
     if(b>=) s2 = (b+)/;
     t = (s2-s1)/;
     ans+=t;
     for(int i =  ;i <  ;i++)
     {
         if(a>=pp[i]) s1 = cc(a,i);
         else s1=;
         if(b>=pp[i]) s2 = cc(b,i);
         else s2=;
         t = (t+s2-s1)/;
         ans+=t;
     }
     while(t)
     {
         t/=;
         ans+=t;
     }
     return ans;
 }
 int main()
 {
     int a,b;
     pp[] = ;
     for(int i = ; i <  ; i++)
     pp[i] = pp[i-]*;
     memset(dp,-,sizeof(dp));
     int t = ;
     while(scanf("%d%d",&a,&b)!=EOF)
     {
         if(a==b)
         printf("0\n");
         else
         printf("%I64d\n",cal(a-,b));
     }
     return ;
 }