UVAlive 5863 Probability Computation 문제 해결

제목.
Tri_integral Summer Training 4
제목:
n자리의 이진수 x는 한 자리가 1 또는 0인 확률을 알고 있으며 xmod q==r의 확률을 구한다.
문제 풀이:
dp[i][j]로 i위와 i위 이전의 디지털 q가 j일 확률을 표시하면 이동할 때 0과 2^i로 이동하고 마지막으로 dp[n][r]가 답이다.
 

//Time:9ms
//Memoty:0KB
//Length:1019B
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <set>
#define MAXN 210
#define INF 1000000007
#define MP(x,y) make_pair(x,y)
#define FI first
#define SE second
using namespace std;
double dp[2][MAXN];
int p[MAXN];
int main()
{
    //freopen("/home/moor/Code/input","r",stdin);
    int r,q,s,n,ncase,step;
    scanf("%d",&ncase);
    while(ncase--)
    {
        scanf("%d",&n);
        for(int i=0;i<n;++i)    scanf("%d",&p[i]);
        scanf("%d%d",&q,&r);
        for(int i=0;i<2;++i)
            for(int j=0;j<q;++j)
                dp[i][j]=0;
        s=0;
        dp[s][0]=1;
        step=1;
        for(int i=n-1;i>=0;--i)
        {
            for(int j=0;j<q;++j)
                dp[!s][j]=0;
            for(int j=0;j<q;++j)
                dp[!s][j]+=dp[s][j]*(100.0-p[i])/100.0,
                dp[!s][(j+step)%q]+=dp[s][j]*p[i]/100.0;
            s=!s;
            step=(step*2)%q;
        }
        printf("%.5f
",dp[s][r]);
    }
    return 0;
}