高精度计算

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;
}
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 157,198评论 4 359
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 66,663评论 1 290
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 106,985评论 0 237
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 43,673评论 0 202
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 51,994评论 3 285
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 40,399评论 1 211
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 31,717评论 2 310
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 30,407评论 0 194
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 34,112评论 1 239
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 30,371评论 2 241
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 31,891评论 1 256
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 28,255评论 2 250
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 32,881评论 3 233
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 26,010评论 0 8
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 26,764评论 0 192
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 35,412评论 2 269
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 35,299评论 2 260

推荐阅读更多精彩内容