一些公式代码（慢慢添加）

求n!的位数

/*原理：log10(n!)=log10(1*2*...*n)=log10(1)+log10(2)+...+log10(n)*/
for(i=1; i<=n; i++)
    sum+=log10((double)i);
sum++; //sum是double型的

欧几里得算法，求最大公因子

EUCLID(int a,int b)
{
    if(b==0) return a;
    else EUCLID(b,a%b);
}

排列组合公式

int pacf(int m,int n) //用的是int型函数
{
    int i,k; //这里的k是最终值，如果数值比较大可以改用64的
    if(n==0||n==m) k=1;
    else
    {
        k=m;
        for(i=m-1; i>m-n; i--)k*=i;
        for(i=2; i<=n; i++)
            k/=i;
    }
    return k;
}

求1..n-1中与n互质的数的个数

int eular(int n)
{
    int ret=1,i;
    for (i=2; i*i<=n; i++)
        if (n%i==0)
        {
            n/=i,ret*=i-1;
            while (n%i==0) n/=i,ret*=i;
        }
    if (n>1)
        ret*=n-1;
    return ret;
}

求π

#include<cmath> //注意头文件
#define pai acos(-1.0)

点到直线的距离

double lenth(point x,point y,point z)   //求点z到直线xy的距离
{
    if(x.x==y.x) return fabs(z.x-x.x);
    double a=(x.y-y.y)/(x.x-y.x);
    double b=x.y-a*x.x;
    double ans=fabs(a*z.x-z.y+b)/sqrt(a*a+1);
    return ans;
}

两点的距离

double juli(point a,point b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}