r=max(r,1);

vector<int> depth(h*w,-2100000000);

vector<PointInfo> pattern;

for (int dx=-r;dx<=r;dx++)

for (int dy=-r;dy<=r;dy++)

if (dx*dx+dy*dy<r*r){

double dz=sqrt(double(r*r-dx*dx-dy*dy));

PointInfo pinfo;

pinfo.x=dx;

pinfo.y=dy;

pinfo.z=dz;

pinfo.r=dz/r;

pinfo.g=dz/r;

pinfo.b=dz/r;

pattern.push_back(pinfo);

}

double zmin=0,zmax=0;

for (int i=0;i<n;i++){

if (i==0){

zmin=xyzs[i*3+2]-r;

zmax=xyzs[i*3+2]+r;

}else{

zmin=min(zmin,double(xyzs[i*3+2]-r));

zmax=max(zmax,double(xyzs[i*3+2]+r));

}

}

for (int i=0;i<n;i++){

int x=xyzs[i*3+0],y=xyzs[i*3+1],z=xyzs[i*3+2];

for (int j=0;j<int(pattern.size());j++){

int x2=x+pattern[j].x;

int y2=y+pattern[j].y;

int z2=z+pattern[j].z;

if (!(x2<0 || x2>=h || y2<0 || y2>=w) && depth[x2*w+y2]<z2){

depth[x2*w+y2]=z2;

double intensity=min(1.0,(z2-zmin)/(zmax-zmin)*0.7+0.3);

show[(x2*w+y2)*3+0]=pattern[j].b*c2[i]*intensity;

show[(x2*w+y2)*3+1]=pattern[j].g*c0[i]*intensity;

show[(x2*w+y2)*3+2]=pattern[j].r*c1[i]*intensity;

}

}

}