ACM小组内部预定函数

 

1.精度计算——大数阶乘

语法：

int result=factorial(int n);

参数：

n：n 的阶乘

返回值：

阶乘结果的位数

 注意：

本程序直接输出n!的结果，需要返回结果请保留long a[]

需要 math.h

源程序：

 

 int factorial(int n)

{

long a[10000];

int i,j,l,c,m=0,w;

 

a[0]=1;

for(i=1;i<=n;i++)

    {

    c=0;

    for(j=0;j<=m;j++)

        {

        a[j]=a[j]*i+c;

        c=a[j]/10000;

        a[j]=a[j]%10000;

    }

    if(c>0) {m++;a[m]=c;}

}

 

w=m*4+log10(a[m])+1;

printf("/n%ld",a[m]);

for(i=m-1;i>=0;i--) printf("%4.4ld",a[i]);

return w;

}

 

2.精度计算——乘法（大数乘小数）

语法：

mult(char c[],char t[],int m);

参数：

c[]：被乘数，用字符串表示，位数不限

t[]：结果，用字符串表示

m：乘数，限定10以内

返回值：

null

注意：

需要 string.h

源程序：

void mult(char c[],char t[],int m)

{

    int i,l,k,flag,add=0;

    char s[100];

    l=strlen(c);

    for (i=0;i<l;i++)

        s[l-i-1]=c[i]-'0';

 

    for (i=0;i<l;i++)

           {

           k=s[i]*m+add;

           if (k>=10) {s[i]=k%10;add=k/10;flag=1;} else {s[i]=k;flag=0;add=0;}

           }

    if (flag) {l=i+1;s[i]=add;} else l=i;

 

    for (i=0;i<l;i++)

        t[l-1-i]=s[i]+'0';

    t[l]='/0';

}

3.精度计算——乘法（大数乘大数）

语法：

mult(char a[],char b[],char s[]);

参数：

a[]：被乘数，用字符串表示，位数不限

b[]：乘数，用字符串表示，位数不限

t[]：结果，用字符串表示

返回值：

null

注意：

空间复杂度为 o(n^2)

需要 string.h

源程序：

void mult(char a[],char b[],char s[])

{

    int i,j,k=0,alen,blen,sum=0,res[65][65]={0},flag=0;

    char result[65];

    alen=strlen(a);blen=strlen(b);

 

    for (i=0;i<alen;i++)

   for (j=0;j<blen;j++) res[i][j]=(a[i]-'0')*(b[j]-'0');

 

    for (i=alen-1;i>=0;i--)

        {

            for (j=blen-1;j>=0;j--) sum=sum+res[i+blen-j-1][j];

            result[k]=sum%10;

            k=k+1;

            sum=sum/10;

        }

 

    for (i=blen-2;i>=0;i--)

        {

            for (j=0;j<=i;j++) sum=sum+res[i-j][j];

            result[k]=sum%10;

            k=k+1;

            sum=sum/10;

        }

    if (sum!=0) {result[k]=sum;k=k+1;}

 

    for (i=0;i<k;i++) result[i]+='0';

    for (i=k-1;i>=0;i--) s[i]=result[k-1-i];

    s[k]='/0';

 

    while(1)

        {

        if (strlen(s)!=strlen(a)&&s[0]=='0')

            strcpy(s,s+1);

        else

            break;

        }

}

 

4.精度计算——加法

语法：

add(char a[],char b[],char s[]);

参数：

a[]：被乘数，用字符串表示，位数不限

b[]：乘数，用字符串表示，位数不限

t[]：结果，用字符串表示

返回值：

 null

注意：

空间复杂度为 o(n^2)

需要 string.h

源程序：

void add(char a[],char b[],char back[])

{

    int i,j,k,up,x,y,z,l;

    char *c;

    if (strlen(a)>strlen(b)) l=strlen(a)+2; else l=strlen(b)+2;

    c=(char *) malloc(l*sizeof(char));

    i=strlen(a)-1;

    j=strlen(b)-1;

    k=0;up=0;

    while(i>=0||j>=0)

        {

            if(i<0) x='0'; else x=a[i];

            if(j<0) y='0'; else y=b[j];

            z=x-'0'+y-'0';

            if(up) z+=1;

            if(z>9) {up=1;z%=10;} else up=0;

            c[k++]=z+'0';

            i--;j--;

        }

    if(up) c[k++]='1';

    i=0;

    c[k]='/0';

    for(k-=1;k>=0;k--)

        back[i++]=c[k];

    back[i]='/0';

}

5.精度计算——减法

语法：

sub(char s1[],char s2[],char t[]);

参数：

s1[]：被减数，用字符串表示，位数不限

s2[]：减数，用字符串表示，位数不限

t[]：结果，用字符串表示

返回值：

 null

注意：

默认s1>=s2，程序未处理负数情况

需要 string.h

源程序：

void sub(char s1[],char s2[],char t[])

{

    int i,l2,l1,k;

    l2=strlen(s2);l1=strlen(s1);

    t[l1]='/0';l1--;

    for (i=l2-1;i>=0;i--,l1--)

        {

        if (s1[l1]-s2[i]>=0)

            t[l1]=s1[l1]-s2[i]+'0';

        else

            {

            t[l1]=10+s1[l1]-s2[i]+'0';

            s1[l1-1]=s1[l1-1]-1;

            }

        }

    k=l1;

    while(s1[k]<0) {s1[k]+=10;s1[k-1]-=1;k--;}

    while(l1>=0) {t[l1]=s1[l1];l1--;}

loop:

    if (t[0]=='0')

        {

        l1=strlen(s1);

        for (i=0;i<l1-1;i++) t[i]=t[i+1];

        t[l1-1]='/0';

        goto loop;

        }

    if (strlen(t)==0) {t[0]='0';t[1]='/0';}

}