CodeForces 441E(Codeforces Round #252 (Div. 2))

思路：dp[i][now][mark][len] i 表示当前第i 次now存的是后8位，mark为第9位为0还是1 len第九位往高位还有几位和第9位相等。 只存后8位的原因：操作只有200次每次都为加法的话后8位可以表示，如果为乘法第八位已知再加上第九位 和往前的长度已知，所以可以表示所有状态。

所存在问题就是 10 1111 1111 此时加上1之后 会变成 11 0000 0000 但这样并处影响结果 如果之后操作都为加法，只有200次，他不可能影响到前面的1， 乘法相当于左移也不会影响。所以前面的1 不可能作为答案。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include <iostream>
using namespace std;
double d[210][260][2][260];
int a[15];
int main() {
int x, k, p, h = 0;
memset(d, 0, sizeof(d));
memset(a, 0, sizeof(a));
scanf("%d%d%d", &x, &k, &p);
while (x) {
a[h++] = x %2;
x /= 2;
}
int now = 0;
for (int i = 0; i < 8; ++i)
if (a[i])
now += (1 << i);

if (h < 9) {
d[0][now][0][0] = 1;
} else {
int cnt = 1;
for (int i = 9; i < h; ++i) {
if (a[i] != a[i - 1])
break;
cnt++;
}
d[0][now][a[8]][cnt] = 1;
}
for (int i = 0; i < k; ++i) {
for (int j = 0; j <= 255; ++j) {
for (int x = 0; x <= 250; ++x) {
if (j != 255) {
d[i + 1][j + 1][0][x] += d[i][j][0][x] * (100 - p) / 100.0;
d[i + 1][j + 1][1][x] += d[i][j][1][x] * (100 - p) / 100.0;
} else {
d[i + 1][0][0][x] += d[i][j][1][x] * (100 - p) / 100.0;
d[i + 1][0][1][1] += d[i][j][0][x] * (100 - p) / 100.0;
}

if (j & (1 << 7)) {
d[i + 1][(j << 1) % 256][1][1] += d[i][j][0][x] * p / 100.0;
d[i + 1][(j << 1) % 256][1][x + 1] += d[i][j][1][x] * p
/ 100.0;
} else {
d[i + 1][(j << 1) % 256][0][1] += d[i][j][1][x] * p / 100.0;
d[i + 1][(j << 1) % 256][0][x + 1] += d[i][j][0][x] * p
/ 100.0;
}
}
}
}
double sum=0;
for(int i=1;i<255;++i)
{
for(int j=0;j<2;++j)
for(int x=0;x<=250;++x)
{
int now=i;
int cnt=0;
while(now%2==0)
{
cnt++;
now/=2;
}
sum+=d[k][i][j][x]*cnt;
}
}
for(int x=0;x<=250;++x)
sum+=d[k][0][1][x]*8;
for(int x=0;x<=250;++x)
sum+=d[k][0][0][x]*(x+8);
printf("%.10lf\n",sum);

return 0;
}

View Code