C语言实现二维泊松碟随机均匀采样算法

#include<stdlib.h>
#include<stdio.h>
#include<time.h>

int main(int argc, char* argv[])
{      int n=50,d=4;
        int i,r,s,a[n],b[n];
        int  c[2][n];
        srand((unsigned)time(NULL));
        i=0;
        while (i<n)
        {
                a[i]=int(n*(rand()/(RAND_MAX+0.5)));
                b[i]=int(n*(rand()/(RAND_MAX+0.5)));
                c[0][i]=a[i];
                c[1][i]=b[i];
                for(s=0;s<i;s++)
                {
                        //      a[s]=int(20*(rand()/(RAND_MAX+0.5)));
                        //     b[s]=int(20*(rand()/(RAND_MAX+0.5)));

                        if ((a[s]-a[i])*(a[s]-a[i])+(b[s]-b[i])*(b[s]-b[i])<d)
                        {
                                i--;
                                break;
                        }

                }
                i++;