高精度练习之除法(codevs_3118)

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
inline int max(int a,int b){return a>b?a:b;}
struct hpc{int a[260],l;};
hpc change(char s[]){
    int len=0,l=strlen(s);
    hpc ans;memset(ans.a,0,sizeof(ans.a));
    for(int i=l-1;i>=0;i-=4){
        int t=0,k=max(i-3,0);
        for(int j=k;j<=i;j++)
            t=t*10+s[j]-'0';
        ans.a[len++]=t;
    }
    ans.l=len;
    return ans;
}
void print(hpc a){
    printf("%d",a.a[a.l-1]);
    for(int i=a.l-2;i>=0;i--)
        printf("%04d",a.a[i]);
    puts("");
}
bool operator <(const hpc a,const hpc b){
    if(a.l<b.l)return 1;
    if(a.l>b.l)return 0;
    for(int i=a.l-1;i>=0;i--){
        if(a.a[i]<b.a[i])return 1;
        if(a.a[i]>b.a[i])return 0;
    }
    return 0;
}
inline hpc operator *(const hpc a,const hpc b){
    hpc ans;int l;
    memset(ans.a,0,sizeof(ans.a));
    if((a.l==1&&a.a[0]==0)||(b.l==1&&b.a[0]==0))
        {ans.l=1;return ans;}
    for(int i=0;i<a.l;i++){
        for(int j=0;j<b.l;j++){
            ans.a[i+j]+=a.a[i]*b.a[j];
            ans.a[i+j+1]+=ans.a[i+j]/10000;
            ans.a[i+j]%=10000;
        }
    }
    l=a.l+b.l-1;
    while(ans.a[l]>0){
        ans.a[l+1]=ans.a[l]/10000;
        ans.a[l]%=10000; l++;
    }
    ans.l=l;
    return ans;
}
inline hpc operator +(const hpc a,const hpc b){
    hpc ans;int l=max(a.l,b.l);
    memset(ans.a,0,sizeof(ans.a));
    for(int i=0;i<l;i++){
        if(i<a.l)ans.a[i]+=a.a[i];
        if(i<b.l)ans.a[i]+=b.a[i];
        ans.a[i+1]=ans.a[i]/10000;ans.a[i]%=10000;
    }
    while(ans.a[l]>0){
        ans.a[l+1]=ans.a[l]/10000;
        ans.a[l]%=10000; l++;
    }
    ans.l=l;
    return ans;
}
inline hpc operator /(const hpc a,const int b){
    hpc ans;int l=a.l;
    memset(ans.a,0,sizeof(ans.a));
    for(int i=a.l-1;i>0;i--){
        ans.a[i-1]=(a.a[i]+ans.a[i])%b*10000;
        ans.a[i]=(a.a[i]+ans.a[i])/b;
    }
    ans.a[0]=(ans.a[0]+a.a[0])/b;
    while(l>1&&ans.a[l-1]==0)l--;
    ans.l=l;return ans;
}
inline hpc operator /(const hpc a,const hpc b){
    hpc l,r=a,mid,t;
    memset(l.a,0,sizeof(l.a));
    memset(t.a,0,sizeof(t.a));
    l.l=t.l=t.a[0]=1;
    while(l+t<r){
        mid=(l+r)/2;
        if(a<mid*b)r=mid;
        else l=mid;
    }
    return l;
}
int main()
{
    char s[510];
    scanf("%s",s);hpc a=change(s);
    scanf("%s",s);hpc b=change(s);
    a=a/b;print(a);
    return 0;
}

本文由Yzyet编写，网址为www.cnblogs.com/Yzyet。非Yzyet同意，禁止转载，侵权者必究。