重载高精度模板(使用方便).docVIP

重载高精度模板(使用方便) #include iostream #include string using namespace std; #define DIGIT 4 #define DEPTH 10000 #define MAX 100 typedef int bignum_t[MAX+1]; int read(bignum_t a,istream is=cin){ char buf[MAX*DIGIT+1],ch; int i,j; memset((void*)a,0,sizeof(bignum_t)); if (!(isbuf)) return 0; for (a[0]=strlen(buf),i=a[0]/2-1;i=0;i--) ch=buf[i],buf[i]=buf[a[0]-1-i],buf[a[0]-1-i]=ch; for (a[0]=(a[0]+DIGIT-1)/DIGIT,j=strlen(buf);ja[0]*DIGIT;buf[j++]=0); for (i=1;i=a[0];i++) for (a[i]=0,j=0;jDIGIT;j++) a[i]=a[i]*10+buf[i*DIGIT-1-j]-0; for (;!a[a[0]]a[0]1;a[0]--); return 1; } void write(const bignum_t a,ostream os=cout){ int i,j; for (osa[i=a[0]],i--;i;i--) for (j=DEPTH/10;j;j/=10) osa[i]/j%10; } int comp(const bignum_t a,const bignum_t b){ int i; if (a[0]!=b[0]) return a[0]-b[0]; for (i=a[0];i;i--) if (a[i]!=b[i]) return a[i]-b[i]; return 0; } int comp(const bignum_t a,const int b){ int c[12]={1}; for (c[1]=b;c[c[0]]=DEPTH;c[c[0]+1]=c[c[0]]/DEPTH,c[c[0]]%=DEPTH,c[0]++); return comp(a,c); } int comp(const bignum_t a,const int c,const int d,const bignum_t b){ int i,t=0,O=-DEPTH*2; if (b[0]-a[0]dc) return 1; for (i=b[0];id;i--){ t=t*DEPTH+a[i-d]*c-b[i]; if (t0) return 1; if (tO) return 0; } for (i=d;i;i--){ t=t*DEPTH-b[i]; if (t0) return 1; if (tO) return 0; } return t0; } void add(bignum_t a,const bignum_t b){ int i; for (i=1;i=b[0];i++) if ((a[i]+=b[i])=DEPTH) a[i]-=DEPTH,a[i+1]++; if (b[0]=a[0]) a[0]=b[0]; else for (;a[i]=DEPTHia[0];a[i]-=DEPTH,i++,a[i]++); a[0]+=(a[a[0]+1]0); } void add(bignum_t a,const int b){ int i=1; for (a[1]+=b;a[i]=DEPTHia[0];a[i+1]+=a[i]/DEPTH,a[i]%=DEPTH,i++); for (;a[a[0]]=DEPTH;a[a[0]+1]=a[a[0]]/DEPTH,a[a[0]]%=DEPTH,a[0]++); } void sub(bignum_t a