double@n,@p,a[1d4];wt(sum[i,2,n+1]((1-p)**(Divisor(i,a)-2)));