BNU 0821 C. Probability Computation

2011 ACM-ICPC Asia Hsinchu Regional Contest

是一个递推式。
s [ i ] [ j ] 表示的是有 i 位（二进制），余数为 j 的概率。
a [ i ] 表示在第 i 位添加 1 ，会为这个数添加的余数。（ % q ）
s [ i ] [ j ] = ( 1 - p [ i ] ) * s [ i - 1 ] [ j ] + p [ i ] * s [ i - 1 ] [ ( j + q - a [ i ] ) % q ] ;
式子里的第一部分好理解，就是有i位的时候，如果最高位为0，则选择乘上（ 1 - p [ i ] ），此时，除去最高位，选择 s [ i - 1 ] [ j ] 即可，因为最高位为 0 ， i 位和 i - 1 位的余数是一样的；
第二部分，当最高位选择为 1 的时候，要乘上 p [ i ] ，此时第 i 位添加 1 ，即为 i - 1 位的时候的余数添加了 a [ i ] ，用（ j + q - a [ i ] ） % q ，即是当最高位为 1 要保证余数为 j 时选择的长度为 i - 1 位的余数。


【说的好乱……

P.S.注意存 p [ i ] 的时候是倒着存的，因为第一个其实是最高位，最后进行计算的。

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
int a[210];
double s[210][110];
double p[210];
int main(){
    int n,q,r;
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%d",&n);
        for(int i=n;i>=1;i--)
            scanf("%lf",&p[i]);
        scanf("%d%d",&q,&r);
        a[1]=1;
        for(int i=2;i<=n;i++){
            a[i]=(2*a[i-1])%q;
           // printf("%d\n",a[i]);
        }
        for(int i=1;i<=n;i++){
            p[i]=1.0*p[i]/100;
            //printf("%.lf\n",p[i]);
        }
        memset(s,0,sizeof(s));
        s[0][0]=1;
        for(int i=1;i<=n;i++){
            for(int j=0;j<q;j++){
                s[i][j]=(1-p[i])*s[i-1][j]+p[i]*s[i-1][(j+q-a[i])%q];
                //printf("%.lf\n",s[i][j]);
            }
        }
        printf("%.5lf\n",s[n][r]);
    }
    return 0;
}