算法请教,有一泊松随机分布(算法给出),要随机的加到512*512个点上,写了一个算法,要执行两分钟,谁能帮看看,简化一下,谢谢了

[zuohongying]

int poisson(double lambda,double probability)//lambda和probability都是定值20和0.25
{
int k=0;
double r1=rand();//r1为0到32768之间的某值

double r=r1/rand_max;// rand_max=32768为r使在0到1之间
double v=1;
int p;
p=::floor(r+probability);//通过probability随机分布概率控制

while (v>=exp(-lambda))
{
v=v*r;
k++;
}
return (k-lambda)*p;
}
///////////////////////////
int main(void)
{
for(int j=0;j<512;j++)
{
for(int i=0;i<512;i++)
{
g=poisson(20,0.25);
B[i][j]=nAxtualx[i][j]+g;//nAxtualx[i][j]的值已知
if(B[i][j]>255)B[i][j]=255;
}
}
return
}



//////////////////
其中泊松分布算法为

double r=r1/rand_max;// rand_max=32768为r使在0到1之间
double v=1;
while (v>=exp(-lambda))
{
v=v*r;
k++;
}
return k;
}

[mathe]

double v=1;
while (v>=exp(-lambda))
{
v=v*r;
k++;
}
return k;
}
等价于
return k=ceil(-lambda/log(v));

[LLnju]

大家帮我看看，分数由你定: http://www.csdn.net/Expert/topic/503/503163.shtm

[msdos]

你的程序写得错误太多，我编译出现许多错误。
测了一下你的程序，在我的机器上只化了不到一秒。
记住：想让别人帮忙就应该把自己得问题和目标描述清楚，另外
发问之前自己先仔细想，不要太依赖别人

[zuohongying]

thank you,let me try agin