[Jobdu] 题目1016：火星A+B

 

题目描述：

    读入两个不超过25位的火星正整数A和B，计算A+B。需要注意的是：在火星上，整数不是单一进制的，第n位的进制就是第n个素数。例如：地球上的10进 制数2，在火星上记为“1,0”，因为火星个位数是2进制的；地球上的10进制数38，在火星上记为“1,1,1,0”，因为火星个位数是2进制的，十位 数是3进制的，百位数是5进制的，千位数是7进制的……

 

输入：

    测试输入包含若干测试用例，每个测试用例占一行，包含两个火星正整数A和B，火星整数的相邻两位数用逗号分隔，A和B之间有一个空格间隔。当A或B为0时输入结束，相应的结果不要输出。

 

输出：

    对每个测试用例输出1行，即火星表示法的A+B的值。

 

样例输入：

1,0 2,1

4,2,0 1,2,0

1 10,6,4,2,1

0 0

 

样例输出：

1,0,1

1,1,1,0

1,0,0,0,0,0

 

跟大数运算加法有点像，对于大数运算第一位都是10进制，在取值时，需要模10，在这里只要模上该位对应的素数就OK了。

先建一个30位的素数表，判断n是否素数时只需检查n是否能被1到sqrt(n)之间的素数整除，这些素数可以从表里取得。若n是素数则将其加到表里。

代码如下：

#include <iostream>

#include <vector>

#include <sstream>

#include <cmath>

#include <algorithm>

using namespace std;



string A,B;

int p[30];

vector<int> a,b,c;



void init(){

    p[0] = 2; p[1] = 3; p[2] = 5;

    int c = 3;

    for(int i=6;c<30;i++){

        bool flag = true;

        for(int j=0;p[j]<=sqrt(i)+1;j++){

            if(i%p[j]==0){

                flag = false;

                break;

            }

        }

        if(flag){

            p[c] = i; c++;

        }

    }

}



void toInt(){

    a.clear(); b.clear();

    for(int i=0;i<A.length();i++){

        if(A[i]==',') A[i]=' ';

    }

    for(int i=0;i<B.length();i++){

        if(B[i]==',') B[i]=' ';

    }

    int v;

    istringstream sina(A);

    while(sina>>v){

        a.push_back(v);

    }

    istringstream sinb(B);

    while(sinb>>v){

        b.push_back(v);

    }

    reverse(a.begin(),a.end());

    reverse(b.begin(),b.end());

}



void getRes(){

    c.clear();

    int la = a.size();

    int lb = b.size();

    int v = 0, carry = 0;

    for(int i=0;i<(la>lb?la:lb);i++){

        if(i>=la){

            v = b[i] + carry;

        }

        else if(i>=lb){

            v = a[i] + carry;

        }

        else{

            v = a[i] + b[i] + carry;

        }

        c.push_back(v%p[i]);

        carry = v/=p[i];

        if(i==(la>lb?la:lb)-1&&carry!=0){

            c.push_back(carry);

        }

    }

    reverse(c.begin(),c.end());

    for(int i=0;i<c.size();i++){

        if(i==c.size()-1){

            cout<<c[i]<<endl;

            break;

        }

        cout<<c[i]<<",";

    }

}



int main(int argc,char* argv[]){

    init();

    while(cin>>A>>B&&A!="0"&&B!="0"){

        toInt();

        getRes();

    }

    return 0;

}