N阶贝塞尔曲线

long fac(int m)
{
	int i;
	long temp=1;
	if(m==0)	return 1;
	else
	{
		for(i=2;i<=m;i++)
			temp=temp*i;
	}
	return temp;
}

double powi(double v, int k)
{
	double temp=1.0;
	if(k==0||v==0)	return 1;
	else	
	{
		for(int i=1;i<=k;i++)
			temp=temp*v;
	}

	return temp;
}
void BezierN(CDC *pDC,CPoint *m_BezierN, int m_numBezierN)
{
	CPen Newpen1(PS_SOLID,1,RGB(0,0,0));
	CPen *pOldPen1=pDC->SelectObject(&Newpen1);

	int x,y,i,j,k=40,a,b;
	double	t,t1,u,v;
	double temp,temp1,temp2,bi;
	t=1.0/k;
	a=m_BezierN[0].x;
	b=m_BezierN[0].y;
	//pDC->MoveTo(a,b);
	for(j=1;j<k;j++ )
	{
		pDC->MoveTo(a,b);
		t1=j*t;
		u=t1;v=1-u;
		x=0;y=0;
		for(i=0;i<=m_numBezierN;i++)
		{
			temp=double(fac(m_numBezierN)/fac(i)/fac(m_numBezierN-i));
			temp1=powi(u,i);
			temp2=powi(v,m_numBezierN-i);
			bi=temp*temp1*temp2;
			x=x+bi*m_BezierN[i].x;
			y=y+bi*m_BezierN[i].y;
		}
		pDC->LineTo(x,y);
		a=x;
		b=y;
	}
	pDC->SelectObject(pOldPen1);
}