标准程序
#include
#include
#include
#include
#include
#include
#include
using namespace std;
double v[6][16][16], p[16][16];
int main() {
for (int i = 0; i < 16; i++)
for (int j = 0; j < 16; j++)
scanf("%lf", &p[i][j]);
memset(v, 0, sizeof(v));
for (int i = 0; i < 16; i++)
v[0][i][0] = 1.0;
for (int i = 1; i <= 4; i++)
for (int j = 0; j < 16 >> i; j++) {
for (int x = 0; x < 1 << (i - 1); x++)
for (int y = 0; y < 1 << (i - 1); y++)
v[i][j][x] += v[i - 1][j << 1][x] * v[i - 1][(j << 1) + 1][y] * p[(j << i) + x][(j << i) + (1 << (i - 1)) + y],
v[i][j][(1 << (i - 1)) + y] += v[i - 1][j << 1][x] * v[i - 1][(j << 1) + 1][y] * (1.0 - p[(j << i) + x][(j << i) + (1 << (i - 1)) + y]);
}
for (int i = 0; i < 16; i++)
printf("%.10f%c", v[4][0][i], i == 15 ? '\n' : ' ');
}