美文网首页
高精度计算

高精度计算

作者: nino天 | 来源:发表于2014-08-22 00:27 被阅读22次

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

    #include <stdio.h>
    #include <string.h>
    #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<MAX_LEN;i++)//加法运算
        {
            an1[i] += an2[i];
            if(an1[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 <stdio.h>
    #include <string.h>
    #include <math.h>
    #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<strlen(Line2);i++)//用an1的一位去和an2的各位相乘
            for(j=0;j<strlen(Line1);j++)
                Result[i+j] += an2[i]*an1[j];
    
        for(i=0;i<MAX_LEN*2;i++)//统一的进位
            if(Result[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<LEN;i++)
        {
            nCarry=0;
            for(j=0;j<LEN-i;j++)
            {
                nTmp=c[i+j]+a[j]*b[i]+nCarry;
                c[i+j]=nTmp%10000;
                nCarry=nTmp/10000;
            }
        }
    

    输出格式为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 <stdio.h>
    #include <memory.h>
    #define LEN 125
    #include <math.h>
    
    
    void Multiply(int *a,int *b)
    {
        int i,j;
        int nCarry;
        int nTmp;
        int c[LEN]={0};
        for(i=0;i<LEN;i++)
        {
            nCarry=0;
            for(j=0;j<LEN-i;j++)
            {
                nTmp=c[i+j]+a[j]*b[i]+nCarry;
                c[i+j]=nTmp%10000;
                nCarry=nTmp/10000;
            }
        }
        memcpy(a,c,LEN*sizeof(int));
    
    }
    
    
    void main()
    {
        int i,p;
        int Pow[LEN]={0},Result[LEN]={0};
        scanf("%d",&p);
    
        printf("%d\n",(int)(p*log10(2))+1);//位数
    
        Pow[0]=2;
        Result[0]=1;
    
        while(p>0)
        {
            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 <stdio.h>
    #include <math.h>
    #include <string.h>
    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;
    }
    

    相关文章

      网友评论

          本文标题:高精度计算

          本文链接:https://www.haomeiwen.com/subject/dfzetttx.html