UVA 11176 Winning Streak

#include <iostream>
#include <stdio.h>
#include <cstring>
#define N 501

using namespace std;
int n;
double p;
double pa

;
double pow
;
void Init()
{
memset(pa,0,sizeof(pa));
pow[0]=1;
for(int i=1;i<N;i++)
{
pow[i] = pow[i-1]*p;
}
}

void Do()
{
for(int i=0;i<=n;i++)
{
for(int j=0; j<=n; j++)
{
if(i-j<=0)
pa[i][j] = 1;
else if(i-j-1==0)
pa[i][j]=pa[i-1][j]-pow[j+1];
else if(i-j-2>=0)
pa[i][j]=pa[i-1][j]-pa[i-j-2][j]*(1-p)*pow[j+1];
}
}

double ans = 0;
for(int i=1;i<=n;i++)
{
ans += (pa
[i]-pa
[i-1])*i;
}
printf("%.6lf\n",ans);
}

int main()
{
while(scanf("%d %lf",&n,&p)!=EOF && n)
{
Init();

Do();
}
return 0;
}