処理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);
}
}
簡単に言うと、プログラムは2つの「シーカー」を送信します。1つは反時計回りに、もう1つは時計回りに移動します。これらのシーカーはそれぞれ、最も近い極を見つけて、鎖が十分に長い場合はそれに弧を描き、それ以外の場合は円を描きます。円弧を描くと、別のシーカーをその極に送信し、プロセスが続行されます。f()各シーカーのプロセスが含まれています。私がこれをもっとゴルフするとすぐに、より詳細な説明が来るでしょう。
![x（[1j、1 + 1j]）](https://i.stack.imgur.com/rTW3q.png)
{0,-.5}何ですか？