Koch 分形，海岸线，雪花

此算法用于生成Koch分形（海岸线，雪花），速度快速，效果绚丽

//支持的初始直线水平角度为60 的倍数， 交换起点与终点坐标可改变生成方向
void Koch(CDC *pDC, int x1, int y1, int x2, int y2, int n)
{
	if (n > 0)
	{
		double k = 1.0*(y2 - y1) / (x2 - x1);
		int x3 = (x2 - x1) / 3;
		int y3 = (y2 - y1) / 3;

		Koch(pDC, x1, y1, x1 + x3, y1 + y3, n - 1);
		Koch(pDC, x2 - x3, y2 - y3, x2, y2, n - 1);
		if (abs(k) <= 1e-6)
		{
			Koch(pDC, x1 + x3, y1, x1 + x3 + x3 / 2, y1 + sqrt(3.0) / 2 * x3, n - 1);
			Koch(pDC, x1 + x3 + x3 / 2, y1 + sqrt(3.0) / 2 * x3, x2 - x3, y1, n - 1);
		}
		else if (k > 0)
		{
			Koch(pDC, x1, y2 - y3, x2 - x3, y2 - y3, n - 1);
			Koch(pDC, x1 + x3, y1 + y3, x1, y2 - y3, n - 1);
		}
		else
		{
			Koch(pDC, x1 + x3, y1 + y3, x2, y1 + y3, n - 1);
			Koch(pDC, x2, y1 + y3, x2 - x3, y2 - y3, n - 1);
		}
	}
	else
	{
		pDC->MoveTo(x1, y1);
		pDC->LineTo(x2, y2);
	}
}

调用：

int k = 4; //递归的深度
	CPen pen, *oldpen;

	pen.CreatePen(PS_SOLID, 1, RGB(255, 0,0 ));
	oldpen = pDC->SelectObject(&pen);
	Koch(pDC, 800, 150, 440, 150, k);
	pen.DeleteObject();

	pen.CreatePen(PS_SOLID, 1, RGB(0, 255, 0));
	pDC->SelectObject(&pen);
	Koch(pDC, 440, 200, 800, 200, k);
	pen.DeleteObject();

	pen.CreatePen(PS_SOLID, 1, RGB(0, 0, 255));
	pDC->SelectObject(&pen);
	Koch(pDC, 200, 100, 300, 275, k);
	Koch(pDC, 400, 100, 200, 100, k);
	Koch(pDC, 300, 275, 400, 100, k);
	pen.DeleteObject();
	pDC->SelectObject(oldpen);

递归4次的效果图如下：