高精度（高精度加、减、乘、除、阶乘）

高精度加法

 代码实现：

# include<stdio.h>
# include<string.h>
int main()
{
    char num1[110],num2[110];
    int a[110]={0},i;
    scanf("%s%s",num1,num2);
    int l1=strlen(num1);
    int l2=strlen(num2);
    for(i=0;i<l1||i<l2;i++)
    {
        if(i<l1)
        a[i]+=num1[l1-i-1]-'0';
        if(i<l2)
        a[i]+=num2[l2-i-1]-'0';
        if(a[i]>=10)
        {
            a[i+1]+=a[i]/10;
            a[i]%=10;
        }
    }
    int len=(l1>l2)?l1:l2;
    if(a[len]>0) len++;
    for(i=len-1;i>=0;i--)
    printf("%d",a[i]);
    return 0; 
}

高精度减法 

代码实现：
 

# include<stdio.h>
# include<string.h>
char num1[999],num2[999];
int len,a[999],b[999],i;
void jianfa()
{
    for(i=0;i<len;i++)
    {
        if(a[i]<b[i])
        {
            a[i+1]--;
            a[i]+=10;
        }
        a[i]-=b[i];
    }
    while(!a[len-1])
    len--;
    for(i=len-1;i>=0;i--)
    printf("%d",a[i]);
}
int main()
{
    scanf("%s%s",num1,num2);
    int l1=strlen(num1);
    int l2=strlen(num2);
    if(l1>l2||(l1==l2&&strcmp(num1,num2)>0))
    {
        len=l1;
        for(i=0;i<l1;i++)
        a[i]=num1[l1-i-1]-'0';
        for(i=0;i<l2;i++)
        b[i]=num2[l2-i-1]-'0';
        jianfa();
    }
    else
    {
        printf("-");
        len=l2;
        for(i=0;i<l1;i++)
        b[i]=num1[l1-i-1]-'0';
        for(i=0;i<l2;i++)
        a[i]=num2[l2-i-1]-'0';
        jianfa();
    }
    return 0;
 } 

高精度乘法

 

代码实现：

# include<stdio.h>
# include<string.h>
int main()
{
    char num1[1010],num2[1010];
    scanf("%s%s",num1,num2);
    int l1,l2;
    l1=strlen(num1);
    l2=strlen(num2);
    int a[l1+10],b[l2+10],i,j;
    for(i=0;i<l1;i++)
    a[i]=num1[l1-i-1]-'0';
    for(i=0;i<l2;i++)
    b[i]=num2[l2-i-1]-'0';
    int c[3000]={0};
    for(i=0;i<l1;i++)
    for(j=0;j<l2;j++)
    c[i+j]+=a[i]*b[j];
    for(i=0;i<l1+l2;i++)
    {
        if(c[i]>=10)
        {
        c[i+1]+=c[i]/10;
        c[i]%=10;
        }
    }
    for(j=l1+l2;j>0;j--)
    {
        if(c[j]!=0)
            break;
    }
    for(i=j;i>=0;i--)
    printf("%d",c[i]);
    printf("\n");
    return 0;
}

 高精度除法

 

代码实现：

#include<iostream>
#include<cstring>
using namespace std;
int a[100],b[100],c[100];
int compare(int a[],int b[])//比较a、b，若a>b为1；若a<b为-1；若a=b为0
{
    int i;
    if(a[0]>b[0])
        return 1;
    if(a[0]<b[0])
        return -1;
    for(i=a[0];i>0;i--)//从高位到低位比较
    {
        if(a[i]>b[i])
            return 1;
        if(a[i]<b[i])
            return -1;
    }
    return 0;
}
void subduction(int a[],int b[])//计算a=a-b
{
    int flag;
    int i;
 
    flag=compare(a,b);
    if(flag==0)//相等
    {
        a[0]=0;
        return;
    }
    if(flag==1)//大于
    {
        for(i=1;i<=a[0];i++)
        {
            if(a[i]<b[i])//若不够向上借位
            {
                a[i+1]--;
                a[i]+=10;
            }
            a[i]-=b[i];
        }
        while(a[0]>0&&a[a[0]]==0)//删除前导0
            a[0]--;
        return;
    }
}
int main()
{
    char str1[100],str2[100];
    int i,j;
    cin>>str1>>str2;
    a[0]=strlen(str1);//a[0]存储串1的位数
    b[0]=strlen(str2);//b[0]存储串2的位数
    for(i=1;i<=a[0];i++)
        a[i]=str1[a[0]-i]-'0';
    for(i=1;i<=b[0];i++)
        b[i]=str2[b[0]-i]-'0';
    int temp[100];
    c[0]=a[0]-b[0]+1;
    for(i=c[0];i>0;i--)
    {
        memset(temp,0,sizeof(temp));
 
        for(j=1;j<=b[0];j++)//从i开始的地方，复制数组b到数组temp
            temp[j+i-1]=b[j];
        temp[0]=b[0]+i-1;
 
        while(compare(a,temp)>=0)//用减法模拟
        {
            c[i]++;
            subduction(a,temp);
        }
    }
    while(c[0]>0&&c[c[0]]==0)//删除前导0
        c[0]--;
    if(c[0]==0)//输出结果
        cout<<0<<endl;
    else
    {
        for(i=c[0];i>0;i--)
            cout<<c[i];
        cout<<endl;
    }
    if(a[0]==0)//输出余数
        cout<<0<<endl;
    else
    {
        for(i=a[0];i>0;i--)
            cout<<a[i];
        cout<<endl;
    }
    return 0;
}

阶乘之和 

 代码实现：
 

#include<stdio.h>
int main()
{
    int i,A[1005]={0},B[1005]={0},n,j;
    scanf("%d", &n);
    A[0]=B[0]=1;
    for (i=2;i<=n;i++){
        for (j=0;j<100;j++)
            B[j]*=i;
        for (j=0;j<100;j++)
            if (B[j]>9){
                B[j+1] += B[j]/10;
                B[j]%=10;
            }
        for (j=0;j<100;j++){
            A[j]+=B[j];
            if (A[j]>9) {
                A[j+1] += A[j]/10;
                A[j]%=10;
            }
        }
    }
    for (i=100;i>=0&&A[i]==0;i--);
    for (j=i;j>=0;j--) printf("%d", A[j]);
    return 0;
}