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;
}