高精度乘法和除法

http://poj.org/problem?id=2389

#include <iostream>
#include<cstring>
#include<cstdio>
#include<cstdlib>
using namespace std;
char s1[305],a[305];
char s2[305],b[305];
int sum[305];
int main()
{
    memset(sum,0,sizeof(sum));
    scanf("%s%s",s1,s2);
    //cout<<s1<<endl<<s2<<endl;
    int len1=strlen(s1);
    int len2=strlen(s2);
    int s,v;
    for(s=0,v=len1-1;s<len1;s++,v--)
        a[s]=s1[v];
    for(s=0,v=len2-1;s<len2;s++,v--)
        b[s]=s2[v];
   // cout<<a<<endl<<b<<endl;
    for(int i=0;i<len1;i++)
        for(int j=0;j<len2;j++)
        sum[i+j]+=(a[i]-'0')*(b[j]-'0');//最后在确定进位的事
    int jw=0;
    for(int i=0;i<304;i++)
    {

        sum[i]=sum[i]+jw;
        jw=sum[i]/10;
        sum[i]%=10;
    }

    int k;

    for(k=304;sum[k]==0;k--);

    for(int i=k;i>=0;i--)
        printf("%d",sum[i]);
    printf("\n");
    return 0;
}

http://bailian.openjudge.cn/practice/2737

除法：

/********************************************
算法思想：
反复做减法，看看从被除数里最多能减去多少个除数，商就是多少;
所以算法核心是写一个大整数的减法函数;
反复调用该函数进行减法操作;

算法步骤：
用数组a表示被除数，数组b表示除数，数组res表示商;
先用被除数a减去除数b得到差的位数k,同时商+1;
再用被除数a减去若干个除数b*(10^k);
不够减了，再减去若干个除数b*(10^(k-1))...b*(10^(k-2))...;
一直减到不够减为止;
每成功减一次，则商的相应位+1;
最后再循环处理商的进位问题;
********************************************/
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<cstdio>
#include<climits>
#include<algorithm>
using namespace std;

const int N=1000;
int a[N],b[N],res[N*2];
char x[N],y[N];
int len1,len2;

void init()
{
    len1=strlen(x);
    len2=strlen(y);
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    memset(res,0,sizeof(res));
    for(int i=len1-1,j=0; i>=0; i--)
    {
        a[j++]=x[i]-'0';
    }
    for(int i=len2-1,j=0; i>=0; i--)
    {
        b[j++]=y[i]-'0';
    }
}

int substraction(int *p1,int *p2,int len1,int len2)//计算长度为len1的大整数减去长度为len2的大整数的结果的长度
{                                                  //减的结果放在数组p1中，不够返回-1,正好返回0
    if(len1<len2)
        return -1;
    bool flag=0;
    if(len1==len2)
    {
        for(int i=len1-1; i>=0; i--)
        {
            if(p1[i]>p2[i])
                flag=1;
            else if(p1[i]<p2[i])
            {
                if(!flag)
                    return -1;
            }
        }
    }

    for(int i=0; i<len1; i++)//减法
    {
        p1[i]-=p2[i];
        if(p1[i]<0)
        {
            p1[i]+=10;
            p1[i+1]--;
        }
    }
    for(int i=len1-1; i>=0; i--)
    {
        if(p1[i])
            return i+1;
    }
    return 0;
}

void output()
{
    for(int i=0; i<N; i++) //进位
    {
        if(res[i]>=10)
        {
            res[i+1]+=res[i]/10;
            res[i]%=10;
        }
    }

    bool flag=0;
    for(int i=N-1; i>=0; i--)//输出
    {
        if(flag)
        {
            printf("%d",res[i]);
        }
        else if(res[i])
        {
            printf("%d",res[i]);
            flag=1;
        }
    }
    if(!flag)
        printf("0");
    printf("\n");
}

void solve()
{
    init();
    if(len1<len2)
    {
        printf("0\n");
        return;
    }

    len1=substraction(a,b,len1,len2);
    if(len1<0)
    {
        printf("0\n");
        return;
    }
    else if(len1==0)
    {
        printf("1\n");
        return;
    }

    res[0]++;//减掉一次了，商+1
    int k=len1-len2;
    if(k<0)//减一次后不能再减了
    {
        output();
        return;
    }
    else if(k>0)//将数组b乘以10的某次幂，使得其长度与数组a相同
    {
        for(int i=len1-1; i>=0; i--)
        {
            if(i>=k)
                b[i]=b[i-k];
            else
                b[i]=0;
        }
    }

    len2=len1;
    for(int j=0; j<=k; j++)//先减去若干个b*(10^k)，不够减了再减去若干个b*(10^(k-1))...
    {
        int temp;
        while((temp=substraction(a,b+j,len1,len2-j))>=0)//一直减到不够减为止
        {
            len1=temp;
            res[k-j]++;//每成功减一次，则商的相应位+1
        }
    }
    output();
}

int main()
{
    freopen("C:\\Users\\Administrator\\Desktop\\kd.txt","r",stdin);
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%s%s",x,y);
        solve();
    }
    return 0;
}