3 815 처리 833 835 876 879 바이트를
불필요한 괄호를 제거하여 @ZacharyT 덕분에 2 바이트를 절약했습니다.
void settings(){size(600,600);}int i,w,x,n;float l,d,t,a,f,g,m,R,U;float[][]N,T;float[]S,p;void s(float[][]t){N=new float[t.length+1][2];N[0][0]=N[0][1]=i=0;for(float[]q:t)N[++i]=q;translate(w=300,w);noFill();pushMatrix();f(N,0,-w,w,1,0);popMatrix();f(N,0,-w,w,0,0);}float p(float a,float b){for(a+=PI*4;a>b;)a-=PI*2;return a;}void f(float[][]P,float x,float y,float L,int c,int I){l=2*PI;d=i=0;S=null;for(;i<P.length;i++){float[]p=P[i];g=atan2(y,x);m=atan2(p[1],p[0]);if(p(f=(c*2-1)*(g-m),0)<l&(t=dist(0,0,p[0],p[1]))<=L&I!=i){l=p(f,0);S=new float[]{g,m};d=t;n=i;}}if(S==null)ellipse(0,0,2*(L-d),2*(L-d));else{arc(0,0,L*2,L*2,p(S[c],S[1-c]),S[1-c]);R=cos(a=S[1]);U=sin(a);translate(d*R,d*U);T=new float[P.length][2];for(int i=0;i<T.length;T[i][1]=P[i][1]-d*U,i++)T[i][0]=P[i][0]-d*R;f(T,(L-d)*R,(L-d)*U,L-d,c,n);}}
이 프로그램을 다음과 같이 실행하십시오 :
void setup() {
s(new float[][]{{0,100},{100,100},{-200,100},{-100,-50}});
}
(함수 s는를받습니다 float[][]). 이것은 기본적으로 테스트 케이스 # 3이지만 창에 맞추기 위해 100을 곱합니다.
몇 가지 참고할 사항 :
- 이 프로그램은 극을 그리지 않습니다
- Processing의 좌표계에서 양의 y 축이 내려 가기 때문에 이미지가 거꾸로 뒤집힌 것처럼 보입니다.
- Processing은 float를 사용하기 때문에 계산이 정확하지 않으므로 이미지에서이를 확인할 수 있습니다. 이러한 부동 소수점 오류가 문제인지 OP에 문의했습니다.
- 창의 크기는 600 x 600 픽셀입니다.
- 스택
pushMatrix()과 popMatrix()연산이 32 개의 행렬 만 보유 할 수 있기 때문에 매우 작은 입력 좌표는 프로그램을 중단시킵니다.
- 개는 (0, -300)에서 시작하고 체인은 300 픽셀 길이에서 시작합니다.
- 아래 이미지는 편의를 위해 축소되었습니다.
위 테스트 케이스의 샘플 출력.
사전 출력을 보려면 translate(w,w);in 함수 바로 다음에이 줄을 추가하십시오 s.
background(-1);scale(1,-1);fill(255,0,0);ellipse(0,0,25,25);fill(0);for(float[]q:N)ellipse(q[0],q[1],25,25);
그리고 이것은 우리에게이 결과를줍니다 :
언 골프 f()와 설명
(디버그 코드도 포함)
void f(float[][]points, float x, float y, float len, int c, int pindex) {
print(asd+++")");
float closest = 2*PI;
float d=0,t;
float[]stuff = null;
int index = 0;
for(int i=0;i<points.length;i++) {
if(pindex != i) {
float[]p = points[i];
float originAngle = atan2(y, x);
float tempAngle = atan2(p[1], p[0]);
//println(x,y,p[0],p[1]);
float diff = c<1?tempAngle-originAngle:originAngle-tempAngle;
println("@\t"+i+"; x=\t"+x+"; y=\t"+y+"; tx=\t"+p[0]+"; ty=\t",p[1], diff, originAngle, tempAngle);
if(p(diff) < closest && (t=dist(0,0,p[0],p[1])) < len) {
println("+1");
closest = p(diff);
stuff = new float[]{originAngle, tempAngle};
d=t;
index = i;
}
}
}
if(stuff == null) {
ellipse(0,0,2*(len-d),2*(len-d));
println("mayday");
} else {
println("d angles",d,p(stuff[c],stuff[1-c],c), stuff[1-c]);
//println(points[0]);
arc(0, 0, len*2, len*2, p(stuff[c],stuff[1-c],c), stuff[1-c]);
float angle = stuff[1];
translate(d*cos(angle), d*sin(angle));
println("Translated", d*cos(angle), d*sin(angle));
println("angle",angle);
float[][]temp=new float[points.length][2];
for(int i=0;i<temp.length;i++){
temp[i][0]=points[i][0]-d*cos(angle);
temp[i][1]=points[i][1]-d*sin(angle);
println(temp[i]);
}
println(d*sin(angle));
pushMatrix();
println();
f(temp, (len-d)*cos(angle), (len-d)*sin(angle), (len-d), c, index);
popMatrix();
//f(temp, (len-d)*cos(angle), (len-d)*sin(angle), (len-d), 0, index);
}
}
간단히 말해, 프로그램은 두 개의 "탐색자"를 전송하는데, 하나는 시계 반대 방향으로, 다른 하나는 시계 방향으로 이동합니다. 이 구도자 각각은 가장 가까운 극점을 찾아 체인이 충분히 길면 원호를 그립니다. 일단 호를 그리면 다른 추적자를 해당 극으로 보내고 프로세스는 계속됩니다. f()각 시커의 프로세스를 포함합니다. 골프를하자마자 더 자세한 설명을하겠습니다.
![x ([1j, 1 + 1j])](https://i.stack.imgur.com/rTW3q.png)
{0,-.5}?