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);
}