码迷,mamicode.com
首页 > 其他好文 > 详细

hdu 1042 N!(大数)

时间:2015-08-14 13:31:11      阅读:151      评论:0      收藏:0      [点我收藏+]

标签:

题意:整数大数乘法

思路:大数模板

技术分享
#include<iostream>
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
using namespace std;

#define MAXN 9999//万进制
#define DLEN 4//4位

class BigNum{
private:
    int a[50000];//可以控制大数位数(500*4)
    int len;//大数长度
public:
    BigNum(){//构造函数
        len=1;
        memset(a,0,sizeof(a));
    }
    BigNum(const int);//将int转化为大数
    BigNum(const char*);//将字符串转化为大数
    BigNum(const BigNum &);//拷贝构造函数
    BigNum &operator=(const BigNum &);//重载复制运算符,大数之间赋值

    BigNum operator+(const BigNum &)const;//大数+大数
    BigNum operator-(const BigNum &)const;//大数-大数
    BigNum operator*(const BigNum &)const;//大数*大数
    BigNum operator/(const int &)const;//大数/int

    BigNum operator^(const int &)const;//幂运算
    int operator%(const int &)const;//取模
    bool operator>(const BigNum &)const;//大数与大数比较
    bool operator>(const int &)const;//大数与int比较

    void print();//输出大数
};

BigNum::BigNum(const int b){//将int转化为大数
    int c,d=b;
    len=0;
    memset(a,0,sizeof(a));
    while(d>MAXN){
        //c=d-(d/(MAXN+1))*(MAXN+1);
        c=d%(MAXN+1);//取出后四位
        d=d/(MAXN+1);//
        a[len++]=c;
    }
    a[len++]=d;
}

BigNum::BigNum(const char *s){//将字符串转化为大数
    int t,k,index,l,i,j;
    memset(a,0,sizeof(a));
    l=strlen(s);
    len=l/DLEN;
    if(l%DLEN)++len;
    index=0;
    for(i=l-1;i>=0;i-=DLEN){
        t=0;
        k=i-DLEN+1;
        if(k<0)k=0;
        for(j=k;j<=i;++j)
            t=t*10+s[j]-0;
        a[index++]=t;
    }
}

BigNum::BigNum(const BigNum &T):len(T.len){//拷贝构造函数
    int i;
    memset(a,0,sizeof(a));
    for(i=0;i<len;++i)
        a[i]=T.a[i];
}

BigNum &BigNum::operator=(const BigNum &n){//重载复制运算符,大数之间赋值
    int i;
    len=n.len;
    memset(a,0,sizeof(a));
    for(i=0;i<len;++i)
        a[i]=n.a[i];
    return *this;
}

BigNum BigNum::operator+(const BigNum &T)const{//大数+大数
    BigNum t(*this);
    int i,big;//位数
    big=T.len>len?T.len:len;
    for(i=0;i<big;++i){
        t.a[i]+=T.a[i];
        if(t.a[i]>MAXN){
            ++t.a[i+1];
            t.a[i]-=MAXN+1;
        }
    }
    if(t.a[big]!=0)t.len=big+1;
    else t.len=big;
    return t;
}

BigNum BigNum::operator-(const BigNum &T)const{//大数-大数
    int i,j,big;
    bool flag;
    BigNum t1,t2;//t1大的,t2小的
    if(*this>T){
        t1=*this;
        t2=T;
        flag=0;//前面的大
    }
    else{
        t1=T;
        t2=*this;
        flag=1;//前面的小
    }
    big=t1.len;
    for(i=0;i<big;++i){
        if(t1.a[i]<t2.a[i]){
            j=i+1;
            while(t1.a[j]==0)++j;
            --t1.a[j--];
            while(j>i)t1.a[j--]+=MAXN;
            t1.a[i]+=MAXN+1-t2.a[i];
        }
        else t1.a[i]-=t2.a[i];
    }
    while(t1.a[t1.len-1]==0&&t1.len>1){
        --t1.len;
        --big;
    }
    if(flag)t1.a[big-1]=-t1.a[big-1];//前面的小,结果为负
    return t1;
}

BigNum BigNum::operator*(const BigNum &T)const{//大数*大数
    BigNum ret;
    int i,j,up;
    int temp,temp1;
    for(i=0;i<len;++i){
        up=0;
        for(j=0;j<T.len;++j){
            temp=a[i]*T.a[j]+ret.a[i+j]+up;
            if(temp>MAXN){
                //temp1=temp-temp/(MAXN+1)*(MAXN+1);
                temp1=temp%(MAXN+1);
                up=temp/(MAXN+1);
                ret.a[i+j]=temp1;
            }
            else{
                up=0;
                ret.a[i+j]=temp;
            }
        }
        if(up!=0)ret.a[i+j]=up;
    }
    ret.len=i+j;
    while(ret.a[ret.len-1]==0&&ret.len>1)--ret.len;
    return ret;
}

BigNum BigNum::operator/(const int &b)const{//大数/int
    BigNum ret;
    int i,down=0;
    for(i=len-1;i>=0;--i){
        ret.a[i]=(a[i]+down*(MAXN+1))/b;
        down=a[i]+down*(MAXN+1)-ret.a[i]*b;
    }
    ret.len=len;
    while(ret.a[ret.len-1]==0&&ret.len>1)--ret.len;
    return ret;
}

BigNum BigNum::operator^(const int &n)const{//幂运算
    BigNum t,ret(1);
    int i;
    if(n<0)exit(-1);
    if(n==0)return 1;
    if(n==1)return *this;
    int m=n;
    while(m>1){
        t=*this;
        for(i=1;i<<1<=m;i<<=1){
            t=t*t;
        }
        m-=i;
        ret=ret*t;
        if(m==1)ret=ret*(*this);
    }
    return ret;
}

int BigNum::operator%(const int &b)const{//取模
    int i,d=0;
    for(i=len-1;i>=0;--i){
        d=((d*(MAXN+1))%b+a[i])%b;
    }
    return d;
}

bool BigNum::operator>(const BigNum &T)const{//大数与大数比较
    int ln;
    if(len>T.len)return true;
    else if(len==T.len){
        ln=len-1;
        while(a[ln]==T.a[ln]&&ln>=0)--ln;
        if(ln>=0&&a[ln]>T.a[ln])return true;
        else return false;
    }
    else return false;
}

bool BigNum::operator>(const int &t)const{//大数与int比较
    BigNum b(t);
    return *this>b;
}

void BigNum::print(){//输出大数
    int i;
    printf("%d",a[len-1]);
    for(i=len-2;i>=0;--i){
        printf("%.4d",a[i]);//%.4d代表4位,不够前面补0
    }
    printf("\n");
}

int main(){
    int n,i;
    BigNum a;
    while(~scanf("%d",&n)){
        if(n==0||n==1)printf("1\n");
        else{
            a=BigNum(1);
            for(i=2;i<=n;++i)
                a=a*BigNum(i);
            a.print();
        }
    }
    return 0;
}
View Code

 

hdu 1042 N!(大数)

标签:

原文地址:http://www.cnblogs.com/bofengyu/p/4729408.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!