Codeforces 441 E. Valera and Number

/*
 Codeforces 441 E. Valera and Number
 题意:给定三个整数,x,k,p.进行k轮操作,每轮有p%的可能性将x乘上2,有(100-p)%的可能性把x加上1,问最后x的二进制表示中末尾有多少个零的期望.
 解法:dp解决,dp[k][i][j]表示进行k轮后,后八位是j并且后八位之前连续的1或者0的长度.(之所以定义为八位是因为考虑到加法操作再怎么影响,也只有后八位,可能影响前面的,但结算过一次之后,就不可能再影响了)
 之后dp的转移就比较好想了,注意转移的方法以及位运算的写法.
 */
#include <cstdio>
#include <cstring>

double dp[201][501][256];

int cal(int i,int j)
{
    int ans = 0;
    while(i % 2 == 0)
    {
        ans ++;
        i /= 2;
		if(ans == 8)
			break;
    }
    if(ans < 8) return ans;
    if(j > 250) return ans + j - 250;
    return ans;
}

int main()
{
    int x, k;
    double p;
    scanf("%d %d %lf", &x, &k, &p);
    p = p / 100.;
    
    //预处理dp
    int loc = x & 255;
    int base;
    if( (x & (1<<8)) > 0)
    {
        base = 1;
        for(int i = 9;;i++)
        {
            if((x&(1<<i)) == 0) break;
            base ++;
        }
        dp[0][base][loc] = 1.;
    }
    else
    {
        if(x >= (1<<8))
        {
            base = 250;
            for(int i = 8;;i++)
            {
                if((x&(1<<i)) > 0) break;
                base ++;
            }
            dp[0][base][loc] = 1.;
        }
        else
        {
            dp[0][0][loc] = 1.;
        }
    }
    //printf("%lf\n",dp[0][0][loc]);
    //dp过程.
    for(int step = 0; step < k; step ++)
    {
        for(int i=0; i<= 500; i++)
        {
            for(int j=0; j<256; j++)
            {
                if(dp[step][i][j] == 0) continue;
                if(i == 0)
                {
                    if((j&128) > 0)
                    {
                        dp[step+1][1][(j<<1)&255] += dp[step][0][j] * p;
                        
                    }
                    else
                    {
                        dp[step+1][0][j<<1] += dp[step][0][j] * p;
                    }
                    if(j < 255)
                    {
                        dp[step+1][0][j+1] += dp[step][0][j] * (1 - p);
                    }
                    else
                    {
                        dp[step+1][1][0] += dp[step][0][j] * (1 - p);
                    }
                    //printf("%lf\n",dp[1][0][2]);
                    continue;
                }
                if(i <= 250)
                {
                    //p的概率乘法
                    if((j & 128) > 0)
                    {
                        dp[step+1][i+1][(j<<1)&255] += dp[step][i][j] * p;
                    }
                    else
                    {
                        dp[step+1][251][(j<<1)&255] += dp[step][i][j] * p;
                    }
                    //1-p的概率加法
                    if(j < 255)
                    {
                        dp[step+1][i][j+1] += dp[step][i][j] * (1 - p);
                    }
                    else
                    {
                        dp[step+1][250 + i][0] += dp[step][i][j] * (1 - p);
                    }
                }
                else
                {
                    //p的概率乘法
                    if((j & 128) > 0)
                    {
                        dp[step+1][1][(j<<1)&255] += dp[step][i][j] * p;
                    }
                    else
                    {
                        dp[step+1][i+1][(j<<1)&255] += dp[step][i][j] * p;
                    }
                    //加法
                    if(j < 255)
                    {
                        dp[step+1][i][j+1] += dp[step][i][j] * (1 - p);
                    }
                    else
                    {
                        dp[step+1][1][0] += dp[step][i][j] * (1 - p);
                    }
                }
            }
        }
        
    }
    //printf("%lf %lf\n",dp[1][0][10],dp[1][0][6]);
    //统计答案
    double ans = 0.;
    for(int i = 0;i <= 500; i++)
    {
        for(int j = 0;j < 256; j++ )
        {
            if(dp[k][i][j] == 0)    continue;
            ans += dp[k][i][j] * cal(j, i);
        }
    }
    printf("%.13lf\n",ans);
    return 0;
}