高精度计算

1.大整数加法:开两队数组,输入为字符串,然后要转换成int数组。最先输入的是char[0]的位置,然而这个应该是int的高位,所以要头尾转换。

#include 
#include 
#define MAX_LEN 250

int an1[MAX_LEN]={0};
int an2[MAX_LEN]={0};
char Line1[MAX_LEN];
char Line2[MAX_LEN];

void main()
{
    scanf("%s",Line1);
    scanf("%s",Line2);
    int i,j;

    j=0;//字符转存成数字,个位放在an1[0]处
    for(i=strlen(Line1)-1;i>=0;i--)
        an1[j++]=Line1[i]-'0';
    
    j=0;
    for(i=strlen(Line2)-1;i>=0;i--)
        an2[j++]=Line2[i]-'0';

    for(i=0;i= 10)
        {
            an1[i] -= 10;
            an1[i+1]++;
        }
    }

    bool flag = false;//跳过多余的0
    for(i=MAX_LEN-1;i>=0;i--)
    {
        if(flag) printf("%d",an1[i]);
        else if(an1[i])
        {
            printf("%d",an1[i]);
            flag = true;
        }
    }
    printf("\n");
}

2.大数乘法:用第三个结果数组去存储,两个乘数每一位相乘的结果,最后统一进行进位操作。

#include 
#include 
#include 
#define MAX_LEN 250

int an1[MAX_LEN]={0};
int an2[MAX_LEN]={0};
char Line1[MAX_LEN];
char Line2[MAX_LEN];
int Result[MAX_LEN*2]={0};

void main()
{
    scanf("%s",Line1);
    scanf("%s",Line2);
    int i,j;

    j=0;//字符转存成数字,个位放在an1[0]处
    for(i=strlen(Line1)-1;i>=0;i--)
        an1[j++]=Line1[i]-'0';
    
    j=0;
    for(i=strlen(Line2)-1;i>=0;i--)
        an2[j++]=Line2[i]-'0';

    for(i=0;i=10)
        {
            Result[i+1] += Result[i]/10;
            Result[i] %= 10;
        }

    bool flag = false;//跳过多余的0
    for(i=MAX_LEN*2-1;i>=0;i--)
    {
        if(flag) printf("%d",Result[i]);
        else if(Result[i])
        {
            printf("%d",Result[i]);
            flag = true;
        }
    }
    printf("\n");

}

2.麦森数:重点是(1)计算2的p次方(2)只留500位(3)为了不超时数组的一个元素用于存放10000进制数的1位

p用二进制表示,可把2p表示成2a0 * 2^(2a1) *... *2(2n-1)(a_n-1)

while(p>0)
    {
        if(p&1)
            Multiply(Result,Pow);
        p>>=1;
        Multiply(Pow,Pow);
    }

大数乘法要把int数组的每一位都存成10000进制的数

for(i=0;i

输出格式为04d,而且遇到情况要拆行输出

for(i=LEN-1;i>=0;i--)
    {
        if(i%25 == 12)
            printf("%02d\n%02d",Result[i]/100,Result[i]%100);
        else
        {
            printf("%04d",Result[i]);
            if(i%25==0)
                printf("\n");
        }
    }
#include 
#include 
#define LEN 125
#include 


void Multiply(int *a,int *b)
{
    int i,j;
    int nCarry;
    int nTmp;
    int c[LEN]={0};
    for(i=0;i0)
    {
        if(p&1)
            Multiply(Result,Pow);
        p>>=1;
        Multiply(Pow,Pow);
    }
    Result[0]--;

    for(i=LEN-1;i>=0;i--)
    {
        if(i%25 == 12)
            printf("%02d\n%02d",Result[i]/100,Result[i]%100);
        else
        {
            printf("%04d",Result[i]);
            if(i%25==0)
                printf("\n");
    
        }
    }
}

3.阶乘的运算

#include 
#include 
#include 
int f[3000]={0};

int main()
{
    int n;
    scanf("%d",&n);
    for(i=2;i<=n;i++)//乘以i
    {
        int c=0;
        for(j=0;j<3000;j++)
        {
            s=f[j]*i+c;
            f[j]=s%10;
            c=s/10;
        }
    }
    for(i=2999;i>=0;i--) if(f[j]) break;
    for(i=j;i>=0;i--) printf("%d",f[i]);
    printf("\n");
    return 0;
}

你可能感兴趣的:(高精度计算)